CC BY
5
2019 Математика и механика № 57
удк 539.3
doi 10.17223/19988621/57/7
В.М. Пестренин, И.В. Пестренина, Л.В. Ландик, Е.Р. Эйсмонт
КОМПОНЕНТЫ НАПРЯЖЕНИЙ НА ЛИНИИ ОСОБЫХ ТОЧЕК СВОБОДНОЙ ОТ НАГРУЗКИ ТОРЦЕВОЙ ПОВЕРХНОСТИ СЛОИСТОГО ЦИЛИНДРА
На основе представления об особой точке деформируемого тела в виде бесконечно малой частицы, получаемой в результате стягивания к ней элементарного объема, изучаются ограничения на параметры состояния в точках края поверхности соединения двухслойного цилиндра. Устанавливаются варианты возможных постановок задач механики деформируемого твердого тела. Выявляются критические сочетания геометрических и материальных параметров, приводящие к сингулярному характеру напряженного состояния в особых точках.
Ключевые слова: особые точки, неклассические задачи, двухслойный цилиндр, температурная нагрузка.
В настоящее время для исследования полей напряжений вблизи особых точек деформируемых тел (вершин клиньев, конусов, многогранников, точек края поверхностей соединения тел и т.п.) применяются два подхода. Первый из них (далее классический или асимптотический) разрабатывается авторами публикаций [1-17 и др.]. Классический подход характерен тем, что особая точка в нем исключается из области поиска решения. Это делается путем помещения в особую точку полюса криволинейной системы координат (полюс не может быть точкой области построения решения, так как в нем отсутствует однозначное соответствие между координатами и точкой тела). Исключение особой точки приводит к неадекватному определению напряжений в ее окрестности, так как при этом не рассматриваются условия непосредственно в особой точке. Кроме того, в данном случае отсутствует механический смысл получаемого асимптотического решения, поскольку невозможно указать элементарный объем, в котором такое решение реализуется. Альтернативный (неклассический) подход к изучению напряжений в особых точках и их окрестностях предлагается в публикациях [18-24]. В этом подходе особая точка считается бесконечно малой частицей среды (тела), получающейся путем стягивания в нее элементарного объема. Такое представление о точке сплошной среды разработано учеными XVIII века (Даниил и Иоганн Бер-нулли, Ж. Л. Д'Аламбер и Л.Эйлер [25]) и используется современными исследователями. Элементарным объемом в рассматриваемом контексте является континуальная модель представительного объема изучаемого тела. Параметры состояния элементарного объема однородны и сохраняют свои значения при стягивании в бесконечно малую частицу. Поэтому параметры состояния в точке сплошной среды (в частности, в особой точке) совпадают с параметрами состояния содержащего ее элементарного объема. Изучение напряженного состояния в особых точках и их окрестностях на основе альтернативного подхода показывает, что необычность (уникальность) этих точек в деформируемых твердых телах проявляется в избыточном количестве (по сравнению с обычными точками поверхности те-
ла) задаваемых в них ограничений [18-24]. Данное обстоятельство обусловливает неклассическую постановку задач механики деформируемого твердого тела, содержащего особые точки. Неклассические (в указанном смысле) задачи рассматривались в работах [18] - однородные плоские клинья, [20, 21]- составные плоские клинья, [22] - составные пространственные ребра, [23] - круговые и составные конусы, [24] - внутренние особые точки плоских элементов конструкций.
В настоящей статье неклассический подход используется для изучения параметров состояния на линии особых точек, расположенной на свободной торцевой поверхности составного цилиндра и представляющей собой ребро, образованное пересечением торцевых поверхностей составляющих элементов. Устанавливаются соответствующие постановки задач механики деформируемого твердого тела (МДТТ). Приводится решение задачи о температурной нагрузке двухслойного полого цилиндра, согласующееся со всеми задаваемыми в особых точках ограничениями.
1. Постановка задачи
Рассматривается двухслойный цилиндр. Его торцевое и осевое сечения представлены на рис. 1 и 2. Принимается, что торцевые поверхности составляющих цилиндров 1 и 2 образуют ребро Г. Считается, что плоскости ф = const являются нормальными плоскостями ребра Г. Углы между касательными к линиям пересечения нормальной плоскости ребра Г и общей образующей цилиндров 1 и 2 обозначаются соответственно а и р (рис. 2).
Рис. 1. Торцевое сечение двухслойного цилиндра Fig. 1. End section of a double-layered cylinder
Рис. 2. Осевое сечение двухслойного цилиндра Fig. 2. Axial section of a double-layered cylinder
Область изменения параметров а и р задается неравенствами:
0 <а <п, 0 <Р<л, 0 <а + р<2п. (1)
Через п , т обозначаются нормали в точках ребра Г к касательным плоскостям поверхностей, образующих ребро. Через п ', т' обозначаются орты, перпендикулярные соответственно направлениям п , т . Приняты обозначения: Ек, Ок, vk, юк (к = 1,2) - модули Юнга, модули сдвига, коэффициенты Пуассона и температурной деформации составляющих цилиндров; ст^), е^) (к = 1,2) - физиче-
ские компоненты напряжений и деформаций; сп, тп1, тп2 - нормальное и касательные напряжения в направлении соответственно n' и Г2 (Т2 — орт, определяющий угловое направление цилиндрической системы координат) на площадке, ориентируемой ортом n ; ст, тт1, тт2 - нормальное и касательные напряжения в направлении соответственно т' и Г2 на площадке, ориентируемой ортом т; AT - приращение температуры. Считается, что на торцевом сечении отсутствует поверхностная нагрузка. К ребру (тА) примыкают элементарные объемы скрепленных цилиндров 1 и 2. В соответствии с принятой концепцией эти элементарные объемы при стягивании к точке А образуют особые точки. Напряжения и деформации в особых точках совпадают с напряжениями и деформациями соответствующих элементарных объемов. Для параметров состояния в рассматриваемых элементарных объемах выполняются следующие условия:
1) равенство нулю поверхностных нагрузок:
с„ = 0, х„1 = 0, х„2 = 0; Cm = 0, ^ = 0, Тт2 = 0; (2)
2) непрерывность напряжений на поверхности соединения составляющих цилиндров
с® =C(2} =C C(1) =C(2) =C C(1) =C(2) =C • (3)
rr rr rr' r) rф r)' rz rz rz ' v >
3) непрерывность деформаций на поверхности соединения составляющих элементов (следует из совпадения поверхностей скрепляемых тел с поверхностью соединения)
P(1) = P(2) p(1) = p(2) _(1) = _(2) zz zz ' фф Ьфф > )z * (4)
Через физические компоненты тензоров напряжений равенства (2) с использованием равенств (3) в цилиндрической системе координат запишутся соотношениями
crr cos2 a- 2crz cos а sin а + С1 sin2 а = 0, (crr -с® )cos а sina + crz (cos2 а - sin2 а) = 0, crr cos2 p + 2crz cos p sin p + с^ sin2 p = 0, ( - Crr) cos p sin p + Crz (cos2 p - sin2 p) = 0;
—c cos а + с® sin а = 0,
(2 (6) Cr) cos P + CVsin P = 0.
Уравнения (5) являются автономной системой линейных однородных уравнений относительно компонент cz^,с^,crr,crz. Определитель этой системы уравнений
A = - sin а sin р sin(а + P) (7)
в области (1) изменения параметров а, р обращается в нуль в точках
а + Р = л. (8)
Равенства (4) с использованием физических уравнений термоупругости и равенств (3) преобразуются к виду
^ (( _ У1СТгг -у1СТфф )ДТ = Е- ( — V ) + ; (9)
Е1 Е2
( - v1аrr ) + Ю1ДТ = ( ^агг ^а® ) + Ю2ДТ ; (10)
Е1 Е2
7 а^ф = Е аС^. (11)
Задача состоит в исследовании возможных решений построенных уравнений (5), (6), (9) - (11) и выявлении в зависимости от геометрических и материальных параметров независимых ограничений на параметры состояния в особых точках.
2. Построение ограничений на параметры состояния в особых точках
Ограничениями на параметры состояния в особых точках являются решения систем уравнений (5) и (6), (11). В этом пункте рассматриваются возможные случаи таких решений. Каждому случаю отвечает независимая постановка задачи, потому что ограничения в каждом случае различаются и количеством и содержанием.
2.1. Случай а + в * п.
Определитель (7) системы уравнений (5) не обращается в нуль. Ее решение
а® =а® =агг = аГГ = 0. (12)
Равенства (9), (10) приводятся к системе уравнений относительно напряжений
»2) АТ, -1 (1) 1
Определитель этой системы
а(1) и а(2) фф фф
Е - Е- афф = К - ®2 )ДТ, - Е афф + Е_ афф = ( -„2 )ДТ. (13)
Е1 Е2 Е1 Е2
А! = -1- (V! -V 2 ) (14)
Е1Е2
обращается в нуль при выполнении равенства у1 = у2 , поэтому возможны случаи
1) (ю1 - ю2)АГ = 0. Система уравнений (13) однородна.
а) если у1 ф V2, ее решение тривиально афф = а)2) = 0.
б) если v1 =v2, уравнения (13) линейно зависимы. Напряжения афф,аф^ф^ связаны равенством
К2афф - Е1афф) = 0. (15)
2) (ю1 - ю2)ЛГ * 0.
а) если v1 ф V2, система уравнения (13) имеет решение
а(1) = I^ )ДТ (1 + V 2 ) Е1 а(2) = I.2 ) ДТ (1 + v1) е2
фф I \ , фф I \ . (16)
( ^ 2 ) ( -v 2 )
б) если v1 = V2, ранг расширенной матрицы уравнений (13) оказывается большим, чем ранг системы, поэтому она несовместна. Напряжения афф, аф2ф при стремлении v1 к v2 в особых точках стремятся к неограниченно большим значениям.
Равенства (6), (11) образуют систему линейных однородных уравнений относительно напряжений стГф, ст®, ст^2-1. Определитель этой системы
Д2 = ^sinрcos a+^sina cos p
r) > zф > zф
1 • n 1 •
—sin p cos а+—
G1 G2
обращается в нуль при условии
Р = 0. (17)
ц ь2
Если условия (17) не выполняются, уравнения (6), (11) имеют тривиальное решение
стгф = 0, ст® = о, ст(,ф = 0. (18)
Если условие (17) выполняется, ранг системы уравнений (6), (11) равен двум, [ряжения ст®, ст$ могут быть выражены через стгф:
с® = Сгф ctg а, с® = -Сгф ctg р. (19)
Zф Гф о Zф Гф
Рассмотренные выше случаи сведены в табл.1.
Таблица 1
Ограничения на компоненты напряжений в особых точках для случая а+
№ Условия на материальные параметры и ДТ Ограничения (дополнительно к равенствам (3)) Количество ограничений
1 2 3 4
1 -Utg a + -^ctg G1 G2 (ю1 -ю2)ДТ = 0, v1 ^v2 ctZZ = 0, CTZ2) = 0, CTrz = 0, СТГГ = 0, Стгф= 0, ctZ— = 0,ст® = 0, стфф = 0, стфф = 0 фф ' фф Двенадцать
2 ~ctg a + -^ctg p*0, G1 G2 (ю1 -ю2)ДТ = 0, v1 = v2 ctZZ> = 0, CTZ2) = a CTrz = ^ СТГГ = 0 Стгф = 0, ctZ— = 0, ст<ф = 0, ЫФ - £1стф2ф) = 0 Одиннадцать
3 ~ctg a + -^ctg p^0, G1 G2 (ю1 -ю2)ДТ ^ 0, v1 ^ v2 ctZZ = 0, ctZ2) = 0, CTrz = 0, CTrr = 0, Ст ф = 0,ст(1ф = 0,ст(ф = 0, Гф 5 Zф 5 Zф 5 ст(1) (a1 -a2 )ДТ (1 +V 2 )E1 фф~ (V1 -V 2 ) , CT(2) (a1 -a2 )ДТ (1 +V1 )E2 фф_ (V1 -V2 ) Двенадцать
4 ~ctg a + -^ctg p^0, G1 G2 (ю1 -ю2)ДТ ^ 0, v1 = v2 Ограничения несовместны —
Окончание табл. 1
1 2 3 4
5 ~ctg а + "~ctg Р = 0, G1 G2 (ю1 -ю2)АГ = 0, v1 ф V2 ст® = 0, ст® = 0, стгс = 0, ст„ = 0, ст!ф = стгф ctg а, = -стгф ctg P, ст(1 ) = 0ст( 2 ) = 0 фф ' фф Одиннадцать
6 ~ctg а + 7~ctg Р = 0, G1 G2 (ю1 -ю2)аГ = 0, v1 = V2 ст® = 0, ст^ = 0, стгс = 0, стгг = 0, ^ф = стгф ctg а, ст?ф = -стгф ctg P, ^фф - Е1стф2ф) = 0 Десять
7 ~ctg а + Р = 0, G1 G2 (ю1 -ю2)аг ф 0, v1 ф v2 ст« = а ст® = 0, = 0, ст„ = (x ^ф = стгф ctg ^ = -стгф ctg P, ст(1) (а1 -а2 )АТ (1 +V 2 )E1 ФФ_ (V -V 2 ) , _(2) (а1 -а2 )АТ (1 +V1 )E2 ФФ_ (V1 -V2 ) Одиннадцать
8 7^ctg а + Р = 0, G1 G2 (ю1 -ю2)аг ф 0, v1 = V2 Ограничения несовместны —
Таким образом, для случая а + рфл возможны шесть различных вариантов формулировки ограничений в особых точках и, следовательно, шесть различных постановок задач МДТТ. Во всех этих вариантах, различающихся условиями для материальных и геометрических параметров, количество ограничений в особых точках превосходит количество обычно принимаемых авторами в подобных задачах ограничений равное девяти (равенства (2), (3)). Это увеличение связано с включением в список ограничений соотношений (4), которые в случае различных материальных параметров составляющих цилиндров не выполняются тождественно.
2.2. Случай а + в = п.
Уравнения (5), (6) принимают вид
стгг cos2 а- 2стгг cos а sin а + ст® sin2 а = 0,
(сг -ст® )cos а sin а + ст„ (cos2 а-sin2 а) = 0,
2 (2) 2 (20) стгг cos а- 2ст rz cos а sin а + ст^ ;srn а= 0,
-(ст® -стгг )cos а sin а + ст rz (cos2 а- sin2 а) = 0;
-стг(р cos а + ст® sin а = 0, -стг(р cos а + ст(г2^ш а = 0. (21)
Ранг матрицы системы уравнений (20) равен трем. Следовательно, три компонента напряжений выражаются через четвертый:
ст® = ctg2 астгг, ст(2) = ctg2 астгг, стгг = ctg астГГ. (22)
Из уравнений (21) следуют зависимости
ст2( = стгф ctgа, ст™ = стгФ ctgа. (23)
Соотношения (9), (10), (11) запишутся равенствами
2iCT(i) v
7 фф Л
1 (,) 1
E СТФФ - Y СТФФ) + F°rr = ( - ®2 )АТ,
-EСТФФ + ТСТФФ) + Н= (Ю1 -Ю2 )АТ, E1 E2
2
G - G i0«»ctg a=o
(24)
(25)
(26)
В этих равенствах обозначено
F = -E-(ctg2 a-Vj)+ -E-(ctg2 a-V2), H = "^](ctg2 a + 1).
Возможны следующие случаи решений системы уравнений (24) и (25)
1) Vj Ф V2. Ранг системы уравнений (24), (25) равен двум. Напряжения
стфф, стф2(р можно выразить через напряжение стгг и приращение температуры:
(1) = (F + Hv 2 ) E°rr + К -Ю2 )АТ (1 + v 2 ) E1 .
стфф=-
(V -V 2 )
(V -V 2 )
(2) = (F + H V1 ) E2CTrr + К -Ю2 )AT (1 + V1 ) E2
стфф =-
(v1 -V 2 )
(v1 -V2 )
(27)
(28)
ctg2 a
2) v1 = V2 Ф-2—. Ранг системы уравнений (24), (25) равен двум. Опреде-
1 + ctg2 a ляется напряжение
К -®2 )ATE1E2
(29)
(30)
rr (E -Е2)[ctg2 a-v(l + ctg2 a)]
и линейная комбинация
_L (1)__1_ (2) = - ctg2 a(ai -a2 )AT
Ei w E2 w ctg2 a-v(l + ctg2 a)' ctg2 a
3) v1 = v2 =--—. Ранг системы уравнений (24), (25) равен единице. Ранг
1 + ctg2 a
расширенной матрицы равен двум, если (ю1 -ю2)АГ Ф 0 . Ограничения на параметры состояния в особых точках несовместны. Когда геометрические и материальные параметры элемента конструкции стремятся к данному сочетанию, напряжения в особых точках неограниченно возрастают. При выполнении условия (ю1 -ю2)АГ = 0 ранги матрицы системы уравнений (24), (25) и ее расширенной матрицы совпадают, уравнения становятся зависимыми. Любое из них определяет
связь между напряжениями ст
(1) фф'
СТфф и CTrr .
Из равенства (26) следует, что при выполнении какого-либо из условий
G1 = G2, ctg a = 0 (31)
на напряжения стгф ограничения не накладываются. Если ни одно из условий (31)
не выполняется, напряжения стгф подчиняются ограничению
стгф = 0. (32)
Результаты данного пункта сведены в табл. 2.
Таблица 2
Ограничения в особых точках для случая a+ß = n
№ Условия на материальные параметры и AT Ограничения (дополнительно к равенствам (1.3)) Количество ограничений
1 2 3 4
стге = ctga^rr, ст® = ctg2aarr, ) = ctg2aCT„.,
1 (G2 -G1)ctga^ 0, (ю1 -ro2)AT ф 0, V1 ф v2 0, = 0, ст% = 0, c(1) (F +Hv 2 )Е1стгг +(a1 -a2 )AT (1+v 2 )E1 ФФ (v1-V 2) (v1-V 2) ' ст(2) (F +Hv1)E2CTrr +(a1 -a 2 )AT (1+V1)E2 ФФ (v1 -V 2 ) (v1 -V 2 ) Одиннадцать
CTrz = ctgaCTrr, СТ® = ctg2aCTrr, CTz? = ctg2«CTrr,
(G2 -G1)ctga^ 0, CTr( = 0, CT« = 0, CT™ = 0,
2 (ю1 -ro2)AT ф 0, ctg2 a V1 =V 2 =v* + 2 1+ctg2 a 1 (1) 1 (2) _ -ctg2 a(«1 -a 2 )AT 77 СТФФ 77 СТФФ = , 2 Л , , 2 V E1 E2 ctg a-v(1+ctg a) (a1 -a2 )ATE1E2 CTrr =(E1 -E2 )(ctg2 a-v(1+ctg2 a)) Одиннадцать
(G2 —G1)ctgaФ 0,
3 (ю1 -ro2)AT ф 0, ctg2 a V1 =V 2 =v = + 2 1+ctg2 a Ограничения несовместим: —
4 (G2 -G1)ctgaФ 0, (ю1 —ro2)AT = 0, V1 ф V2 CTrz = ctgaCTrr , CT<ZZ = ctg2aCTrr , CT2 = ctg2«CTrr , CT Ф= 0,CT(1) = 0,CT(() = 0, r ф 5 Zф 5 Zф 5 CT(1) (F +Hv2 )E1CTrr CT(2) (F +Hv1)E2CTrr ФФ_ (v1 —v 2) , ФФ (v1 -v2) Одиннадцать
(G2 -G1)ctgaФ 0, CTrz = ctgaCTrr , CTZZ' = ctg2aCTrr , CT2 = ctg2«CTrr ,
5 (ю1 —ro2)AT = 0, ctg2 a V1 =V 2 =v* + 2 1+ctg2 a CTr( = 0, ctZ( = 0, ctZ( = 0, 1 CT() -1ct(2 ) = 0 CT = 0 E ФФ E ФФ r -C1 2 Одиннадцать
(G2 -G1)ctgaФ 0, CTrz = ctgaCTrr , CTZZ' = ctg2aCTrr , CTZ2) = ctg2«CTrr ,
6 (ю1 —ro2)AT = 0, ctg2 a V1 =V 2 =v = + 2 1+ctg2 a Стгф = 0, ctZ( = 0, ctZ( = 0, —E ^Фф +1 ^Ф2^ +H CT rr = 0 -C1 2 Десять
Окончание табл. 2
1 2 3 4
7 (G2 -G1)ctga=0, (ю1 -ю2)ат ф 0, V1 ф v2 arz = ctgа^гг, ^ = ctg2aarr, ) = rtg2«^, =стгфctgа, ^ =стирctgа, ст(1 ) (F +Hv 2 ,(a1 -а2 )AT (1+V 2 )E1 ФФ (v1-V 2) (v1-V 2) ' ст(2) = (F +Hv1)E2CTrr ,(a1 -а2)AT(1+v1)E2 ФФ (v1 -V 2 ) (v1 -V 2 ) Десять
8 (G2 -G1)ctga=0, (ю1 -c»2)AT ф 0, ctg2 а V1 =V 2 =v* + 2 1+ctg2 а rz = ^g^rr , ^ = ctg2, CTZ2) = ctg2«CTrr , ^ф =стгф ^^ =стгф Ctg^ 1 (1) 1 (2) -ctg2 а(а1 -а 2 )AT Е%ф Е^фф ctg2а-v(l+ctg2а), (а1 -а2 )ATE1E2 CTrr =(E1 -E2 )(cctg2 а-v(l+ctg2 а)) Десять
9 (G2 -G1)ctga=0, (ю1 -ю2)аГ ф 0, аФп/2 ctg2 а V1 =V 2 =v = + 2 1+ctg2 а Ограничения несовместим: —
10 (G2 -G1)ctgа=0, (ю1 -ю2)аГ = 0, V1 ф V2 rz = ^g^rr , ^ = ctg2астгг , CTZ2) = ctg2«CTrr , ^ =стгфCtga, ^ =стгфCtga, _(1) ( +H v2 )E1CTrr ст(2) (F +Hv1)E2CTrr фф_ (v1 -v 2) , фф (v1 -v2) Десять
11 (G2 -G1)ctgа=0, (ю1 -ю2)аГ = 0, CTrz = ^g^rr , ^ = ctg2астгг , CTZ2) = ^g^rr , ^ =стгфCtga, ^ =стгфCtga Десять
ctg2 а V1 =V 2 =v* + 2 1+ctg2 а 1 ст ф ф -1-CT 2 ) = 0 E фф E фф -C1 2 CT rr = 0
12 (G2-G1)=0, аФп/2 (ю1 -ю2)аГ = 0, ctg2 а V1 =V 2 =V= + 2 1+ctg2 а Соединяемые элементы имеют одинаковые материальные параметры. Постановка упругой задачи совпадает с постановкой в классическом случае
13 G2 ф G1, (ю1 -ю2)аГ = 0, v1 =v 2 =v а=п/2 ст = 0, ст(1) = 0, ст(2) = 0, rzzzzz ст(1 ) = 0 ст( 2 ) = 0 --ст^ ) +1-стф2 ) + Hст = 0 ^ф ' ^ф ' р фф р фф rr E1 E2 Девять
Из табл. 2 видно, что для случая а + р = п возможны одиннадцать совместных вариантов формулировки ограничений в особых точках. В варианте 12 материальные параметры О и V соединяемых элементов совпадают, поэтому он отвечает нагружению однослойного цилиндра (при АТ = 0 различие коэффициентов ю1 и ю2 не сказывается на решении задачи). Равенства (24) - (26) выполняются тождественно. Задача становится классической. Во всех остальных вариантах количество ограничений в особых точках превосходит количество ограничений, вытекающих из равенств (2), (3), принимаемых во внимание в классическом подходе.
3. Пример. Температурное нагружение двухслойного полого цилиндра
Рассматривается температурное нагружение двухслойного полого цилиндра (рис. 1 и 2) при следующих геометрических и материальных параметрах: Я = 0.1 м, Н1 = 0.005 м, И2 = 0.001 м, Н = 0.3 м, Е1 = 100 ГПа, v1 = 0.3, v2 = 0.28, ю1 = 0.510-5 1/град, ю2 = 0.810-5 1/град, АТ = 100°С, а = р = л/2. Жесткость внешнего цилиндра Е2 варьируется (60-146 ГПа). В табл. 2. (случай 1) выписаны ограничения, накладываемые на физические компоненты напряжений. Решение задачи строится численно-аналитическим итерационным методом, описанным в работе [19]. Характерный размер конечно-элементной сетки вблизи точек ребра был равен 0.510-5 м. Результаты вычислений приводятся на рис. 3 - 4.
На рис. 3 показаны радиальные напряжения вблизи линии особых точек в торцевом сечении. Интересно, что знак этих напряжений зависит от жесткости цилиндров. При повышении температуры в точках края поверхности соединения напряжения растягивающие, если внутренний цилиндр более жесткий, и сжимающие, - если менее жесткий. Это обстоятельство можно использовать, к примеру, в задачах защиты цилиндрических конструкций от коррозии.
<rr, МПа 200
100
0 -100 -200
0.104995 0.104998 0.105000 0.105003 r, м
Рис. 3. Напряжения <rr в точках торцевого сечения: сплошная линия во внутреннем цилиндре, штриховая - во внешнем. 1 - E2 = 60 ГПа; 2 - E2 = 146 ГПа Fig. 3. Stress values <rr at the points of end section:
the solid line corresponds to the inner cylinder; the dashed, to the outer cylinder. E2 = (1) 60 and (2) 146 GPa
2 \ y'\
/V \
\ Д I \ t L
v !
1
На рис. 4, а показаны напряжения стгг в составляющих цилиндрах вблизи линии особых точек в направлении общей образующей. Видно, что эти напряжения в окрестности особой точки имеют значительную величину, но равны нулю в самой точке, что согласуется с соотношениями, представленными в табл. 2.
стгг, МПа 300 200 100
0 -100
астфф, МПа
40 30 20 10
1 а
~ i Л\
____- A Z \
2 V
1 b
2 \
\\
\\
0
0.029995
0.029997
0.029999
Рис. 4. Напряжения в окрестности особой точки вдоль образующей поверхности соединения цилиндров: 1 - внутренний цилиндр; 2 - внешний цилиндр, E2 = 60 ГПа: а - azz; b - невязка выполнения равенств (27), (28). Сплошная линия - итерационное решение; штриховая линия - ANSYS-решение Fig. 4. Stress values in the vicinity of singular point along a generator of cylinders' joint surface: 1, inner cylinder and 2, outer cylinder at E2 = 60 GPa: a) azz and b) discrepancy in the equalities (27) and (28). The solid line indicates the iterative solution; the dashed line, the ANSYS-solution
z, м
На рис. 4, Ь показана умноженная на разность ^ -V2) невязка выполнения равенств (27), (28). На линии особых точек торцевого сечения эти равенства должны выполняться, поэтому невязка в этих точках должна обращаться в нуль. Видно, что при приближении к линии особых точек невязка по равенствам (27), (28) в решении задачи стремится к нулю. В А№У8-решении (классический подход), построенном на той же конечно элементной сетке, ограничения (27), (28) на компоненты напряжений в особой точке не выполняются. Вне малой окрестности особой точки классическое и неклассическое решения совпадают.
Заключение
Представление о точке сплошной среды в виде бесконечно малой частицы, получаемой в результате стягивания к ней элементарного объема, распространяется на особые точки деформируемых тел. Предлагаемый подход позволяет сформулировать в особых точках задаваемые условия (ограничения) на параметры состояния. Показано, что количество таких ограничений в рассматриваемой задаче определяется материальными и геометрическими параметрами элемента конструкции и, как правило, превышает количество граничных условий в обычной точке поверхности тела. Это обстоятельство обусловливает неклассическую постановку задач МДТТ. Установлены сочетания материальных и геометрических параметров, указывающие на критическое состояние элемента конструкции в том смысле, что при выполнении таких сочетаний напряжения в особых точках неограниченно возрастают. С использованием итерационного численно-аналитического метода [19] исследовано напряженное состояние двухслойного цилиндра при температурной нагрузке. Показано, что полученное решение на линии особых точек торцевой поверхности согласуется со всеми задаваемыми в этих точках ограничениями.
Приведенные в работе результаты найдут применение в исследовании кромочных эффектов слоистых конструкций, механике композитов, изучении концентрации напряжений в точках края поверхностей клеевых соединений, и т.п.
ЛИТЕРАТУРА
1. Уфлянд Я.С. Интегральные преобразования в задачах теории упругости. М.; Л.: Изд-во АН СССР, 1967. 402 c.
2. Bogy D.B. Two Edge-bonded Elastic Wedges of Different Materials and Wedge Angles under Surface Tractions // Trans. ASME. Ser. E. 1971. V. 38. No. 2. P. 87-96. DOI: 10.1115/ 1.3408786.
3. Чобанян К.С. Напряжения в составных упругих телах. Ереван: Изд-во АН АрмССР, 1987. 338 с.
4. Sinclear G.B. Stress singularities in classical elasticity - I: Removal, interpretation and analysis // App. Mech. Rev. 2004. V. 57. No. 4. P. 251-297.
5. Sinclear G.B. Stress singularities in classical elasticity - II: Asymptotic identification // App. Mech. Rev. 2004. V. 57. No. 4. P. 385-439.
6. Barut A., Guven I., Madenci E. Analysis of singular stress fields at junctions of multiple dissimilar materials under mechanical and thermal loading // Int. J. of Solid and Structures. -2001. V. 38. No. 50-51. P. 9077-9109.
7. Шемякин Е.И. О краевых задачах теории упругости для областей с угловыми точками (плоская деформация) // ДАН. - 1996. - Т. 347. № 3. С. 342-345.
8. Задоян М.А. Прочность соединения составных плит // Механика твердого тела. - 2003. № 1. - С.111-122.
9. Paggi M., Carpintery A. On the stress singularities at multimaterial interfaces and related analogies with fluid dynamics and diffusion // Mech. Rev., 61 (2008) 020801. http://dx.doi.org/10.1115/1.2885134.
10. Xu L.R., Kuai H., Sengupta S. Dissimilar material joints with and without free-edge stress singularities: Part I. A Biologically Inspired Design // Experimental mechanics. 2004. V. 44. No. 6. P. 608-615.
11. Xu L.R., Kuai H., Sengupta S. Dissimilar material joints with and without free-edge stress singularities: Part II. An integrated numerical analysis // Experimental mechanics. 2004. V. 44. No. 6. P. 616-621.
12. Коваленко М.Д., Меньшова И.В., Шуляковская Т.Д. Разложения по функциям Фадля -Папковича. Примеры решений в полуполосе // Изв. РАН. МТТ. 2013. № 5. C. 136-158.
13. Luangarpa C., Koguchi H. Analysis of singular stresses at a vertex and along a singular line in three-dimensional bonded joints using a conservative integral // Europ. J. of Mechanics -A/Solids. 2016. V. 60. P. 208-216.
14. Кондратьев В.А. Краевые задачи для эллиптических уравнений в областях с коническими или угловыми точками // Тр. Моск. матем. об-ва. 1967. Т. 16. С. 209-292.
15. Мазья В.Г., Пламеневский Б.А. О коэффициентах в асимптотике решений эллиптических краевых задач вблизи ребра // ДАН СССР. 1976. Т. 229. № 1. С. 33-36.
16. Партон В.З., Перлин П.И. Методы математической теории упругости. М.: Наука, 1981. 688 с.
17. Kozlov V.A., Mazya V.G., Schwab C. On singularities of solutions of the displacement problem of linear elasticity near the vertex of a cone // Arch. Ration. Mech. Anal. 1992. No. 119. P. 197-227.
18. Пестренин В.М., Пестренина И.В., Ландик Л.В. Нестандартные задачи для однородных элементов конструкций с особенностями в виде клиньев в условиях плоской задачи // Вестник Томского государственного университета. Математика и механика. 2014. № 1(27). С. 95-109.
19. Пестренин В.М., Пестренина И.В., Ландик Л.В. Итерационный конечно-элементный алгоритм исследования напряженного состояния элементов конструкций с особыми точками и его реализация // Вестник Пермского национального исследовательского политехнического университета. Механика. 2015. № 4. С. 171-187. DOI: 10.15593/ perm.mech/2015.4.11
20. Pestrenin V. M., Pestrenina I. V., Landik L.V. Stress state at the vertex of a composite wedge, one side of which slides without friction along a rigid surface // Latin American Journal of Solids and Structures, 2017. V. 14. No. 11. P. 2067-2088. http://dx.doi.org/10.1590/1679-78253826.
21. Пестренин В.М., Пестренина И.В., Ландик Л.В. Напряженно-деформированное состояние в окрестности вершины составного клина с жестко защемленными образующими // Вестник Пермского национального исследовательского политехнического университета. Механика. 2016. № 3. С. 131-147. DOI: 10.15593/perm.mech/2016.3.09.
22. Pestrenin V.M., Pestrenina I.V., Landik L.V. Nonstandart problems for structural elements with spatial composite ribs // Mechanics of Composite Materials. 2015. V. 51. No. 4. P. 489-504.
23. Пестренин В.М., Пестренина И.В., Ландик Л.В. Ограничения на параметры напряженного состояния в вершине кругового конуса // Вестник Томского государственного университета. Математика и механика. 2018. № 52. С. 89-101.
24. Pestrenin V.M., PestreninaI.V. Constraints on stress components at the internal singular point of an elastic compound structure // Mechanics of Composite Materials. 2017. V. 53. No. 1. P. 107-116. DOI: 10.1007/s11029-017-9644-1. http://link.springer.com/article/10.1007/ s11029-017-9644-1
25. Ковалев Б.Д. Формирование Эйлеровой гидродинамики // Исследования по истории механики. Сборник статей. М.: Наука, 1983. С. 146 -167.
Статья поступила 05.09.2018 г.
Pestrenin V.M., Pestrenina I.V., Landik L.V., Eismont E.R.(2019) STRESS COMPONENTS NEAR THE SINGULAR POINTS AT A LOAD-FREE END SECTION OF A LAYERED CYLINDER. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. 57. pp. 85-99
DOI 10.17223/19988621/57/7
Keywords: singular points, non-classical problems, double-layered cylinder, temperature load.
In this paper, special points of a composite deformable body (points on the edge of joint surface) are considered as ordinary points of continuous medium representing infinitely small particles obtained by contracting elementary volumes towards these points. Such an approach supposes the singular points to be located in the solution area and makes it possible to formulate
the conditions (restrictions) for state parameters in these points. It is shown that the number of restrictions on the line of singular points on the end face of double-layered cylinder is defined by material and geometric parameters of the structural element, and, generally, exceeds the number of boundary conditions specified at the ordinary points of body surface. This fact specifies a non-classical formulation of the problem of deformable solid mechanics. In this work, various statements of the solid mechanics problems are developed for a considered construction element. The critical combinations of geometric and material parameters leading to a singular nature of the stress state at the singular points are revealed. A solution to the problem of double-layered cylinder under temperature loading is obtained using the numerical-analytical iterative method. The results presented could be applied in mechanics of composites and during the study of edge effects in the layered structures or stress concentration at the points of surface edge of glue joints.
PESTRENIN Valery Mikhaylovich (Candidate of Physics and Mathematics, Perm State National Research University, Perm, Russian Federation). E-mail: PestreninVM@mail.ru
PESTRENINA Irena Vladimirovna (Candidate of Technical Science, Perm State National Research University, Perm, Russian Federation). E-mail: IPestrenina@gmail.com
LANDIK Lidiya Vladimirovna (Perm State National Research University, Perm, Russian Federation). E-mail: LidiaLandik@gmail.com
EISMONT Elena Romanovna (Perm State National Research University, Perm, Russian Federation). E-mail: eis-elena@yandex.ru
REFERENCES
1. Uflyand Ya.S. (1967) Integral'nye preobrazovaniya v zadachakh teorii uprugosti [Integral transformations in problems of the theory of elasticity]. Moscow; Leningrad: AN SSSR Publ.
2. Bogy D.B. (1971) Two edge-bonded elastic wedges of different materials and wedge angles under surface tractions. Trans. ASME. Ser.E. 38(2). pp. 87-96. D01:10.1115/1.3408786.
3. Chobanyan K.S. (1987) Napryazheniya v sostavnykh uprugikh telakh [Stresses in the composite elastic bodies]. Yerevan: AN ArmSSR Publ.
4. Sinclear G.B. (2004) Stress singularities in classical elasticity - I: Removal, interpretation and analysis. App. Mech. Rev. 57(4). pp. 251-297.
5. Sinclear G.B. (2004) Stress singularities in classical elasticity - II: Asymptotic identification. App. Mech. Rev. 57(4). pp. 385-439.
6. Barut A., Guven I., Madenci E. (2001) Analysis of singular stress fields at junctions of multiple dissimilar materials under mechanical and thermal loading. Int. J. of Solid and Structures. 38(50-51). pp. 9077-9109.
7. Shemyakin E.I. (1996) O kraevykh zadachakh teorii uprugosti dlya oblastey s uglovymi tochkami (ploskaya deformatsiya) [0n the boundary value problems of the theory of elasticity for domains with angular points (plane deformation)]. Dokl. Akad. Nauk. 347(3). pp. 342-345.
8. Zadoyan M.A. (2003) Prochnost' soedineniya sostavnykh plit [Bonding strength of the composite plates]. Mechanics of Solids. 1. pp. 111-122.
9. Paggi M., Carpintery A. (2008) On the stress singularities at multimaterial interfaces and related analogies with fluid dynamics and diffusion. Mech. Rev. 61(020801). DOI: 10.1115/1.2885134.
10. Xu L.R., Kuai H., Sengupta S. (2004) Dissimilar material joints with and without free-edge stress singularities: Part I. A Biologically Inspired Design. Experimental mechanics. 44(6). pp. 608-615.
11. Xu L.R., Kuai H., Sengupta S. (2004) Dissimilar material joints with and without free-edge stress singularities: Part II. An integrated numerical analysis. Experimental mechanics. 44(6). pp. 616-621.
12. Kovalenko M.D., Men'shova I.V., Shulyakovskaya T.D. (2013) Razlozheniya po funktsiyam Fadlya-Papkovicha. Primery resheniy v polupolose [Expansions in Fadle-Papkovich. Examples of solutions in a half-strip]. Izv. Akad. NaukМТТ. 5. pp.136-158.
13. Luangarpa C., Koguchi H. (2016) Analysis of singular stresses at a vertex and along a singular line in three-dimensional bonded joints using a conservative integral. Europ. J. of Mechanics - A/Solids. 60. pp. 208-216. DOI: 10.1016/j.euromechsol.2016.08.002.
14. Kondrat'ev V.A. (1967) Kraevye zadachi dlya ellipticheskikh uravneniy v oblastyakh s konicheskimi ili uglovymi tochkami [Boundary value problems for elliptic equations in the domains with conical or angular points]. Proceedings of Moscow Mathematical Society. 16. pp. 209-292.
15. Maz'ya V.G., Plamenevsky B.A. (1976) O koeffitsientakh v asimptotike resheniy ellipticheskikh kraevykh zadach vblizi rebra [On coefficients in the asymptotics of solutions to the elliptic boundary value problems in the vicinity of edge]. Dokl. Akad. Nauk SSSR. 229(1). pp. 33-36.
16. Parton V.Z., Perlin P.I. (1984) Mathematical Methods of the Theory of Elasticity. Moscow: Mir.
17. Kozlov V.A., Mazya V.G., Schwab C. (1992) On singularities of solutions of the displacement problem of linear elasticity near the vertex of a cone. Arch. Ration. Mech. Anal. 119. pp. 197-227.
18. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2014) Nestandartnye zadachi dlya odnorodnykh elementov konstruktsiy s osobennostyami v vide klin'ev v usloviyakh ploskoy zadachi [Nonstandard problems of homogeneous structural elements with wedge shape features in the plane case]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika -Tomsk State University Journal of Mathematics and Mechanics. 27(1). pp. 95-109.
19. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2015) Iteratsionnyy konechno-elementnyy al-goritm issledovaniya napryazhennogo sostoyaniya elementov konstruktsiy s osobymi tochkami i ego realizatsiya [Iterative finite element algorithm and its implementation for stress state of structural elements with singular points]. Vestnik Permskogo natsional'nogo issledovatel'skogo politekhnicheskogo universiteta. Mekhanika - PNRPUMechanics Bulletin. 4. pp. 171-187. DOI: 10.15593/perm.mech/2015.4.11.
20. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2017) Stress state at the vertex of a composite wedge, one side of which slides without friction along a rigid surface. Latin American Journal of Solids and Structures. 14(11). pp. 2067-2088. DOI: 10.1590/1679-78253826.
21. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2016) Napryazhenno-deformirovannoe sostoy-anie v okrestnosti vershiny sostavnogo klina s zhestko zashchemlennymi obrazuyushchimi [Stress-strain state near the wedge top with rigidly fastened sides]. Vestnik Permskogo nat-sional'nogo issledovatel'skogo politekhnicheskogo universiteta. Mekhanika - PNRPU Mechanics Bulletin. 3. pp. 131-147. DOI: 10.15593/perm.mech/2016.3.09.
22. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2015) Nonstandard problems for structural elements with spatial composite ribs. Mechanics of Composite Materials. 51( 4). pp. 489-504. DOI: 10.1007/s11029-015-9520-9.
23. Pestrenin V.M., Pestrenina I.V., Landik L.V. (2018) Ogranicheniya na parametry napryaz-hennogo sostoyaniya v vershine krugovogo konusa [Restrictions on stress components in the top of round cone]. Vestnik Tomskogo Gosudarstvennogo Universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 52. pp. 89-101. DOI: 10.17223/19988621/56/9.
24. Pestrenin V.M., Pestrenina I.V. (2017) Constraints on stress components at the internal singular point of an elastic compound structure. Mechanics of Composite Materials. 53(1). pp. 107-116. DOI: 10.1007 / s11029-017-9644-1.
25. Kovalev B.D. (1983) Formirovanie Eylerovoy gidrodinamiki [Formation of the Euler hydrodynamics]. Issledovaniya po istorii mekhaniki. Sbornik statey [Research of the history of mechanics. Collected papers]. Мoscow: Nauka. pp. 146 -167.
Received: May 5, 2018
