Разработка алгоритмов численного исследования конструкций из неоднородной среды методом конечных элементов Текст научной статьи по специальности «Строительство и архитектура»

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

УДК 001.891.573; 519.6

А. С. Васильев, Н. А. Тарануха

РАЗРАБОТКА АЛГОРИТМОВ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ КОНСТРУКЦИЙ ИЗ НЕОДНОРОДНОЙ СРЕДЫ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ

В статье предложена новая математическая модель исследования прочности конструкций из неоднородной среды. Использован метод конечных элементов. Разработаны численные процедуры для реализации математической модели. Разработана новая матрица жёсткости для композитного конечного элемента. Матрица жёсткости нелинейная. Матрица жёсткости учитывает изменение механических характеристик материалов в составе конечного элемента. Включение различных материалов в конечный элемент осуществляется на основе идеи комбинирования объёмов. Количество различных компонентов в составе неоднородного материала может быть достаточно большим. Разработан программный комплекс для численного исследования напряжённо-деформированного состояния композитных конструкций. Приведены алгоритмы работы программ. Алгоритмы построены в виде блок-схем. Выполнено тестирование и практическая апробация математической модели и программного комплекса. Сформулирована новая методика исследования композитных материалов.

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

Ключевые слова: матрица жёсткости; неоднородная среда; композит; алгоритм; блок-схема; математическая модель; критерий прочности; предельное состояние.

Введение

Неоднородные материалы, к которым относится большое число композитных материалов, активно применяются в различных современ-

Васильев Алексей Сергеевич — старший преподаватель (Приамурский государственный университет имени Шолом-Алейхема, Биробиджан); e-mail: [email protected].

Тарануха Николай Алексеевич — доктор технических наук, профессор, заведующий кафедрой кораблестроения (Комсомольский-на-Амуре государственный технический университет, Комсомольск-на-Амуре); e-mail: [email protected].

© Васильев А. С., Тарануха Н. А., 2016

78

ных отраслях промышленности. В частности, активно используются композиты из нестандартизированного материала [6; 7; 8]. Большая распространённость композитов требует изучения их свойств. Эти свойства влияют на жёсткость и несущую способность конструкций на разных этапах их работы под нагрузкой. Расчёт таких конструкций целесообразно выполнять методом конечных элементов с использованием шаго-во-итерационных процедур на базе общих принципов механики деформируемого твёрдого тела и численных методов решения физически-нелинейных задач.

Существует много программных продуктов, позволяющих моделировать и производить расчёт композитных конструкций, в частности, из железобетона. К ним относятся МЛБТКЛМ, ЛМБУБ, ЛИРА. Данные программные продукты выполняют расчёты методом конечных элементов с использованием современных теорий прочности, соответствующих специфике работы конкретного материала.

Аппроксимационные зависимости для диаграмм деформирования материалов

Как известно, зависимость между деформациями и напряжениями может быть учтена при помощи диаграмм деформирования материалов. При этом в данной работе говорится о кратковременном нагружении. Рассматриваются зависимости, которые связывают относительные деформации с напряжениями (еЬ — оЬ) при одноосном сжатии и растяжении. Следует отметить, что иллюстрирующим материалом, характеризующим работу композитов, может быть железобетон. Это наиболее распространённый и изученный композитный материал.

а)

О;

СЬ0=1

Gm

Рис. 1. Диаграммы деформирования материалов в составе композита:

а — связующего (бетона), б — армирующих элементов (арматуры)

Диаграмма деформирования связующего показана на рисунке 1а. Прочность и деформативность армирующих элементов характеризуется диаграммой на рисунке 1б. При этом согласно работе [4]:

79

Л; _ <T.„-R.

Е,

£ ,=-!L, Е - W

s

Здесь os, у — временное сопротивление, И8 —расчётное сопротивление, Б8 — модуль упругости армируюшдх элементов, es,y = 0.02.

Матрица жёсткости нелинейного композитного конечного элемента

Композитный конечный элемент представляет собой совместную работу материала матрицы и «размазанной» по объёму конечного элемента арматуры при сохрании свойства сплошности материала. В работах [9; 10] авторами использовалась матрица жёсткости на основе идеи комбинирования площадей материалов в поперечном сечении конечного элемента. В данной работе для вычисления упругих характеристик композитного КЭ применим комбинирование объёмов материалов, из которых он состоит:

