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

УДК 593.3

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

© 2012 А.В. Зайцев, А.В. Кислицын, А.В. Кутергин, А.А. Фукалов

Пермский национальный исследовательский политехнический университет

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

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

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

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

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

Зайцев Алексей Вячеславович, кандидат физико-математических наук, доцент кафедры «Механика композиционных материалов и конструкций». E-mail: [email protected]

Кислицын Антон Владимирович, аспирант Кутергин Алексей Владимирович, магистрант Фукалов Антон Александрович, аспирант

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

Определение эффективных упругих модулей квазиизотропной матрицы на основе полиэфирной смолы Виналкид 550 (модуль Юнга и коэффициент Пуассона - Е=5,3 ГПа и и=0,25), наполненной мраморной мукой (модули Юнга в плоскости изотропии и в нормальном ей направлении £=55,0 ГПа и Е = 23,0 ГПа и коэффициенты Пуассона 1)=0.29 и и = 0,32) до концентрации 26%, проводилось на основе численного решения краевых задач в конечно-элементном пакете ANSYS 12.0 для независимых реализаций представительных объемов, характерные размеры которых определялись при помощи метода, предложенного авторами [1-4], а сами пространственные случайные структуры были сгенерированы при помощи алгоритмов [5-8]. Полученные результаты показали, что наполненная матрица статистически изотропна, имеет эффективный модуль Юнга Е=8,276 ГПа и коэффициент Пуассона и=0,24, которые попадают в ограничения Фойгта-Рейса и находятся вблизи нижней границы оценки Хашина-Штрихмана. Для проверки этих результатов были также спрогнозированы свойства полимербетонов аналитически в рамках полидисперсных моделей механики композитов [9, 10]. Отклонение численных значений эффективного модуля объемного сжатия матрицы, наполненной мраморной мукой, от аналитической оценки не превышает 5%.

Вычислительные эксперименты в пакете ANSYS 12.0 позволили также определить эффективные упругие модули полимербетона (модули Юнга и коэффициенты Пуассона в плоскости

изотропии Е=29,20 ГПа, и=0,19 и в нормальном к ней направлении ЕЕ = 27,30 ГПа и v = 0,21), матрица которого - полиэфирная смола Виналкид 550, наполненная мраморной мукой, а армирующая фаза - сферические трансверсально-изотропные частицы диабаза, образующие пространственную гексагональную структуру с объемной долей минеральных включений 60%.

Контейнеры для хранения высокоагрессивных веществ являются тяжелыми составными цилиндрами, изготовленными из анизотропных материалов. Составные части этих конструкций ограничены цилиндрическими поверхностями радиусов а и b (a<b), а поверхность контакта находится на расстоянии c от оси симметрии. Будем считать, что материал обоих цилиндров линейно-упругий однородный, с постоянной плотностью по всему объему, обладает цилиндрической ортотропией. Наиболее распространенными видами нагрузки являются статическое или изменяющееся во времени внешнее и/или внутреннее давление и собственный вес. В простейшем случае, пренебрегая влиянием торцевых нагрузок, контейнер может быть представлен в виде бесконечного составного горизонтального ортотропного цилиндра, жестко закрепленного по внешней поверхности и находящегося в равновесии под действием массовых сил.

В силу бесконечной протяженности цилиндра функции перемещений, деформаций и напряжений не зависят от осевой координаты z, цилиндрической ортогональной системы координат r,0 и z, вследствие чего уравнения равновесия могут быть записаны в следующем виде:

дг

1 дт( к) _(к) (к) + 1 дтг8 + агг ст99 = р( к)

г д8

cos 8

1 да!

.(к)

г д8

- +

г

дт(к)

^ + 2

дг

(к)

тг 8 =_р( к)

г

(1)

sin 8

Здесь р(к) - удельный вес материла. Все константы и функции будем обозначать верхними индексами 1 и 2 для внутреннего и вешнего цилиндров соответственно. Определяющие соотношения

а

(к) = K(k )Jk) + K(к )F(k)

— к11 гг + Kl "> 500

а

(к )= к (к к )■

12 гг

а

(к )= E( к)

(Лк )а(к)

ГТ гг

7(к)

12

- K(к)£(к)

K22 588

ц(к )а(к)^

Нz8 а88

Р(к)

т(к )= С(к Пк)

тг8 = 4-8 1 г8

(2)

для тяжелого составного ортотропного осесиммет-ричного тела, находящегося в условиях плоской деформации, содержат коэффициенты

m Е(к) Лк ) _ Ег.

