Русская Википедия:Изгиб пластин

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

Файл:BendingCircularPlate.png
Изгиб круглой пластины с закреплёнными краями под действием пореречной силы. Левая половина пластины показывает деформированную форму, а правая половина — недеформированную. Этот расчет произведён с использованием программы Ansys.

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

Изгиб пластин в теории Кирхгофа — Лява

Шаблон:Main

Файл:PlateForcesMoments upd.png
Силы и моменты действующие на плоской пластины для рассчёта равновесия элементарного объёма

Определения

Для тонкой прямоугольной пластины толщиной <math>H</math>, модулем Юнга <math>E</math> и коэффициентом Пуассона <math>\nu</math>, можно определить упругие параметры в терминах прогиба пластины <math>w</math>.

В декартовой системе координат жёсткость при изгибе определяется

<math>
   D = \frac{EH^3}{12\left(1-\nu^2\right)}.

</math>

Моменты

Изгибные моменты на единицу длины задаютсяШаблон:Sfn

<math>
   M_{x} = -D \left( \frac{\partial^2 w}{\partial x^2} + \nu \frac{\partial^2 w}{\partial y^2} \right),

</math>

<math>
   M_{y} = -D \left( \nu \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} \right).

</math>

Крутящий момент на единицу длины определяется

<math>
   M_{xy} = -D \left( 1 - \nu \right) \frac{\partial^2 w}{\partial x\partial y}.

</math>

Силы

Сдвиговые силы на единицу длины определяются выражениемШаблон:Sfn

<math>
   Q_{x} = -D \frac{\partial}{\partial x} \left( \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} \right),

</math>

<math>
   Q_{y} = -D \frac{\partial}{\partial y} \left( \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} \right).

</math>

Напряжения

Компоненты изгибных напряжений определяются выражением

<math>
   \sigma_{x} = -\frac{12Dz}{H^3} \left( \frac{\partial^2 w}{\partial x^2} + \nu \frac{\partial^2 w}{\partial y^2} \right),

</math>

<math>
   \sigma_{y} = -\frac{12Dz}{H^3} \left( \nu \frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} \right).

</math> Напряжение сдвига задается

<math>
   \tau_{xy} = -\frac{12Dz}{H^3} \left(1-\nu\right) \frac{\partial^2 w}{\partial x\partial y}.

</math>

Деформации

Изгибающие деформации в теории для малых отклонений определяются

<math>
   \epsilon_{x} = \frac{\partial u}{\partial x} = -z\frac{\partial^2 w}{\partial x^2},

</math>

<math>
   \epsilon_{y} = \frac{\partial v}{\partial y} = -z\frac{\partial^2 w}{\partial y^2}.

</math>

Деформации сдвига в теории для малых отклонений задаются

<math>
   \gamma_{xy} = \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} = -2z\frac{\partial^2 w}{\partial x\partial y}.

</math>

В теории для больших отклонений пластин рассматривают деформации мембран в виде

<math>
   \epsilon_{x} = \frac{\partial u}{\partial x} + \frac{1}{2}\left(\frac{\partial w}{\partial x}\right)^2,

</math>

<math>
   \epsilon_{y} = \frac{\partial v}{\partial y} + \frac{1}{2}\left(\frac{\partial w}{\partial y}\right)^2,

</math>

<math>
   \gamma_{xy} = \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} + \frac{\partial w}{\partial x} \frac{\partial w}{\partial y}.

</math>

Прогибы

Эти прогибы определяются

<math>
   u = -z\frac{\partial w}{\partial x},

</math>

<math>
   v = -z\frac{\partial w}{\partial y}.

</math>

Вывод

В теории пластин Кирхгофа — Лява система определяющих уравнений состоит из[1]

<math>
    N_{\alpha\beta,\alpha} = 0 
</math>

и

<math>
   M_{\alpha\beta,\alpha\beta} - q = 0.
</math>

Или в развёрнутой (координатной) форме

<math>
     \cfrac{\partial N_{11}}{\partial x_1} + \cfrac{\partial N_{21}}{\partial x_2} = 0 ~;~~
     \cfrac{\partial N_{12}}{\partial x_1} + \cfrac{\partial N_{22}}{\partial x_2} = 0,
</math>

и

<math>
     \cfrac{\partial^2 M_{11}}{\partial x_1^2} + 2\cfrac{\partial^2 M_{12}}{\partial x_1 \partial x_2} +
     \cfrac{\partial^2 M_{22}}{\partial x_2^2} = q.

</math> где <math>q(x)</math> приложенная поперечная нагрузка на единицу площади, а толщина плиты равна <math>H=2h</math>, напряжения <math>\sigma_{ij}</math>, и

<math>
  N_{\alpha\beta} := \int_{-h}^h \sigma_{\alpha\beta}~dx_3 ~;~~
  M_{\alpha\beta} := \int_{-h}^h x_3~\sigma_{\alpha\beta}~dx_3~.

</math> Величина <math>N</math> имеет размерность единицы силы на единицу длины. Величина <math>M</math> имеет размерность единицы момента на единицу длины.

Для изотропных, однородных пластин с модулем Юнга <math>E</math> и коэффициентом Пуассона <math>\nu</math> эти уравнения сводятся к[2]

