Осреднение нелинейных задач в механике композитов Текст научной статьи по специальности «Физика»

УДК 511

ОСРЕДНЕНИЕ НЕЛИНЕЙНЫХ ЗАДАЧ В МЕХАНИКЕ КОМПОЗИТОВ

С. В. Шешенин1, М. И. Савенкова2

Предлагается развитие метода осреднения для решения физически нелинейных задач о равновесии слоистых пластин или пластин из функционально-градиентных материалов. Согласно методу осреднения, решением задачи является суперпозиция решения глобальной задачи во всей области и решения локальной задачи для представительной области, например ячейки периодичности. Для нелинейной задачи суперпозиция неверна, что осложняет применение метода в случае нелинейности. Выходом может служить процедура объединения метода осреднения и метода линеаризации при решении краевой (или вариационной) задачи. Определяющие соотношения в механике деформируемого твердого тела можно рассматривать как уравнения относительно скоростей или дифференциалов напряжений и деформаций по времени или параметру нагружения. В том случае, если они линейны относительно скоростей, можно применить процедуру метода осреднения. В статье такой подход демонстрируется на примере симметричной слоистой пластины, изгибающейся под воздействием равномерно распределенной нагрузки, изменяющейся во времени.

Ключевые слова: метод осреднения, эффективные модули, пластичность, теория деформаций, изгиб, композит, слоистая пластина, линеаризация, метод Эйлера, нелинейность.

In this paper we consider the method of homogenization that helps to solve the nonlinear equilibrium problems for layered plates and FGM plates. The homogenization method uses the superposition idea to combine the solution on the global domain with the fluctuations obtained as the solutions to local problems on the representative volume element domain. Superposition is not correct in the case of any nonlinearity. Nevertheless, it is quite possible to use the linearization of a nonlinear problem by differentiating the variational equation with respect to time or with respect to a loading parameter and by considering the constitutive equations formulated in terms of strain and stress rates. These are equations of hypo-elastic materials. As an example, a symmetric laminated plate subjected to bending is studied when a load uniformly depends on time.

Key words: homogenization method, effective moduli, plasticity, deformation theory, bending, composite, laminated plate, linearization, Euler method, nonlinearity.

Метод осреднения [1, 2] является математически строгим подходом для решения задач деформирования элементов конструкций из композитов, например слоистых пластин или пластин из функционально-градиентных материалов [3]. Прежде всего это касается линейных задач. Комбинированием решений локальных задач с решением глобальной задачи для материала с осредненными свойствами метод осреднения позволяет вычислить не только эффективные свойства, но и все характеристики напряженно-деформированного состояния. Наибольшее развитие он получил для периодически неоднородных сред, так как в этом случае представительной областью является любая периодическая ячейка структуры. Эффективность решения линейной задачи основывается на том, что достаточно решать локальные задачи только на одной ячейке.

В случае нелинейности средние по представительным областям свойства зависят от медленных координат [4]. Поэтому для определения эффективных свойств как функций медленных координат требуется решать локальные задачи для каждой представительной области. Во избежание этого используется параллельное решение локальных задач на многопроцессорных системах. Ниже описывается реализация процедуры комбинирования метода осреднения и метода линеаризации при решении задачи о равновесии слоистой пластины. В случае, когда напряженно-деформированное состояние в глобальной задаче является однородным, локальные задачи не зависят от медленных координат, и, как в линейном случае, достаточно решить их для одной представительной области. При этом решение локальных задач

1 Шешенин Сергей Владимирович — доктор физ.-мат. наук, проф. каф. механики композитов мех.-мат. ф-та МГУ, e-mail: sergey. sheshenin@mail. ru.

2 Савенкова Маргарита Ивановна — асп. каф. механики композитов мех.-мат. ф-та МГУ, e-mail: madgista@gmail.com.

можно рассматривать как численный эксперимент, заменяющий реальный эксперимент по определению эффективных материальных функций осредненного нелинейного материала. Такой подход применялся, например, в [5, 6].

Используем определяющее соотношение, линейное относительно дифференциалов по параметру на-гружения напряжений и деформаций в рамках геометрически линейной теории деформирования:

= Сцк1(§)ёк1-

Здесь и далее точка обозначает дифференцирование по параметру нагружения. Уравнения равновесия с граничными условиями также формулируются в терминах скоростей перемещений

( {С^ы(и)ик,г),э = 0,

\tii\j, = и0, и)ик, щ 2 =

где £! — кинетическая поверхность, ^2 — статическая поверхность, и0 — граничные скорости перемещений, £0 — поверхностные силы, щ — компоненты вектора нормали к ^

В качестве примера используем хорошо известную деформационную теорию изотропной пластичности [7], для которой краевая задача записывается в виде (1), причем при активном нагружении

