Деформирование и разрушение композитов
..pdfРис. 1. Расчетная схема.
Варьирование в (1) осуществляется по Н и по составляющим вектора U. Причем / 3, входящие в комплексы
# /(2 /3), по U не варьируются. Вариа ционные уравнения для (1) записы ваются в виде
х [(/з - |
1) - |
аА Д -1 6ffdV, = О, |
|
|
|
|
|
|
|
|
||||
L |
|
|
2/3 |
|
|
|
|
|
|
|
|
|
|
|
[ < [ у V h |
Q + (* 1 + * S) [ ( / , - |
1) - |
аА •£-] /.O '-*]• • 8а - |
|||||||||||
|
|
|
— p°7f8U > dK„ — j |
PcbUdS0=0. |
|
(2) |
||||||||
Из последних выражений легко показать, что уравнением |
||||||||||||||
Эйлера для |
(1) будет уравнение |
равновесия, а естественными гра |
||||||||||||
ничными условиями — граничные |
условия |
в напряжениях. Отме |
||||||||||||
тим, что редукция функционала |
(1) |
|
на |
линейные задачи непос |
||||||||||
редственно дает принцип Геррманна2. |
|
|
за редким исклю |
|||||||||||
Реализация функционала |
(1) предполагает, |
|||||||||||||
чением, |
его |
численное |
решение. |
В |
|
качестве метода численного |
||||||||
решения используется метод конечного элемента |
(МКЭ). |
В соот |
||||||||||||
ветствии с процедурой |
МКЭ составляющие вектора перемещения |
|||||||||||||
и функция |
среднего давления |
Н в элементе |
аппроксимируются |
|||||||||||
функциями |
формы |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Up= y KUNp. |
|
|
|
|
|
Я = < |
|
|
(3) |
|||
Здесь N, |
М — номера узлов элемента; IJK . VMK — ко- и контр А — |
|||||||||||||
вариантные |
составляющие |
вектора |
перемещений; |
выбраны |
||||||||||
квадратичными функциями координат, |
— линейными. |
значения, |
||||||||||||
Учитывая, что в (3) |
варьируются |
только |
узловые |
|||||||||||
и в силу |
произвольности |
Шр, |
|
ЬН уравнения для элемента за |
||||||||||
пишутся |
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
f [2о</ (A?iN+B?iKMNUMK)- |
р » ^ ^ ] dV, - |
j |
P > MdS»=0, |
|||||||||||
Vo |
|
|
|
|
|
|
|
|
|
|
|
Л |
|
|
2 Г е р р м а н н |
Л. Р. Вариационный принцип для |
уравнении |
упругости |
|||||||||||
несжимаемых |
и почти несжимаемых |
материалов.— Ракетная техника и кос |
||||||||||||
монавтика, 1965, № |
10, с. 139—144. |
|
|
|
|
|
|
|
|
|
111
где
Ayw= g f 'i V . ; - 'V f /.
в,ут =§№мЛ . / - УмЛГ/К+V » ,.r« - 'W ^ rfs.
Суммирование выражений (4) по всем элементам приводит к раз решающей системе нелинейных уравнений.
Для линеаризации системы (4) перемещения, функция сред него давления и нагрузки записываются в приращениях;
и м к ^ , ) м к + ги к 1 Я " = Я М+Т1«
PN = PN+S N,
причем величины £м/с и rjM считаются малыми и квадратами их можно пренебречь.
112
Разложение величин, входящих в (4) и зависящих от UMK и Нм, в ряд по \ мк и r|M, и сохранение только линейных членов приводит к системе линейных уравнений относительно прираще ний перемещений и функции riM.
Для решения задачи о растяжении цилиндра из слабосжимаемого и несжимаемого материала с осевыми шаровыми вклю чениями принята следующая схема: жесткие шары радиуса Rlt помещенные в цилиндр радиуса R2, получают осевое смещение w\ первоначальное расстояние между шарами /; материал цилиндра может быть несжимаемым или слабосжимаемым (см. рис. 1).
Согласно изложенному, |
задача решается в приращениях: задан |
||||||||
ное жесткое смещение шаров |
разбивается |
на п шагов, причем п |
|||||||
выбирается |
таким |
образом, |
чтобы |
на |
каждом шаге квадратами |
||||
деформаций, |
вызываемых растяжением |
шаров, можно было пре |
|||||||
небречь. |
|
|
|
|
|
результаты для цилиндра |
|||
На рис. 2 представлены некоторые |
|||||||||
из несжимаемого |
(а = 0) |
и |
слабосжимаемого |
(а = 0,078) матери |
|||||
алов |
типа |
неогука {Кл = 0,5 |
кг/см2, |
/С2 = 0) |
с геометрическими |
||||
характеристиками |
R2 = 1 |
см, |
R±= 0,625 см, |
/ = 0,75 см. На каж |
|||||
дом шаге задается смещение wn = 0,005 см, |
при этом приращения |
||||||||
деформаций |
не превышают |
1,5%. На рис. 2 а, б представлены |
|||||||
линии |
равных осевых напряжений |
на 50-м шаге, при этом мак |
|||||||
симальные |
деформации |
достигают |
100%. |
На |
рис. 2в показано |
||||
распределение среднего напряжения о* в |
цилиндре при а = 0 и |
||||||||
осевом растяжении |
ш = 0,5 см. |
|
|
|
|
АКАДЕМИЯ НАУК СССР |
УРАЛЬСКИЙ НАУЧНЫЙ ЦЕНТР |
|
ДЕФОРМИРОВАНИЕ И РАЗРУШЕНИЕ КОМПОЗИТОВ |
1985 |
Е. А. КОЛЧАНОВА
ПРИМЕНЕНИЕ РЕШЕНИЯ ОБРАТНОЙ КРАЕВОЙ ЗАДАЧИ ТЕРМОУПРУГОСТИ К СИНТЕЗУ
ПОПЕРЕЧНЫХ СЕЧЕНИЙ ЦИЛИНДРИЧЕСКИХ ТЕЛ
В некоторых отраслях промышленности нашли применение составные цилиндрические конструкции с поперечным сечением сложной формы. При проектировании их поперечных сечений (в виде двухсвязных областей, ограниченных снаружи окруж ностью, а изнутри — кривой сложной формы) оптимальный с точки зрения прочности вариант конструктивного оформления выбирается путем применения так называемого вариантного оптимального проектирования. В случае использования такого подхода рассчитывается напряженно-деформированное состоя ние ряда конструкций заданной формы и размеров с тем, чтобы из класса допустимых найти вариант, удовлетворяющий, наряду с другими, требованиям прочности. При этом оценка НДС по перечных сечений при действии на конструкцию силовой и тем пературной нагрузок осуществляется на основе решения прямой краевой задачи теории термоупругости в плоской постановке.
В данной работе предлагается более эффективный метод ре шения указанной задачи, позволяющий без расчета НДС опре делить множество вариантов конструктивного оформления по перечных сечений, удовлетворяющих требованиям прочности, а также заданной площади, и выбрать из них наиболее оптималь ный, с точки зрения других требований.
Возможности разработанного метода проиллюстрированы на примере решения задачи для цилиндра, скрепленного с оболоч кой, который под действием постоянного температурного пере пада находится в условиях плоской деформации. Относительно оболочки сделано допущение, что она абсолютно жесткая.
1. Предлагаемый метод проектирования конструкций, удов летворяющих указанным требованиям, основан на решении обратной краевой задачи плоской теории упругости (термо упругости) [1, 2] для поперечных сечений конструкций изучае мого класса. Решение осуществляли в' естественной криволиней ной системе координат р, ft.
Принятое допущение относительно оболочки позволяет ре-
114
шение обратной краевой задачи термоупругости свести к реше нию соответствующей задачи теории упругости при нулевых значениях нормальной и касательной компонент тензора напря жений на внутренней границе сечения:
®»=0, Т„ = 0 |
(1) |
и заданном векторе перемещений — на внешней:
|
|
un = AaATb, vn = 0. |
(2) |
|
Здесь |
<т„, тл, ип, |
vn— нормальная и касательная компоненты |
||
тензора напряжений и вектора перемещений; п — нормаль к гра |
||||
нице; Аа = ах — а0— разность |
коэффициентов |
линейного темпера |
||
турного |
расширения |
наполнителя и материала оболочки; АТ = |
||
= Тг — Т0— разность |
между |
температурой |
полимеризации Тх |
и окружающей среды Г0; Ь— радиус окружности, являющейся внешней границей сечения.
Кроме краевых условий, исходными при решении данной задачи были: величина предельного отрывного напряжения апЬ между слоями наполнитель—оболочка и предельного напряже ния аь (деформации гь) на растяжение для наполнителя, а также
формула |
связи |
его |
с |
компонентами тензора напряжений; пред |
||||||||||||||
полагаемая |
форма |
поперечного |
сечения |
(число |
концентраторов |
|||||||||||||
и их соотношение |
размеров), соответствующая |
закону |
распреде |
|||||||||||||||
ления напряжения |
|
|
|
(деформации ед=(1 — vz)a^/E) по внут |
||||||||||||||
ренней |
границе |
сечения. |
Если |
в |
область |
q — симметричная, то |
||||||||||||
достаточно |
задать |
этот |
закон |
одном |
симметричном элементе |
|||||||||||||
0 <д<сг/<7. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Максимум этого напряжения рассчитывали на основании фор |
||||||||||||||||||
мулы |
|
для |
аь (деформации |
еь), |
с учетом (1) и соотношения сг2 = |
|||||||||||||
= va^,. |
При |
этом |
необходимо иметь в виду, что количество кон |
|||||||||||||||
центратов |
|
в |
поперечном сечении |
равно числу |
максимумов функ |
|||||||||||||
ции /(О). Характер зависимости |
последней от угла О |
обусловли |
||||||||||||||||
вает |
размеры и форму |
концентраторов, а также их относительное |
||||||||||||||||
местоположение. |
решения |
задачи |
находится ряд вариантов кон |
|||||||||||||||
В результате |
||||||||||||||||||
структивного оформления поперечного сечения |
изделия, |
удовлет |
||||||||||||||||
воряющего |
требованию |
прочности. |
Затем из элементарного рас |
|||||||||||||||
чета |
определяется |
конфигурация |
поперечного |
сечения |
заданной |
|||||||||||||
площади. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2. Разрешающая система уравнений в терминах теории функ |
||||||||||||||||||
ций |
комплексного |
переменного с применением конформных ото |
||||||||||||||||
бражений |
[3] |
включает: |
|
|
|
|
уравнений |
в |
криволинейной |
|||||||||
1) |
|
систему |
функциональных |
|||||||||||||||
системе координат |
р, |
О: |
|
|
|
|
|
|
|
|
|
|||||||
[х*Ф (ра)ц-Ф (per)] со' (per) — а~2 [со (ра) Ф' (ра)+Чг (ра) ю' (рег)]= |
||||||||||||||||||
|
|
|
|
) 0 |
|
при |
х* = 1, |
p = pj; |
|
|
|
|
*(3) |
|||||
|
|
|
|
J —2рАаА7’а)' (рог) при |
х* = —х, |
р = 1. |
|
(4) |
115
(здесь |
ф (£), |
Ч1, (£), |
о)(£) — комплексные |
потенциалы |
и функция |
|||||
конформного |
отображения, |
голоморфные |
в |
области |
P i< |C |< l; |
|||||
;p= £,/2(l + v), |
Е — модуль |
Юнга, |
v — коэффициент |
Пуассона, |
||||||
х = 3 — 4v; а = е^\ ^= per; |
последнее |
соотношение есть деформа |
||||||||
ционный вариант геометрических граничных условий [3]); |
= |
|||||||||
2) |
уравнение |
связи между |
заданной |
компонентой |
||||||
тензора |
напряжений |
на |
внутренней |
границе области |
(р = рх) |
|||||
с одним из комплексных потенциалов: |
|
|
|
|
|
|||||
|
|
2[Ф(р1(г)+ф(р1ог)]=/(в); |
|
|
(5) |
3)соотношение между этим комплексным потенциалом и пре
дельным отрывным напряжением опЬ = в рь на |
линии контакта |
(р=1), имеющее место при 0 = 0: |
|
х-±1[Ф(1)+Ф.(1)]-рДаДГ=сГрб; |
(6) |
4) уравнение |
|
х+1 [Ф (СУ) — Ф (а)] = 0, |
(7) |
которое является следствием уравнения (4) и представляет собой условие равенства нулю касательного напряжения на линии кон
такта |
(р= 1). |
|
поставленной задачи осуществляли методом рядов, |
||||||
3. |
Решение |
||||||||
согласно которому комплексные потенциалы и конформно-ото- |
|||||||||
бражающая функция |
аппроксимируются рядами |
Лорана |
в |
||||||
форме |
[4]: |
|
|
|
|
|
|
|
|
|
Ф(С)= |
2 |
0*£‘+ |
2 |
в_»р‘Е-‘, |
( 8) |
|||
|
|
k —q n , |
к — q n , |
|
|
|
|
||
|
|
л=0 |
|
л= 1 |
|
Р.<|£1<1 |
|
||
|
|
|
|
|
|
|
|
||
|
ЧЧ£)= |
2 |
|
|
2 Ь_„р‘Е-‘- 2, |
(9) |
|||
|
|
к—q n , |
|
k= q n , |
|
|
|
||
|
|
п=0 |
|
|
п = |
1 |
|
|
|
|
ш(0= |
2 ад*+ 2 с_»р*Е-‘, Pl<it|<i. |
( |
||||||
|
|
k = q n + 1, |
|
k= qn— l, |
|
|
|||
|
|
л—0 |
|
|
л=1 |
|
|
|
|
Коэффициенты |
ряда (8) определяются по следующим |
формулам: |
|||||||
|
|
|
|
|
|
n/q |
|
|
|
|
|
|
|
a°= f |
\ |
П Щ М ' |
(11) |
||
|
|
|
|
|
n/q |
|
|
|
|
|
a_h = d - - l — |
\ /(О)cosMrfft, |
( 12) |
||||||
|
|
|
l+Pi* о |
|
|
|
|
||
|
|
|
|
|
Oft |
Р |
—ft» |
( 13) |
|
которые получены |
из |
условий |
(5), |
(7). |
|
|
116
Величина |
радиуса рх находится из уравнения (6), которое |
с учетом (13) |
принимает вид |
|
(14) |
При ограничении в сумме одним членом ряда, что является достаточным в том случае, когда величина a_2 ( J первого отбро
шенного члена мала по сравнению с величиной, стоящей в пра вой части, радиус рх = р* вычисляется по формуле
где |
Р* = е с, |
(15) |
У+^АаАГ |
|
|
1 |
|
|
a _ qq [ |
х+ 1 |
|
Проверку точности расчетов и дальнейшее уточнение значе ния (15) можно проводить наиболее простым способом. Для этого
значение р= р! |
(15) |
необходимо |
подставить |
в (6) и вычислить |
|||||
погрешность, |
с |
которой |
выполняется |
это |
условие. |
Если она |
|||
укладывается |
в |
пределы |
погрешности |
для |
а рь, то определение |
||||
рх можно |
считать |
законченным. |
Если |
расчетное |
напряжение |
||||
^ ( 0 = 0) |
по |
формуле (6) окажется меньше |
а рь, то значение р! |
можно увеличить на бесконечно малую величину Apx; если боль ше, то pi необходимо уменьшить и снова сделать проверку и т. д.
Решение этой задачи можно осуществлять также методом по следовательных приближений, используя в качестве нулевого приближения р —р* (15).
Следует иметь в виду, что при анализе исходных данных для решения поставленной задачи остался незафиксированным пара метр q. При выбранном множестве значений рх>р1 он может быть произвольным.
Таким образом, после того, как найдены входные параметры задачи, можно перейти к ее решению — определению функции конформного отображения из системы (3), (4).
В результате соответствующих алгебраических операций [3, 4] исключается комплексный потенциал Ф (С) (3). Полученная си стема линейных алгебраических уравнений относительно неиз вестных коэффициентов ряда (10) в общем случае является одно
родной. |
Ее |
можно привести к неоднородной в том случае, если |
||||
одна |
из границ у проектируемой конструкции поперечного сече |
|||||
ния |
представляет собой окружность. В таком случае один коэф |
|||||
фициент |
в (10) |
(в рассматриваемом случае Сх) можно [4, |
5] счи |
|||
тать известным и равным величине радиуса окружности. |
методом |
|||||
Решение |
|
неоднородной |
системы разыскивается |
|||
Гаусса. |
В |
результате находится искомая функция (10), |
с помо |
|||
щью |
которой |
описывается |
поперечное сечение изделия, |
удовлет |
воряющее требованиям прочности.
117
Рис. 1. Вид поперечного сечения составно го цилиндра.
4.
подходе в результате решения за дачи для оболочки с наполнителем определены форма и размеры кон центраторов, при которых прочность поперечного сечения обеспечена.
|
|
|
|
При |
этом не учитываются те участ |
|||||
вые выступы), |
которые |
ки |
поперечного |
сечения |
(межщеле |
|||||
находятся |
вдали от |
концентраторов |
и |
|||||||
не влияют |
на |
прочность (рис. |
1, |
заштрихованная |
часть). |
|
||||
Такие участки, |
как |
показывают практические |
расчеты НДС |
|||||||
конструкций |
изучаемого |
класса |
и |
специальные |
исследования, |
|||||
проведенные |
нами, |
имеются в |
тех конструкциях |
изделий, |
где |
жесткость поперечного сечения оболочки больше соответствую щей жесткости поперечного сечения цилиндрического тела.
В связи со сказанным выше при построении поперечных сечений заданной площади к конфигурации поперечного сечения,
полученной |
из |
решения задачи, следует пристраивать соответст |
||||
вующим |
образом |
участки необходимой площади, не влияющие |
||||
на |
прочность поперечного сечения. |
|||||
|
Для |
простоты |
расчетов площади искомых участков попереч |
|||
ное сечение |
целесообразно представить состоящим из двух фигур: |
|||||
кругового кольца |
pib< |£ |<6 и перфорированной части. Тогда |
|||||
площадь S* |
последней будет равна |
|||||
|
|
|
|
|
|
S* = S — S°, |
где |
через |
S |
и |
S 0 |
обозначены, соответственно, общая площадь |
|
поперечного |
сечения |
и площадь кольца, равная |
S° = nb2{1- р2) .
Таким образом, в результате решения задачи при одном законе распределения о'^= Д^) получается счетное множество
сечений, отличающихся величиной радиуса pt и числом осей сим-
Рис. 2. Конструкция попереч ного сечения, соответствующая закону распределения напря жения Стф с наличием одного
(а) и двух (б) максимумов в секторе О с д с л / д (а = 4).
118
Рис. 3. Характер распределе ния напряжений по внутренней границе сечения, соответствую щий одному (а) и двум (б) кон центраторам в симметричном
элементе.
0,Z 0,if О,В О,В |
0,2 0,if О, В 0,2 |
V-’/ (зг/$ \ град
метрик q. Все конструкции этого множества удовлетворяют тре бованиям прочности. Теперь из этого множества можно выбирать оптимальный вариант конструкции с точки зрения других тре бований.
Такой подход, как это может показаться на первый взгляд, нарушает традиционный путь проектирования конструкций. Одна ко это не так, поскольку удовлетворяющая условиям прочности конструкция (выбранная из данного семейства) может быть допроектирована уже традиционным путем. Следует также отме тить, что предлагаемый метод проектирования обладает большим преимуществом: для выбранной конструкции уже известна ото бражающая функция, которая позволяет сделать проверочный расчет НДС на основании решения [3—5] прямой задачи.
На основе предложенного и разработанного метода с помощью программы для ЭЦВМ был просчитан ряд вариантов, отличаю щихся законом распределения напряжения сг^ (функции /(ft))
по внутренней границе области.
Вкачестве примеров на рис. 2 приведены конфигурации об ластей, у которых максимальное напряжение (в некоторых точках внутреннего контура) одинаково (рис. 3), а распределение напря жения (T^ по контуру области различное.
Взаключение следует отметить, что построенные области
описываются функцией |
конформного |
отображения, |
содержащей |
|||||||||
12 членов ряда. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Л И Т Е Р А Т У Р А 1234 |
|
|
|
|
|
||
|
1. К о л ч а н о в а |
Е. |
А., |
О с и п а н о в |
А. А. Решение |
обратной |
краевой |
|||||
задачи плоской теории |
упругости.— В кн.: |
Всесоюзный |
симпозиум |
по |
тео |
|||||||
рии |
механической |
переработки полимерных |
материалов: |
[Тез. |
докл. |
Все |
||||||
союзна симпозиума]. Пермь, 1976, с. 70. |
А. А. Об одном методе решения об |
|||||||||||
|
2. К о л ч а н о в а |
Е. А., О с и п а н о в |
||||||||||
ратной краевой задачи плоской теории |
упругости.— В |
сб.: |
Вопросы |
тео |
||||||||
рии упругости и вязкоупругости. Свердловск: УНЦ АН СССР, |
1978, с. 73—76. |
|||||||||||
ции |
3. К о л ч а н о в а |
Е. А., |
Ш у м о в с к а я |
А. А. Определение |
концентра |
|||||||
напряжений |
в |
составной двухсвязной |
области.— Прикл. |
мех., |
1974, |
|||||||
т. 10, вып. 10, с. 120—126. |
Х а с а н о в а |
Г. |
3. О решении некоторых |
задач |
||||||||
|
4. К о л ч а н о в а |
Е. А., |
119
плоской теории упругости для сложных двухсвязных областей.— В сб.: Прикладные вопросы теории упругости и вязкоупругости. Пермь: Перм.
политехи, |
ин-т, 1971, с. |
11—17. |
|
5. К о л ч а н о в а |
Е. А., |
М и л я к о в а Л. В. Алгоритм построения кон |
|
формного |
отображения для |
сложной замкнутой двухсвязной области.— |
В кн.: Методы решения задач теории упругости и вязкоупругости. Сверд ловск: УНЦ АН СССР, 1974, с. 36—41.
6. М у с |
х е л и ш в и л и Н. И. |
Некоторые основные задачи математиче |
ской теории |
упругости. М.: Наука, |
1966. 707 с. |