<math>
  \nabla^2\nabla^2 w = -\cfrac{q}{D} ~;~~ D := \cfrac{2h^3E}{3(1-\nu^2)} = \cfrac{H^3E}{12(1-\nu^2)}
</math>

где <math>w(x_1,x_2)</math> прогиб средней поверхности пластины.

Малые прогибы тонких прямоугольных пластин

Малые прогибы тонких прямоугольных пластин описываются уравнением тонкой пластины Жермен — Лагранжа

<math>
  \cfrac{\partial^4 w}{\partial x^4}
  + 2\cfrac{\partial^4 w}{\partial x^2\partial y^2}
  + \cfrac{\partial^4 w}{\partial y^4}
  = \cfrac{q}{D}.
</math>

Это уравнение было впервые выведено Лагранжем в декабре 1811 г. который исправил доклад Софи Жермен.

Большой прогиб тонких прямоугольных пластин

Большой прогиб тонких прямоугольных пластин описывается уравнениями для пластины Феппля — фон Кармана

<math>
  \cfrac{\partial^4 F}{\partial x^4} +
  2\cfrac{\partial^4 F}{\partial x^2\partial y^2} +
  \cfrac{\partial^4 F}{\partial y^4} =
  E\left[\left(\cfrac{\partial^2 w}{\partial x \partial y}\right)^2 -
  \cfrac{\partial^2 w}{\partial x^2}
  \cfrac{\partial^2 w}{\partial y^2}\right],

</math>

<math>
  \cfrac{\partial^4 w}{\partial x^4} +
  2\cfrac{\partial^4 w}{\partial x^2\partial y^2} +
  \cfrac{\partial^4 w}{\partial y^4} =
  \cfrac{q}{D} + \cfrac{H}{D}\left(
  \cfrac{\partial^2 F}{\partial y^2}\cfrac{\partial^2 w}{\partial x^2} +
  \cfrac{\partial^2 F}{\partial x^2}\cfrac{\partial^2 w}{\partial y^2} -
  2\cfrac{\partial^2 F}{\partial x \partial y}\cfrac{\partial^2 w}{\partial x \partial y}
  \right),

</math> где <math>F</math> функция напряжения.

Круглые пластины Кирхгофа-Лява

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

Основное уравнение в безкоординатной форме имеет вид

<math>
 \nabla^2 \nabla^2 w = -\frac{q}{D} \,.

</math> В цилиндрических координатах <math>(r, \theta, z)</math>,

<math>
 \nabla^2 w \equiv \frac{1}{r}\frac{\partial }{\partial r}\left(r \frac{\partial w}{\partial r}\right) + 
     \frac{1}{r^2}\frac{\partial^2 w}{\partial \theta^2} + \frac{\partial^2 w}{\partial z^2} \,.

</math> Для симметрично нагруженных круглых пластин, где изгиб зависит от только радиуса <math> w = w(r)</math> получим

<math>
 \nabla^2 w \equiv \frac{1}{r}\cfrac{d }{d r}\left(r \cfrac{d w}{d r}\right) \,.

</math> Следовательно, основное уравнение приобретёт вид обыкновенного дифференциального уравненияШаблон:Sfn

<math>
  \frac{1}{r}\cfrac{d }{d r}\left[r \cfrac{d }{d r}\left\{\frac{1}{r}\cfrac{d }{d r}\left(r \cfrac{d w}{d r}\right)\right\}\right] = -\frac{q}{D}\,.

</math> Если <math>q</math> и <math>D</math> постоянны, то прямое интегрирование основного уравнения имеет решение

<math>

w(r) = -\frac{qr^4}{64 D} + C_1\ln r + \cfrac{C_2 r^2}{2} + \cfrac{C_3r^2}{4}(2\ln r - 1) + C_4 </math>

где <math>C_i</math> константы интегрирования. Наклон отклоняющей поверхности равен

<math>
 \phi(r) = \cfrac{d w}{d r} = -\frac{qr^3}{16D} + \frac{C_1}{r} + C_2 r + C_3 r \ln r \,.

</math> Для круглой пластины требование конечности прогиба и крутизны прогиба при <math>r = 0</math> подразумевает, что <math>C_1 = 0</math>. Однако, <math>C_3</math> не обязательно равняется 0, так как правый предел <math>r \ln r\,</math> существует по мере приближения к началу координат <math>r = 0</math>.

Закрепленные края

Для круглой пластины (радиуса a) с зажатыми краями <math>w(a) = 0</math> и <math>\phi(a) = 0</math> на краю пластины. Подставляя эти граничные условия в общее решение получаемШаблон:Sfn

<math>

w(r) = -\frac{q}{64 D} (a^2 -r^2)^2 \quad \text{and} \quad \phi(r) = \frac{qr}{16 D}(a^2-r^2) \,. </math>

Смещения пластины в плоскости равны

<math>
  u_r(r) = -z\phi(r) \quad \text{and} \quad u_\theta(r) = 0 \,.

</math> Плоские деформации в пластине равны

<math>
 \varepsilon_{rr} = \cfrac{d u_r}{d r} = -\frac{qz}{16D}(a^2-3r^2) ~,~~
 \varepsilon_{\theta\theta} = \frac{u_r}{r} = -\frac{qz}{16D}(a^2-r^2) ~,~~
 \varepsilon_{r\theta} = 0 \,.

</math> Напряжения в плоскости пластины равны

