Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Деформирование и разрушение композитов

..pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
6.86 Mб
Скачать

Рис. 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

см,

= 0,625 см,

/ = 0,75 см. На каж­

дом шаге задается смещение wn = 0,005 см,

при этом приращения

деформаций

не превышают

1,5%. На рис. 2 а, б представлены

линии

равных осевых напряжений

на 50-м шаге, при этом мак­

симальные

деформации

достигают

100%.

На

рис. показано

распределение среднего напряжения о* в

цилиндре при а = 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 с.