р = Е(.ех\У1+Е(.ех)2У2 + ... + Е(.ехХУя (2)

кУь ТЛ , ТЛ , , ТЛ

Mhyb ~~

(3)

У1+¥2+... + К

где Е1, Е2, ... Еп — соответственно нелинейные модули упругости материалов в составе композитного КЭ; У1, У2, ...Уп — соответственно объёмы каждого из материалов в составе композитного КЭ (где п может быть достаточно большим), (^)х, (^2)х и (^п)х — соответственно коэффициенты Пуассона для каждого материала в композитном КЭ. Подставляя (2) и (3) в стандартную матрицу упругости, получим матрицу упругости для композитного нелинейного объёмного конечного элемента [0]ьуь, сочетающего совместную работу различных сред.

ДуЬ =

^hyb (1 ("hyb)

(1 + ^X1-2/^)

1

Ayb

Мпр

АуЬ АуЬ

1-/V 1-/V

/'hyb

W'hyb

М\ф

1

1 Aiyb 0

0

0

0 0 0

0 0 0

!-2/'hyb 2(1" /'hyb) 0

0

0 0 0 0

1-2/V 2(1"/^) 0

(4)

1-2 fi,

hyb

Именно с матрицей [0]ьуь, связаны особенности данной работы.

Конечный элемент в форме прямоугольного параллелепипеда содержит восемь узловых точек (рис. 2). Матрица деформаций [В] содер-

0

0

0

0

0

80

жиг восемь блоков [в] = [|в]0) [в]С2) ■■■ [в]'4 ■■■ [я]{3)Стандартный блок её равен [4]:

¿(1 + М)(1 + С,С)

Bf =[Ф]С\(Ш)=1~

о

о

а

0 b О О

0

о

b 0

М+ММ+М

ЬП Чу! Л ЬкЬ a

(5)

При этом каждый из блоков матрицы жёсткости вычислен по фор-

муле:

(6)

Здесь [^ — матрица жёсткости 1-го КЭ; k и ] — номера узлов КЭ. При этом каждый элемент матрицы жёсткости [Щкц состоит из подматрицы третьего порядка:

(7)

Вводя переменные k = 1:8 и ] = 1:8 после взятия определённого интеграла (8), имея в виду постоянство [0]ьуь по объёму элемента, получим:

1 /'hyb

1-Х

2

1

1

1

1

1

1

3

32(l+/ihb)(l-2/;hb)

....(8)

K

US-

ab 3

Zklj +Иъyb^jflk -Vjbk)

уа{Е{ех\Г1+Е{ех)2Г2... + Е{£х)1,Г1,) "' 32(1 + ^X1-2/^)

Остальные элементах матрицы жесткости определяются аналогичным образом. Здесь п = 1+3(]-1); т = 1+3(к-1). В результате получаем матрицу жёсткости КЭ размером 24*24.

В качестве метода решения нелинейных задач был использован метод последовательных приближений [3; 4]:

c

K

l+n.l-m

1

2

81

где n — номер этапа приближения. Алгоритмы, реализующие данную математическую модель на основе МКЭ, использованы в разработанных авторами программах: «Composit» [1], «Strength» [2], «Ultimate State», «Программа для расчёта конструкций из композитных материалов». Данные алгоритмы предназначены для расчёта конструкций из композитных материалов при простом напряжённом состоянии, т. к. использует диаграммы деформирования материала при одноосном растяжении и сжатии. При последовательном увеличении нагрузки на конструкцию происходит постепенное её разрушение вследствие изменения механических характеристик материалов. Характеристики при этом меняются согласно диаграммам е—о.

Алгоритм расчета конструкций из композитных материалов.

Расчёт вёлся с помощью алгоритмов программ Composit [1], Strength [2] и Ultimate State, основанных на языке Matlab Программы Composit и Strength, зарегистрированы в Реестре программ для ЭВМ и имеют свидетельства об официальной регистрации.

Сборка глобальной матрицы жёсткости выглядит следующим образом (рис. 2), где z — множество i-ых КЭ, принадлежащих промежутку 1:n.; W — матрица индексов конструкции, состоящая из n строк и p столбцов; а, в = 1, ..., p; p — количество перемещений в КЭ. Для данной задачи n = 4, так как при стыковке элементов суммируются матрицы жёсткости КЭ бетона, композитные КЭ бетона и арматуры, КЭ арматуры без уёета бетона, КЭ, не имеющие в своём составе материала, когда матрица упругости D = 0.

Алгоритм работы программы представлен рисунком 3.

Здесь R2 — расчётное сопротивление второго материала (арматуры), V1 и V2 — объёмы соответственно первого и второго материалов в конечном элементе.

Сравнение расчёта автора с экспериментальными исследованиями статически определимых балок из композитного материала.

Рассмотрим железобетонную статически определимую балку прямоугольного сечения без предварительного напряжения [5]. Нагрузка в виде сосредоточенных сил прикладывалась симметрично ступенями с шагом ДР = 5 кН. Критерием исчерпания прочности железобетонных образцов являлось появление текучести в армирующих элементах (арматуре).

Расчётная схема, геометрические размеры и схема армирования представлены на рисунке 4.

Исходные характеристики материалов соответствуют опытным значениям и принимались равными: связующее (бетон) — Rb = 18.5 МПа, Rbt = 1.6 МПа, Еь = 30000 МПа, ц = 0.2; армирующие элементы (арматура): Rs = 510 МПа, Es = 190000 МПа.

(9)

82

_±_

Узловые перемещения ¡-го элемента

с—■

Рис. 2. Схема формирования и расчета системы МКЭ

Предварительные расчёты по СНиП: высота сжатой зоны х = 72.9 мм, разрушающая нагрузка Р = 129.44 кН, первые трещины появились при 17.67 кН. Прогиб к моменту разрушения 7.01 мм. Разрушение конечно-элементной модели произошло вследствие преобладающего действия изгибающего момента в середине пролёта балки.

Значения прогиба характеризуются максимальными перемещениями в середине балки под действием вертикальной нагрузки. На рисунке 5а показано изменение максимальных перемещений балки на различных этапах нагружения. Приводится сравнение расчётов в ПК COMSOL [5], ПК ЛИРА, автора и результатов эксперимента. Разрушающая нагрузка в расчёте автора — 140 кН. Прогиб при этом составил 8,6 мм.

83

Да

Конец

Рис. 3. Блок-схема алгоритма работы программы

Рис. 4. Модель балки: а — расчётная схема железобетонной балки, б — схема разбивки на конечные элементы

84

а)