<math>
 \sigma_{rr} = \frac{E}{1-\nu^2}\left[\varepsilon_{rr} + \nu\varepsilon_{\theta\theta}\right] ~;~~
 \sigma_{\theta\theta} = \frac{E}{1-\nu^2}\left[\varepsilon_{\theta\theta} + \nu\varepsilon_{rr}\right] ~;~~
 \sigma_{r\theta} = 0 \,.

</math> Для плиты толщиной <math>2h</math>, жесткость на изгиб <math>D = 2Eh^3/[3(1-\nu^2)]</math> и

<math>

\begin{align} \sigma_{rr} &= -\frac{3qz}{32h^3}\left[(1+\nu)a^2-(3+\nu)r^2\right] \\ \sigma_{\theta\theta} &= -\frac{3qz}{32h^3}\left[(1+\nu)a^2-(1+3\nu)r^2\right]\\ \sigma_{r\theta} &= 0 \,. \end{align} </math>

Результирующие моменты (изгибные моменты) равны

<math>
 M_{rr} = -\frac{q}{16}\left[(1+\nu)a^2-(3+\nu)r^2\right] ~;~~
 M_{\theta\theta} = -\frac{q}{16}\left[(1+\nu)a^2-(1+3\nu)r^2\right] ~;~~
 M_{r\theta} = 0 \,.

</math> Максимальное радиальное напряжение при <math>z = h</math> и <math>r = a</math>:

<math>
 \left.\sigma_{rr}\right|_{z=h,r=a} = \frac{3qa^2}{16h^2} = \frac{3qa^2}{4H^2}

</math> где <math>H := 2h</math>. Изгибающие моменты на границе и в центре пластины равныШаблон:Sfn

<math>
 \left.M_{rr}\right|_{r=a} = \frac{qa^2}{8} ~,~~
 \left.M_{\theta\theta}\right|_{r=a} = \frac{\nu qa^2}{8} ~,~~
 \left.M_{rr}\right|_{r=0} = \left.M_{\theta\theta}\right|_{r=0} = -\frac{(1+\nu) qa^2}{16} \,.

</math>

Круговая пластина нагруженная силой зависящей от радиуса

Шаблон:Sfn

Прямоугольные пластины Кирхгофа-Лява

Файл:RectangularPlateBending.svg
Изгиб прямоугольной пластины под действием распределенной силы <math>q</math> на единицу площади.

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

Синусоидальная нагрузка

Предположим, что нагрузка имеет видШаблон:Sfn

<math>
  q(x,y) = q_0 \sin\frac{\pi x}{a}\sin\frac{\pi y}{b} \,.

</math> Здесь <math>q_0</math> амплитуда, <math>a</math> ширина пластины в направлении <math>x</math> и <math>b</math> ширина пластины в направлении <math>y</math>.

Поскольку пластина просто поддерживается на краях, то смещение <math>w(x,y)</math> на краях пластины равно нулю, и изгибающий момент <math>M_{xx}</math> также равен нулю на границах <math>x=0</math> и <math>x=a</math>, <math>M_{yy}</math> равен нулю на границах <math>y=0</math> и <math>y=b</math>.

При этих граничные условиях и решение уравнения для пластины имеет видШаблон:Sfn

<math>
  w(x,y) = \frac{q_0}{\pi^4 D}\,\left(\frac{1}{a^2}+\frac{1}{b^2}\right)^{-2}\,\sin\frac{\pi x}{a}\sin\frac{\pi y}{b} \,.

</math> Где D жесткость на изгиб

<math>

D=\frac{Et^3}{12(1-\nu^2)}. </math> Analogous to flexural stiffness EI.[3] Напряжения и деформации в пластине можно рассчитать, если знаем смещение.

При общей нагрузки в виде

<math>
  q(x,y) = q_0 \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} 

</math> где <math>m</math> и <math>n</math> целые, получим решениеШаблон:Sfn

<math> \text{(1)} \qquad

w(x,y) = \frac{q_0}{\pi^4 D}\,\left(\frac{m^2}{a^2}+\frac{n^2}{b^2}\right)^{-2}\,\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} \,. </math>

Решение Навье

Уравнение для двумерного тригонометрического ряда

Определяем общую нагрузку <math>q(x,y)</math> в видеШаблон:Sfn

<math>
  q(x,y) = \sum_{m=1}^{\infty} \sum_{n=1}^\infty a_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}

</math> где <math>a_{mn}</math> коэффициент Фурье, определяемый формулойШаблон:Sfn

<math>
  a_{mn} = \frac{4}{ab}\int_0^b \int_0^a q(x,y)\sin\frac{m\pi x}{a}\sin\frac{n\pi y}{b}\,\text{d}x\text{d}y

</math>. Таким образом, классическое уравнение прямоугольной пластины для малых прогибов принимает следующий вид:

<math>
  \cfrac{\partial^4 w}{\partial x^4} + 2\cfrac{\partial^4 w}{\partial x^2\partial y^2} + \cfrac{\partial^4 w}{\partial y^4} =
  \cfrac{1}{D} \sum_{m=1}^{\infty} \sum_{n=1}^\infty a_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}

</math>

Свободно опёртая пластина с общей нагрузкой

Предполагаем решение <math>w(x,y)</math> вида

