open all | close all

3.2.2.1. Пример расчета цилиндрических оболочек

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

Параметры оболочек

D1=20cm, H1=15cm, t1=2mm; D2=20cm, H2=15cm, t2=1mm; E=2.1 106kg/cm2.

При генерации оболочек размер конечного элемента выбирался таким образом, чтобы на одной длине полуволны теоретического решения задачи располагалось не менее 3х-4х конечных элементов. В связи с этим, число делений по окружности первой оболочки (t=2mm) было принято 96. второй (t=1mm) - 128.

Оболочка в недеформированном состоянии.

Форма потери устойчивости "идеальной" оболочки (на примере оболочки t=2mm).

Форма потери устойчивости для цилиндрической оболочки всегда получалась несимметричной относительно диаметральной плоскости. Возможно это связано с тем, что первое собственное значение всегда являлось кратным (попытка построить сопряженную форму потери устойчивости не предпринималась).

Результаты расчета

Статический расчет
Главные напряжения s1 оболочки t=2mm с погибью глубиной 0,5mm от нагрузки 0.5 Pкр идеальной оболочки

Расчет по деформированной схеме
Главные напряжения s1 оболочки t=2mm с погибью глубиной 0,5mm от нагрузки 0.5 Pкр идеальной оболочки

Ниже показаны графики изменения критического параметра для обеих оболочек в зависимости от относительной глубины погиби (глубина погиби отнесена к толщине оболочки). Значение критического параметра 1 - соответствует теоретическому значению линейной теории - Pкр = 0,605E ( t / R ). (Поиск критического параметра во всех случаях производился с точностью до 1%. Точность установки продольного напряженного состояния в нелинейных расчетах - до 1%).

Относительная глубина погиби