б)

P-100kH I^SOkH

p=lokH р-юкн

J 15°

\ \

Расчет A.A. Попова и др. в \\ * ПК COMSOL Multiphysics \

Расчет автора | /

Напряжения, МПэ

Рис. 5. Сравнение результатов расчёта: a — сравнение изменений прогибов под нагрузкой; б — сравнение некоторых эпюр нормальных напряжений в поперечном сечении балки (расчёта автора с наибольшими значениями аналогичных эпюр, полученными в ПК COMSOL Multiphysics)

При потере несущей способности нормальные напряжения в верхних слоях балки составили —16.3 МПа. На рисунке 5б изображены эпюры напряжений в плоскости симметрии балки на различных этапах на-гружения и их сравнение с эпюрами, полученными автором.

В таблицах 1 и 2 приводятся сопоставление результатов расчётов: прогибов под нагрузкой в момент разрушения и предельной нагрузки на конструкцию.

Таблица 1

Максимальный прогиб Отклонение результатов расчётов от COMSOL Multiphysics 4.3 b

№ COMSOLM ultiphysics 4.3 b Расчёт автора Эксперимент СНиП Расчёт ПК ЛИРА Расчёт автора Эксперимент СНиП Расчёт ПК ЛИРА

А, мм А, мм А, мм А, мм А, мм % % % %

1 7.28 8.6 4.15 7.01 8.35 18.13 42.93 3.85 14.70

2 7.28 8.6 5.32 7.01 8.35 18.13 26.92 3.85 14.70

3 7.28 8.6 6.1 7.01 8.35 18.13 16.21 3.85 14.70

4 7.28 8.6 6.72 7.01 8.35 18.13 7.69 3.85 14.70

Таблица 2

Разрушающая нагрузка Отклонение результатов расчётов от COMSOL Multiphysics 4.3 b

№ COMSOLM ultiphysics 4.3 b Расчёт автора Эксперимент СНиП Расчёт ПК ЛИРА Расчёт автора Эксперимент СНиП Расчёт ПК ЛИРА

P, кН P, кН P, кН P, кН P, кН % % % %

1 130 140 106 129.44 150 -7.69 18.4 0.43 15.38

2 130 140 121 129.44 150 -7.69 6,92 0.43 15.38

3 130 140 133 129.44 150 -7.69 -2.31 0.43 15.38

4 130 140 135 129.44 150 -7.69 -3.84 0.43 15.38

85

Из таблиц видно, что погрешность расчёта авторов от расчёта в программе COMSOL Multiphysics 4.3 b составляет около 18 % и 8 % для максимального прогиба и разрушающей нагрузки соответственно. В то же время общая погрешность результатов экспериментов от расчёта COMSOL составляет 23.44 % и 4.79 % для максимального прогиба и разрушающей нагрузки. Наблюдается удовлетворительное совпадение расчетных и опытных данных.