<math>
  w(x,y) = \sum_{m=1}^{\infty} \sum_{n=1}^\infty w_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}

</math> Частные дифференциалы этой функции даются выражениями

<math>
  \cfrac{\partial^4 w}{\partial x^4} = \sum_{m=1}^{\infty} \sum_{n=1}^\infty
  \left(\frac{m \pi}{a}\right)^4 w_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}.

</math>

<math>
  \cfrac{\partial^4 w}{\partial x^2 \partial y^2} = \sum_{m=1}^{\infty} \sum_{n=1}^\infty
  \left(\frac{m \pi}{a}\right)^2 \left(\frac{n \pi}{b}\right)^2 w_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b},

</math>

<math>
  \cfrac{\partial^4 w}{\partial y^4} = \sum_{m=1}^{\infty} \sum_{n=1}^\infty
  \left(\frac{n \pi}{b}\right)^4 w_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}.

</math> Подставляя эти выражения в уравнение для пластины, получим

<math>
  \sum_{m=1}^{\infty} \sum_{n=1}^\infty
  \left( \left(\frac{m \pi}{a}\right)^2 + \left(\frac{n \pi}{b}\right)^2 \right)^2
  w_{mn}\sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} =
  \sum_{m=1}^{\infty} \sum_{n=1}^\infty \cfrac{a_{mn}}{D} \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b}

</math> Приравнивая два ряда получим для коэффициентов

<math>
  \left( \left(\frac{m \pi}{a}\right)^2 + \left(\frac{n \pi}{b}\right)^2 \right)^2 w_{mn} = \cfrac{a_{mn}}{D}

</math> или при перестановки получим

<math>
  w_{mn} = \frac{1}{\pi^4 D}\frac{a_{mn}}{\left(\frac{m^2}{a^2} + \frac{n^2}{b^2}\right)^2}

</math> Прогиб свободно опертой пластины (на углах) при общей нагрузке задаётся выражениемШаблон:Sfn

<math>

w(x,y) = \frac{1}{\pi^4 D} \sum_{m=1}^\infty \sum_{n=1}^\infty \frac{a_{mn}}{\left(\frac{m^2}{a^2} + \frac{n^2}{b^2}\right)^2} \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} </math>

Свободно опёртая пластина с постоянной нагрузкой

Шаблон:Multiple image

Для равномерно распределенной нагрузки имеем

<math>
  q(x,y) = q_0

</math> Таким образом, соответствующий коэффициент Фурье определяется выражением

<math>
  a_{mn} = \frac{4}{ab} \int_0^a \int_0^b q_0\sin\frac{m\pi x}{a}\sin\frac{n\pi y}{b}\,\text{d}x\text{d}y

</math>. Вычисляя двойной интеграл, имеем

<math>
  a_{mn} = \frac{4q_0}{\pi^2 mn}(1 - \cos m\pi)(1 - \cos n\pi)

</math>, или в другом виде кусочно-заданной функции

<math>
  a_{mn} = \begin{cases}
           \cfrac{16q_0}{\pi^2 mn} & m~\text{and}~n~\text{odd} \\
           0 & m~\text{or}~n~\text{even}
           \end{cases} 

</math>

Прогиб свободно опертой пластины (с условиями на углах) с равномерно распределенной нагрузкой определяется выражением

<math>

w(x,y) = \frac{16 q_0}{\pi^6 D} \sum_{m=1,3,5,...}^\infty \sum_{n=1,3,5,...}^\infty \frac{1}{mn\left(\frac{m^2}{a^2} + \frac{n^2}{b^2}\right)^2} \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} </math>

Изгибающие моменты на единицу длины в пластине определяются выражением

<math>

M_{x} = \frac{16 q_0}{\pi^4} \sum_{m=1,3,5,...}^\infty \sum_{n=1,3,5,...}^\infty \frac{\frac{m^2}{a^2} + \nu\frac{n^2}{b^2}} {mn\left(\frac{m^2}{a^2} + \frac{n^2}{b^2}\right)^2} \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} </math>

<math>

M_{y} = \frac{16 q_0}{\pi^4} \sum_{m=1,3,5,...}^\infty \sum_{n=1,3,5,...}^\infty \frac{\frac{n^2}{b^2} + \nu\frac{m^2}{a^2}} {mn\left(\frac{m^2}{a^2} + \frac{n^2}{b^2}\right)^2} \sin\frac{m \pi x}{a}\sin\frac{n \pi y}{b} </math>

Решение Леви

Другой подход был предложен Леви[4] в 1899 году. В этом случае мы начинаем с предполагаемой формы смещения и пытаемся подогнать параметры так, чтобы выполнялись определяющее уравнение и граничные условия. Цель — найти рашения основного уравнения <math>\nabla^2 \nabla^2 w = q/D</math><math>Y_m(y)</math> такие, что они удовлетворяют граничным условиям при <math>y = 0</math> и <math>y = b</math>.

Предположим, чтоШаблон:Sfn

<math>
  w(x,y) = \sum_{m=1}^\infty Y_m(y) \sin \frac{m\pi x}{a} \,.

</math> Для пластины, которая свободно опирается краями при <math>x=0</math> и <math>x=a</math>, граничные условия: <math>w=0</math> и <math>M_{xx}=0</math>. Обратите внимание, что на этих краях нет изменений смещения, что означает <math>\partial w/\partial y = 0</math> и <math>\partial^2 w/\partial y^2 = 0</math>, сводя, таким образом, моментное граничное условие к эквивалентному выражению <math>\partial^2 w/\partial x^2 = 0</math>.

