employee
Moskva, Moscow, Russian Federation
student
Moskva, Moscow, Russian Federation
GRNTI 67.11 Строительные конструкции
BBK 385 Строительные конструкции
The calculation is based on the approximation of the initial differential equation of the bent of a thin isotropic plate by due to generalized equations of the finite difference method (FDM) as to the second derivatives of the deflection function. These equations allow to solve a problem taking into account the discontinuities of a required function, its first-order derivative and the right-hand side of a primitive differential equation within integration domain. Solving the system of equations. We can obtain the values of the second derivatives of the required function at each calculated point of mesh. Using known values, we can move to bent moments, which simplifies the solution. The article describes the algorithm for solving problems using the proposed methodology. Examples of calculations with various boundary conditions and load with the minimum number of partitions. The results are compared with the decision of S.P. Timoshenko in the sequences. Such approach can be used as methodological recommendations, also for carrying out verification calculations in the design of structures.
bent plate, refined, isotropic, the required function, discontinuities of the function, computational solution, generalized equations of finite difference method
Разрешающее дифференциальное уравнение изгиба тонкой пластинки [1] имеет вид:
(1)
где – прогиб;
– интенсивность распределенной нагрузки;
– цилиндрическая жесткость.
Перейдем к безразмерным параметрам и понизим порядок уравнения (1), введя обозначения
(2)
получим:
(3)
где
- интенсивность нагрузки в какой-либо точке,
– длина одной из сторон плиты.
Внутренние усилия также запишем в безразмерном виде:
(4)
Приведем разностную аппроксимацию дифференциального уравнения (3), используя обобщенное уравнение метода конечных разностей МКР (2.1.17) [2] при ,
с заменой
на
.
(5)
(6)
Здесь и
– величины разрывов безразмерной поперечной силы в направлениях
соответственно в точке, расположенной на бесконечно малом расстоянии от точки
.
Уравнения записываем на квадратной сетке с шагом h. Фрагмент сетки на которой строится численное решение показан на рис. 1.
Алгоритм расчета сводится к следующему: для определения и
необходимо совместное решение систем уравнений (5) и (6) с учетом граничных условий.
При шарнирном опирании на контуре , поэтому для решения достаточно уравнений (5) и (6).
|
|
Рис. 1. Шаблон с расчётными точками Рис. 2. Шаблон для записи граничных условий
Получим уравнение описывающее граничные условия при жестком защемлении стороны , параллельной оси
. Граничные условия можно записать в виде уравнения метода последовательных аппроксимаций используя (3.1.6) [2]:
(7)
и выражение , при
, где
– точка края, получим
(8)
Для других краев плиты эти уравнения записываются по аналогии с заменой на
. Фрагмент сетки приведен на рис. 2.
Для вычисления вторых производных воспользуемся известными формулами метода конечных разностей:
(9)
Формула для следует из (9) с заменой
на
.
Приведем решение нескольких задач.
1. Рассмотрим квадратную шарнирно опёртую по контуру плиту, нагруженную распределенными моментами по краям . Расчетная схема изображена на рис.3. Принимаем
. При такой нагрузке момент на контуре
.
найдем из зависимостей (4) -
Для решения необходимо записать уравнения (5) и (6) для расчетных точек поля с учетом симметрии при
.
Решая систему получим: Из (4) определим:
Результаты, приведенные в работе С.П. Тимошенко [1]:
. Погрешность по моментам
– 9,7%,
– 6,7%.
|
Рис. 3. Расчетная схема к задаче 1
2. Квадратная пластинка, три края которой свободно оперты и один защемлен, загружена равномерно распределенной нагрузкой . Граничные условия как на рис 4.
Для решения этой задачи наряду с уравнениями (5) и (6) составленных для точки поля 22, необходимо записать уравнение (8) для точки контура 23, учитывающее граничные условия. Значение прогиба в точке 22, в уравнении (8), запишем с использованием известного уравнения метода конечных разностей (9).
(10)
Решая (10) найдем: С учётом (4)
Прогиб
Решение в рядах [1] для этой задачи дает:
Погрешность по прогибам -
по моменту в заделке –
.
Рис. 4. Расчетная схема к задаче 3
|
Рис. 5. Расчетная схема к задаче 4
|
3. Рассмотрим плиту, представленную на рис. 4 под действием гидростатической нагрузки. Уравнения (5) и (6) записываются аналогично предыдущему примеру, при этом ,5. Уравнение (8) записываем с учетом (9) полагая
. В результате получим:
С учётом (4) при
–
Прогиб
Решение [1]:
Погрешность по прогибам –
по моменту в заделке -
.
4. Пластинка, жестко закрепленная двумя противоположными сторонами и свободно опертая двумя другими под действием равномерно распределенной нагрузки по всей поверхности (рис. 5). Для решения записываем уравнения (5) и (6) для точки 22, уравнение (8), с учетом (9) для точки 23. При шаге , получим следующие результаты:
С учётом (4) при
-
Прогиб
Решение [1]:
Погрешность по прогибу -
по моменту в заделке -
.
Таблица 1
Результаты расчетов задач на сетке с шагом и
|
|
|
|
|
|
задача №1 |
1/2 |
0,0312 |
0,422 |
0,422 |
- |
1/4 |
0,0356 |
0,281 |
0,369 |
- |
|
решение [1] |
0,0368 |
0,256 |
0,394 |
- |
|
задача №2 |
1/2 |
0,0026 |
0,027 |
0,027 |
-0,083 |
1/4 |
0,0027 |
0,035 |
0,031 |
-0,0838 |
|
решение [1] |
0,0028 |
0,039 |
0,034 |
-0,084 |
|
задача №3 |
1/2 |
0,0011 |
0,012 |
0,012 |
-0,051 |
1/4 |
0,0013 |
0,015 |
0,017 |
-0,043 |
|
решение [1] |
0,0013 |
0,016 |
0,019 |
-0,048 |
|
задача №4 |
1/2 |
0,0018 |
0,019 |
0,019 |
-0,067 |
1/4 |
0,00187 |
0,022 |
0,028 |
-0,069 |
|
решение [1] |
0,0019 |
0,024 |
0,033 |
-0,069 |
Решен ряд тестовых задач. Из анализа приведенных результатов следует, что решения с помощью изложенного алгоритма достоверны, погрешность при минимальном числе разбиений плиты на элементы не значительная. Нагрузка может быть любого типа, в том числе сосредоточенная с использованием подхода, показанного в [5], [8], [9], [11]. Исходя из этого, данный метод может быть использован для проведения расчетов подобных задач и для задач с другими вариантами нагрузок и граничных условий. Работа может представлять интерес с методической точки зрения, а также в связи с возможностью определения непосредственно моментов. В заключении отметим, что существует обширная литература, посвященная расчету изгибаемых плит. Расчет методом конечных разностей в традиционной форме представлен в работах [3], [4], [12]. Вариационные методы использованы в [1], [7], [10]. А также решение методом конечных элементов реализовано в вычислительных комплексах [6], [13].
1. Timoshenko S.P. Voinovskiy-Kriger S. Plate and shell translation from English. M, Science, 1966, 635 p.
2. Gabbasov R.F., Gabbasov A.R. Filatov V.V. Numerical construction of discontinuous solutions to the problems of structural mechanics. M.: ASV, 2008, 280 p.
3. Kiselev V.A. Calculation of plates. M: Stroiizdat, 1973, 151 p.
4. Vainberg D.V. Handbook of strength, stability and oscillations of wafers. K.: Budevilnik, 1973. 488 p.
5. Mikhailov B.A. Plates and shells with discontinuous parameters. L.: Publishing house of the Leningrad university, 1980, 196 p.
6. Zolotov A.B., Akimov P.A., Sidorov V.A. Mozgaleva M.L. Numerical and analytical methods for calculation building structures. M. : ASV, 2009. 336 p.
7. Vardanyan G.S., Andreev V.I., Atarov N.M., Gorshkov A.A. Resistance of materials with the fundamentals of the theory of elasticity and plasticity. M.: Infra-M; izdanie 2-e ispr. i dop. 2011. 638 p.
8. Gabbasov R. F., Filotov V.V., Ovarova N.B., Mansour A.M. Dissection method applications for complex shaped membranes and plates. Procedia engineering, 2016, pp. 444-449.
9. Uvarova N.B., Paramonov E.E. Application of the generalized equations of finite difference method to the calculation of bent plates of local and discontinuous loadings. Bulletin of BSTU named after V.G. Shokhov, 2018, no.1, pp. 56-59.
10. Birger I.A., Panovko Ya.G. Strenght, stability, vibrations. Reference book in 3 volumes. Volume 1. M : Mashinostroenie, 1968, 831 p.
11. Gabbasov R.F., An Hoang Tuan, An Nguen Hoang Comparison of the results of calculating thin bending plates using generalized equations of finite difference methods and successive approximations. Industrial and civil constrictions, 2014, no. 1, pp. 62-64.
12. Zhilkin V.A. Calculation of articulated rectangular plates by the finite difference method in MathCAD. APK Russian, 2017, vol. 24, no.1, pp. 119-129.
13. Karpilovskiy V.S., Kriskunov V.S., Malarenko E.Z., Perelmuter A.V., Perelmuter M.F., Fialko S.U. SCAD Office. Version 21. Computer complex SCAD++. «SCAD Soft», 2015. 848 p.