D У

)_ Ег Л (к)..(к)\

K11 ^TÖT Ц _H8z Н z8 )

т Е(к)

(к) _ Ег

Dк

(к)

K(к) = к12 =

Dк)

Е

D

(н((8) + Н(гк )Н((к8))

K (к ) = Е^(1 _Н(к ук))

K22 - (к) (1 Нгг Нгг )

D(к) = 1 _2^«' _мГk9)н9í) _H8'Z,HZS) ^W?

определяемые модулями Юнга Ег(к), Ее(к) и Е/к) в направлениях г,е и г, коэффициентами Пуассона

цге , и ц[кг^, поперечным модулем сдвига

(к) г8 .

^8.

Последовательная подстановка геометрических соотношений Коши

5(к ) =

ди(к) дг

5(к ) = ь88 =

) 1 ду(к)

■ + -

г г д8

(к ) = 1 v^

1г8 г д8 дг г

(3)

и(2 >1 = v( 2 )|

= 0

г=b

а

(1)1

= 0

г(1)

г8

= 0

г=a

и(1) I ' = и (2)| v(1 )| = v(2 )|

1г=С 1г=С г=С 1г

а^ )| гг 1г=С = а(2 )| ' гг 1г=С т(1 )| тг 8 1г=С = т(2)| = тг 8 1г

(4)

не нарушают симметрию задачи, что позволяет применить метод, основанный на разложении компонент вектора перемещений по окружной и радиальной координате в тригонометрические и обобщенные степенные ряды [11-14]. Разложение по окружной координате в тригонометрические ряды

в определяющие (2), а затем полученного результата в уравнения равновесия (1) позволяет записать неоднородную систему дифференциальных уравнений Ламе в частных производных, которые не приводятся ввиду громоздкости. Граничные условия на жестко закрепленной внешней (т=Ъ), свободной от нагрузок внутренней (г=а) и контактной (г=с) поверхностях

(к)_„( к)

и - и.

+ и(к W 8 + Х"

2 и

(к),

cos i 8

„(к) ).

= Vf W 8 + Х;

■2 V

,(к ).

sin i 8

позволяет понизить размерность задачи и перейти к последовательности систем обыкновенных дифференциальных уравнений, однородных при ¡ф1. Вид граничных условий (3) позволяет сделать вывод о том, что все системы при имеют единственное тривиальное решение: пг- —0 и у/к)=0, а совместное использование метода вариации произвольных постоянных и разложения вектора

г

перемещений по радиальной координате в обобщенные степенные ряды позволяют получить точное аналитическое решение задачи о равновесии

толстостенного тяжелого составного ортотропного цилиндра с жестко закрепленной внешней поверхностью [13]

^Ца^* + а(к)г"Р* + а(к)1п г + а(к) + Б(к)г2) сов 9 у(к) = (а(к)а(к)А + а(к)а(к)г~Рк - а(к)1п г - а« + Л(к)г'2) в1п 9

(5)

Здесь

а( к ) =

рк(кк2 +1) + р2 +1 р2 - р2 -1

(к) Рк(к'+1)-Рк2 -1

а2 ) =--^-?-

2 р2 - р2 -1 ,

Рк = >/пк + Рк2 - (2 + кк) +1

^к) К*к)- 1)р(к)

ю( к )Ок)

ю3 Ог9

,(к) в(к)р(к)

о

(к) г9

Л к )0 к )+®( к )К(к )

к) = ю3 0г9 + юг9 К11

К к)(»(к)»(к)

•'2 ю4

®1к) = 4 - т2 - п1, ®2к)= 2Чк + тк - П

ю3к)= 3-р2 Ю(4к)= 2к2 + р2 + 3

V

ок)

ог9

К1к)

V

К

(к) 22

К ^

(к) 12 (к)

Рк =

К(к)

К22 и- _

Ог9

К

(к)

21

О

( к) г9

Постоянные интегрирования а1(^), а2(к), а3(к) и а4(к), входящие в общее решение (5) неоднородной системы обыкновенных дифференциальных уравнений при /=1 определяются первым членом разложения граничных условий (3) в тригонометрический ряд по окружной координате и не приводятся ввиду громоздкости. Подстановка уравнений (5) последовательно в геометрические (4) и определяющие (2) соотношения позволяет записать выражения для напряжений. Обратим внимание на то, что в частном случае для изотропного тяжелого горизонтального составного цилиндра новое решение совпадает с результатами, приведенным в работе [11], а в случае, когда можно пренебречь вкладом массовых сил - совпадает с решением