Моменты на краях

Рассмотрим случай чисто моментной нагрузки. В этом случае <math>q = 0</math> и функция <math>w(x,y)</math> должна удовлетворять уравнению <math>\nabla^2 \nabla^2 w = 0</math>. в В прямоугольных декартовых координатах основное уравнение выражается как

<math>
  \frac{\partial^4 w}{\partial x^4} + 2 \frac{\partial^4 w}{\partial x^2\partial y^2}
  + \frac{\partial^4 w}{\partial y^4}  = 0 \,.

</math> Подставляем выражение для <math>w(x,y)</math> в основное уравнение что приводит кШаблон:Sfn

<math>
 \sum_{m=1}^\infty \left[\left(\frac{m\pi}{a}\right)^4 Y_m \sin\frac{m\pi x}{a}
  - 2\left(\frac{m\pi}{a}\right)^2 \cfrac{d^2 Y_m}{d y^2} \sin\frac{m\pi x}{a}
  + \frac{d^4Y_m}{dy^4} \sin\frac{m\pi x}{a}\right] = 0 

</math> или

<math>
  \frac{d^4Y_m}{dy^4}  - 2 \frac{m^2\pi^2}{a^2} \cfrac{d^2Y_m}{dy^2} + \frac{m^4\pi^4}{a^4} Y_m = 0 \,.

</math> Это обыкновенное дифференциальное уравнение, имеющее общее решениеШаблон:Sfn

<math>
 Y_m = A_m \cosh\frac{m\pi y}{a} + B_m\frac{m\pi y}{a} \cosh\frac{m\pi y}{a} + 
 C_m \sinh\frac{m\pi y}{a} + D_m\frac{m\pi y}{a} \sinh\frac{m\pi y}{a}  

</math> где <math>A_m, B_m, C_m, D_m</math> константы, которые можно определить из граничных условий. Следовательно, изгибное решение имеет вид

<math>

w(x,y) = \sum_{m=1}^\infty \left[ \left(A_m + B_m\frac{m\pi y}{a}\right) \cosh\frac{m\pi y}{a} + \left(C_m + D_m\frac{m\pi y}{a}\right) \sinh\frac{m\pi y}{a} \right] \sin \frac{m\pi x}{a} \,. </math>

Выберем систему координат так, чтобы границы пластины находились на краях при <math>x = 0</math> и <math>x = a</math>, при <math>y = \pm b/2</math>. Тогда граничные условия на моменты при <math>y = \pm b/2</math>

<math>
 w = 0 \,, -D\frac{\partial^2 w}{\partial y^2}\Bigr|_{y=b/2} = f_1(x) \,, 
 -D\frac{\partial^2 w}{\partial y^2}\Bigr|_{y=-b/2} = f_2(x)  

</math> где <math>f_1(x), f_2(x)</math> известные функции. Решение можно найти, используя эти граничные условия. Можно показать, что для симметричного случая, когда

<math>
  M_{yy}\Bigr|_{y=-b/2} = M_{yy}\Bigr|_{y=b/2}

</math> и

<math>
  f_1(x) = f_2(x) = \sum_{m=1}^\infty E_m\sin\frac{m\pi x}{a}

</math> получимШаблон:Sfn

<math>

w(x,y) = \frac{a^2}{2\pi^2 D}\sum_{m=1}^\infty \frac{E_m}{m^2\cosh\alpha_m}\, \sin\frac{m\pi x}{a}\, \left(\alpha_m \tanh\alpha_m \cosh\frac{m\pi y}{a} - \frac{m\pi y}{a}\sinh\frac{m\pi y}{a}\right) </math>

где

<math>
 \alpha_m = \frac{m\pi b}{2a} \,.

</math> Аналогично для антисимметричного случая, когда

<math>
  M_{yy}\Bigr|_{y=-b/2} = -M_{yy}\Bigr|_{y=b/2}

</math> получимШаблон:Sfn

<math>

w(x,y) = \frac{a^2}{2\pi^2 D}\sum_{m=1}^\infty \frac{E_m}{m^2\sinh\alpha_m}\, \sin\frac{m\pi x}{a}\, \left(\alpha_m \coth\alpha_m \sinh\frac{m\pi y}{a} - \frac{m\pi y}{a}\cosh\frac{m\pi y}{a}\right) \,. </math>

Используя симметричные и антисимметричные решения, можно составить более общие решения.

Опёртая пластина с равномерно распределенной нагрузкой

Для равномерно распределенной нагрузки

<math>
  q(x,y) = q_0

</math>

Отклонение опёртой пластины с центром при <math>\left(\frac{a}{2}, 0\right)</math> с равномерно распределенной нагрузкой определяется выражениемШаблон:Sfn

<math>