Заключение

Предложенные алгоритмы и математическую модель возможно применять:

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

— для исследования напряжённо-деформированного состояния, установления перемещений, напряжений и деформаций на каждом этапе нагружения конструкции. Предельное состояние конструкции может наступать при разрушении связующего в поперечном сечении конструкции;

— для установления опасных зон и зон разрушения конструкции на различных этапах нагружения, включая предельные;

— для расчётов композитов, имеющих в своем составе инородные включения, газовые пузыри, технологические отверстия. При этом можно учитывать любое количество составляющих внутри композита в пределах конечного элемента.

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

1. Васильев А. С., Тарануха Н. А. Программное обеспечение для расчёта перемещений в конструкциях из композитных материалов методом конечных элементов (ПО «Composit»): свидетельство о регистрации программы для ЭВМ. № 2015610761. 2015.

2. Васильев А. С., Тарануха Н. А. Программное обеспечение для исследования напряжённо-деформированного состояния объектов на различных стадиях нагружения (ПО «Strength»): свидетельство о регистрации программы для ЭВМ. № 2015610762. 2015.

3. Клованич С. Ф., Мироненко И. Н. Метод конечных элементов в механике железобетона. Одесса, 2007. 111 с.

4. Клованич С. Ф., Безушко Д. И. Метод конечных элементов в расчётах пространственных железобетонных конструкций. Одесса: Изд-во ОНМУ, 2009. 89 с.

5. Попов А. А., Хатунцев А. А., Шашков И. Г., Кочетков А. В. Пространственный деформационный нелинейный расчёт железобетонных изгибаемых конструкций методом конечных элементов / / Науковедение: интернет журнал. 2013. № 5. URL: http: / / cyberleninka.ru/article/n/prostranstvennyy-deformatsionnyy-

86

ne1meynyy-raschet-zhe1ezobetonnyh-izgibaemyh-konstruktsiy-metodom-копесЬпуЬ-е1ешеп1о¥

6. Тарануха Н. А., Жеребко К. В., Петрова А. Н., Петров М. Р. Математическая модель шарнирной стержневой системы с большими перемещениями узлов / / Известия высших учебных заведений. Строительство. 2003. № 3. С. 12 — 18.

7. Тарануха Н. А., Петрова А. Н., Любушкина Н. Н. Механика морских динамических систем с большими деформациями из нестандартизированного материала / / Морские интеллектуальные технологии. 2010. № 3(9). С. 56 — 59.

8. Тарануха Н. А., Петрова А. Н., Любушкина Н. Н. Колебания динамических систем с большими деформациями из нестандартизированного материала // Ученые записки КнАГТу. 2010. № Ш-1(3). С. 4 — 11.

9. Тарануха Н. А., Васильев А. С. Алгоритмы и модели при численном проектировании композитных сред на заданные характеристики для морских сооружений / / Ученые записки КнАГТУ: Науки о природе и технике. 2015. № 1-1(21). С. 81—86.

10. Тарануха Н. А., Васильев А. С. Численное исследование предельной несущей способности конструкций из композитных материалов / / Морские интеллектуальные технологии: Кораблестроение, информатика, вычислительная техника и управление. 2015. № 3(29). Т. 2. С. 27 — 32.

* * *

Vasiliev Alexei S., Taranukha Nikolay A.

DEVELOPMENT OF ALGORITHMS OF NUMERICAL INVESTIGATION OF STRUCTURES MADE OF INHOMOGENEOUS MEDIA BY FINITE ELEMENT METHOD

(Sholom-Aleichem Priamursky State University, Birobidzhan; Komsomolsk-on-Amur State Technical University, Komsomolsk-on-Amur) The paper proposes a new mathematical model for researching strength of structures of inhomogeneous medium. The method of finite elements is used. Numerical procedures for the implementation of the mathematical model are designed. A new stiffness matrix composite finite element is developed. The stiffness matrix is nonlinear. The stiffness matrix takes changes of the mechanical properties of materials as a part of the finite element. The inclusion of different materials in the final element is based on the idea of combining the volumes. The quantity of different components in inhomogeneous material can be quite large. A software package is developed for the numerical study of stress-strain state of composite structures. The work algorithms of programs are shown. The algorithms are shown in a block diagram form. Testing and practical testing of mathematical models and software are performed. A new method of researching inhomogeneous materials is formulated.

