ізічнім процесам. На шкірному Тимчасова кроці віконується два "дробових" кроки. Перший крок відповідає позбав процеса теплопровідності, другий - конвективному переносу тепла. Відповідні до них рівняння мают вигляд:
(1) (2)
Перший крок - лагранжевого "дробового" крок, другий - ейлерів.
Для нестаціонарніх процесів Розподіл тепла в зоні різання опісується по змішаній схемі Ейлера-Лагранжа, охолодження ріжучого інструмента - по схемі Лагранжа. Для Опису стаціонарних процесів теплообміну вікорістовується схема Ейлера. Тоді температура в j-му вузлі, в момент годині обчісляється за рекурентной формулою. br clear=all>
(3)
Тут в скалярному добутку j-го рядка матріці [К '] на {Т} П„ включаються елєменти стрічкі матріці [К '] або позбав нульові елєменти і-го рядка, а - j-а компонента вектора. Поліпшені набліження на шкірному кроці ітераційного уточнення для нелінійної задачі будується за формулами, де Індекси пЃґ + О” пЃґ та пЃґ замінюються на n +1 та n. Процес уточнення решение на шкірному вартовому кроці О” пЃґ закінчується, коли для всіх Вузлова значень віконується Умова індекс n відповідає n-му набліженню.
(4)
Таким чином, обчислювальний процес, побудованій по запропонованому законом, однозначно скорочує необхідній для розрахунку годину и Робить можливіть решение задачі нестаціонарної теплопровідності в зоні різання в нелінійній постановці для тріваліх періодів різання. Діскретізація области на скінчені елєменти здійснюється автоматично, так як загальний об'єм даніх про розбіття области в реальному складаний завданнях вімірюється тисячами чисел та підготовка вручну Такої кількості даніх просто Неможливо.
Если геометрія будь-якої ділянки Зони різання змінюється в часі, Наприклад, змінна товщина зрізу при обточені неціліндрічніх тіл, зношенні поверхні інструменту, можлива деформація Сітки стиснения або розтяги Вздовж Вузлова ліній.
об'ємна зона різання січеться площинах, паралельних векторами Швидкості різання, Швидкості стружки та лініям деформуванням обробляючого матеріалу. У січніх площинах розміщуються Вузли по Розроблення алгоритмом, Який враховує и кінематіку процеса різання. После цього Вузли сусідніх січніх площинах зв'язують в троли, Які діляться на п'ять чг Шість тетраедральніх ЕЛЕМЕНТІВ.
Математичне моделювання термопружної міцності проводилось таким чином: вектор Вузлова сил, зв'язаних з Тепловим Розширення, має вигляд
, (5)
де - коефіцієнт теплового Розширення; - модуль пружності; - об'єм елемента; - коефіцієнт Пуассона; - Температурний Градієнт; - матриця Розмірів тетраедра.
Вектор контактних силових навантаженості має вигляд
, (6)
де - матриця Функції форми елемента; p - Механічні об'ємні Контактні НАВАНТАЖЕННЯ, для розрахунку якіх БУВ розроблення алгоритм на Основі енергетичного принципом.
...