\begin{align} &w(x,y) = \frac{q_0 a^4}{D} \sum_{m=1,3,5,...}^\infty \left( A_m\cosh\frac{m\pi y}{a} + B_m\frac{m\pi y}{a}\sinh\frac{m\pi y}{a} + G_m\right) \sin\frac{m\pi x}{a}\\\\ &\begin{align} \text{where}\quad &A_m = -\frac{2\left(\alpha _m\tanh\alpha _m + 2\right)}{\pi^5 m^5 \cosh\alpha _m}\\ &B_m = \frac{2}{\pi^5 m^5 \cosh\alpha _m}\\ &G_m = \frac{4}{\pi^5 m^5}\\\\ \quad &\alpha _m = \frac{m\pi b}{2a} \end{align} \end{align} </math>

Изгибающие моменты на единицу длины в пластине определяются выражениями

<math>

M_x = -q_0\pi^2 a^2\sum_{m=1,3,5,...}^\infty m^2\left( \left(\left(\nu -1\right)A_m + 2\nu B_m\right)\cosh\frac{m\pi y}{a} + \left(\nu -1\right)B_m\frac{m\pi y}{a}\sinh\frac{m\pi y}{a} - G_m\right) \sin\frac{m\pi x}{a} </math>

<math>

M_y = -q_0\pi^2 a^2\sum_{m=1,3,5,...}^\infty m^2\left( \left(\left(1-\nu\right)A_m + 2B_m\right)\cosh\frac{m\pi y}{a} + \left(1-\nu\right)B_m\frac{m\pi y}{a}\sinh\frac{m\pi y}{a} - \nu G_m\right) \sin\frac{m\pi x}{a} </math>

Равномерная и симметричная моментная нагрузка

Для частного случая, когда нагрузка симметрична и момент однороден, при <math>y=\pm b/2</math>,

<math>
  M_{yy} = f_1(x) = \frac{4M_0}{\pi}\sum_{m=1}^\infty \frac{1}{2m-1}\,\sin\frac{(2m-1)\pi x}{a} \,.

</math>

Шаблон:Multiple image

Результирующий изгиб равен

<math>

\begin{align} & w(x,y) = \frac{2M_0 a^2}{\pi^3 D}\sum_{m=1}^\infty \frac{1}{(2m-1)^3\cosh\alpha_m}\sin\frac{(2m-1)\pi x}{a} \times\\ & ~~ \left[ \alpha_m\,\tanh\alpha_m\cosh\frac{(2m-1)\pi y}{a} -\frac{(2m-1)\pi y}{a} \sinh\frac{(2m-1)\pi y}{a}\right] \end{align} </math>

где

<math>
 \alpha_m = \frac{\pi (2m-1)b}{2a} \,. 

</math> Изгибающие моменты и поперечные силы, соответствующие смещению <math>w</math> находятся по формулам

<math>
 \begin{align}
    M_{xx} & = -D\left(\frac{\partial^2 w}{\partial x^2}+\nu\,\frac{\partial^2 w}{\partial y^2}\right) \\
           & = \frac{2M_0(1-\nu)}{\pi}\sum_{m=1}^\infty\frac{1}{(2m-1)\cosh\alpha_m}\,\times \\
           & ~    \sin\frac{(2m-1)\pi x}{a} \,\times \\
           & ~   \left[
                 -\frac{(2m-1)\pi y}{a}\sinh\frac{(2m-1)\pi y}{a} + \right. \\
           & \qquad \qquad \qquad \qquad
             \left. \left\{\frac{2\nu}{1-\nu} + \alpha_m\tanh\alpha_m\right\}\cosh\frac{(2m-1)\pi y}{a}
               \right] \\
    M_{xy} & = (1-\nu)D\frac{\partial^2 w}{\partial x \partial y} \\
           & = -\frac{2M_0(1-\nu)}{\pi}\sum_{m=1}^\infty\frac{1}{(2m-1)
                   \cosh\alpha_m}\,\times \\
           & ~  \cos\frac{(2m-1)\pi x}{a} \, \times \\
           & ~  \left[\frac{(2m-1)\pi y}{a}\cosh\frac{(2m-1)\pi y}{a} + \right. \\
           & \qquad \qquad \qquad \qquad
             \left. (1-\alpha_m\tanh\alpha_m)\sinh\frac{(2m-1)\pi y}{a}\right] \\
    Q_{zx} & = \frac{\partial M_{xx}}{\partial x}-\frac{\partial M_{xy}}{\partial y} \\
           & = \frac{4M_0}{a}\sum_{m=1}^\infty \frac{1}{\cosh\alpha_m}\,\times \\
           & ~    \cos\frac{(2m-1)\pi x}{a}\cosh\frac{(2m-1)\pi y}{a}\,.
 \end{align}

</math> Напряжения

<math>
 \sigma_{xx} = \frac{12z}{h^3}\,M_{xx} \quad \text{and} \quad 
 \sigma_{zx} = \frac{1}{\kappa h}\,Q_{zx}\left(1 - \frac{4z^2}{h^2}\right)\,. 

</math>

Изгиб цилиндрической пластины

Цилиндрический изгиб возникает, когда прямоугольная пластина имеющая размеры <math>a \times b \times h</math>, где <math>a \ll b</math> и малую толщину <math>h</math>, подвергается равномерной распределенной нагрузке, перпендикулярной плоскости пластины. Такая пластина имеет форму поверхности цилиндра.

С помощью методов Навье и Леви также можно найти решения для свободно опёртых пластин при цилиндрическом изгибе с различным количеством незакреплённых краёвШаблон:Sfn.

Изгиб толстых пластин Миндлина