The report examines the structure of the inhomogeneous medium. We consider the stress-strain state of structures of inhomogeneous medium in their limiting states. The practical calculations of reinforced concrete beams are performed. We get specific results for specific designs. The comparison of the results on the proposed mathematical model with experimental results of other authors is made. The comparison with the calculations in other software products, and with the results of calculations by the method of normative documents is made. General conclusions and recommendations on the application of the proposed algorithms, the mathematical model and stiffness matrix are formulated.

Keywords: stiffness matrix; inhomogeneous medium; composite; algorithm; block diagram; mathematical model; strength criterion; limit state.

References

1. Vasil'ev A. S., Taranukha N. A. Programmnoe obespechenie dlya rascheta pereme-shcheniy v konstruktsiyakh iz kompozitnykh materialov metodom konechnykh ele-mentov

87

(PO «Composit»): svidetel'stvo o registratsii programmy dlya EVM (Software for the calculation displacements in structures made of composite materials using finite element method (PO «Composit»): The certificate of registration of a computer program), no. 2015610761, 2015.

2. Vasil'ev A. S., Taranukha N. A. Programmnoe obespechenie dlya issledovaniya napryazhenno-deformirovannogo sostoyaniya ob"ektov na razlichnykh stadiyakh nagruzheniya (PO «Strength»): svidetel'stvo o registratsii programmy dlya EVM (Software for the study stress-strain state of objects at various stages of loading (PO «Strength»): The certificate of registration of a computer program), no. 2015610762, 2015.

3. Klovanich S. F., Mironenko I. N. Metod konechnykh elementov v mekhanike zhele-zobetona (Finite Element Method in the mechanics of reinforced concrete), Odessa, 2007. 111 p.

4. Klovanich S. F., Bezushko D. I. Metod konechnykh elementov v raschetakh prostranstvennykh zhelezobetonnykh konstruktsiy (Finite Element Method nonlinear analysis of spatial reinforced concrete structures), Odessa, ONMU Publ., 2009. 89 p.

5. Popov A. A., Khatuntsev A. A., Shashkov I. G., Kochetkov A. V. Spatial deformation nonlinear calculation of reinforced concrete bent designs by the method of final element [Prostranstvennyy deformatsionnyy nelineynyy raschet zhelezobetonnykh izgibaemykh konst-ruktsiy metodom konechnykh elementov], Naukovedenie: internet journal, 2013, no. 5. Available at: http://cyberleninka.ru/ article/n/prostranstvennyy-deformatsionnyy-nelineynyy-raschet-zhelezobetonnyh-izgibaemyh-konstruktsiy-metodom-konechnyh-elementov

6. Taranukha N. A., Zherebko K. V., Petrova A. N., Petrov M. R. Mathematical model of the hinged rod system with large displacements of nodes [Matematicheskaya model' sharnirnoy sterzhnevoy sistemy s bol'shimi peremeshcheniyami uzlov], Izvestiya vysshikh uchebnykh zavedeniy. Stroitel'stvo, 2003, no. 3, pp. 12 — 18.

7. Taranukha N. A., Petrova A. N., Lyubushkina N. N. The mechanics of marine dynamical systems with large deformations of non-standard material [Mekhanika morskikh dinamiche-skikh sistem s bol'shimi deformatsiyami iz nestandartizirovannogo materiala], Morskie intellektual'nye tekhnologii, 2010, no. 3(9), pp. 56—59.

8. Taranukha N. A., Petrova A. N., Lyubushkina N. N. The Determining characteristics of the liquid non-standard material of the dynamic system of the elastic system [Kolebaniya dinamicheskikh sis-tem s bol'shimi deformatsiyami iz nestandartizirovannogo materiala], Uchenye zapiski KnAGTU, 2010, no. III-1(3), pp. 4—11.

9. Taranukha N. A., Vasil'ev A. S. Algorithms and models for numerical designing composite media to the specified characteristics for marine structures [Algoritmy i modeli pri chislennom proektirovanii kompozitnykh sred na zadannye kharakteristiki dlya morskikh sooruzheniy], Uchenye zapiski KnAGTU: Nauki o prirode i tekhnike, 2015, no. I-1(21), pp. 81—86.

10. Taranukha N. A., Vasil'ev A. S. Numerical research of the limit bearing capacity of constructions from composite materials [Chislennoe issledovanie predel'noy nesushchey sposobnosti konstruktsiy iz kompozitnykh materialov], Morskie intellektual'nye tekhnologii: Korablestroenie, informatika, vychislitel'naya tekhnika i upravlenie, 2015, no. 3(29), vol. 2, pp. 27—32.

•Jc -Jc -Jc

88