Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом
The initial-boundary problem of thermomechanics for a bimetallic spherical element under homogenous non-stationary electromagnetic action is formulated. The azimuthal component of the magnetic field strength vector, temperature, and the radial component of the displacement vector were selected as ke...
Збережено в:
| Дата: | 2023 |
|---|---|
| Автори: | , , , |
| Формат: | Стаття |
| Мова: | Українська |
| Опубліковано: |
Інститут прикладних проблем механіки і математики ім. Я. С. Підстригача НАН України
2023
|
| Теми: | |
| Онлайн доступ: | https://www.fmmit.lviv.ua/index.php/fmmit/article/view/242 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Physico-mathematical modeling and informational technologies |
| Завантажити файл: | |
Репозитарії
Physico-mathematical modeling and informational technologies| _version_ | 1867479610729431040 |
|---|---|
| author | Мусій, Роман Мельник, Наталія Орищин, Оксана Шиндер, Валентин |
| author_facet | Мусій, Роман Мельник, Наталія Орищин, Оксана Шиндер, Валентин |
| author_institution_txt_mv | [
{
"author": "Роман Мусій",
"institution": "Національний університет „Львівська політехніка”"
},
{
"author": "Наталія Мельник",
"institution": null
},
{
"author": "Оксана Орищин",
"institution": null
},
{
"author": "Валентин Шиндер",
"institution": null
}
] |
| author_sort | Мусій, Роман |
| baseUrl_str | http://www.fmmit.lviv.ua/index.php/fmmit/oai |
| collection | OJS |
| datestamp_date | 2023-06-26T23:17:44Z |
| description | The initial-boundary problem of thermomechanics for a bimetallic spherical element under homogenous non-stationary electromagnetic action is formulated. The azimuthal component of the magnetic field strength vector, temperature, and the radial component of the displacement vector were selected as key functions. To find the key functions, the approximation of their distributions in the constituent layers of the spherical element by quadratic polynomials in the radial variable is proposed. This approximation makes it possible to accurately satisfy the specified boundary conditions both on the surfaces of this element and on the surface of the connection of the component layers, where ideal electromagnetic, thermal and mechanical contacts take place. The coefficients of the polynomials approximating the key functions in the component layers of the spherical element are given by a linear combination of functions describing all given boundary conditions and integral characteristics of the key functions on both component layers. As a result, the original initial-boundary value problems for the defining functions are reduced to Cauchy problems for their integral (summed over the package of layers) characteristics. With the use of Laplace transformation in time, the general solutions of these problems under homogeneous non-stationary electromagnetic action are recorded. The change in time of radial and azimuthal stresses and stress intensity was numerically analyzed by technological processing with an electromagnetic impulse, as well as the performance and properties of the contact connection of the bimetallic spherical element were investigated. |
| doi_str_mv | 10.15407/fmmit2022.34-35.014 |
| first_indexed | 2026-06-09T01:09:01Z |
| format | Article |
| fulltext |
14
УДК 539.3
doi.org/10.15407/fmmit2022.34-35.014
Математичне моделювання і аналіз термомеханічної
поведінки біметалевого сферичного елемента при
технологічній обробці електромагнітним імпульсом
Роман Мусій1, Наталія Мельник2, Оксана Орищин3, Валентин Шиндер4
1 д. ф.-м. н., професор, Національний університет „Львівська політехніка”, вул. С. Бандери, 12, Львів, 79013,
e-mail: musiy@lp.edu.ua
2 к. ф.-м. н., доцент, Національний університет „Львівська політехніка”, вул. С. Бандери, 12, Львів, 79013,
e-mail: melnyk.n.pz@gmail.com
3 к. ф.-м. н., доцент, Національний університет „Львівська політехніка”, вул. С. Бандери, 12, Львів, 79013,
e-mail: oksana.h.oryshchyn@lpnu.ua
3 к. ф.-м. н., доцент, Національний університет „Львівська політехніка”, вул. С. Бандери, 12, Львів, 79013,
e-mail: valentyn.k.shynder@lpnu.ua
Сформульовано початково-крайову задачу термомеханіки для біметалевого сферичного
елемента за однорідної нестаціонарної електромагнітної дії. За ключові функції вибрано
азимутальну компоненту вектора напруженості магнітного поля, температуру та радіальну
компоненту вектора переміщень. Для знаходження ключових функцій запропоновано
апроксимацію їх розподілів в складових шарах сферичного елемента квадратичними поліномами
по радіальній змінній. Дана апроксимація дає змогу точно задовольнити задані граничні умови
як на поверхнях даного елемента, так і на поверхні з’єднання складових шарів, на яких мають
місце ідеальні електромагнітний, тепловий і механічний контакти. Коефіцієнти поліномів, що
апроксимують ключові функції у складових шарах сферичного елемента подаються лінійною
комбінацією функцій, що описують всі задані граничні умови та інтегральних характеристик
ключових функцій по обох складових шарах. У результаті вихідні початково-крайові задачі на
визначальні функції зведено до задач Коші на їх інтегральні (сумарні по пакету шарів)
характеристики. З використанням перетворень Лапласа за часом записано загальні розв’язки
цих задач за однорідної нестаціонарної електромагнітної дії. Числово проаналізовано за
технологічної обробки електромагнітним імпульсом зміну в часі радіальних і азимутальних
напружень та інтенсивності напружень, а також досліджено роботоздатність і
властивості контактного з’єднання біметалевого сферичного елемента.
Ключові слова. Біметалевий сферичний елемент, електромагнітна обробка,
електромагнітний імпульс, термомеханічна поведінка, роботоздатність.
Вступ. У сучасних механізмах і пристроях широко використовують кусково-
однорідні електропровідні структурні елементи, які перебувають в умовах
багатофакторних навантажень [1]. Одним з таких елементів, особливо в
електромагнітних перехідниках, є біметалевий сферичний елемент, який в процесі
виготовлення та експлуатації зазнає впливу імпульсних електромагнітних полів
(ЕМП) [2-4]. Від збереження його роботоздатності залежить надійність роботи
відповідних пристроїв, які містять біметалеві сферичні елементи [2, 5].
В літературі відомі дослідження термонапруженого стану та несучої здатності
однорідних тіл канонічної форми за нестаціонарних силових і теплових дій [6-9]. У
монографії [10] вивчено термопружний стан однорідних електропровідних тіл за
індукційного нагріву усталеними та квазіусталеними ЕМП. У працях [11, 12]
вивчено термонапружения стан однорідних електропровідних тіл канонічної форм за
ISSN 1816-1545 Фізико-математичне моделювання та інформаційні технології
2022, вип. 34-35, 14-23
15
дії електромагнітного імпульсу (ЕМІ) та електропровідних тіл з плоскопаралельними
межами за дії імпульсних ЕМП з модуляцією амплітуди. У статтях [13, 14]
розглянуто термонапружений стан біметалевого шару за імпульсної
електромагнітної з модуляцією амплітуди та за дії електромагнітного імпульсу.
Проте недостатньо вивчено термомеханічну поведінку біметалевого сферичного
елемента за дії ЕМІ.
В даній роботі здійснено математичне моделювання термомеханічної
поведінки біметалевого сферичного елемента за дії однорідного нестаціонарного
ЕМП. Числово проаналізовано за дії ЕМІ зміну в часі радіальних і азимутальних
напружень та інтенсивності напружень, а також досліджено несучу здатність і
властивості контактного з’єднання біметалевого сферичного елемента.
1. Формулювання задачі та розрахункова модель
Розглядається біметалевий сферичний елемент з внутрішнім 0rr і зовнішнім
2rr радіусами, віднесений до сферичної системи координат ,,r , центр O якої
співпадає з центром елемента. Матеріали складових шарів елемента однорідні
ізотропні неферомагнітні, а їх фізико-механічні характеристики постійні. Елемент
перебуває за дії нестаціонарного ЕМП, заданого значеннями азимутальної
компоненти nH
вектора напруженості магнітного поля 0;;0 nn HH
на його
внутрішній і зовнішній поверхнях, які теплоізольовані від довкілля і вільні від
зовнішнього силового поверхневого навантаження. На поверхні 1rr з’єднання
складових шарів елемента виконуються умови ідеальних електромагнітного,
теплового і механічного контактів [10, 15]. Всі визначальні функції задачі –
компонента
trH n , вектора nH
, температура
trT n , та відмінна від нуля
радіальна компонента
tru n
r
, вектора переміщень
0;0;n
r
n uu
- є функціями
радіальної координати r та часу t .
Розрахункова схема визначення термонапруженого стану елемента за дії
імпульсного ЕМП складається з трьох етапів [10, 15]. На першому етапі зі
співвідношень Максвелла визначається нестаціонарне ЕМП, яке описується
вектором nH
, та зумовлені ним тепло Джоуля nQ і пондеромоторні сили nF
в
кожному n - му 2,1n складовому шарі за врахування заданих граничних і
контактних умов. На другому етапі з рівняння теплопровідності, в якому тепло
Джоуля nQ є об’ємним джерелом тепла, визначається нестаціонарне температурне
поле nT за конвективного теплообміну з довкіллям і заданих контактних умов. На
третьому етапі з рівнянь динамічної термопружності за знайденими температурою
nT та пондеромоторною силою nF
визначаємо вектор переміщень
nu
та тензор
динамічних напружень
n̂ за заданих граничних і контактних умов.
Для оцінки несучої здатності елемента визначаємо інтенсивності
n
i в n - му
шарі сумарних напружень
Роман Мусій, Наталія Мельник, Оксана Орищин, Валентин Шиндер Математичне моделювання і аналіз
термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці
електромагнітним імпульсом
16
Fn
jj
Qn
jj
n
jj ,,rj , (1)
які обчислюється за формулою [16]
2/))ˆ()ˆ(3( )(2
1
)(
2
)( nnn
i II , (2)
і порівнюємо їх величину з межею пружності
n
d матеріалу n - го шару. Тут
Qn
jj і Fn
jj - напруження, зумовлені теплом Джоуля nQ і пондеромоторною
силою nF
, )(ˆ n
jI )2,1( j – j-й інваріант тензора напружень
n̂ .
Біметалевий сферичний елемент зберігає несучу здатність, якщо в обох його
шарах виконується умова [15]
n
d
n
i . (3)
Відповідно властивості контактного з’єднання зберігаються за виконання на
поверхні 1rr з’єднання складових шарів елемента умови
M
n
i
, (4)
де M – межа міцності контактного з’єднання.
Згідно вибраної розрахункової схеми розглянемо послідовно визначення ЕМП,
температури та термонапруженого стану біметалевого сферичного елемента.
2. Вихідні рівняння та співвідношення
2.1. Визначення ЕМП. Азимутальна компонента
trH n , вектора nH
в кожному
n -му складовому шарі елемента визначається з рівняння [11]
0
2
2
2
t
H
r
H
rr
H n
nn
nn
(5)
за граничних умов на внутрішній і зовнішній поверхнях елемента
tHtrH ,
0
1 , tHtrH ,
2
2 (6)
та умов ідеального електромагнітного контакту на поверхні з’єднання шарів
елемента
trHtrH ,,
1
2
1
1
,
1
1
2
1
2
1
1
1
1
1
,
,
,
,
r
trH
r
trH
k
r
trH
r
trH
, (7)
а також початкових умов
00, rH n
. (8)
ISSN 1816-1545 Фізико-математичне моделювання та інформаційні технології
2022, вип. 34-35, 14-23
17
Тут
n
- коефіцієнт електропровідності,
n
- магнітна проникливість матеріалу n -го
шару елемента,
21
k , t – час.
За знайденою функцією trH n , питомі густини тепла Джоуля trQ n , та
пондеромоторної сили 0;0;,trFF
n
r
n
визначаються за формулами
2
1
r
H
r
H
Q
nn
n
n
,
n
nn
H
r
H
r
H
F n
n
r
. (9)
2.2. Визначення температурного поля. На другому етапі за знайденою питомою
густиною тепла Джоуля trQ n , визначаємо розподіл температури trT n , в
кожному n -му шарі елемента з рівняння теплопровідності [7,11]
n
nnnn Q
t
T
r
T
rr
T
1
2
n
2
2
, (10)
за крайових умов теплоізоляції на внутрішній і зовнішній поверхнях елемента
0
, 0
1
r
trT
,
0=
, 2
2
r
trT
(11)
та умов ідеального теплового контакту на поверхні з’єднання його складових шарів
trTtrT ,, 1
2
1
1 ,
r
trT
k
r
trT
,
, 1
2
1
1
, (12)
а також за початкових умов
00, rT n . (13)
Тут n
,
n
- коефіцієнти температуро- і теплопровідності матеріалу n -го шару
елемента,
12
k .
2.3. Визначення термонапруженого стану. Для визначення відмінних від нуля
радіальної n
rr , азимутальної n
та меридіанної n
компонент тензора
динамічних напружень
trn ,̂ в кожному n - му шарі елемента за вихідну
вибираємо систему рівнянь центрально-симетричної задачі термопружності в
переміщеннях. Тоді радіальна компонента tru n
r
, вектора переміщень в n - му шарі
елемента визначається з рівняння [7, 9]
12
2
2
2
222
2
t
u
c
u
rr
u
rr
u n
r
n
n
r
n
r
n
r
n
r
nn
nn
n
n
n
n
F
Er
T
1
211
1
1
(14)
за крайових умов на внутрішній і зовнішній поверхнях елемента
Роман Мусій, Наталія Мельник, Оксана Орищин, Валентин Шиндер Математичне моделювання і аналіз
термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці
електромагнітним імпульсом
18
trTtru
rr
tru
r
r ,
1
1
,
2
1
,
0
1
1
1
10
1
01
10
1
, (15)
trTtru
rr
tru
r
r ,
1
1
,
2
1
,
2
2
2
2
22
2
22
22
2
,
умов ідеального механічного контакту на поверхні з’єднання її шарів
trutru r ,, 1
2
1
1
, (16)
,,1,
2
,
1
,1,
2
,
1
1
2
221
2
1
1
1
2
2
1
1
111
1
1
1
1
1
1
trTtru
rr
tru
k
trTtru
rr
tru
r
r
E
r
r
а також за початкових умов
00, ru
n
r
,
0
0,
r
ru
n
r . (17)
Тут 21
1
2111
nnnnnn
Ec - швидкість пружної хвилі розширення в
n -му шарі елемента при його ізотермічній деформації; n
, n
- коефіцієнти
лінійного теплового розширення та Пуассона, n
E - модуль Юнга, n
- густина
матеріалу n -го шару, 12
EE
E
k . За знайденою компонентою tru n
r
, вектора
переміщень nu
визначаємо компоненти
n
jj
тензора напружень n̂ в n - му шарі
елемента за формулами [7,9]
n
nn
n
rn
n
r
n
nn
n
rr
Tu
rr
u
1
2
1
)21(1
E
n , (18)
n
nn
n
r
n
r
n
nn
nn T
r
u
r
u
1
)21(1
E
n .
За відомими сумарними переміщеннями n
r
u і напруженнями
n
jj
записуємо
інтенсивності напружень n
i
в n - му шарі елемента, які будуть
nn
rr
n
i . (19)
3. Методика побудови розв’язків початково-крайових задач.
Для розв’язування сформульованих початково-крайових задач електродинаміки,
теплопровідності і термопружності стосовно визначальних функцій
ISSN 1816-1545 Фізико-математичне моделювання та інформаційні технології
2022, вип. 34-35, 14-23
19
trutrTtrH n
r
nnn ,,,,, )()()()(
апроксимуємо їх розподіли за радіальною
змінною r в n - му шарі елемента квадратичними многочленами [13-15]
2
0
)()(
,
i
in
i
n
rtatr . (20)
Коефіцієнти ta n
i
)( апроксимаційних поліномів (20) визначаємо через інтегральні
(сумарні по пакету шарів елемента) характеристики t
s
ключових функцій
trn ,)(
2
1
1
,
n
r
r
lsn
s
n
n
drtrt , 2,1s (21)
та задані крайові умови на внутрішній і зовнішній поверхнях і на поверхні з’єднання
шарів елемента. Тут 0l для визначення trH n ,)(
, trT n ,)( і 1l для визначення
tru n
r
,)( .
Для знаходження інтегральних характеристик t
s
вихідні рівняння (5), (10),
(14) інтегруємо відповідно до (21) за врахування виразів (20). В результаті вихідні
початково-крайові задачі на ключові функції зведено до відповідних задач Коші на
інтегральні характеристики цих функцій, які описуються системами рівнянь
,
872615
2
432211
1
tHdtHdtHdtHd
dt
tdH
tHdtHdtHdtHd
dt
tdH
(22)
,*
22413
2
*
12211
1
tQTdTd
dt
dT
tQTdTd
dt
dT
TT
TT
(23)
.*
22
*
71
*
62
2
2
*
12
*
21
*
12
1
2
tWudud
dt
ud
tWudud
dt
ud
rr
r
rr
r
(24)
Тут коефіцієнти
81
d , Td
41
, *
41
d визначаються через геометричні параметри та
фізико-механічні характеристики складових шарів елемента, tQ
s
* , tW
s
* ( 2,1s ) –
праві частини рівнянь (10) і (14), проінтегровані відповідно до (21) з урахуванням
виразів
)(nT і )(n
r
F на поверхні з’єднання шарів.
Роман Мусій, Наталія Мельник, Оксана Орищин, Валентин Шиндер Математичне моделювання і аналіз
термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці
електромагнітним імпульсом
20
Системи (22)-(24) розв’язуються з використанням інтегрального перетворення
Лапласа за початкових умов, що відповідають заданим початковим умовам (8), (13),
(17) на ключові функції.
Розв’язки задач Коші (22)-(24) за однорідної нестаціонарної електромагнітної
дії подаємо у вигляді згорток функцій, що описують задані крайові умови та
однорідні розв’язки цих задач.
Отримуємо вирази компоненти trH n ,
t
t
k
p
ss
k
n
is
si
n deHkAHkAatrH
0
21
2
1=
)(
2
1=
2
0=
)( =,
in
i
n
i
rtHatHa
)(
4
)(
3 , (25)
температури trT n ,)(
k
k s m
t
tpQ
s
Q
s
n
ks
n rdeWmBWmBbtrT m
2
0
2
1
2
1 0
2211
)()( ,
(26)
та складників tru Qn
r
, і tru Fn
r
, радіальної компоненти tru n
r
, вектора
переміщень
nu
i
i s k
t
t
k
p
ss
n
is
n
r
rdeWkRWkRctru
2
0=
2
=1
4
1 0
*
2
*
2
*
1
*
1
)(=,
. (27)
За відомими виразами (27) функцій tru n
r
, зі співвідношень (18) визначаємо
радіальні n
rr
, азимутальні n
і меридіанні n
компоненти тензора напружень та
інтенсивності напружень )(n
i
у біметалевому сферичному елементі.
4. Числовий аналіз за дії електромагнітного імпульсу
Дія ЕМІ математично описується функцією tH
у вигляді [11, 12]
tt
eeHktH 21
00
. (28)
Тут 0k - нормувальний множник, 1 і 2 – параметри, що характеризують
відповідно часи фронтів наростання incr
t і спадання
decrt ЕМІ тривалості it , 0H -
максимальне значення напруженості магнітного поля в імпульсі.
Підставляючи вираз (28) у формули (25)-(27), записуємо розв’язки задачі
термомеханіки для розглядуваного біметалевого сферичного елемента за дії ЕМІ.
Проаналізовано термонапружений стан і несучу здатність біметалевого
сферичного елемента з радіусами поверхонь ммr 8
0
, ммr 9
1
, ммr 10
2
,
ISSN 1816-1545 Фізико-математичне моделювання та інформаційні технології
2022, вип. 34-35, 14-23
21
складові шари якого виготовлені із нержавної сталі і міді. Характеристики ЕМІ
приймались рівними 410it
с, 310 с, i
t9,6
1
, 12
2 , 4
0
k , МПа
d
300 –
для сталі і МПа
d
70 – для міді, МПа
M
250230 [17].
На рис. 1–3 показано зміну в часі складників радіальних n
rr
і азимутальних
n
напружень та інтенсивностей напружень )(n
i
у сталевому (1) і мідному (2)
шарах елемента на поверхні їх з’єднання.
Рис. 1. Зміна в часі радіальних
напружень на поверхні з’єднання
сталевого і мідного шарів.
Рис. 2. Зміна в часі азимутальних напружень
на поверхні з’єднання сталевого і мідного
шарів
На рис.4 подано залежності максимальних значень
n
i max інтенсивностей
напружень в n-му шарі від величини 0H за різних тривалостей ЕМІ.
Рис. 3. Зміна в часі інтенсивностей
напружень на поверхні з’єднання
сталевого і мідного шарів.
Рис. 4. Залежність максимальних значень
n
i max
інтенсивностей напружень від величини 0
H
напруженості магнітного поля
Роман Мусій, Наталія Мельник, Оксана Орищин, Валентин Шиндер Математичне моделювання і аналіз
термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці
електромагнітним імпульсом
22
З аналізу наведених на рис. 4 залежностей отримано, що роботоздатність
розглядуваного біметалевого сферичного елемента за його технологічної обробки
ЕМІ втрачається при 41,5 10it
с. Відповідно при 41,5 10it
c також втрачаються
і властивості контактного з’єднання.
5. Висновки
Запропоновано методику побудови загальних розв’язків початково-крайових задач
електродинаміки, теплопровідності і термопружності для біметалевого сферичного
елемента. Методика використовує апроксимацію всіх ключових функцій
(азимутальної компоненти вектора напруженості магнітного поля, температури та
радіальної компоненти вектора переміщень) в кожному складовому шарі елемента
по радіальній змінній квадратичними поліномами. У результаті вихідні початково-
крайові задачі на визначальні функції зведено до задач Коші на їх інтегральні
характеристики. Розв’язки цих задач за однорідної нестаціонарної електромагнітної
дії знайдено з допомогою перетворення Лапласа у вигляді згорток функцій, що
описують задані крайові умови та однорідні розв’язки. Чисельно проаналізовано
радіальні і азимутальні напруження та інтенсивності напружень за технологічної
обробки біметалевого сферичного елемента з допомогою ЕМІ.
Виявлено закономірності термомеханічної поведінки розглядуваного елемента
за дії ЕМІ і встановлено граничні значення параметрів ЕМІ, за перевищення яких
біметалевий сферичний елемент втрачає роботоздатність та властивості контактного
з’єднання.
Література
1. F. Matthews, R. Rawlings, Composite materials: engineering and science. London: Chapman&Hall, 1994,
pp. 20–46.
2. Батыгин Ю.В., Лавинский В.И., Хименко Л.Т. Импульсные магнитные поля для прогрессивных
технологий – Харьков: МОСТ – Торнадо, 2003. – 288 с.
3. G.A.Shneerson, M.I.Dolotenko, S.I.Krivosheev, Strong and Superstrong Pulsed Magnetic Fields
Generation. De Gruyter, 2010.
4. F. Herlach. Pulsed Magnets, Rep. Prog. Phys. 62 (6), 1999, pp. 859–920.
5. Suneet Singh, Prashant K. Jain, Rizwan-uddin. Analytical Solution for Three-Dimensional, Unsteady Heat
Conduction in a Multilayer Sphere J. Heat Transfer. Oct 2016, 138(10): 101301 (11 p.)
6. И.А. Мотовиловец, В.И. Козлов. Механика связанных полей в элементах конструкций, т.1:
Термоупругость. – К.: Наук. думка, 1987. – 264 с.
7. Подстригач Я.С., Коляно Ю.М. Обобщенная термомеханика. - Киев: Наукова думка, 1976. - 312 с.
8. Лурье А.И. Теория упругости. – М.: Наука, 1970. – 940с.
9. Грибанов В.Ф., Паничкин Н.Г. Связанные и динамические задачи термоупругости.-
М.:Машиностроение, 1984- 184 с.
10. Я.С. Подстригач, Я.И. Бурак, А.Р. Гачкевич, Л.В. Чернявская Термоупругость электропроводных
тел. – К.: Наук. думка, 1977. – 248 с.
11. Мусій Р. С. Динамічні задачі термомеханіки електропровідних тіл канонічної форми. – Львів: Вид-
во «Растр-7», 2010. – 216 с.
12. Гачкевич О. Р., Мусій Р. С., Тарлаковський Д.В. Термомеханіка неферомагнетних електропровідних
тіл за дії імпульсних електромагнетних полів з модуляцією амплітуди. – Львів: СПОЛОМ, 2011. –
216 с.
13. Hachkevych O.R., Musij R., Melnyk N., Szymczak J. Rearing capability of bimetallic plate at effect of
pulse electromagnetic field with amplitude modulation // Optimization of manufacturing processes and
javascript:;
javascript:;
javascript:;
ISSN 1816-1545 Фізико-математичне моделювання та інформаційні технології
2022, вип. 34-35, 14-23
23
work environ-ment/ Ed/ by M.Gajek/. Studia i monografie z. 276. Ch. 6. Opole: Oficyna Wydawnicza
Politechniki Opolskiej. 2010. – Pp. 57-70.
14. Musii R., Melnyk N., Dmytruk V. Mathematical modeling and study of thermoelastic behavior of a
bimetallic layer with plane-parallel boundaries under the action of electromagnetic pulses // Journal of
thermal stresses, 2021. Vol. 44, no. 10, P. 1275–1288.
15. R.Musii, N.Melnyk, V.Dmytruk, Thermoelastic processes analyzer for piecewise homogeneous conductive
structures subjected to pulsed electromagnetic action, Journal of Thermal Stresses, 41 (9), 2018, pp. 1125–
1135.
16. Ионов В.Н., Огибалов П.М. Напряжения в телах при импульсивном нагружении. – М.: Высш. школа,
1975. – 463 с.
17. Баженов В.Г., Петров М.В. О применении
магнитоимпульсного способа деформирования для исследования вязкопластических характеристик
материалов. // Прикладные проблемы прочности и пластичности. Методы решения задач упругости
и пластичности. Всесоюз. межвуз. сб. Горьк. ун-та, – 1980. – С. 18-25.
Mathematical modeling and analysis of the thermomechanical
behavior of a bimetallic spherical element during technological
processing with an electromagnetic impulse
Roman Musij, Nataliya Melnyk, Oksana Oryshchyn, Valentyn Shynder
The initial-boundary problem of thermomechanics for a bimetallic spherical element under
homogenous non-stationary electromagnetic action is formulated. The azimuthal component of the
magnetic field strength vector, temperature, and the radial component of the displacement vector were
selected as key functions. To find the key functions, the approximation of their distributions in the
constituent layers of the spherical element by quadratic polynomials in the radial variable is proposed.
This approximation makes it possible to accurately satisfy the specified boundary conditions both on
the surfaces of this element and on the surface of the connection of the component layers, where ideal
electromagnetic, thermal and mechanical contacts take place. The coefficients of the polynomials
approximating the key functions in the component layers of the spherical element are given by a linear
combination of functions describing all given boundary conditions and integral characteristics of the
key functions on both component layers. As a result, the original initial-boundary value problems for
the defining functions are reduced to Cauchy problems for their integral (summed over the package of
layers) characteristics. With the use of Laplace transformation in time, the general solutions of these
problems under homogeneous non-stationary electromagnetic action are recorded. The change in time
of radial and azimuthal stresses and stress intensity was numerically analyzed by technological
processing with an electromagnetic impulse, as well as the performance and properties of the contact
connection of the bimetallic spherical element were investigated.
Отримано 23.06.22.
|
| id | oai:ojs2.www.fmmit.lviv.ua:article-242 |
| institution | Physico-mathematical modeling and informational technologies |
| keywords_txt_mv | keywords |
| language | Ukrainian |
| last_indexed | 2026-06-09T01:09:01Z |
| publishDate | 2023 |
| publisher | Інститут прикладних проблем механіки і математики ім. Я. С. Підстригача НАН України |
| record_format | ojs |
| resource_txt_mv | wwwfmmitlvivua/cc/a83c5111aa7da3743cd1233b4f5eb3cc.pdf |
| spelling | oai:ojs2.www.fmmit.lviv.ua:article-2422023-06-26T23:17:44Z Mathematical modeling and analysis of the thermomechanical behavior of a bimetallic spherical element during technological processing with an electromagnetic impulse Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом Мусій, Роман Мельник, Наталія Орищин, Оксана Шиндер, Валентин Біметалевий сферичний елемент, електромагнітна обробка, електромагнітний імпульс, термомеханічна поведінка, роботоздатність. The initial-boundary problem of thermomechanics for a bimetallic spherical element under homogenous non-stationary electromagnetic action is formulated. The azimuthal component of the magnetic field strength vector, temperature, and the radial component of the displacement vector were selected as key functions. To find the key functions, the approximation of their distributions in the constituent layers of the spherical element by quadratic polynomials in the radial variable is proposed. This approximation makes it possible to accurately satisfy the specified boundary conditions both on the surfaces of this element and on the surface of the connection of the component layers, where ideal electromagnetic, thermal and mechanical contacts take place. The coefficients of the polynomials approximating the key functions in the component layers of the spherical element are given by a linear combination of functions describing all given boundary conditions and integral characteristics of the key functions on both component layers. As a result, the original initial-boundary value problems for the defining functions are reduced to Cauchy problems for their integral (summed over the package of layers) characteristics. With the use of Laplace transformation in time, the general solutions of these problems under homogeneous non-stationary electromagnetic action are recorded. The change in time of radial and azimuthal stresses and stress intensity was numerically analyzed by technological processing with an electromagnetic impulse, as well as the performance and properties of the contact connection of the bimetallic spherical element were investigated. Сформульовано початково-крайову задачу термомеханіки для біметалевого сферичного елемента за однорідної нестаціонарної електромагнітної дії. За ключові функції вибрано азимутальну компоненту вектора напруженості магнітного поля, температуру та радіальну компоненту вектора переміщень. Для знаходження ключових функцій запропоновано апроксимацію їх розподілів в складових шарах сферичного елемента квадратичними поліномами по радіальній змінній. Дана апроксимація дає змогу точно задовольнити задані граничні умови як на поверхнях даного елемента, так і на поверхні з’єднання складових шарів, на яких мають місце ідеальні електромагнітний, тепловий і механічний контакти. Коефіцієнти поліномів, що апроксимують ключові функції у складових шарах сферичного елемента подаються лінійною комбінацією функцій, що описують всі задані граничні умови та інтегральних характеристик ключових функцій по обох складових шарах. У результаті вихідні початково-крайові задачі на визначальні функції зведено до задач Коші на їх інтегральні (сумарні по пакету шарів) характеристики. З використанням перетворень Лапласа за часом записано загальні розв’язки цих задач за однорідної нестаціонарної електромагнітної дії. Числово проаналізовано за технологічної обробки електромагнітним імпульсом зміну в часі радіальних і азимутальних напружень та інтенсивності напружень, а також досліджено роботоздатність і властивості контактного з’єднання біметалевого сферичного елемента. Інститут прикладних проблем механіки і математики ім. Я. С. Підстригача НАН України 2023-03-09 Article Article application/pdf https://www.fmmit.lviv.ua/index.php/fmmit/article/view/242 10.15407/fmmit2022.34-35.014 PHYSICO-MATHEMATICAL MODELLING AND INFORMATIONAL TECHNOLOGIES; No. 34-35 (2022): PHYSICO-MATHEMATICAL MODELLING AND INFORMATIONAL TECHNOLOGIES; 14-23 ФІЗИКО-МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ ТА ІНФОРМАЦІЙНІ ТЕХНОЛОГІЇ; № 34-35 (2022): ФІЗИКО-МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ ТА ІНФОРМАЦІЙНІ ТЕХНОЛОГІЇ; 14-23 2617-5258 1816-1545 10.15407/fmmit2022.34-35 uk https://www.fmmit.lviv.ua/index.php/fmmit/article/view/242/219 |
| spellingShingle | Біметалевий сферичний елемент електромагнітна обробка електромагнітний імпульс термомеханічна поведінка роботоздатність. Мусій, Роман Мельник, Наталія Орищин, Оксана Шиндер, Валентин Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title | Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title_alt | Mathematical modeling and analysis of the thermomechanical behavior of a bimetallic spherical element during technological processing with an electromagnetic impulse |
| title_full | Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title_fullStr | Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title_full_unstemmed | Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title_short | Математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| title_sort | математичне моделювання і аналіз термомеханічної поведінки біметалевого сферичного елемента при технологічній обробці електромагнітним імпульсом |
| topic | Біметалевий сферичний елемент електромагнітна обробка електромагнітний імпульс термомеханічна поведінка роботоздатність. |
| topic_facet | Біметалевий сферичний елемент електромагнітна обробка електромагнітний імпульс термомеханічна поведінка роботоздатність. |
| url | https://www.fmmit.lviv.ua/index.php/fmmit/article/view/242 |
| work_keys_str_mv | AT musíjroman mathematicalmodelingandanalysisofthethermomechanicalbehaviorofabimetallicsphericalelementduringtechnologicalprocessingwithanelectromagneticimpulse AT melʹniknatalíâ mathematicalmodelingandanalysisofthethermomechanicalbehaviorofabimetallicsphericalelementduringtechnologicalprocessingwithanelectromagneticimpulse AT oriŝinoksana mathematicalmodelingandanalysisofthethermomechanicalbehaviorofabimetallicsphericalelementduringtechnologicalprocessingwithanelectromagneticimpulse AT šindervalentin mathematicalmodelingandanalysisofthethermomechanicalbehaviorofabimetallicsphericalelementduringtechnologicalprocessingwithanelectromagneticimpulse AT musíjroman matematičnemodelûvannâíanalíztermomehaníčnoípovedínkibímetalevogosferičnogoelementapritehnologíčníjobrobcíelektromagnítnimímpulʹsom AT melʹniknatalíâ matematičnemodelûvannâíanalíztermomehaníčnoípovedínkibímetalevogosferičnogoelementapritehnologíčníjobrobcíelektromagnítnimímpulʹsom AT oriŝinoksana matematičnemodelûvannâíanalíztermomehaníčnoípovedínkibímetalevogosferičnogoelementapritehnologíčníjobrobcíelektromagnítnimímpulʹsom AT šindervalentin matematičnemodelûvannâíanalíztermomehaníčnoípovedínkibímetalevogosferičnogoelementapritehnologíčníjobrobcíelektromagnítnimímpulʹsom |