Для толстых пластин необходимо учитывать влияние сдвиговых напряжений по толщине на ориентацию нормали к средней поверхности после деформации. Теория Миндлина предлагает единый подход к нахождению деформации и напряжений в таких пластинах. Решения теории Миндлина можно получить из эквивалентных решений Кирхгофа-Лява с использованием канонических соотношений[5].

Основные уравнения

Канонические уравнения для изотропных толстых пластин можно записать в виде[5]

<math>
 \begin{align}
   & \nabla^2 \left(\mathcal{M} - \frac{\mathcal{B}}{1+\nu}\,q\right) = -q \\
   & \kappa G h\left(\nabla^2 w + \frac{\mathcal{M}}{D}\right) =
     -\left(1 - \cfrac{\mathcal{B} c^2}{1+\nu}\right)q \\
   & \nabla^2 \left(\frac{\partial \varphi_1}{\partial x_2} - \frac{\partial \varphi_2}{\partial x_1}\right)
     = c^2\left(\frac{\partial \varphi_1}{\partial x_2} - \frac{\partial \varphi_2}{\partial x_1}\right)
 \end{align}

</math> где <math>q</math> приложенная поперечная нагрузка, <math>G</math> модуль сдвига, <math>D = Eh^3/[12(1-\nu^2)]</math> жесткость на изгиб, <math>h</math> толщина пластины, <math>c^2 = 2\kappa G h/[D(1-\nu)]</math>, <math>\kappa</math> коэффициент поправки сдвигового напряжения, <math>E</math> модуль Юнга, <math>\nu</math> коэффициент Пуассона и

<math>
    \mathcal{M}  = D\left[\mathcal{A}\left(\frac{\partial \varphi_1}{\partial x_1} + \frac{\partial \varphi_2}{\partial x_2}\right)
   - (1-\mathcal{A})\nabla^2 w\right] + \frac{2q}{1-\nu^2}\mathcal{B}  \,.

</math> Согласно теории Миндлина <math>w</math> поперечное смещение средней поверхности пластины, а величины <math>\varphi_1</math> и <math>\varphi_2</math> соответственные повороты нормали к средней поверхности относительно <math>x_2</math> и <math>x_1</math>-осей. Канонические параметры этой теории <math>\mathcal{A} = 1</math> и <math>\mathcal{B} = 0</math>. Коэффициент поправки сдвигового напряжения <math>\kappa</math> обычно принимают за <math>5/6</math>.

Решения основных уравнений можно найти, если знать соответствующие решения Кирхгофа-Лява с помощью соотношений

<math>
 \begin{align}
   w & = w^K + \frac{\mathcal{M}^K}{\kappa G h}\left(1 - \frac{\mathcal{B} c^2}{2}\right)
       - \Phi + \Psi \\
   \varphi_1 & = - \frac{\partial w^K}{\partial x_1}
    - \frac{1}{\kappa G h}\left(1 - \frac{1}{\mathcal{A}} - \frac{\mathcal{B} c^2}{2}\right)Q_1^K
    + \frac{\partial }{\partial x_1}\left(\frac{D}{\kappa G h \mathcal{A}}\nabla^2 \Phi + \Phi - \Psi\right)
    + \frac{1}{c^2}\frac{\partial \Omega}{\partial x_2} \\
   \varphi_2 & = - \frac{\partial w^K}{\partial x_2}
    - \frac{1}{\kappa G h}\left(1 - \frac{1}{\mathcal{A}} - \frac{\mathcal{B} c^2}{2}\right)Q_2^K
    + \frac{\partial }{\partial x_2}\left(\frac{D}{\kappa G h \mathcal{A}}\nabla^2 \Phi + \Phi - \Psi\right)
    + \frac{1}{c^2}\frac{\partial \Omega}{\partial x_1}
 \end{align}

</math> где <math>w^K</math> это смещение, предсказанное для пластины Кирхгофа-Лява, <math>\Phi</math> бигармоническая функция такая, что <math>\nabla^2 \nabla^2 \Phi = 0</math>, <math>\Psi</math> функция, удовлетворяющая уравнению Лапласа, <math>\nabla^2 \Psi = 0</math> и

<math>
 \begin{align}
   \mathcal{M} & = \mathcal{M}^K + \frac{\mathcal{B}}{1+\nu}\,q + D \nabla^2 \Phi ~;~~ \mathcal{M}^K := -D\nabla^2 w^K \\
   Q_1^K & = -D\frac{\partial }{\partial x_1}\left(\nabla^2 w^K\right) ~,~~
   Q_2^K = -D\frac{\partial }{\partial x_2}\left(\nabla^2 w^K\right) \\
   \Omega & = \frac{\partial \varphi_1}{\partial x_2} - \frac{\partial \varphi_2}{\partial x_1} ~,~~ \nabla^2 \Omega = c^2\Omega \,.
 \end{align}

</math>

Свободно опёртые прямоугольные пластины

Для свободно опертых пластин сумма моментов Маркуса равна нулю

<math>
  \mathcal{M} = \frac{1}{1+\nu}(M_{11}+M_{22}) = D\left(\frac{\partial \varphi_1}{\partial x_1}+\frac{\partial \varphi_2}{\partial x_2}\right) = 0 \,.