Сг]ш{е) = А(е)5ц5к1 + ¡л(е) + 5й5^к)

еи (2)

ц(е) = С(1 - и(е)), Х(е) = К - 2/Ъц(е), ец = 1/Ъ05ц, еи = л/ёТе,

где в = екк, и(е) — функция пластичности, К — модуль объемного упругого растяжения, С — касательный модуль сдвига в области упругости материала, Л — постоянная Ламе. При разгрузке

С^ы = Л5^5кь + С (5ц~6]1 + 5ц5ук), Л = К - 2/3 С. (3)

Поскольку при пошаговом решении задачи для композиционного материала в момент гт перемещение и известно, то определяющее уравнение можно переписать в виде

ёц = С%ы (£, х) еы, СЦы (£, х) = (е (£, х, гт)),

где х — глобальные (медленные) координаты, заданные в материале в целом; £ = х/е — локальные (быстрые) координаты, заданные на ячейке периодичности.

Далее рассматривается применение метода осреднения для решения задачи (1)—(3) на примере задачи изгиба слоистой пластины толщиной Н с постоянным прямоугольным поперечным сечением под воздействием равномерно распределенной нагрузки р(Ь). Решение задачи (1) ищется, согласно методу

осреднения, в первом приближении относительно скорости и в момент гт в виде и = ие + иЬ [6], где

• е ■ Ь

ие — скорость перемещения, возникающего вследствие растяжения-сжатия в плоскости пластины; иЬ —

скорость перемещения вследствие изгиба:

ие (С, х, гт) = Щ (х, гт) + бЩ^ (С, х, гт) Ур(х, ^ , г = 1, 2, 3;

и} (С, х,гт) = -6^,1 (х,гт)+ е2N'ЬpQ (С, х,гт) W,PQ (х,1т), I = 1, 2; (4)

