Размещено на http://www.allbest.ru/

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

А.Ш. Кусяков

Аннотация

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

Ключевые слова: композит; пластина; ребра жесткости; проектирование.

© Кусяков А.Ш., 2013Рассматривается подкрепленная композитная пластинка длиной a и шириной b, находящаяся под действием сжимающих нагрузок (рис. 1).

Размещено на http://www.allbest.ru/

Рис. 1. Подкрепленная пластинка под действием сжимающих нагрузок

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

Введем следующие обозначения:

- толщина многослойной пластинки;

и - относительные содержания продольных и поперечных слоев соответственно;

- шаг, высота, площадь поперечного сечения и момент инерции продольных ребер соответственно;

- шаг, высота, площадь поперечного сечения и момент инерции поперечных ребер соответственно;

- плотности материалов пластинки, продольных и поперечных ребер соответственно;

- модули упругости вдоль и поперек волокон материала пластинки соответственно;

- коэффициенты Пуассона материала пластинки;

- модуль сдвига материала пластинки;

- предел прочности материала слоя полотна пластины при сжатии вдоль волокон;

- модули упругости материалов продольных и поперечных ребер;

- предел прочности материала продольных ребер при сжатии.

В наиболее общем виде задача оптимального проектирования подкрепленной пластинки формулируется следующим образом: найти минимум массы подкрепленной пластины при условии сохранения устойчивости и прочности.

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

, (1)

Где

(2)

Здесь - прогиб пластины; - координаты (рис. 2); - заданная сжимающая нагрузка; - компоненты матрицы изгибной жесткости (i, j = 1, 2, 6).

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

, (i, j = 1, 2, 6). (3)

Здесь , , и - компоненты матрицы жесткости многослойного пакета:

(4)

Компоненты матрицы жесткости слоя в главных осях , , и связанны с упругими постоянными следующими зависимостями:

(5)

Для шарнирно опертой по всему контуру пластинки решение основного разрешающего уравнения можно представить в виде

, (6)

Где

, . (7)

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

. (8)

Критическая нагрузка определяется путем минимизации выражения (8) по параметрам волнообразования (7). Очевидно, что в выражении (8) надо положить , то есть

.

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

(9)

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

 (10)

Отсюда найдем:

. (11)

Подставив выражение (11) в формулу (9), получим выражение для критической нагрузки [2]

. (12)

Известно, что формула вида (12) может быть использована для исследования как общей, так и местной форм потери устойчивости. В случае общей формы потери устойчивости изгибные жесткости ортотропной пластины следует заменить на "приведенные" изгибные жесткости подкрепленной пластины:

;

, (13)

Где

; (14)

. (15)

Здесь и - расстояния между срединной поверхностью пластины и центрами тяжести продольных и поперечных ребер соответственно (эксцентриситеты).

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

. (16)

Здесь - редукционный коэффициент, вычисляемый по формуле

. (17)

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

, , (18)

Где

, . (19)

Здесь и - распределенные нагрузки, воспринимаемые полотном пластинки и продольными ребрами соответственно.

В качестве параметров оптимизации подкрепленной пластинки следующие величины:

- толщина многослойной пластинки;

- относительное содержание продольных слоев;

- относительное содержание поперечных слоев;

- условная "толщина" продольных ребер

;

- высота продольных ребер;

- шаг продольных ребер;

- условная "толщина" поперечных ребер

;

- высота поперечных ребер;

- шаг поперечных ребер.

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

(20)

при наличии ограничений

; (21)

, ; (22)

, ; (23)

, . (24)

Здесь равенство (21) - структурное ограничение; неравенства (22) - ограничения на высоты продольных и поперечных ребер ( и - максимально допустимые по условиям технологии значения высот продольных и поперечных ребер); неравенства (23) - ограничения по общей и местной формам потери устойчивости; неравенства (24) - ограничения по прочности для полотна пластинки и продольных ребер. Для решения данной задачи можно воспользоваться любым известным методом нелинейного математического программирования [5].

В результате многочисленных расчетов было установлено, что независимо от величины сжимающей нагрузки

1) происходит полное вырождение поперечных слоев полотна пластинки и поперечных ребер;

2) всегда активны ограничения по прочности.

На основании полученных результатов можно рекомендовать следующий алгоритм проектирования подкрепленной пластины (в предположении, что ребра и полотно пластины изготовлены из одного и того же материала):

1. Полагаем

, , ,

, , . (25)

Таким образом, остаются только три искомых параметра: - толщина многослойной пластинки, - условная "толщина" продольных ребер, - шаг продольных ребер.

2. Из условия прочности вычисляем полную условную "толщину" подкрепленной пластины

. (26)

3. Полагая , находим критическую нагрузку для гладкой пластинки . Проверяем выполнение условия устойчивости

. (27)

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

4. Если условие устойчивости, полученное на предыдущем шаге, не выполняется, решаем уравнение относительно безразмерного параметра

, (28)

Где

. (29)

Здесь - критическая нагрузка (12), соответствующая общей форме потери устойчивости подкрепленной пластинки. Изгибные жесткости подкрепленной пластины, входящие в выражение для критической нагрузки, выражаются через величины и следующим образом:

,

, ,

(30)

Для решения уравнения (28) можно воспользоваться, например, методом половинного деления.

5. По найденным значениям параметра и полной условной "толщины" вычислим величины и

, . (31)

6. Шаг продольных ребер находим из условия сохранения местной устойчивости

. (32)

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

, ,

, ; (33)

.

Подставив выражения (33) и (34) в уравнение (33), после несложных преобразований получим явное выражение для нахождения шага продольных ребер:

. (35)

После получения оптимального варианта конструкции рекомендуется произвести уточненный прямой расчет, используя какую- либо доступную систему инженерного анализа, например систему ANSYS [1, 3]. Библиотека конечных элементов этой системы содержит, в частности, оболочечный элемент SHELL63 и балочный элемент BEAM3, которые могут быть использованы для расчета подкрепленных пластин. алгоритм композитный пластина

Список литературы

1. Басов К.А. ANSYS: Справочник пользователя. М.: ДМК Пресс, 2005. 640 с.

2. Вольмир А.С. Устойчивость упругих систем. М.: Физматгиз, 1963. 880 с.

3. Кусяков А.Ш. Конечно-элементное моделирование в среде ANSYS / Перм. гос. ун-т. Пермь, 2007. 150 с.

4. Тетерс Г.А., Рикардс Р.Б., Нарусберг В.Л. Оптимизация оболочек из слоистых композитов. Рига: Зинатне, 1978. 240 с.

5. Химмельблау Д. Прикладное нелинейное программирование. М.: Мир, 1975. 534 с.

Размещено на Allbest.ru