</math> В этом случае функции <math>\Phi</math>, <math>\Psi</math>, <math>\Omega</math> равны нулю, а решение Миндлина связано с соответствующим решением Кирхгофа соотношением

<math>
 w = w^K + \frac{\mathcal{M}^K}{\kappa G h} \,.

</math>

Изгиб консольно-закреплённых пластин Рейсснера-Штейна

Теория Рейсснера-Штейна для консольных пластин[6] приводит к следующим связанным обыкновенным дифференциальным уравнениям для консольной пластины с сосредоточенной нагрузкой на торце <math>q_x(y)</math> в точке <math>x=a</math>.

<math>
 \begin{align}
   & bD \frac{\mathrm{d}^4w_x}{\mathrm{d}x^4}  = 0 \\
   & \frac{b^3D}{12}\,\frac{\mathrm{d}^4\theta_x}{\mathrm{d}x^4} - 2bD(1-\nu)\cfrac{d^2 \theta_x}{d x^2} = 0
 \end{align}

</math> и граничных условиях в точке <math>x=a</math>

<math>
 \begin{align}
 & bD\cfrac{d^3 w_x}{d x^3} + q_{x1} = 0 \quad,\quad
  \frac{b^3D}{12}\cfrac{d^3 \theta_x}{d x^3} -2bD(1-\nu)\cfrac{d \theta_x}{d x} + q_{x2} = 0 \\
 & bD\cfrac{d^2 w_x}{d x^2} = 0 \quad,\quad  \frac{b^3D}{12}\cfrac{d^2 \theta_x}{d x^2} = 0 \,.
 \end{align}

</math> Решение этой системы двух ОДУ дает

<math>
 \begin{align}
   w_x(x) & = \frac{q_{x1}}{6bD}\,(3ax^2 -x^3) \\
   \theta_x(x) & = \frac{q_{x2}}{2bD(1-\nu)}\left[x - \frac{1}{\nu_b}\,
      \left(\frac{\sinh(\nu_b a)}{\cosh[\nu_b (x-a)]} + \tanh[\nu_b(x-a)]\right)\right]
 \end{align}

</math> где <math>\nu_b = \sqrt{24(1-\nu)}/b</math>. Изгибные моменты и поперечные силы, соответствующие смещению <math>w = w_x + y\theta_x</math>

<math>
 \begin{align}
    M_{xx} & = -D\left(\frac{\partial^2 w}{\partial x^2}+\nu\,\frac{\partial^2 w}{\partial y^2}\right) \\
           & = q_{x1}\left(\frac{x-a}{b}\right) - \left[\frac{3yq_{x2}}{b^3\nu_b\cosh^3[\nu_b(x-a)]}\right]
               \times \\
           & \quad \left[6\sinh(\nu_b a) - \sinh[\nu_b(2x-a)] + 
                 \sinh[\nu_b(2x-3a)] + 8\sinh[\nu_b(x-a)]\right] \\
    M_{xy} & = (1-\nu)D\frac{\partial^2 w}{\partial x \partial y} \\
           & = \frac{q_{x2}}{2b}\left[1 - 
               \frac{2+\cosh[\nu_b(x-2a)] - \cosh[\nu_b x]}{2\cosh^2[\nu_b(x-a)]}\right] \\
    Q_{zx} & = \frac{\partial M_{xx}}{\partial x}-\frac{\partial M_{xy}}{\partial y} \\
           & = \frac{q_{x1}}{b} - \left(\frac{3yq_{x2}}{2b^3\cosh^4[\nu_b(x-a)]}\right)\times 
               \left[32 + \cosh[\nu_b(3x-2a)] - \cosh[\nu_b(3x-4a)]\right. \\
           & \qquad \left. - 16\cosh[2\nu_b(x-a)] +
                23\cosh[\nu_b(x-2a)] - 23\cosh(\nu_b x)\right]\,.
 \end{align}

</math> Напряжения

<math>
 \sigma_{xx} = \frac{12z}{h^3}\,M_{xx} \quad \text{and} \quad 
 \sigma_{zx} = \frac{1}{\kappa h}\,Q_{zx}\left(1 - \frac{4z^2}{h^2}\right)\,. 

</math> Если приложенная нагрузка на краю постоянна, мы восстанавливаем решения для балки при сосредоточенной торцевой нагрузке. Если приложенная нагрузка — линейная функция <math>y</math>, то

<math>
 q_{x1} = \int_{-b/2}^{b/2}q_0\left(\frac{1}{2} - \frac{y}{b}\right)\,\text{d}y = \frac{bq_0}{2} ~;~~
 q_{x2} = \int_{-b/2}^{b/2}yq_0\left(\frac{1}{2} - \frac{y}{b}\right)\,\text{d}y = -\frac{b^2q_0}{12} \,.

</math>

Ссылки

Шаблон:Примечания

Литература

Шаблон:Книга

  1. Reddy, J. N., 2007, Theory and analysis of elastic plates and shells, CRC Press, Taylor and Francis.
  2. Timoshenko, S. and Woinowsky-Krieger, S., (1959), Theory of plates and shells, McGraw-Hill New York.
  3. Cook, R. D. et al., 2002, Concepts and applications of finite element analysis, John Wiley & Sons
  4. Lévy, M., 1899, Comptes rendues, vol. 129, pp. 535—539
  5. 5,0 5,1 Шаблон:Статья
  6. Шаблон:Статья