иЬ3 (С, х, гт) = ^ (х, гт) + е2^^ (С, х, ^ W,pQ (х, ^ ,

где Р = 1, 2; ( = 1, 2; х = (х1,х2); С = СЬ = ХЬ/е; VI(х,гт) (I = 1, 2, УЬ = 0) — компоненты гладкой составляющей скорости в плоскости пластины; й(х, гт) — скорость прогиба; е — малый параметр, являющийся отношением толщины Н пластины к ее характерному размеру в плане; С — единственная быстрая координата.

Вообще говоря, метод осреднения позволяет удовлетворить уравнениям равновесия в первом приближении, а граничным условиям — в нулевом. Поэтому граничные условия приходится формулировать относительно средних скоростей перемещений VI и -ш. При этом возможно рассматривать граничные условия разных типов, например жесткое закрепление (из = 0, ^ = 0) или шарнирное закрепление (ь) = 0,

ди2 '

Будем считать, что проведено обезразмеривание всех соотношений, тогда характерный размер пластины в ее плоскости равен единице, а е = Н. Поэтому будем использовать символ Н для обозначения толщины пластины и малого параметра.

Локальные функции х,Ьт), NlbpQ(С, х,Ьт) находятся из решения локальных задач, получа-

емых путем подстановки соотношений (4) в (1). Для слоистой пластины, как и в линейном случае [8], локальные задачи в момент Ьт имеют вид

P3PQ (С, х,Ьт), = 0, QiзpQ (С, х,гт). = 0, С е (-1/2,1/2);

(5)

PiзpQ (С, х,гт) = о, QiзpQ (С, х,гт ) = 0, С = ±1/2,

где использованы обозначения:

PijPQ (С, х, Ьт) = PijPQ = CijPQ + Cijk3NkPQ>Z,

QijPQ (С, х, Ьт) = QíjPQ = CijPQ С + ■

Существенное отличие от линейного случая состоит в том, что эти задачи зависят от х как от параметра. Для слоистой пластины или пластины из функционально-градиентных материалов локальные задачи (5) решаются аналитически и имеют нетривиальные решения для функций N£pQ(С, х,Ьт)

и N3bpQ(С, х,Ьт). При численной реализации локальные задачи решаются для каждого узла двумерной сетки в плоскости пластины.

В итоге эффективные жесткости [7, 8] как функции медленных координат определяются по формулам

1/2 1/2

лт^ (х) = н | рт^ (с, х) йС, вт^ (х) = Н2 | ^т^ (С, х) йс,

-1/2 -1/2

1/2 1/2

(6)

в^ (х) = Н2 ! СРГ^ (С, х) йС, БТ^ (х) = Н3 I ^Т^ (С, х) йС. -1/2 -1/2

Тогда осредненная система уравнений равновесия примет вид

(ЛТШЬ (х) Vт (х),ь ),з - (ВЪКЬ (х) (х),кь ),з = 0, {В^Гкь (х) Vт (х),ь) ,и + (ВТшь (х) ЪТ (х),кь ),и = Г-

(7)

Процедура решения в приращениях типа явного метода Эйлера состоит в том, что на шаге т решение ит = и(С, х,Ьт) предполагается известным. Тогда с помощью формул (2), (3) определяются модули жесткости CТpQ, затем решаются локальные задачи (5) и вычисляются эффективные жесткости (6), после чего решаются глобальные уравнения (7) с соответствующими граничными условиями и получаются перемещения на шаге т + 1: = + V'АЬ, ът+1 = Ът + ЪитАЬ, и затем алгоритм повторяется для

следующего шага нагружения.

Неявный метод Эйлера заключается в том, что вычисленные с помощью явного метода Эйлера перемещения считаются только нулевым приближением: угГ+1'° = )Т+1, ът+1'0 = Ът+1, поэтому производится уточнение модулей. По формулам (2), (3) вычисляются = CijPQ(um'+1'0), далее решаются локальные задачи (5) и с помощью (6) вычисляются эффективные жесткости АТ/+Кь, ВТш 'ь ,

БТт+1'0 ГЧТ+1'0 , т+1'1 т , -т+1'1 д ,

икь , Бик'ь , после чего находятся уточненные перемещения на шаге т + 1: = )Т + АЬ,

Ът+1'1 = ът + Ъ]т+1'1АЬ. Таким образом, для итерации с номером в имеем у'^^^1 = у' + у'^'^1 АЬ, Ът+1'3+1 = ът + Ъ]т+1'3+1 АЬ. После некоторого числа итераций функции )Т+1'8+1, -ът+1'3+1 принимаются за окончательные значения перемещений Ът+1 на шаге т + 1.

Изложенные явный и неявный методы были реализованы в виде программ; осуществлено решение задачи цилиндрического изгиба шарнирно опертой, слоистой симметричной пластины бесконечной ширины —то < Х2 < с постоянным прямоугольным поперечным сечением 0 ^ Х1 ^ I, -0,5 ^ С ^ 0,5 под действием нагрузки р(Ь), равномерно распределенной по верхней границе пластины С = 0,5 и линейно изменяющейся во времени. Отношение Н/1 было выбрано равным 0,1. Пластина пятислойная,

каждый из слоев имеет высоту 0,2 и является элементом одного из трех множеств по координате £: А\ = (0,5; -0,3) и (0,3; 0,5), А2 = (-0,3; -0,1) и (0,1; 0,3), Аз = (-0,1; 0,1). Численное решение задачи искалось для упругопластического материала с билинейной диаграммой а3(е3) (рис. 1). Параметры материалов приведены в таблице, где а8 — безразмерный предел текучести, О — модуль сдвига упругой зоны деформации, ОТАМ — касательный модуль сдвига в пластической зоне, V — коэффициент Пуассона.

Вследствие того что приложенная к пластине сила не зависит от координаты х2, эффективные жесткости и прогиб также от нее не зависят, поэтому далее будем обозначать х1 = х. Поскольку рассматриваемый изгиб пластины цилиндрический, в формулах (6), (7) имеем Р = Q = 1.

Кроме того, рассматриваемая пластина симметрична, при ее изгибе функции Q^JPQ являются антисимметричными по £, поэтому соотношения (6) для изгибных жесткостей принимают вид

1/2

^п (х) = ь3/ -№л! (е,х) -1/2

Уравнения равновесия и граничные условия следующие:

Рис. 1. Билинейная диаграмма упруго-пластического материала

Параметры материала Слои, составляющие пластину

А1 А2 Аз

Т., 1,2 0,7 1

С 0,5 0,32 0,56

ОГАМ 0,05 0,032 0,056

V 0,32 0,34 0,2

(^п (х) ад(х)) ,и=рт, й,т(0) = адт(1) = 0, й,т (0) = й,т (о = 0, Тогда решения локальных задач в нулевом приближении дают

« 1/2

й0(х) = 0.

N311 (х,е)= / у

Сзззз (X, у)

^у + СОП81,

-1/2 -1/2 и для каждого значения х деформации и напряжения имеют вид

.(0)т

£ С3311 1£ _ 0

4 Сзззз (■*,£) 4

4Г = -и^ ,

'I3

= 0,

-й(0)т = и Кт • т ;

= Qm7PQй,mQ ,

40)т = 0.

На рис. 2 представлена вычисленная зависимость напряжения стц от вертикальной координаты £ £ (-1/2, 1/2) в поперечном сечении х = 1/2. По горизонтальной оси отложены значения координаты £ с шагом 0,1, по вертикальной оси — значения напряжений в пластине. Для построения графиков было рассмотрено пять значений параметра нагруже-ния, которым соответствуют кривые 1-5. Ввиду симметричности пластины график напряжений симметричен относительно начала координат, а наличие разрывов — прямое следствие неоднородности пластины. Каждый непрерывный отрезок графика представляет собой зависимость напряжения стц от вертикальной координаты £ в соответствующем Рис. 2 График зависимости напряжения ап от вертикаль-слое. Зона пластичности возникает между тре- ной координаты £ для центрального поперечного сечения тьим и четвертым значением параметра на- пластины х = 1/2

гружения, что хорошо заметно на рисунке: зависимость а(£), соответствующая четвертому значению, уже не является линейной на каждом отрезке. В момент возникновения пластичности отношение силы в рассматриваемый момент к общей силе, изгибающей пластину, составляет 0,637.

Таким образом, в работе показано, что метод осреднения может быть обобщен на нелинейные задачи с использованием метода линеаризации (или прослеживания по параметру нагружения). Основная возникающая при этом сложность — решение большого количества локальных задач — в настоящее время

может быть сведена к минимуму благодаря возможности проведения вычислений на многопроцессорных компьютерных системах. Такой подход позволяет значительно сократить расчетное время вследствие высокой степени параллелизма при решении локальных задач. Кроме того, ввиду того что все рассмотренные в статье соотношения записывались в терминах скоростей напряжений и деформаций, возможно использование не только теории малых деформаций, выбранной в качестве примера, но и других теорий, например теории течения, не вносящих никаких изменений в предлагаемый алгоритм.

СПИСОК ЛИТЕРАТУРЫ

1. Победря Б.Е. Механика композиционных материалов. М.: Изд-во МГУ, 1984.

2. Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. М.: Наука, 1984.

3. Hui-Shen Shen. Functionally graded materials: nonlinear analysis of plates and shells. Boca Raton: CRC Press, 2009.

4. Шешенин С.В., Fu M., Ивлева Е.А. Об осреднении периодических в плане пластин ^ Тр. Междунар. конф. "Теория и практика расчета зданий, сооружений и элементов конструкций. Аналитические и численные методы". М.: Изд-во МГСУ, 2008. 148-158.

5. Муравлева Л.В., Шешенин С.В. Эффективные свойства железобетонных плит при упругопластических деформациях У У Вестн. Моск. ун-та. Матем. Механ. 2004. № 3. 62-65.

6. Муравлева Л.В. Эффективные свойства ортотропных композитов при упругопластических деформациях УУ Упругость и неупругость: Мат-лы Междунар. науч. симп. по проблемам механики деформируемых тел, посвященного 95-летию со дня рождения А. А. Ильюшина. M.: URSS, 2006. 371-378.

7. Ильюшин А.А. Пластичность. Ч. I. М.: ОГИЗ, 1948.

8. Шешенин С.В. Применение метода осреднения к пластинам, периодическим в плане У У Вестн. Моск. ун-та. Матем. Механ. 2006. № 1. 47-51.

Поступила в редакцию 14.11.2011

УДК 532.528+539.374

СХЛОПЫВАНИЕ ОДИНОЧНОГО ПУЗЫРЬКА В ЖИДКОСТИ ШВЕДОВА-БИНГАМА ПОД ДЕЙСТВИЕМ ПУЛЬСИРУЮЩЕГО ВНЕШНЕГО ДАВЛЕНИЯ

А. А. Шевчук1

Проведено численно-аналитическое исследование поведения сферического пузырька в идеальной жидкости, вязкой среде и несжимаемой вязкопластической среде с пределом текучести при действующем на достаточном удалении от поверхности пузырька давлении, изменяющемся со временем по периодическому закону. Исследованы и классифицированы различные режимы схлопывания. Найдены критические значения характеризующих поведение системы ключевых параметров, одним из которых является безразмерная частота внешнего давления.

Ключевые слова: сферическая полость, пузырек, динамика, задача Рэлея, задача За-бабахина, сжатие, вязкопластическая среда, предел текучести.

Numerical and analytical study of the behavior of a spherical bubble in an ideal fluid, a viscous medium and an incompressible viscoplastic medium with a yield strength under time-varying periodic pressure at a sufficient distance from the surface of the bubble is carried out. Different modes of collapse are examined and classified. The critical values of the key parameters that characterize the behavior of the system (one of them is the dimensionless frequency of external pressure) are found.

Key words: spherical cavity, bubble, dynamics, Rayleigh problem, Zababakhin problem, compression, viscoplastic medium, yield limit.

1. Среда без предела текучести. Пусть все пространство заполнено идеальной несжимаемой жидкостью с плотностью р. На бесконечности действует пульсирующее давление p= pa cos Uot [1, 2].

1 Шевчук Александр Александрович — асп. каф. механики композитов мех.-мат. ф-та МГУ, e-mail: rockie@yandex.ru.