CC BY
99
ДИНАМИКА, ПРОЧНОСТЬ, НАДЕЖНОСТЬ \
УДК 539.3
ЧИСЛЕННЫЙ АНАЛИЗ БОЛЬШИХ ПЛОСКИХ ДЕФОРМАЦИЙ АРОЧНОГО АМОРТИЗАТОРА
Представлены определяющие соотношения инкрементальной теории расчета больших плоских деформаций резиновой детали арочного амортизатора. Для описания упругости резины применена модель сжимаемого неогуковского материала. Приведены результаты расчета нагрузочных характеристик амортизаторов и полей напряжений в резиновой детали. Проанализировано влияние малой сжимаемости резины на расчетные характеристики амортизаторов.
E-mail: [email protected]
Ключевые слова: резиновые амортизаторы, большие деформации, нагрузочные характеристики, численные методы.
Арочные амортизаторы (рис. 1) применяются для виброизоляции оборудования на судах и плавсредствах различного назначения. Для обеспечения требуемых характеристик амортизаторов необходимо располагать методом расчета резиновых деталей амортизаторов при больших деформациях.
В работах [1, 2] расчеты арочных амортизаторов выполнялись с помощью метода Л.В. Канторовича с предсказанием закона изменения перемещений точек резиновой полосы по одному из направлений. Резина рассматривалась как неогуковский материал, т.е. несжимаемый
A. Е. Белкин, В. В. Семенов
B. К. Семенов
Рис. 1. Арочный амортизатор
материал с упругим потенциалом Трелоара [3]. Модель несжимаемого материала обеспечивает достаточную точность, если напряженное состояние резиновой детали характеризуется невысокими значениями компонент шарового тензора напряжений при существенных компонентах девиатора. В низких амортизаторах вследствие стеснения поперечных деформаций большой объем резины может находиться в условиях всестороннего сжатия. Поэтому предположение об абсолютной несжимаемости резины может привести к значительной погрешности при вычислении жесткости амортизатора.
Для детального анализа поведения амортизатора под нагрузкой, включающего в себя расчет напряжений и деформаций с учетом малой сжимаемости резины, необходимо применение метода конечных элементов (МКЭ).
Рассмотрим большие деформации резиновой детали амортизатора, происходящие в плоскости декартовых координат х^х2 (см. рис. 1). Для расчета используем инкрементальную теорию конечно-элементного анализа больших деформаций [4-7]. Все компоненты исследуемого напряженно-деформированного состояния считаем зависящими от параметра нагружения £ (условного времени). Координаты произвольной точки резиновой полосы амортизатора в недефор-мированном состоянии, т.е. в момент времени 0, обозначим как 0х^ (г = 1, 2,3), координаты той же материальной точки в деформированном состоянии в момент времени £ — как ьХг. Недеформированное состояние детали принимаем за отсчетную конфигурацию, в метрике которой формулируются все уравнения задачи. При плоской деформации 4Х1(°Х1, 0х2), 4Х2(°Х1, 0х2), 1 х3 = 0х3.
Для описания деформированного состояния введем градиент вектора места 4Е, компоненты которого образуют матрицу
[% ] =
д *хг д °х,
'д * xi д *xi
д°Х! д°Х2
д * Х2 д *Х2
д°Х1 д°Х2
0 0
0
1
Через градиент места выразим тензор меры деформаций Коши-Грина
ЬС = 4Ет •1Е и тензор деформаций Грина-Лагранжа
*е = ^(4Ет • *Е - I), 2
где I — единичный тензор, точкой обозначено скалярное произведение тензоров.
0
В качестве основного уравнения задачи равновесия рассмотрим уравнение принципа виртуальных работ, записанное в отсчетной конфигурации:
: 5*е = (1)
0 V
где * Б — второй тензор напряжений Пиолы-Кирхгофа; 0 V — объем тела в недеформированном состоянии; 5íW(exí) — виртуальная работа внешних сил; символ : обозначает двойное скалярное произведение тензоров.
В соответствии с идеей инкрементального метода предполагается, что на малом шаге нагружения £ + Д£ малыми являются приращения компонент тензора напряжений и тензора деформаций
ДБ = *+д*8 - Де = *+А*е — *е; связь между ними может быть линеаризована:
ДБ = Е : Де,
где Е — тензор секущих модулей упругости для рассматриваемого шага нагружения.
В первом приближении секущие модули упругости заменяются ка-
д * Б
сательными модулями Е « *Е = ——.
д*е
Для шага нагружения Д£ проводится линеаризация уравнения (1) в окрестности предположительно известного состояния £ по отношению к малым приращениям перемещений Дм = *+л*х — (г = 1, 2). Линеаризованное уравнение виртуальных работ имеет вид [4]
JfE : Ae) : 8Aed°V + j 'S : SA^dOV =
0у 0 V
= 5 í+AíWexí — / : 5Де <ТУ, (2)
0V
где Де, Дп — соответственно линейная и квадратичная относительно приращений перемещений части тензора Де = Де + Дп с компонентами
Де%3 = -^(Диг,3 + Дм Ц + ЬПи<гДии<3 + ЬПи<3 Дщ,г),
ДПгз = ^Ди^Ди^.
Отметим, что уравнение (2) составлено для того случая нагружения, когда внешние силы не зависят от перемещений, и поэтому не возникает задачи линеаризации выражения виртуальной работы внешних
мации г+дг£(1), а затем с помощью уравнении состояния и напряжения
сил 8г+дгЖежг. Именно такое нагружение рассматривается при расчете амортизатора.
Уравнение (2) используется для нахождения первого приближения Аи(1) к точным значениям приращений перемещении на шаге нагружения. Добавляя полученные приращения Аи(1) к уже накопленным перемещениям гщ, получаем г+дги(1). По найденным перемещениям, используя геометрические соотношения, вычисляют деформации г+дгг
г+дгс(1) ^ .
Ошибка такого вычисления (невязка уравнения виртуальных работ) составляет
Я _ г+дгБ(1) : 8г+дг£(1) (V.
0У
Для уменьшения невязки после каждого шага нагружения осуществляется итерационный процесс, исполняемый до тех пор, пока разница между работами внешних и внутренних сил не окажется меньше заданного малого значения.
Полагая, что на итерации к + 1
г+дги(^+1) _ г+дги(к) + Аи(к+1)
г+дг £(й+1) _ г+дг£(*0 + А е(к+1) + А п^+Ч г+дг8(^+1) _ г+дг§(й) + г+дгЕ(^) : (Ае(к+1) + А^+1)),
(Аи(к+1), Ае(к+1), Ап(к+1) — малые уточнения перемещений и деформаций), приходим к уравнению метода Ньютона
у (г+дгЕ(к) : Ае(к+1)) : 8Ае^^У + / г+дгБ(к) : 8Ап^СУ =
= 8г+дгЖ(ежг) - I г+дгБ(к) : 8*+дге(к)СУ (3)
0У
После достижения равновесия по найденным напряжениям Пиолы-Кирхгофа определяются истинные напряжения в материале г+дгр (напряжения Коши) с помощью следующего преобразования:
г+дгр _ г+дгр _ г+дгф _ г+дгрт/г+дг J
где _ $+д*У/($У _ det[í+ДíF] — мера объемной деформации
материала.
В настоящей работе в качестве модели упругого поведения резины рассматривается модель сжимаемого неогуковского материала, для
которого упругий потенциал имеет вид [6]
^(*/1с,* 7) = 2С(*/1С — 3) + 4л(*72 — 1) — (1л + с) 1п *7, (4)
где в, Л — постоянные упругости материала; * 11с — первый инвариант тензора меры деформаций *С.
Исходя из упругого потенциала (4), можно получить выражение для тензора напряжений Пиолы-Кирхгофа
д^ 1 = 2д^ = С(1 —* С-1) + -Л(*72 — 1)*С-1, (5)
д С 2
где *С-1 — тензор, обратный к тензору меры деформаций.
Запишем соотношения упругости (5) в осях главных деформаций в компонентной форме:
31 = в (1 — Л-2) + 1л (72 — 1) Л-2 (1 ^ 2 ^ 3), (6) 2
где Л1, Л2, Л3 — кратности удлинений; 7 = Л1Л2Л3; левый верхний знак £, указывающий значение параметра нагружения, опущен.
Истинные главные напряжения выражаются через напряжения Пиолы-Кирхгофа (6)
а1 = А"31 = С7-1(Л1 — 1) + ^Л(7 — 7-1) (1 ^ 2 ^ 3). (7)
Л2Л3 2
При малых деформациях соотношения упругости (7) переходят в закон Гука, постоянные упругости материала в, Л принимают смысл постоянных Ламе. Чтобы охарактеризовать жесткость материала при объемной деформации, введем параметр в = Л/2С. Отметим, что в пределах малых деформаций, подчиняющихся закону Гука, значение параметра в связано со значением коэффициента Пуассона V соотношением в = v/(1 — 2v).
Для иллюстрации выбранного закона связи напряжений и деформаций рассмотрим примеры простейших типов напряженного состояния. В случае одноосного напряженного состояния, исходя из условия = 0з = 0, можно получить сначала кратности удлинений в поперечных направлениях
1 Г^1 + 4в (в +1)Л? — 1
\2 _ \2 _ Л2 _ Л2 _ 2вЛ2
а затем зависимость а1(Л1). На рис.2,а показаны графики этой зависимости при значениях параметра в = 4,5 (коэффициент Пуассона V = 0,45) и в = 49,5 (V = 0,495). Здесь же приведены графики изменения условного напряжения, т.е. напряжения, отнесенного к начальной площади. Напряжения выражены в долях параметра упругости в. Как и следовало ожидать, при одноосном напряженном состоянии влияние
-1
-2
-4
- /3 = 49,5 /3 = 4,5
/г У /' / // - / / ß = 49 ' > = 4,5 ,5
-1
-2
0,5
1,0
1,5
-4
Ii S
Л gl // Ht /// fil //// tß' = 49,5
' H'ß -hl i r l'l 1 Il 1 1 i4'5 ,
0,5
1,0
1,5
0\/G 16
12
-4
-12 -16
- "So II V.4°
/ £ = 4,5
1 1 1 1
0,96 0,98 1,0 1,02 1,04 Ai
Рис. 2. Диаграммы растяжения-сжатия сжимаемого неогуковского материала:
а — одноосное состояние; б — равномерное двухосное состояние; в — равномерное трехосное состояние; сплошными линиями показаны истинные напряжения, штриховыми — условные напряжения
параметра ß, т.е. сжимаемости материала, на диаграмму растяжения-сжатия не существенно.
В случае двухосного растяжения-сжатия (рис. 2, б) влияние сжимаемости также мало, поскольку деформация в направлении третьей оси не стеснена. Однако при трехосном напряженном состоянии ситуация меняется. Как видно из графиков на рис. 2, в, теперь малая объемная деформация становится значимым фактором в законе упругого поведения материала.
Вернемся к плоскому деформированному состоянию амортизатора. Для этого состояния соотношения упругости (5) в компонентной форме записываются следующим образом:
tSil = G (1 - 4C22/ J2) + h (1 - 1/ J2)tC22,
= G (1 - *СП J2) + 2л (1 - 1/íJ2) *СШ ^12 = ^12/ J2 - 2л(1 -1/^)^12,
íSзз = 2л(^2 -1).
Дифференцируя выражения для напряжений по компонентам тензора деформаций, найдем касательные модули упругости
^111 = ^+Л)СC22/íJ2)2, ^22 = ^2211 = Л+(2G+Л)(íCl2/íJ2)2, 'Е1112 = 'Е1121 = 'Е1211 = 'Е2111 = -(2G + Л)(ÍC12ÍC22/Í J4),
^2222 = (2G + Л)(íCll/íJ2)2, 'Е2212 = 'Е2221 = 'Е1222 = 'Е2122 = - (2G + Л)(í C12í C11/í J4),
t тр _t тр _t тр _t TP _
e1212 — e1221 — e2112 — E2121 —
= (G + 2л) (tCiitC22 + tC122)/tJ4 - ^A,
из которых составим симметричную матрицу касательных модулей
['е ] =
ttt Eiiii E1122 E
1111 e1122 e1112 tt 2211 e2222 e2212 tt 1211 e1222 e1212
Переходя к матричной записи уравнений, принятой в МКЭ, представим связь приращений напряжений и деформаций в виде
где } = { ^11 ÍS22 %2 }, {е} = { ^11 ^22 2*^ } -векторы-столбцы.
Из линеаризованного уравнения виртуальных работ (2) в результате конечно-элементной дискретизации следуют уравнения МКЭ для шага нагружения [4]:
(ВД + рк^]) Д{д} = {*+д*Р}еа* - ГРЫ, (8)
где Д{д} — искомый вектор приращений узловых перемещений; К^], — матрицы тангенциальной жесткости конструкции при линейных (Де) и нелинейных (Дп) деформациях соответственно; }еХ и {Р— векторы внешних и внутренних узловых сил.
Независимо от типа используемого конечного элемента общий вид матриц, входящих в уравнение (8), таков:
М = / №]ТЕV,
0у
ЗД = 1 [Б^]т ^[Б^ь 0 V,
0 у
{Р}ш = | ГБь]т{S0V,
0у
где ] — матрица связи линейных деформаций { Де11 Де22 2Де12 } с узловыми перемещениями Д{д}; ] — матрица связи производных перемещений { Ди1,1 Ди1,2 Ди2,1 Ди2,2 } с узловыми перемещениями Д{д };
4Sи 0 0
4 Sl2 4 S22 0 0 о 0 4 ^ 4 Sl2 0 0 t Sl2 t S22
— матрица начальных напряжений.
Вид матриц ], S] может быть пояснен следующим представлением удельной работы начальных напряжений:
48 : 5Дп = 4^2Д^МДи^) =
= 5{ ДиМ Ди1,2 Ди2,1 Ди2,2 }т S]{ Ди1,1 Ди1,2 Д^д Д^,2 } =
= 5 Д{д }т [Б^ь ]т ^[Б^ь ]Д{? }.
На каждом шаге нагруже-ния приращения перемещений Д{?} уточняются итерационно по методу Ньютона:
([¿+д* Кь ](к) +
+ ГД К^ь ](к) )Д {д }(к+1) =
= { Р}еХ -{ Р}гп4,
Д{д }(к+1) =
= Д{д }(к) + Д {д }(к+1).
В настоящей работе расчеты выполнялись с помощью изопараметрического четырех-узлового конечного элемента с билинейной аппроксимацией перемещений. Далее приведены некоторые результаты расчетов.
f S ] —
Рис. 3. Нагрузочные характеристики амортизаторов Сплошные линии соответствуют Ь/Н = 0,3; штриховые — Ь/Н = 0,4)
Геометрия резинового элемента амортизатора в исходном состоянии определяется двумя параметрами: углом наклона граней элементов к горизонтальному основанию а и относительной толщиной Ь/Н. На рис. 3 показаны нагрузочные характеристики для амортизаторов с параметрами а = 60°, Ь/Н = 0,3; 0,4. Из рис.3 следует, что учет сжимаемости материала заметно влияет на значения предельных сил. При параметрах жесткости материала в = 4,5 и 49,5 различие между предельными силами составляет 21 и 25 % для амортизаторов с относительными толщинами Ь/Н = 0,3 и 0,4 соответственно. В рассмотренном диапазоне изменения в жесткость амортизатора на линейном участке нагрузочной характеристики изменяется на 15 % и более.
На рис. 4, а и б показаны конфигурация нагруженного амортизатора и законы распределения по толщине одного из главных истинных напряжений — напряжения, направленного вдоль деформированной
А/Я = 0,2 А/Я =0,3
б
Рис. 4. Распределение напряжений в резиновом элементе амортизатора при различных значениях параметров:
а — b/H = 0, 4; ß = 4, 5; б — b/H = 0, 4; ß = 49, 5; напряжения отнесены к модулю упругости G
полосы. Из трех действующих главных напряжений это напряжение является наибольшим по абсолютной величине. Эпюры напряжений построены в координатных сечениях 0x2 = const. Приведенные результаты свидетельствуют о весьма существенном влиянии сжимаемости резины на значения напряжений, что, конечно, является ожидаемым при режиме заданных смещений торцов амортизатора. В целом можно утверждать, что учет малой сжимаемости резины вносит заметные изменения в результаты расчета амортизаторов.
СПИСОК ЛИТЕРАТУРЫ
1. Клюева Г. П. Исследование напряженного и деформированного состояния резиновых амортизаторов арочного типа // Расчеты на прочность. - 1977. -Вып. 18. -C. 59-72.
2. Б и д е р м а н В. Л., Ж и с л и н А. Я. Нелинейные характеристики резиноме-таллических упругих элементов // Расчеты на прочность. - 1979. - Вып. 20. -C. 107-122.
3.Трелоар Л. Физика упругости каучука / Пер. с англ. под ред. Е.В. Кувшинского. - М.: ИИЛ, 1953. - 240 с.
4. Bathe K. J. Finite element procedures. - Prentice-Hall, Englewood Cliffs, 1996. - 1036 p.
5. Васидзу К. Вариационные методы в теории упругости и пластичности / Пер. с англ. под ред. Н.В. Баничука. - М.: Мир, 1987. - 542 с.
6. Wriggers P. Nonlinear finite element methods. - Springer-Verlag Berlin Heidelberg, 2008. - 559 p.
7. B o n e t J., Wood R. D. Nonlinear continuum mechanics for finite element analysis. - Cambridge University Press, 1997. - 245 p.
Статья поступила в редакцию 15.12.2010
Александр Ефимович Белкин родился в 1951г., окончил в 1974 г. МВТУ им. Н.Э. Баумана. Д-р техн. наук, профессор кафедры "Прикладная механика" МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в области численных методов расчета конструкций, механики пневматических шин.
A.Ye. Belkin (b. 1951) graduated from the Bauman Moscow Higher Technical School in 1974. D. Sc. (Eng.), professor of "Applied Mechanics" department of the Bauman Moscow State Technical University. Author of more than 60 publications in the field of numerical methods of design of constructions, mechanics of pneumatic tires.
Виктор Владимирович Семенов родился в 1985 г., окончил в 2008 г. МГТУ им. Н.Э. Баумана. Аспирант кафедры "Прикладная механика" МГТУ им. Н.Э. Баумана. Специализируется в области численных методов расчета конструкций.
V.V. Semenov (b. 1985) graduated from the Bauman Moscow State Technical University in 2008. Post-graduate of "Applied Mechanics" department of the Bauman Moscow State Technical University. Specializes in the field of numerical methods of design of constructions.
Владимир Константинович Семенов родился в 1987 г. Студент МГТУ им. Н.Э. Баумана, инженер научно-технического центра "Интайр". Специализируется в области механики пневматических шин.
V.K. Semenov (b. 1987) — student of the Bauman Moscow State Technical University. Engineer of Scientific and Technical Center "Intair". Specializes in the field of mechanics of pneumatic tires.
