CC BY
34
УДК 531
РАСЧЕТ ДЕФОРМИРОВАНИЯ ГЕОМЕТРИЧЕСКИ НЕЛИНЕЙНОЙ ПЛОСКОЙ КОНСТРУКЦИИ ИЗ ИДЕАЛЬНО ПЛАСТИЧНОГО МАТЕРИАЛА
А.О. Щербакова1
Предложенная расчетная модель позволяет, зная нагрузки, приложенные к плоской конструкции из идеально пластичного материала, найти соответствующие смещения ее точек и, наоборот, по заданным смещениям найти приложенные к конструкции силы. Модель включает в себя геометрические, статические и физические соотношения и основывается на методе конечных элементов, где конечные элементы представляют треугольные симплексы.
Ключевые слова: геометрическая нелинейность, идеально пластичный материал, чистый изгиб, деформирование.
Введение
В настоящее время в механике сплошной среды имеется значительное количество различных мер напряжений (тензор напряжений Коши [1-3], тензоры Пиолы-Кирхгофа первого и второго рода [2, 3]) и деформаций (тензоры Коши-Грина, Грина-Лагранжа, Альманси, Генки и т.д.) [15]. Такое разнообразие, конечно, дает определенную свободу в создании новых и усовершенствовании уже имеющихся моделей расчета конечных деформаций конструкций. Однако в большинстве случаев понимание физического смысла соотношений, используемых в этих моделях, затруднено, особенно, когда речь идет о конструкциях из неупругого материала, где кроме геометрической нелинейности имеет место физическая. Целью данной работы является разработка такой модели расчета конечных деформаций неупругой конструкции (в качестве примера выбран идеально пластичный материал), которая была бы основана на тензорах, наиболее понятных с точки зрения физического смысла их координат, а именно: на тензоре напряжений Коши о и тензоре инженерной деформации е, наследуемом из геометрически линейного подхода.
Основные соотношения расчетной модели
Предложенная модель основана на методе конечных элементов, где элементы представляют треугольные симплексы. В рамках каждого элемента поле смещений является однородным и линейным, что обуславливает постоянство полей напряжений и деформаций. Однородность деформированного состояния элемента означает наличие линейной связи между начальными векторами волокон l0 этого элемента и векторами его волокон l в деформированном состоянии:
l = F • lo. (1)
Оператор этой связи - тензор дисторсии F, координаты которого определяются векторами смещений узлов элемента. Выражение для дисторсии конструкции, состоящей из нескольких (например, k) конечных элементов, удобно записывать в матричной форме:
[F] = [L] [q] + [/1]. (2)
В этом выражении [F] - это столбец, содержащий координаты тензоров дисторсии элементов, перечисленные последовательно в виде блоков
[F] = [[F1 F2 F3 F4] 1 [F1 F2 F3 F4]2 ••• [F1 F2 F3 F4]k]T;
[/1] - столбец, содержащий блоки координат единичного тензора
[/1] = [[1 0 0 1] [1 0 0 1] ... [1 0 0 1]]T;
а [q] - столбец, содержащий смещения узлов конструкции
[q] = [q1 q2 ... qm]T,
где число m - степень свободы конструкции. Матрица [L] размерностью 4k*m определяет связь между смещениями [q] и дисторсией [F]; для каждой конструкции [L] постоянна и ее координаты не зависят от нагрузок или смещений.
1 Щербакова Алла Олеговна - кандидат технических наук, доцент, кафедра Прикладной механики, динамики и прочности машин, Южно-Уральский государственный университет.
E-mail: AllaScherbakova@list.ru
Механика
Тензор дисторсии представляет скалярное произведение тензора жесткого поворота Я и тензоров упругого ие и неупругого ир растяжения:
^ = Я ■ ие ■ ир, ие = I + ее, ир = I + ер, (3)
здесь тензоры £е и £р - это тензоры упругой и неупругой деформации соответственно, а I - единичный тензор. В матричной форме тензоры ие, ир, £е и ер записываются в виде столбцов (по аналогии с дисторсией) с учетом симметрии. Например [£е] - это столбец размерностью 3к, который содержит компоненты тензоров упругих деформаций (£ и е3 - продольные деформации, £2 -сдвиговая) к элементов конструкции:
И = [[£1 £2 £3]1 [£1 £2 £3]2 ••• [£1 £2 £3]к]Т
Тензор Я определяется единственным параметром - углом жесткого поворота р, отложенным против часовой стрелки
Я = (е1е1+е2е2) 008 р + (е2е1-е1е2) Ыр, а угол рпри известном тензоре ^ и, например, ир определяется следующим образом:
р= аГ018[(^зир3-^4ир2+^1ир2-^2ир1) / (^1 ^3-^2 ^2+^2-^)].
Тензор упругой деформации, который понадобится позже для расчета напряжений и усилий, согласно выражениям (3) имеет вид
£е= и - I = ЯT ■ ^ ■ ир-1 - I.
Для конструкции это выражение записывается в матричной форме:
[£е] = [В] [*] - И. (4)
Здесь матрица [В] размерностью 3к*4к определяет связь между дисторсией и упругими деформациями конструкции в зависимости от смещений ее узлов и накопленной пластической деформации, а столбец [!2] содержит координаты единичного тензора и имеет размерность 3к:
и = [[1 0 1] [1 0 1] ... [1 0 1]]Т.
Свойства материала описываются моделью идеальной пластичности. Согласно этой модели при упругом деформировании интенсивность напряжений о (длина вектора напряжений в девиа-торном пространстве) меньше радиуса поверхности текучести Яс (найденного по пределу текучести при растяжении оТ и, следовательно, равного величине л/273 от), а при пластическом деформировании о = Яс; при этом упругие деформации связаны с напряжениями законом Гука:
о = С-- Ее,
где С - тензор констант упругости. Для конструкции закон Гука приобретает матричный вид
[О] = [С] [ее], (5)
где столбец [0] размерностью 3к содержит компоненты тензоров напряжений Коши (о1 и а3 -нормальные напряжения, о2 - касательное) элементов конструкции:
[о] = [[01 02 О3] 1 [01 02 Оз]2 ... [01 02 Оз]к]Т, а матрица [С] размерностью 3к*3к содержит константы упругости материала (Е - модуль упругости и ц - коэффициент Пуассона).
Плоская конструкция, состоящая из набора треугольных элементов, может быть представлена в виде наложенного на нее «корсета» - фермы, стержни которой шарнирно связанны в узлах и определяют сетку конечных элементов. Они могут растягиваться или сжиматься, но не могут изгибаться, что обеспечивает однородность напряженно-деформированного состояния элементов. Тензор напряжений элемента с помощью граничных условий позволяет найти распределенные силы, действующие на его сторонах, а наличие корсета - заменить систему этих распределенных сил сосредоточенными силами, приложенными к узлам корсета. Для произвольной плоской конструкции выражение, связывающее напряжения элементов и узловые силы имеет следующий вид:
[0] = №], (6)
где столбец [0] содержит координаты сил, приложенных к узлам конструкции:
[б] = [в 1 в 2 - в -]Т
(по аналогии со смещениями [д]), а матрица [Н] связывает координаты тензоров напряжений элементов конструкции с координатами сил, приложенных к узлам конструкции в ее деформированном состоянии. Очевидно, что матрица [Н] зависит от смещений, так как равновесие элемента рассматривается в деформированном состоянии.
Щербакова А.О. Расчет деформирования геометрически нелинейной
плоской конструкции из идеально пластичного материала
При кинематическом нагружении конструкции (то есть при известных смещениях узлов) расчет ведется шагами по времени, при этом на каждом шаге используются пластические деформации, найденные на предыдущем шаге расчета.
Схема расчета при кинематическом нагружении
1. Номер шага нагружения п = 1.
2. Начальные пластические растяжения [ир*] = [!2].
3. Смещения на текущем шаге [дп].
4. Функция уточнения пластических растяжений: [ирп] = Адп, ир*):
a. Дисторсии [^] = [£][дп]+ [Il];
b. Приближенные упругие деформации и напряжения:
[£е*] = [В(д, Ч*)т - [І2], [о*] = №е*];
о. Приближенный вектор напряжения [5] в девиаторном пространстве и его модуль 0!* (интенсивность напряжений). Для конечного элемента с координатами тензора напряжений 01*, 02* и 03* получим
51^^(01*+0з*), 52=-^(01*-0з*), 5з=л/2 02*, 0г* = (5^ + 522 + 5з2)Ш;
Л/6 л/2
d. Если в каком-либо элементе конструкции 01* > Яо, тогда укорачиваем вектор 5 до размера Яо, затем уточняем координаты тензора напряжений 0;
e. Координаты тензоров упругих растяжений
и = [С]1 [0] + И;
£ Координаты тензоров пластических растяжений
[ирп] = №п, ие)][^1
(здесь матрица [П] выражает связь между пластическими растяжениями и дисторсией элементов конструкции);
g. Если 0* < Яо, тогда не меняем координаты тензоров пластических растяжений
[ирп] = [ир*].
5. Координаты тензоров пластических растяжений после поворота элемента как жесткого целого [Ррп]. Например, для одного элемента
[Урп] = [Яп]Т [ирп] [Яп].
6. Координаты тензоров упругих растяжений после поворота элемента как жесткого целого
[^еп] = [В(Чп, Урп)] [*].
7. Вычисляем напряжения [0п] = [С] ([Уеп] - [I2]).
8. [вп] = №)] [0п].
9. [ир*]= [ирп], п = п + 1, возвращаемся к п. 3.
Расчет силового нагружения конструкции также ведется шагами по времени, при этом приходится использовать численные методы решения системы нелинейных уравнений
[д] = КТ1 ([в] - [Ь]),
где
[К] = [Н][С][В][Ь], [Ь] = ннсшвш - [І2]).
Верификация модели
С целью верификации предложенной расчетной модели в работе были сопоставлены величины изгибающих моментов М*, возникающих в прямоугольной балке при ее чистом изгибе, вычисленные с помощью предложенной модели, а также моменты
М = Ек3/(12р), М2 = (0т к3 / 4) - 0т3 р2/(3Е2), вычисленные с использованием точных решений, известных из курса сопротивления материалов. Здесь М1 - это момент, вычисленный при упругом деформировании балки, а М2 - за пределами упругости, Е - модуль упругости материала, р - радиус кривизны балки, к - высота балки. В расчетах использовали следующие исходные данные: материал балки - изотропный, идеально пластичный (Е = 2105 МПа, 0Т = 10 ГПа, коэффициент Пуассона л для простоты принят равным нулю), к = 8 м, длина балки I = 15 м, толщина балки равна 1 м, напряженное состояние - плоское, высота конечного элемента - 0,5 м.
На рис. 1 приведена схема балки в начальном состоянии (с разбивкой на конечные элементы) и в конечном, а на рис. 2 показаны зависимости изгибающих моментов от а - угла поворота пра-
Механика
вого торца балки. Расчет показал, что отличие между предложенным решением и соответствующим точным решением не превышает 1 %.
У -ж-
Q.
V
м
<-
Рис. 2. Зависимости изгибающих моментов от угла а Рис. 1. Схема чистого изгиба балки при у = 0 M1 и М2 - точные решения;
M* - решение с помощью предложенной модели;
М3 - решение с помощью пакета ANSYS
Для сравнения на рис. 2 показано изменение изгибающего момента M3, возникающего при чистом изгибе балки, вычисленное с помощью пакета прикладных программ ANSYS (в расчете использовали квадратные элементы plane182 соответствующего размера). Геометрические параметры конструкции и свойства материала использовали те же, что и в расчете с помощью предложенной модели. Расчет показал не только количественное, но и качественное различие между решением рассмотренной задачи в пакете ANSYS и точным решением: с увеличением а значение момента М3 не стабилизируется на величине 16• 1010 Нм, а продолжает расти. Это объясняется использованием в пакете ANSYS тензора логарифмической деформации, который при повороте главных осей конечных элементов в результате деформирования, теряет физический смысл [6-8].
Заключение
Предложена расчетная модель позволяющая, зная нагрузки, приложенные к плоской конструкции из идеально пластичного материала, найти соответствующие смещения ее точек и, наоборот, по заданным смещениям найти приложенные к конструкции силы. Модель включает в себя геометрические, статические и физические соотношения и основывается на методе конечных элементов. В качестве меры напряжений используется тензор напряжений Коши, а меру деформаций представляет тензор инженерной деформации.
Для верификации предложенной модели выполнен расчет деформирования балки при чистом изгибе. Сопоставление результатов расчета с соответствующим точным решением, известным из курса сопротивления материалов, позволяет сделать вывод об адекватности предложенной модели (различие между величинами изгибающих моментов, полученными двумя методами, составляет не более 1 %).
Литература
1. Truesdell, K. A first course in rational continuum mechanics / K. Truesdell. Baltimore: The Johns Hopkins University, 1972. - 264 p.
2. Chadwick, P. Continuum mechanics: concise theory and problems (2 ed.) / P. Chadwick. - Dover publications, 1998. - 193 р.
3. Belytschko, T. Nonlinear finite elements for continua and structures / T. Belytschko, W.K. Lin, B. Moran. - New York: John Wiley and sons, 2000. - 660 p.
4. Mase, G. Theory and problems of continuum mechanics / G. Mase. - New-York: McGraw-Hill Book Company, 1970. - 221 p.
5. An objective time-integration procedure for isotropic rate-independent rate-dependent elastic-plastic constitutive equations / G.G. Weber, A.M. Lush, T.A. Zavaliangos, L. Anand // International journal of plasticity. - 1990. - Vol. 6. - P. 701-744.
Щербакова А.О.
Расчет деформирования геометрически нелинейной плоской конструкции из идеально пластичного материала
6. Садаков, О.С. Конечные деформации в механике деформируемого твердого тела / О.С. Садаков // Вестник ЮУрГУ. Серия «Математика, физика, химия». - 2005. - Вып. 6. -№6(46). - С. 114-121.
7. Садаков, О.С. К расчетам напряженно-деформированного состояния конструкций в геометрически нелинейной постановке/ О.С. Садаков // Труды Международной конференции «Сне-жинск и наука. - 2003. Современные проблемы атомной науки и техники». - Снежинск: СГФТА, 2003. - С. 73-74.
8. Буслаева, О.С. Скаляр и тензор логарифмической деформации / О.С. Буслаева,
О.С. Садаков, А.А. Шапиро // Научно-технические ведомости СПбГТУ. - 2003. - № 3(33). -С. 125-129.
Поступила в редакцию 16 июня 2013 г.
Bulletin of the South Ural State University Series “Mathematics. Mechanics. Physics” ________________2014, vol. 6, no. 2, pp. 35-39
CALCULATION OF DEFORMATION FOR GEOMETRICALLY NONLINEAR PLANE STRUCTURE MADE OF PERFECT PLASTIC MATERIAL
A.O. Scherbakova
A computational model given helps to find displacements of a structure, made of a perfect plastic material, using loads applied to this structure. On the other hand it can find loads applied to the structure using displacements given. The model includes geometrical, statical and physical relations and it is based on the finite element method where finite elements represent triangular simplexes.
Keywords: geometrical nonlinearity, perfect plastic material, pure bending, deformation.
References
1. Truesdell, K. A first course in rational continuum mechanics. Baltimore: The Johns Hopkins University, 1972. - 264 p.
2. Chadwick P. Continuum mechanics: concise theory and problems (2 ed.). Dover publications, 1998. 193 p.
3. Belytschko T., Lin W.K., Moran B. Nonlinear finite elements for continua and structures. New York, John Wiley and sons, 2000. ISBN 0-471-98773-5. 660 p.
4. Mase G. Theory and problems of continuum mechanics. New-York: McGraw-Hill Book Company, 1970. 221 p.
5. Weber G.G., Lush A.M., Zavaliangos T.A., Anand L. An objective time-integration procedure for isotropic rate-independent rate-dependent elastic-plastic constitutive equations. International journal of plasticity. 1990. Vol. 6. pp. 701-744.
6. Sadakov O.S. Konechnye deformatsii v mekhanike deformiruemogo tverdogo tela (Finite deformations in mechanics of deformable solids). Vestnik YuUrGU. Seriya «Matematika, fizika, khimiya». 2005. Issue 6. no. 6(46). pp. 114-121. (in Russ.).
7. Sadakov O.S. K raschetam napryazhenno-deformirovannogo sostoyaniya konstruktsiy v geomet-richeski nelineynoy postanovke (To the calculation of tensed deformable state of the structure in geometrically nonlinear performance). Trudy Mezhdunarodnoy konferentsii «Snezhinsk i nauka - 2003. Sovremennye problemy atomnoy nauki i tekhniki» (Proceedings of the International Conference "Snezhinsk and Science - 2003. Modern Problems of Nuclear Science and Engineering"). Snezhinsk: SGFTA, 2003. pp. 73-74. (in Russ.).
8. Buslaeva O.S., Sadakov O.S., Shapiro A.A. Nauchno-tekhnicheskie vedomosti SPbGTU. 2003. no. 3(33). pp. 125-129. (in Russ.).
Received 16 June 2013
1 Scherbakova Alla Olegovna is Cand. Sc. (Engineering), Associate Professor, Applied Mechanics, Dynamics and Strength of Machines Department, South Ural state university.
E-mail: AllaScherbakova@list.ru
