Simulation of the orthotropic plate response behavior during dynamic load.pdf Рассматривается задача о напряженно-деформированном состоянии (НДС) элемента конструкции, представляющего собой пластину, закрепленную по ее контуру. Точки (или линии) закрепления находятся на жестком основании, причем это основание получает импульс ускорения известной величины и формы в направлении, перпендикулярном плоскости пластины. Если принять, что пластина является абсолютно жесткой (известная из теоретической механики модель абсолютно твердого тела), то все ее точки будут иметь те же ускорения, скорости и перемещения (параметры движения), что и заданные на контуре величины. В действительности за счет деформирования пластины в разных ее точках эти параметры движения будут различными. Интерес представляют их экстремальные значения, в частности параметры ускорения, поскольку именно они определяют массовые силы (силы инерции) как в самой пластине, так и в любых элементах, связанных с ее поверхностью. Эти силы могут вызвать нарушение целостности пластины или элементов. НДС ортотропной пластины при действии динамической нагрузки предлагается оценивать с использованием средств вычислительной механики на основе соотношений теории пластин, т.е. использовать для анализа конструкции математическое моделирование в сочетании с численной реализацией. Изотропные и ортотропные модели пластин Рассматривается реакция пластины на ударную нагрузку [1, 2], вызывающую ускорение пластины вдоль нормали к ней. Закон изменения нагрузки во времени задается кусочно-линейной функцией, в частности, в виде трапеции. Такой выбор 112 И.В. Щербаков, Б.А. Люкшин определяется тем, что при ограниченном числе параметров, определяющих форму импульса, можно получить хорошую аппроксимацию реальной нагрузки. Выбирая систему координат, связанную с жестким основанием, можно исключить из рассмотрения общую для всех точек пластины постоянную составляющую перемещений. Это важно при численной реализации с использованием конечно-разностных соотношений, так как перемещения, определяющие деформации, много меньше этой составляющей, а разности относительно близких величин считаются с большими ошибками за счет округления. Поскольку эта система в силу постановки задачи не является инерциальной (закон движения основания - и соответственно начала координат - задается в виде импульса ускорения), в уравнения движения необходимо ввести массовые силы, отражающие характер системы отсчета. Для случая, когда материал пластины является изотропным, в простейшем линейном случае уравнение движения можно записать в виде [3, 4] D V4 w +ρd-W + ρG(t) = 0, (1) h ∂t2 D Eh^ Е где D = - так называемая цилиндрическая жесткость пластины; Е - модуль упругости материала пластины, h - ее толщина; w(-,y,t) - прогибы пластины, т.е. перемещения вдоль нормали; -, y - декартовы координаты вдоль сторон пластины; t - время; G - заданное ускорение, тогда ρG(t) - массовые силы, возникающие 4 ∂4∂4∂4 за счет выбора системы отсчета; V4 =-- + 2-;-- +---- дифференциальный ∂-4 ∂-^ ∂y2 ∂y4 4* оператор четвертого порядка. Во многих случаях при изготовлении конструкций или их элементов рациональным является вариант, когда материал обладает анизотропией деформационно-прочностных свойств [5, 6]. Для упрочненных волокнами полимерных композиций анизотропия может регулироваться относительно простыми средствами, например, путем варьирования направлений укладки и соотношений между количествами подкрепляющих волокон. В работе рассматривается случай достаточно распространенных ортотропных материалов. В этом случае уравнение движения по форме несколько усложняется и принимает вид ∂4w∂4w∂4w ∂2w (2) D,, -+ 2Di. -■-“ + D.2 -"Г + ρhG(t) + ρh -- = 0 ; 11 ∂-4 ∂-2∂y2 ∂y^ ∂t2 ah3 Dnm = (n, m = 1,2), (3) где (3) - цилиндрическая жесткость пластины, в данном случае различная для разных направлений в силу разных модулей упругости материала в двух взаимно перпендикулярных направлениях, а коэффициенты апт для ортотропного материала определяются по соотношениям =E2ν1 , a... = , ai2 =- 1-ν1ν2 22 1-ν1ν2 12 1-ν1ν2 В дополнение к ранее введенным обозначениям принимается, что E1, Е2 - модули Юнга вдоль сторон прямоугольной пластины; ν1, ν2 - коэффициенты Пуассона. =E1=E2 (4) a11 ---, a22 =------ Моделирование поведения отклика ортотропной пластины 113 Начальные условия принимаются в виде равенства нулю перемещений и скоростей перемещений при t = 0: ∂w(x,y,0) w (x, y,0) = 0, - _ = 0 (5) ∂t Для описания закрепления пластины на ее контуре были применены два типа граничных условий - шарнирное и жесткое закрепление. Для шарнирного закрепления равны нулю перемещения и вторые производные от них, последнее означает отсутствие моментов на шарнирах. При так называемом жестком защемлении равны нулю перемещения и первая производная, что характеризует запреты на смещения и повороты кромок пластины. Условия для свободных границ записываются в виде равенств второй и третьей производных от перемещений по нормали к контуру пластины нулю. Это означает отсутствие приложенных моментов и перерезывающих сил на этих границах. В известном смысле жесткая и шарнирная заделки берут в вилку реальное закрепление на контуре. Для того чтобы модель описывала реальное закрепление, необходимо в областях крепления задать определенную степень свободы, величина которой подбирается экспериментальным путем. Таким образом, в каждой точке пластины на ее контуре задаются условия для свободных границ, для шарнирного или жесткого закрепления: d2W = 0 fd^ = 0І d3WW = 0 I -w = 0 I - свободные границы; l∂y3 J (6) ∂x^^ ∂y J ∂x^ 0 w 0 ґ w w = 0 -- = 0 I -- ∂x ^∂y 0 - шарнирное закрепление; (7) ∂w „ I ∂w ] (8) w = 0 - = 0 I - = 0 I - жесткое закрепление. ∂X k∂y J Поскольку в уравнениях движения производные по пространству имеют четвертый порядок, количество условий на границах (по два в каждой точке контура пластины) замыкает математическую модель. Производные по времени в уравнениях движения (1) или (2) имеют второй порядок, и каждое из этих уравнений можно привести к системе двух уравнений первого порядка путем введения функции скоростей .,(X, y,t) = dw^d-y,il ∂t (9) Расчеты проводились с использованием конечно-разностных соотношений. Для этого вводилась двумерная сетка Xj = i×∆x,yj∙ = j×∆y, i = 0, 1,_ , Nx, j = 0, 1,_, Ny, ∆x = a/Nx, ∆y = b/Ny, Nx, Ny - номера последних узлов вдоль направлений осей x, y, параллельных сторонам пластины, а и b - размеры пластины в плане. Число узлов по осям ОХ и OY составляет 24 и 14 соответственно. Проверка сеточной сходимости показала, что при этом погрешность расчетов вполне приемлема. Следует учесть, что сама модель пластины предполагает наличие погрешности порядка h/a, где а - минимальный из размеров пластины в плане. Далее вводятся соответствующие сеточные функции перемещений wij(t) = w(xi, yj, t) и скоростей перемещений w1ij (t) = w1(xi, yj, t). Дифференциальные операторы в уравнениях движения заменялись разностными. Для сеточных функций получается система обыкновенных дифференциальных уравнений, которая 114 И.В. Щербаков, Б.А. Люкшин решается методом Рунге - Кутты второго порядка. В результате получаются поля перемещений (прогибов), скоростей и ускорений точек пластины как функции времени. Поскольку аппроксимация производных возможна лишь для внутренних точек, получаемая система алгебраических уравнений относительно сеточных функций не замкнута. Для решения этой системы необходимо дополнить ее конечно-разностными соотношениями, аппроксимирующими граничные условия в каждой точке контура вида (6), или (7), или (8). Моделирование отклика пластинки при внешнем воздействии Нагрузка в виде зависимости ускорения точек (линий) крепления пластины показана на рис. 1. Задача решается для случая Е1 = Е2 = 22000 МПа [7]. Результаты приведены для узла с координатами по осям ОХ и OY - 12 и 7 соответственно. Первоначально нагрузка растет от нуля до 2 мс, далее уровень остается постоянным. Как видно по графику, отклик на пластине так же растет, как и нагрузка, до 2 мс, а после того как внешнее воздействие приобретает статический характер, отклик ускорений на плате начинает уменьшаться. Рис. 1. Внешнее воздействие и отклик, возникший на плате Fig. 1. External impact and a response arising on the plate Приведены результаты расчетов для случая, когда в качестве материала был выбран стеклотекстолит, это распространенный материал с ортотропными свойствами. Размеры пластины a×b×h составляют соответственно 0.180×0.080×0.002 м. Две стороны пластины (длиной 0.18 м) закреплены шарнирно, а две другие свободны. Модуль Е2 отвечает направлению вдоль коротких свободных сторон. Форма выпучивания пластины (рис. 2) напоминает арку, края которой закреплены, максимальные прогибы наблюдаются в ее середине. При этом на одной стороне пластины возникают сжимающие, на другой - растягивающие напряжения, связанные с модулем Е2. Рис. 2. Форма колебания пластины, условие закрепления: вдоль оси х - шарнир, вдоль оси у - свободные границы Fig. 2. Plate oscillation shape and fixing condition: a hinge along the x axis and free boundaries along the y axis На рис. 3 приведены зависимости максимальных прогибов пластинки от времени для случая ортотропной пластинки. Непрерывной сплошной линией показаны результаты для изотропного материала пластинки, другие кривые отвечают разным вариантам соотношений между модулями упругости вдоль сторон прямоугольной пластины, соответствующие обозначения показаны на рисунке. Физически результаты закономерны: изменение модуля упругости Е1 не оказывает Рис. 3. Зависимости прогибов пластинки в ее центре от времени для ортотропной пластинки при различных соотношениях между упругими модулями вдоль сторон (осей ортотропии), условия закрепления см. рис. 2 Fig. 3. Deflections of the plate center as functions of time for an orthotropic plate with different ratios between elastic modules along the sides (orthotropy axes); for fixing conditions, see Fig. 2 116 И.В. Щербаков, Б.А. Люкшин качественного влияния на результат, поскольку вдоль более длинной стороны пластина практически не деформируется. Что касается модуля упругости Е2, то с его увеличением цилиндрическая жесткость пластинки в соответствующем направлении увеличивается, а с уменьшением соответственно убывает, что приводит к противоположному изменению величин прогибов. Результаты натурных испытаний Для получения исходных данных по математическому моделированию и доказательства наличия ортотропных свойств материала стеклотестолит проведен эксперимент на разрывной машине Instron 5582 [8]. Размеры образцов составляли: 110.00×5.50×0.45 мм. Зона захвата тисков 25 мм с каждой стороны, рабочая зона -60 мм. Образцы вырезались в двух взаимно перпендикулярных направлениях, что позволяет оценить соответствующие свойства ортотропного материала. Результаты механических испытаний приведены на рис. 4. Рис. 4. Диаграммы деформаций, полученные для материала стеклотекстолит при одноосных испытаниях вдоль осей ортотропии Fig. 4. Stress-strain diagrams obtained for a fiberglass laminate material during uniaxial testing along the orthotropy axes Первые три кривые, изображенные на диаграмме деформации, соответствуют образцам, полученным вдоль оси OX, 4 и 5 - вдоль оси OY. Результаты испытаний подтверждают наличие ортотропных свойств материала стеклотекстолит, что видно по различию характеру отклика материала на нагружение. Модули упругости составляют (15 ± 0.1) ГПа при воздействии вдоль оси OX и (11 ± 0.1) ГПа -вдоль OY. Эти результаты позволяют задать свойства материалов при математическом моделировании. Моделирование поведения отклика ортотропной пластины 117 Заключение Ортотропия деформационно-прочностных свойств материала пластины, в частности, наличие разных модулей упругости, оказывает заметное влияние на характер НДС, а предложенная и реализованная модель позволяет получить количественные оценки этого влияния. Созданная модель и полученные с ее применением результаты основаны на данных физических экспериментов, полученных в лабораторных условиях. Разработанная и реализованная модель применима для оценки параметров НДС в реальных условиях эксплуатации, которые нельзя полностью и точно воспроизвести в лабораторных условиях.

Скачать электронную версию публикации