Расчет светопрозрачных элементов заполнения проемов из поликарбоната
Божанов Павел Валерьевич
кандидат технических наук, генеральный директор ООО «Инженерный центр промышленного проектирования», [email protected]
Широкое применение конструктивных элементов из поликарбоната требуют теоретически обоснованных методов расчета, учитывающих дилатационные свойства рассматриваемого материала. С одной стороны, склонность к дилата-ции у поликарбоната особенно ярко проявляется при высоком уровне напряжений, при котором происходит образование пластических деформаций, а с другой стороны, зависит от вида напряженного состояния. В связи с этим, особенную практическую ценность для прикладных инженерных расчетов имеют решения задач по определению предельного состояния конструкций при сложных видах напряженного состояния, ярким примером которого может служить изгиб тонких пластинок. Именно к указанным типам конструкций относятся элементы светопрозрачного заполнения проемов. В настоящей статье решена задача по определению предельного состояния тонких квадратных пластин, находящихся под воздействием поперечной нагрузки. Были рассмотрены два вида опирания плиты: подвижное защемление контура и свободное опирание. Расчет проводился по двум вариантам расчета. В первом варианте при развитии пластических деформаций учитывались дилатационные свойства поликарбоната. Во втором случае использовались классические подходы, не учитывающие специфических свойств рассматриваемого материала. Получены картины развития пластичности по поверхности пластин и распространение по глубине сечения при увеличении уровня нагрузки. Построены графики зависимости прогибов центральной точки от величины интенсивности поперечной нагрузки. Проведено сравнение полученных результатов и установлено, что игнорирование склонности поликарбоната к дилатации и пластическому разрыхлению приводит к значительной разнице в значениях нагрузок, соответствующих предельному состоянию и не может быть допустимым обстоятельством при выполнении инженерных расчетов конструкций, выполненных из рассматриваемого материала.
Ключевые слова: поликарбонат, пластины, предельное состояние, изгиб, дилатирующие материалы, полимеры
Светопрозрачное заполнение проемов, с расчетной точки зрения, выполняется в виде тонких пластин. Как известно, несущая способность пластин не ограничивается переходом материала из упругого состояния в пластическое в какой - либо отдельно взятой точке, а определяется исходя из значения внешнего усилия, при котором возможно образование пластического шарнира. Таким образом, задача по определению несущей способности рассматриваемых элементов должна ставиться с учетом развития упруго - пластических деформаций. Экспериментальные данные [1] свидетельствуют о том, что образование пластических деформаций в поликарбонате зависит от вида напряженного состояния.
Для определения исходной точки, при которой в материале начинают возникать пластические деформации, воспользуемся подходом, изложенным в работах [3,4]. Условие пластичности для дилатирующих разносопротивляю-щихся материалов предлагается сформулировать в достаточно общем для указанного класса материалов виде
^(о у) = т/(£)=кт (1)
где /(£) - функция вида напряженного состояния;
£ = о / £ о - октаэдрическое нормальное нормированное напряжение;
о = о у ■ 5у /3 - среднее напряжение,
где 5у - символ Кронекера;
8 у = о у - 5у ■ о - девиатор напряжений;
т = ^ 8 у ■ 8 у / 3 - касательное напряжение на октаэдрической площадке;
.2
8о =л/о2 + Т - модуль вектора полного напряжения на октаэдрической площадке;
кт = V2/3^т£, где т$ - предел текучести при чистом сдвиге.
X X
о
го А с.
X
го т
о
ю 3
м о
О)
о
см
О!
О Ш
т
X
3
<
т О X X
(2)
Функция /(4), входящая в выражение (1) определяется из экспериментальных зависимостей деформационных характеристик поведения материала, полученных при сложных видах напряженного состояния. В работе [5] была проведена обработка экспериментальных диаграмм напряженно - деформированного состояния образцов из поликарбоната, полученная при сдвиге и растяжении под различными значениями гидростатического давления и предложена аппроксимация функции вида напряженного состояния в форме (2). Значение константы кт принималось равным 33,0 МПа.
/(4) = 1 + 0,3774, при -1 <4< 0] / (4) = 1 - 0,3394, при 0 <4< 1 Обозначим основные гипотезы и предпосылки, принятые при решении изгиба пластин из поликарбоната за пределом упругости при конечных прогибах:
Используются обычные положения технической теории изгиба тонких пластин - гипотеза плоских сечений и гипотеза плоского напряженного состояния.
Диаграмма напряжений - деформаций материала обладает ярко выраженной площадкой текучести, что позволяет применить модель идеального упруго - пластического тела. Нагружение считается простым. Задача по изгибу пластин решается в нелинейной постановке.
Задача по исследованию напряженно - деформирования состояния пластин в нелинейной постановке сводится к решению системы двух нелинейных дифференциальных уравнений [6]:
д ф дф „ - + —+ 2-
дх1
д2т дх12
дх2 д® дх22
+ 2
д 4ф дх^дх^
д2т дх1дх2
= Е
д 2т
д т д т
дх1дх2) дх1 дх2
(3)
к | ц + д2ф д2т + д2ф д2т 2 д2ф д2т Б1 к дх2 дх2 дх2 дх^ дх1дх2 дх1дх2
где Е - модуль упругости материала, д - интенсивность распределенной нагрузки, Л - толщина пластины,
,3
Б =-
Ек
3(1 -И2)
И - коэффициент Пуассона.
Первое выражение в системе (3) является уравнением совместности или неразрывности деформаций, записанное через функции напряжений ф(х^, х2) и прогибовт(х^, х2) , второе - является уравнением равновесия, записанное
цилиндрическая жесткость,
через изгибающие моменты Му , а также функции напряжений и прогибов пластины. Причем функция напряжений непосредственно связана с продольными усилиями Ыу следующими зависимостями:
Ып = к ^
11 дх2
N22 = к ЦФ
Ы„ = -к
д 2Ф
дхлдх2
Система нелинейных уравнений (3) определяет работу пластины в упругой стадии работы материала.
Рассмотрим состояние пластины, когда в материале появляются пластические деформации. Очевидно, что в начале подобного процесса, пластичность будет возникать только с одной стороны пластины, либо в нижней, либо в верхних зонах, как показано на рис. 1.а. Данное утверждение основывается на двух обстоятельствах. Во - первых, при больших прогибах в серединной плоскости возникает растягивающее продольное усилие, которое с одной стороны увеличивает растягивающие напряжения от изгиба и соответственно уменьшает сжимающие изгибные напряжения. Во - вторых, введение функции вида напряженного состояния приводит к зависимости появления пластичности в верхней и нижней зонах от соотношения компонент напряжений в указанных зонах, которые будут различаться между собой, что приведет, в конечном счете, к разнице пределов текучести для одной и той же точки на плоскости пластины, но с разными координатами по её толщине. Следует обратить особое внимание, что последнее замечание подтверждается экспериментальными данными, которые были приведены в начале статьи.
С увеличением уровня нагрузки пластичность будет развиваться по глубине пластины и, в определенный момент, компоненты напряженного состояния могут достигнуть величин, при которых образуется текучесть в противоположной зоне сечения от первоначально зафиксированного положения пластичности в рассматриваемом сечении. В таком случае, в сечении пластины будет наблюдаться двусторонняя пластичность, показанная на рис. 1.б. На рис.1 введены следующие обозначения: Ау - напряжения вызывающие текучесть в нижней зоне; Ву - напряжения вызывающие текучесть в верхней зоне; ау , Ьу- координаты зон пластичности отсчитываемые от серединной поверхности по
2
толщине пластины; с у - координата нейтральной оси сечения, в которой напряжения принимают нулевые значения и отсчитываемая от серединной поверхности по толщине пластины. Следует отметить, что зоны распространения пластичности по толщине пластины ау, Ьу
одинаковы для каждого определенного сечения, поскольку они формируются комплексным соотношением напряжений о у, а не каждым отдельно взятым напряжением. Вследствие этого будет удобно опустить в дальнейших математических выкладках индексы для указанных величин, приняв для них общие обозначения а, Ь. С другой стороны, зоны пластичности являются функциями координат поверхности пластины, поскольку определяются в конкретной точке.
Очевидно, что с ростом нагрузки будет происходить ослабление поперечного сечения пластин за счет развития пластичности, при принятой модели идеального упруго - пластического тела, и соответственно будет уменьшаться значение цилиндрической жесткости пластины. При принятых ранее обозначениях и, в соответствие с рис.1, значение цилиндрической жесткости, ослабленного пластичностю сечения пластины -8 , определиться, как показано в работе [7] следующим выражением:
8 =
Е(а3 -Ь3) 3(1 V)
(4)
Поскольку зоны пластичности, помимо зависимости от уровня нагрузки, также связаны с координатами расположения точки на поверхности пластины, то цилиндрическая жесткость пластины, ослабленная текучестью, тоже будет являться функцией системы координат поверхности пластины. Определение аналитической формы записи рассматриваемой функции, если и возможно чисто теоретически, то практического смысла не имеет, поскольку систему нелинейных дифференцированных уравнений (3)
предлагается вначале линеаризовывать численным методом последовательных нагруже-ний, как рекомендовано в работе [8], а затем решать линеаризованную систему методом конечных разностей. При описанном ранее подходе, значение функции 8 будет определено в каждом узле сетки пластины численно, что позволит определить интегральные характеристики напряженно - деформированного сечения пластины, без поиска аналитических определений рассматриваемой функции.
Система двух нелинейных дифференциальных уравнений (3), с учетом сделанных замечаний и обозначений, принятых в формуле (4), после несложных преобразований, в соответствие с работой [7] запишется в следующей форме:
д 4р д4р д V - + —^ + 2- ^
дх1 дх2 дх дх2
-Е
( д2,
д 2® д 2®
\удх1дх2) дх12 дх2
„ , д4® д4® д4® 8-I—г +-т + 2
дх1 дх2
дх^дх^
д V д2® д V д2® „
= -I д + —— + - 2-
д2р д2®
(5)
дх22 дх12 дх12 дх22 дх1дх2 дх1дх2
Линеаризация системы двух нелинейных дифференциальных уравнений (5), как указывалось ранее, проводилась методом последовательных нагружений, как рекомендовано в работе [8]. Линеаризованная система уравнений (5) будет иметь следующую форму записи:
V 46р+У2Х 6® = 0
У^ 6р+ 8 - V46® - V2д6® = -5д где
У4 =-Ё1+Х. + 2- д4
(6)
дх1
дх24 дх12дх^
у2
УХ = Х22 2
дх1
д2 2 ' Х11 ~ 2 + 2Х12 дх2
дх1дх2
д2 д2 д2
V д = + + 2 N12 —
дх1 дх2
дх1 дх2
6р, 6® - приращения функции напряжений и прогибов, соответствующих малому увеличению внешней нагрузки 6 д;
8, N22^ Ни, Хц, Х22, Х12 - начальные характеристики системы, соответствующие определенному уровню внешней нагрузки.
Решение системы линеаризованных уравнений (6) проводилось численным методом конечных разностей [9].
Расчет проводился по двум вариантам: в первом варианте в качестве условия предельного состояния использовалось уравнение (1), а аппроксимация функции вида напряженного со-
X X
о
го А с.
X
го т
о
ю 3
ю о
<е
®
2
д
О)
о
сч
fO OI
О Ш
m
X
3
<
m О X X
стояния принималась по выражению (2); во втором варианте расчета в качестве условия предельного состояния применялась классическое уравнение Треска - Мизеса - Генки [10], не учитывающее зависимость деформационных характеристик от вида напряженного состояния. Условие пластичности (2) легко приводится к классическому условию Треска - Мизеса - Генки [10], при f (4) = 1 (const).
Расчет выполнялся для двух видов опира-ния: в первом случае принималась подвижная заделка контура пластины; во втором принималось шарнирно подвижное опирание контура. Пластина принималась квадратной в плане, толщиной h равной 4,0 см, длина и ширина пластины - 80,0 см, модуль упругости - 2485 МПа, коэффициент поперечной деформации - 0.39. Поверхность пластин была покрыта сеткой 20x20. В силу симметрии рассчитывалась четвертая часть пластины.
Результаты расчета квадратной пластины с подвижным защемлением контура. Пластичность, при расчете пластины для первого варианта расчета возникает впервые в верхней зоне пластины в точке Х (рис.2.а). Нагрузка, соответствующая появлению пластичности составила 1.03 МПа. Предельная нагрузка, соответствующая образованию пластического шарнира, составила 8.14 МПа. Полученные картины распространения пластичности по верхней поверхности пластины при нагрузке 3.2 МПа приведены на рис. 2.а, по нижней поверхности на рис. 2.б. На рис. 2.в и 2.г. приведены аналогичные картины распространения пластичности при нагрузке 5.0 МПа. На рис. 3.а и 3.б приведены картины распространения пластичности по глубине пластины поперечного сечения при указанных ранее значениях нагрузки соответственно.
Пластичность, при расчете пластины для второго варианта расчета (при f (4) = 1 ) возникает в верхней зоне пластины в точке Х. Нагрузка, соответствующая появлению текучести составила 0.70 МПа. Предельная нагрузка, соответствующая образованию пластического шарнира, составила 4.85 МПа.
Сравнительная диаграмма прогибов центральной точки приведена на рис. 4, где сплошной линией показаны прогибы полученные при расчете пластины по условию (1), пунктиром -при расчете по классическому варианту расчета.
Результаты расчета квадратной пластины с шарнирно подвижным опиранием контура.
Пластичность, при расчете пластины для первого варианта расчета возникает впервые в нижней зоне пластины при нагрузке 0.76 МПа. Предельная нагрузка, соответствующая образованию пластического шарнира, составила 6.56
МПа. Полученные картины распространения пластичности по нижней поверхности пластины при нагрузках 1.8 МПа и 3,5 МПа приведены на рис. 5.а, и 5.б. соответственно. На рис. 6.а и 6.б приведены картины распространения пластичности по глубине пластины поперечного сечения при указанных ранее значениях нагрузки.
Пластичность, при расчете пластины для второго варианта расчета (при /(4) = 1 ) возникает в нижней зоне пластины при нагрузке 0.72 МПа. Предельная нагрузка, соответствующая образованию пластического шарнира, составила 3.64 МПа.
Сравнительная диаграмма прогибов центральной точки приведена на рис. 7, где сплошной линией показаны прогибы полученные при расчете пластины по условию (1), пунктиром -при расчете по классическому варианту расчета.
Рис. 5
Проводя сравнительный анализ, выполненный по полученным результатам расчетов, можно сделать следующие выводы:
Пластические деформации при расчете по условию (1) возникают в пластине при больших нагрузках, чем при расчете по классическому критерию Треска - Мизеса - Генки [10]. Для пластины с подвижным защемлением контура такое отличие составляет 32%, для пластины с шарнирно подвижным опиранием - 6%.
Пластический шарнир при расчете при расчете по условию (1) возникают в пластине при больших нагрузках, чем при расчете по классическому критерию Треска - Мизеса - Генки [10]. Для пластины с подвижным защемлением контура нагрузки различаются на 40%, для пластины с шарнирно подвижным опиранием - на 45%.
Литература
1. Айнбиндер С. Б., Алксне К. И., Тюнина Э.Л., Лака М. Г. Свойства полимеров при высоких давлениях. - М.: Химия, 1973. - 192с.
2. Silano A.A., Bhateja S.K., Pae K.D. Effects of Hydrostatic Pressure on the Mechanical Behavior of Polymers: Polyurethane, Polyoxymethylene, and
Branched Polyethylene //International Journal of Polymeric Materials. - 1974, vol. 3, N 2, p. 117-131.
3. Трещёв А.А. Зависимость предельного состояний конструкционных материалов от вида напряженного состояния //Известия Вузов. -Строительство. - №10. - 1999. - С. 13 - 18.
4. Трещев А.А., Божанов П.В. Об ограниченности применения классических подходов при исследовании пластического изгиба пластин из дилатирующих материалов// Проблемы машиностроения и автоматизации. - 2004. - №2. - С. 53 - 57.
5. Божанов П.В., Захарченко В.А., Трещев А.А. Упруго - пластический изгиб прямоугольных пластин из изотропных дилатирующих материалов при больших прогибах.// Строительная механика и расчет сооружений. - 2008. - №6. -С. 55 - 60.
6. Тимошенко С.П., Войновский - Кригер С. Пластинки и оболочки. - М.: Физматгиз, 1963. -647 с.
7. Божанов П.В., Трещев А.А. Определение прочностных критериев при возникновении пластических деформаций в поликарбонате// Инновации и инвестиции. - 2018. - №12. - С. 323 -326.
8. Петров В.В. Метод последовательных на-гружений в нелинейной теории пластин и оболочек. - Саратов: СГУ. 1975. - 119 с.
9. Варвак П.М. Развитие и приложение метода сеток к расчету пластинок. К.: Изд-во АН УССР.- 1957. - 339 с.
10. Гольденблат И. И., Копнов В. А. Критерии прочности и пластичности конструкционных материалов. - М.: Машиностроение, 1968. -191 с.
Calculation of translucent elements for filling openings
made of polycarbonate Bozhanov P.V.
LLC «Engineering Center for Industrial Design» The widespread use of polycarbonate structural elements requires theoretically grounded calculation methods that take into account the dilatation properties of the material in question. On the one hand, the tendency to dilation in polycarbonate is especially pronounced with a high level of stress, at which plastic deformations occur, and on the other hand, it depends on the type of stress state. In this regard, of particular practical value for applied engineering calculations are the solutions aimed to determine the limiting state of structures for complex types of stress state, an example of which is the bending of thin plates. These types of structures include elements for translucent openings filling. This article solves the problem of determining the limiting state of thin square plates under the influence of a transverse load. Two types of plate support have been considered: moving clamping of the contour and free bearing. The calculation was carried in two variants. In the first variant, in the development of plastic deformations, the dilatation properties of polycarbonate were taken into account. In the second case, classic approaches were used that did not consider the specific properties of the material in question. The picture of the development of plasticity over the surface of the plates and the distribution along the depth of the section with increasing load level has obtained. The graphs of the deflection of the central point versus the intensity of
X X О го А С.
X
го m
о
ю 3
м о
to
the transverse load have been plotted. The obtained results were compared and it was found that ignoring the tendency of polycarbonate to dilation and plastic loosening leads to a significant difference in the values of the loads corresponding to the limit state and cannot be a permissible circumstance when performing engineering calculations of structures made of the material under consideration. Keywords: polycarbonate, plates, limit state, bending, dilating
materials, polymers References
1. Ainbinder S. B., Alksne K. I., Tyunina E. L., Laka M. G. Prop-
erties of polymers at high pressures. - M .: Chemistry, 1973. - 192s.
2. Silano A.A., Bhateja S.K., Pae K.D. Hydrostatic Effects of
Polymers: Polyurethane, Polyoxymethylene, and Branched Polyethylene // International Journal of Polymeric Materials. -1974, vol. 3, N 2, p. 117-131.
3. Treschev A.A. Dependence of limiting states of structural
materials on the type of stress state // News of universities. -Building. - №10. - 1999. - pp. 13-18.
4. Treschev A.A., Bozhanov P.V. On the limitations of the use of
classical approaches in the study of plastic bending plates from dilating materials // Problems of mechanical engineering and automation. - 2004. - №2. - pp. 53 - 57.
5. Bozhanov P.V., Zakharchenko V.A., Treshchev A.A. Elastic -
plastic bending of rectangular plates of isotropic dilating materials with large deflections. / / Structural mechanics and calculation of structures. - 2008. - №6. - S. 55 - 60.
6. Timoshenko S.P., Voinovsky - Krieger S. Records and shells.
- M .: Fizmatgiz, 1963. -647 p.
7. Bozhanov P.V., Treshchev A.A. Determination of strength criteria for the occurrence of plastic deformations in polycarbonate // Innovations and investments. - 2018. - №12. - p. 323 - 326.
8. Petrov V.V. The method of successive loads in the nonlinear
theory of plates and shells. - Saratov: SSU. 1975. - 119 p.
9. Varvak P.M. Development and application of the grid method
to the calculation of plastinoc. K .: Publishing House of the Academy of Sciences of the Ukrainian SSR.- 1957. - 339 p.
10. Goldenblat I.I., Kopnov V.A. Criteria of strength and plasticity of structural materials. - M .: Mashinostroenie, 1968. -191 p.