[15].

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

(к)

осевом и радиальном направлении, от сдвигов в поперечных плоскостях.

Будем рассматривать контейнер, представляющий собой жестко закрепленное по внешней поверхности цилиндр с геометрическими параметрами: а=0,3 м, с=0,5 м и 6=0,6 м, состоящий из внутренней полимербетонной (Ё,г(1"|=Ё'е(1)=29,2 ГПа,

£г(1)=27,3 ГПа,

Ое(1)=11,5

К

ГПа, ц(ге = 0,19, (1) = 0,21 и у(1)=25 кН/м3) и внешней стеклопла-

стиковой (Ег -11,8 ГПа, Ее(2)=56,8 ГПа, £г12=15,7 ГПа, 0^С2)=9,8 ГПа, ^ = 0,06, ц,^ = 0,38,

= 0,21 и у(2)=19 кН/м3) составных частей. На рис. 1 и 2 представлены распределения ненулевых компонент тензора напряжений в поперечных сечениях контейнеров, вызванные массовыми силами и внутренним давлением р=0,4 МПа, вдоль меридиональной и обезразмеренной радиальной Р = (г -Р1)/ (Р2 -Р1) координаты.

Вклад массовых сил в распределение напряжений для конструкции контейнера заданной геометрии не превышает 5%. Касательные напряжения, возникающие в поперечных сечениях, определяются исключительно массовыми силами, в результате чего их величина не зависит от заданного равномерного давления (рис. 1). Радиальные напряжения по абсолютной величине убывают от нагруженной давлением внутренней к закрепленной внешней поверхности. Окружные и осевые напряжения имеют скачок на поверхности контакта.

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

к

к

Рис. 1. Распределение напряжений (МПа) на нагруженной давлением р=0,4 МПа внутренней а™, закрепленной внешней аи контактной ас поверхности контейнера

Наибольший вклад в распределение напряжений в поперечных сечениях контейнера, вызванных сочетанием нескольких факторов, вносит контактное давление рС =1,96 МПа (рис. 2). Наибольшие по абсолютной величине значения радиальных, окружных и осевых напряжений достигаются на верхнем и нижнем сводах цилиндрического контейнера и, следовательно, эти точки являются наиболее опасными с точки зрения возможности

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

Рис. 2. Распределение напряжений (МПа) на нагруженной давлением р=0,4 МПа внутренней , закрепленной внешней а^ и контактной аС поверхности контейнера, с учетом контактного давления

Выводы: несмотря на то, что вклад массо- 5%, полученное аналитическое решение позволяет вых сил в напряженное состояние контейнера из учесть еще и неоднородные реакции основания [14], стеклопластика и полимербетона не превышает а также любые другие нагрузки, симметричные

относительно вертикальной диаметральной плоскости, которые можно разложить в тригонометрические ряды.

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

1. Зайцев, А.В. Статистическое описание структуры двухфазных волокнистых композитов / А.В. Зайцев, А.В. Лукин, Н.В. Трефилов // Вестник ПНИПУ. Механика. 2002. № 10. С. 52-62.

2. Зайцев, А.В. Корреляционные функции случайных структур двухфазных волокнистых композитов, синтезированных методом статистических испытаний / А.В. Зайцев, А.В. Лукин, Н.В. Трефилов // Вестник ПГТУ. Аэрокосмическая техника. 2002. Вып. 14. С. 3-7.

3. Зайцев, А.В. Закономерности случайных полей структуры двухфазных композитов с эллиптическими в поперечном сечении волокнами / А.В. Зайцев, М.С. Кос-тоусов, А.В. Лукин // Вестник ПГТУ. Аэрокосмическая техника. 2003. № 15. С. 3-10.

4. Зайцев, А.В. Закономерности скейлинга в случайных полях структуры однонаправленно армированных волокнистых композитов / А.В. Зайцев, А.В. Лукин, Н.В. Трефилов // Физическая мезомеханика. Спецвыпуск. Ч. 1. 2004. Т. 7. С. 58-61.

5. Зайцев, А.В. Случайные структуры двухфазных композитов: синтез, закономерности, новая оценка характерных размеров представительных объемов / А.В. Зайцев, А.В. Лукин, А.А. Ташкинов, Н.В. Трефилов // Вестник ПНИПУ. Механика. 2004. № 12. С. 30-44.

6. Зайцев, А.В. Компьютерный синтез, закономерности случайных структур и моделирование процессов разрушения волокнистых композитов при продольном сдвиге / А.В. Зайцев, А.В. Лукин, Н.В. Трефилов // Физическая мезомеханика. 2004. Т. 7, № 5. С. 73-79.

7. Зайцев, А.В. Новый метод построения моментных функций второго порядка случайной структуры по-лимербетонов / А.В. Зайцев, Я.К. Покатаев // Вестник ПНИПУ. Механика. 2007. № 15. С. 28-45.

8. Зайцев, А.В. Общие закономерности структуры, случайных полей напряжений и деформаций в волокнистых и дисперсно--упрочненных композитах / А.В. Зайцев, А.В. Кислицын, В.С. Кокшаров // Вестник ННГУ им. Н.И. Лобачевского. 2011. № 4. Ч. 4. С. 1485-1487.

STRESS DISTRIBUTION IN A CROSS SECTION OF CONTAINER FOR LONG-TERM STORAGE OF HIGH-CORROSIVE MEDIA

© 2012 A.V. Zaitsev, A.V. Kislitsyn, A.V. Kutergin, A.A. Fukalov Perm National Research Polytechnical University

On the basis of the predicted effective elastic moduli for polymer concrete and obtained new exact analytical solutions to problem on equilibrium state of thick-walled heavy combined horizontal anisotropic cylinder, which is fixed on external lateral surface, the stress states of containers for long-term storage of high-corrosive solid and liquid media were investigated. External cylinders of containers were usually produced by fiber-glass plastics and coaxial internal one by polymer concrete. Influence of gravity forces and contact pressure on distribution of stresses in cross-sections is analysed, and an initial strength analysis of containers was cared out on the basis of multicriteria approach.

Key words: container, fiber-glass, heavy combined horizontal anisotropic cylinders, composite, mineral aggregates, multiscale model, analytical solutions, multicriteria analysis

Aleksey Zaitsev, Candidate of Physics and Mathematics, Associate Professor at the Department "Mechanics of Composite Materials and Constructions". E-mail: [email protected]; Anton Kislitsyn, Post-graduate Student; Aleksey Kutergin, Master Student; Anton Fukalov, Post-graduate Student

9. Зайцев, А.В. Эффективные упругие модули двухфазных композитов на основе термореактивных смол с минеральным наполнителем / А.В. Зайцев, А.В. Новго-родова, Д.И. Федоров // Вестник ПГТУ. Аэрокосмическая техника. 2005. № 25. С. 62-68.

10. Зайцев, А.В. Эффективные модули объемного сжатия дисперсно--упрочненных композитов со сплошными и полыми анизотропными сферическими включениями / А.В. Зайцев, А.А. Фукалов // Вестник ПНИПУ. Механика. 2010. №4. С. 46-54.

11. Кузнецов, Г.Б. Упругость, вязкоупругость и длительная прочность цилиндрических и сферических тел. -М.: Наука, 1979. 112 с.

12. Зайцев, А.В. Упругое равновесие тяжелой трансвер-сально-изотропной толстостенной сферы с жестко закрепленной внутренней поверхностью / А.В. Зайцев, А.А. Фукалов // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2010. №5(21). С. 85-95.

13. Зайцев, А.В. Упругое равновесие тяжелого горизонтального толстостенного ортотропного цилиндра, находящегося под действием неравномерно распределенного бокового давления / А.В. Зайцев, А.В. Кутер-гин // Вестник ПНИПУ. Механика. 2010. №4. С. 36-45.

14. Фукалов, A.A. Точные аналитические решения задач о равновесии упругих анизотропных тяжелых тел с центральной и осевой симметрией и их приложения / А.А. Фукалов, А.В. Кутергин // Вестник ННГУ им. Н.И. Лобачевского. 2011. № 4. Ч. 4. С. 1831-1833.

15. Зайцев, А.В. Об одном решении задачи Ламе для составного протяженного элемента конструкции, состоящего из посаженных с натягом толстостенного трансверсально-изотропного внешнего цилиндра на соосный изотропный внутренний / А.В. Зайцев, А.В. Кислицын // Вестник СамГТУ. Физико-математические науки. 2007. № 1(14). С. 164-167.

16. Вильдеман, В.Э. Механика неупругого деформирования и разрушения композиционных материалов / В.Э. Вильдеман, Ю.В. Соколкин, А.А. Ташкинов. - М.: Наука, 1997. 288 с.

17. Цыплаков, О.Г. Научные основы технологии композиционно-волокнистых материалов. Ч. 1. - Пермь: Пермское кн. изд-во, 1974. 316 с.