Свободные колебания квадратной консольной пластины
Цель: Модальный анализ квадратной консольной пластины.
Файл с исходными данными: 5.5.spr
Формулировка задачи: Определить собственные формы и частоты колебаний ω квадратной консольной пластины с плотностью материала ρ.
Ссылки: И. А. Биргер, Я. Г. Пановко, Прочность, устойчивость, колебания, Справочник в трех томах, Том 3, Москва, Машиностроение, 1968, стр. 382.
Исходные данные:
E = 2.06·108 кПа | - модуль упругости; |
ν = 0.3 | - коэффициент Пуассона; |
ρ = 7.85 т/м3 | - плотность материала; |
h = 0.01 м | - толщина пластины; |
a1 = 1.0 м | - размер длинной стороны пластины (вдоль оси X общей системы координат); |
a2 = 1.0 м | - размер короткой стороны пластины (вдоль оси Y общей системы координат). |
Конечноэлементная модель: Расчетная схема – балочный ростверк / плита, 400 элементов плиты типа 20. Сетка конечных элементов разбита по длинам сторон пластины (вдоль осей X, Y общей системы координат) с шагом 0.05 м. Обеспечение граничных условий достигается за счет наложения связей по направлению степеней свободы Z, UX, UY для одной из кромок, расположенных вдоль оси Y общей системы координат. Распределенная масса задается преобразованием статической нагрузки от собственного веса пластины ow = γ∙h, где γ = ρ∙g = 77.01 кН/м3. Количество узлов в расчетной схеме – 441. Определение собственных форм и частот выполнено методом Ланцоша. При расчете используется согласованная матрица масс.
Результаты решения в SCAD
Расчетная схема
1-я собственная форма колебаний
2-я собственная форма колебаний
3-я собственная форма колебаний
4-я собственная форма колебаний
5-я собственная форма колебаний
Сравнение решений:
Собственные частоты колебаний ω, рад / с
Форма колебаний |
Узловые линии |
Теория |
SCAD |
Отклонения, % |
---|---|---|---|---|
1 |
54.2 |
53.8 |
0.71 |
|
2 |
132.5 |
131.9 |
0.43 |
|
3 |
332.4 |
330.0 |
0.72 |
|
4 |
425.7 |
421.8 |
0.91 |
|
5 |
483.2 |
480.3 |
0.61 |
Замечания: При аналитическом решении собственные частоты колебаний ω квадратной консольной пластины с плотностью материала ρ могут быть вычислены по следующей формуле, полученной на основе решения по методу Рэлея-Ритца:
\[\omega =\frac{\omega_{m}^{\ast } }{a_{1}^{2}}\cdot \left( {\frac{D}{\rho \cdot h}} \right)^{\frac{1}{2}},\quad где \quad при \quad\frac{a_{2} }{a_{1} }=1\quad: \] \[ \omega_{1}^{\ast } =3.494, \quad \omega_{2}^{\ast } =8.547, \quad \omega_{3}^{\ast } =21.44, \quad \omega_{2}^{\ast } =27.46, \quad \omega_{2}^{\ast } =31.17, \quad D=\frac{E\cdot h^{3}}{12\cdot \left( {1-\nu^{2}} \right)} \]