Теория Кирхгофа — Лява

Когда деформации в пластине бесконечно малы и повороты нормалей средней поверхности меньше 10° соотношения деформация-смещение (то есть используется разложение до первого порядка малости) принимают вид

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}\left({\frac {\partial u_{\alpha }}{\partial x_{\beta }}}+{\frac {\partial u_{\beta }}{\partial x_{\alpha }}}\right)\equiv {\frac {1}{2}}(u_{\alpha ,\beta }+u_{\beta ,\alpha })\\\varepsilon _{\alpha 3}&={\frac {1}{2}}\left({\frac {\partial u_{\alpha }}{\partial x_{3}}}+{\frac {\partial u_{3}}{\partial x_{\alpha }}}\right)\equiv {\frac {1}{2}}(u_{\alpha ,3}+u_{3,\alpha })\\\varepsilon _{33}&={\frac {\partial u_{3}}{\partial x_{3}}}\equiv u_{3,3}\end{aligned}}}$

где ${\displaystyle \beta =1,2}$ как и ${\displaystyle \alpha }$ .

Используя кинематические предположения, получим

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\tfrac {1}{2}}(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0})-x_{3}~w_{,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}&=-w_{,\alpha }^{0}+w_{,\alpha }^{0}=0\\\varepsilon _{33}&=0\end{aligned}}}$

Следовательно, ненулевые деформации возникают только в плоскостях.

Уравнения равновесия пластины выводятся из принципа виртуальной работы. Для тонкой пластины при квазистатической поперечной нагрузке ${\displaystyle q(x)}$ эти уравнения имеют вид

${\displaystyle {\begin{aligned}&{\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\\&{\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\end{aligned}}}$

где толщина пластины ${\displaystyle 2h}$. В индексной записи

${\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\quad \quad N_{\alpha \beta }:=\int _{-h}^{h}\sigma _{\alpha \beta }~dx_{3}\\M_{\alpha \beta ,\alpha \beta }-q&=0\quad \quad M_{\alpha \beta }:=\int _{-h}^{h}x_{3}~\sigma _{\alpha \beta }~dx_{3}\end{aligned}}}$

где ${\displaystyle \sigma _{\alpha \beta }}$ — механические напряжения.

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

Bending moments and normal stresses	Torques and shear stresses

Вывод уравнения равновесия для малых вращений

Для ситуации, когда напряжения и вращения пластины малы, вариация внутренняя энергия записывается в виде ${\displaystyle {\begin{aligned}\delta U&=\int _{\Omega ^{0}}\int _{-h}^{h}{\boldsymbol {\sigma }}:\delta {\boldsymbol {\epsilon }}~dx_{3}~d\Omega =\int _{\Omega ^{0}}\int _{-h}^{h}\sigma _{\alpha \beta }~\delta \varepsilon _{\alpha \beta }~dx_{3}~d\Omega \\&=\int _{\Omega ^{0}}\int _{-h}^{h}\left[{\frac {1}{2}}~\sigma _{\alpha \beta }~(\delta u_{\alpha ,\beta }^{0}+\delta u_{\beta ,\alpha }^{0})-x_{3}~\sigma _{\alpha \beta }~\delta w_{,\alpha \beta }^{0}\right]~dx_{3}~d\Omega \\&=\int _{\Omega ^{0}}\left[{\frac {1}{2}}~N_{\alpha \beta }~(\delta u_{\alpha ,\beta }^{0}+\delta u_{\beta ,\alpha }^{0})-M_{\alpha \beta }~\delta w_{,\alpha \beta }^{0}\right]~d\Omega \end{aligned}}}$ где толщина пластины ${\displaystyle 2h}$ и усилия и моменты определяются как ${\displaystyle N_{\alpha \beta }:=\int _{-h}^{h}\sigma _{\alpha \beta }~dx_{3}~;~~M_{\alpha \beta }:=\int _{-h}^{h}x_{3}~\sigma _{\alpha \beta }~dx_{3}}$ Интегрирование по частям приводит к ${\displaystyle {\begin{aligned}\delta U&=\int _{\Omega ^{0}}\left[-{\frac {1}{2}}~(N_{\alpha \beta ,\beta }~\delta u_{\alpha }^{0}+N_{\alpha \beta ,\alpha }~\delta u_{\beta }^{0})+M_{\alpha \beta ,\beta }~\delta w_{,\alpha }^{0}\right]~d\Omega \\&+\int _{\Gamma ^{0}}\left[{\frac {1}{2}}~(n_{\beta }~N_{\alpha \beta }~\delta u_{\alpha }^{0}+n_{\alpha }~N_{\alpha \beta }~\delta u_{\beta }^{0})-n_{\beta }~M_{\alpha \beta }~\delta w_{,\alpha }^{0}\right]~d\Gamma \end{aligned}}}$ Симметрия тентора напряжений подразумевает, что ${\displaystyle N_{\alpha \beta }=N_{\beta \alpha }}$. Отсюда ${\displaystyle \delta U=\int _{\Omega ^{0}}\left[-N_{\alpha \beta ,\alpha }~\delta u_{\beta }^{0}+M_{\alpha \beta ,\beta }~\delta w_{,\alpha }^{0}\right]~d\Omega +\int _{\Gamma ^{0}}\left[n_{\alpha }~N_{\alpha \beta }~\delta u_{\beta }^{0}-n_{\beta }~M_{\alpha \beta }~\delta w_{,\alpha }^{0}\right]~d\Gamma }$ Проинтегрируем по частям ещё раз ${\displaystyle \delta U=\int _{\Omega ^{0}}\left[-N_{\alpha \beta ,\alpha }~\delta u_{\beta }^{0}-M_{\alpha \beta ,\beta \alpha }~\delta w^{0}\right]~d\Omega +\int _{\Gamma ^{0}}\left[n_{\alpha }~N_{\alpha \beta }~\delta u_{\beta }^{0}+n_{\alpha }~M_{\alpha \beta ,\beta }~\delta w^{0}-n_{\beta }~M_{\alpha \beta }~\delta w_{,\alpha }^{0}\right]~d\Gamma }$ В случае, отсутствия внешних сил, принцип виртуальной работы подразумевает, что вариация ${\displaystyle \delta U=0}$. Уравнения равновесия для пластины задаются ${\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }&=0\end{aligned}}}$ Если пластина испытывает внешнюю распределенную нагрузку ${\displaystyle q(x)}$, которая направлена по нормали к срединной плоскости и напрвлена вдоль направления ${\displaystyle x_{3}}$, то внешняя виртуальная работа из-за нагрузки ${\displaystyle \delta V_{\mathrm {ext} }=\int _{\Omega ^{0}}q~\delta w^{0}~d\Omega }$ Принцип виртуальной работы приводит к уравнениям равновесия ${\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }-q&=0\end{aligned}}}$

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

${\displaystyle {\begin{aligned}n_{\alpha }~N_{\alpha \beta }&\quad \mathrm {or} \quad u_{\beta }^{0}\\n_{\alpha }~M_{\alpha \beta ,\beta }&\quad \mathrm {or} \quad w^{0}\\n_{\beta }~M_{\alpha \beta }&\quad \mathrm {or} \quad w_{,\alpha }^{0}\end{aligned}}}$

Обратите внимание, что ${\displaystyle n_{\alpha }~M_{\alpha \beta ,\beta }}$ — это эффективная сила сдвига.

Соотношения между напряжениями и деформациями для линейной упругой среды имеют вид

${\displaystyle {\begin{aligned}\sigma _{\alpha \beta }&=C_{\alpha \beta \gamma \theta }~\varepsilon _{\gamma \theta }\\\sigma _{\alpha 3}&=C_{\alpha 3\gamma \theta }~\varepsilon _{\gamma \theta }\\\sigma _{33}&=C_{33\gamma \theta }~\varepsilon _{\gamma \theta }\end{aligned}}}$

поскольку ${\displaystyle \sigma _{\alpha 3}}$, а также ${\displaystyle \sigma _{33}}$ не фигурируют в уравнениях равновесия, то неявно предполагается, что эти величины не влияют на баланс импульса и ими можно пренебречь. Остальные соотношения напряжение-деформация в матричной имеют вид

${\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}={\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}}$

Тогда

${\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}=\int _{-h}^{h}{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}dx_{3}=\left\{\int _{-h}^{h}{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}u_{1,1}^{0}\\u_{2,2}^{0}\\{\frac {1}{2}}~(u_{1,2}^{0}+u_{2,1}^{0})\end{bmatrix}}}$

и

${\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=\int _{-h}^{h}x_{3}~{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}dx_{3}=-\left\{\int _{-h}^{h}x_{3}^{2}~{\begin{bmatrix}C_{11}&C_{12}&C_{13}\\C_{12}&C_{22}&C_{23}\\C_{13}&C_{23}&C_{33}\end{bmatrix}}~dx_{3}\right\}{\begin{bmatrix}w_{,11}^{0}\\w_{,22}^{0}\\w_{,12}^{0}\end{bmatrix}}}$

Жесткости — это величины

${\displaystyle A_{\alpha \beta }:=\int _{-h}^{h}C_{\alpha \beta }~dx_{3}}$

Изгибные жесткости (также называемая жесткостью на изгиб) — это величины

${\displaystyle D_{\alpha \beta }:=\int _{-h}^{h}x_{3}^{2}~C_{\alpha \beta }~dx_{3}}$

Основные предположения Кирхгофа — Лява приводят к нулевым поперечным силам. Тогда для определения поперечных сил в тонких пластинах Кирхгофа — Лява должны использоваться уравнения равновесия пластины. Для изотропных пластин эти уравнения приводят к выражению

${\displaystyle Q_{\alpha }=-D{\frac {\partial }{\partial x_{\alpha }}}(\nabla ^{2}w^{0})\,.}$

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

${\displaystyle Q_{\alpha }={\mathcal {M}}_{,\alpha }}$

где

${\displaystyle {\mathcal {M}}:=-D\nabla ^{2}w^{0}\,.}$

Если повороты нормалей к срединной поверхности находятся в диапазоне от 10${\displaystyle ^{\circ }}$ до 15${\displaystyle ^{\circ }}$, то зависимости деформации от смещения можно аппроксимировать как

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\tfrac {1}{2}}(u_{\alpha ,\beta }+u_{\beta ,\alpha }+u_{3,\alpha }~u_{3,\beta })\\\varepsilon _{\alpha 3}&={\tfrac {1}{2}}(u_{\alpha ,3}+u_{3,\alpha })\\\varepsilon _{33}&=u_{3,3}\end{aligned}}}$

Тогда кинематические допущения теории Кирхгофа — Лява приводят к классической теории пластин с деформациями фон Кармана.

${\displaystyle {\begin{aligned}\varepsilon _{\alpha \beta }&={\frac {1}{2}}(u_{\alpha ,\beta }^{0}+u_{\beta ,\alpha }^{0}+w_{,\alpha }^{0}~w_{,\beta }^{0})-x_{3}~w_{,\alpha \beta }^{0}\\\varepsilon _{\alpha 3}&=-w_{,\alpha }^{0}+w_{,\alpha }^{0}=0\\\varepsilon _{33}&=0\end{aligned}}}$

Эта теория нелинейна из-за квадратичных членов в соотношениях деформация-перемещение.

Если соотношения деформация-перемещение принимают форму фон Кармана, то уравнения равновесия перепишутся в виде

${\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\\M_{\alpha \beta ,\alpha \beta }+[N_{\alpha \beta }~w_{,\beta }^{0}]_{,\alpha }-q&=0\end{aligned}}}$

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

${\displaystyle {\frac {\partial ^{4}w^{0}}{\partial x_{1}^{4}}}+2{\frac {\partial ^{4}w^{0}}{\partial x_{1}^{2}\partial x_{2}^{2}}}+{\frac {\partial ^{4}w^{0}}{\partial x_{2}^{4}}}=0\,.}$

Здесь мы предположили, что смещения в плоскости не зависят от ${\displaystyle x_{1}}$ и ${\displaystyle x_{2}}$. В индексной записи

${\displaystyle w_{,1111}^{0}+2~w_{,1212}^{0}+w_{,2222}^{0}=0}$

и в прямой записи

${\displaystyle \nabla ^{2}\nabla ^{2}w=0}$

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

${\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}w_{,11}^{0}\\w_{,22}^{0}\\w_{,12}^{0}\end{bmatrix}}}$

Вывод уравнений равновесия для чистого изгиба

Для изотропных, однородных пластин под действием чистого изгиба основные уравнения ${\displaystyle {\begin{aligned}N_{\alpha \beta ,\alpha }&=0\implies N_{11,1}+N_{21,2}=0~,~~N_{12,1}+N_{22,2}=0\\M_{\alpha \beta ,\alpha \beta }&=0\implies M_{11,11}+2M_{12,12}+M_{22,22}=0\end{aligned}}}$ и соотношения напряжения-деформации ${\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{12}\end{bmatrix}}={\cfrac {E}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{12}\end{bmatrix}}}$ Тогда ${\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}={\cfrac {2hE}{(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}u_{1,1}^{0}\\u_{2,2}^{0}\\{\frac {1}{2}}~(u_{1,2}^{0}+u_{2,1}^{0})\end{bmatrix}}}$ и ${\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}w_{,11}^{0}\\w_{,22}^{0}\\w_{,12}^{0}\end{bmatrix}}}$ Дифференцирование приводит к ${\displaystyle {\begin{aligned}N_{11,1}&={\cfrac {2hE}{(1-\nu ^{2})}}\left(u_{1,11}^{0}+\nu ~u_{2,21}^{0}\right)~;~~N_{22,2}={\cfrac {2hE}{(1-\nu ^{2})}}\left(\nu ~u_{1,12}^{0}+u_{2,22}^{0}\right)\\N_{12,1}&={\cfrac {hE(1-\nu )}{(1-\nu ^{2})}}\left(u_{1,21}^{0}+u_{2,11}^{0}\right)~;~~N_{12,2}={\cfrac {hE(1-\nu )}{(1-\nu ^{2})}}\left(u_{1,22}^{0}+u_{2,12}^{0}\right)\end{aligned}}}$ и ${\displaystyle {\begin{aligned}M_{11,11}&=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}\left(w_{,1111}^{0}+\nu ~w_{,2211}^{0}\right)\\M_{22,22}&=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}\left(\nu ~w_{,1122}^{0}+w_{,2222}^{0}\right)\\M_{12,12}&=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}(1-\nu )~w_{,1212}^{0}\end{aligned}}}$ Подставим результат в основные уравнение, получим ${\displaystyle {\begin{aligned}&u_{1,11}^{0}+\nu ~u_{2,21}^{0}+{\tfrac {1}{2}}(1-\nu )\left(u_{1,22}^{0}+u_{2,12}^{0}\right)=0\\&\nu ~u_{1,12}^{0}+u_{2,22}^{0}+{\tfrac {1}{2}}(1-\nu )\left(u_{1,21}^{0}+u_{2,11}^{0}\right)=0\\&w_{,1111}^{0}+\nu ~w_{,2211}^{0}+2(1-\nu )~w_{,1212}^{0}+\nu ~w_{,1122}^{0}+w_{,2222}^{0}=0\end{aligned}}}$ Поскольку порядок дифференциации не имеет значения, то ${\displaystyle u_{1,12}^{0}=u_{1,21}^{0}}$, ${\displaystyle u_{2,21}^{0}=u_{2,12}^{0}}$, и${\displaystyle w_{,2211}^{0}=w_{,1212}^{0}=w_{,1122}^{0}}$. Отсюда ${\displaystyle {\begin{aligned}&u_{1,11}^{0}+{\tfrac {1}{2}}(1-\nu )~u_{1,22}^{0}+{\tfrac {1}{2}}(1+\nu )~u_{2,12}^{0}=0\\&u_{2,22}^{0}+{\tfrac {1}{2}}(1-\nu )~u_{2,11}^{0}+{\tfrac {1}{2}}(1+\nu )~u_{1,12}^{0}=0\\&w_{,1111}^{0}+2~w_{,1212}^{0}+w_{,2222}^{0}=0\end{aligned}}}$ В прямой тензорной нотации, основное уравнение для пластины ${\displaystyle \nabla ^{2}\nabla ^{2}w=0}$ где мы предположили, что перемещения ${\displaystyle u_{1}^{0},u_{2}^{0}}$ постоянны.

Если распределенная поперечная нагрузка ${\displaystyle -q(x)}$ применяется к пластине, то определяющее уравнение ${\displaystyle M_{\alpha \beta ,\alpha \beta }=-q}$ . Следуя процедуре, из предыдущего раздела, получаем[3]

${\displaystyle \nabla ^{2}\nabla ^{2}w={\cfrac {q}{D}}~;~~D:={\cfrac {2h^{3}E}{3(1-\nu ^{2})}}}$

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

${\displaystyle w_{,1111}^{0}+2\,w_{,1212}^{0}+w_{,2222}^{0}=-{\cfrac {q}{D}}}$

а в цилиндрических координатах принимает вид (для круглой пластинки с аксиально-симметричной нагрузкой)

${\displaystyle {\frac {1}{r}}{\cfrac {d}{dr}}\left[r{\cfrac {d}{dr}}\left\{{\frac {1}{r}}{\cfrac {d}{dr}}\left(r{\cfrac {dw}{dr}}\right)\right\}\right]=-{\frac {q}{D}}\,.}$

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

Вывод уравнений равновесия для поперечной нагрузки

Для поперечно нагруженной пластины без аксиальных деформаций, основное уравнение примет вид ${\displaystyle M_{\alpha \beta ,\alpha \beta }=q\implies M_{11,11}+2M_{12,12}+M_{22,22}=q}$ где ${\displaystyle q}$ распределенная поперечная нагрузка (на единицу площади). Замена выражений на производные ${\displaystyle M_{\alpha \beta }}$ в основном уравнении приводит к ${\displaystyle -{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}\left[w_{,1111}^{0}+2\,w_{,1212}^{0}+w_{,2222}^{0}\right]=q\,.}$ Используя для изгибной жёсткости выражение ${\displaystyle D:={\cfrac {2h^{3}E}{3(1-\nu ^{2})}}}$ запишем основное уравнение в виде ${\displaystyle \nabla ^{2}\nabla ^{2}w=-{\frac {q}{D}}\,.}$ В цилиндрических координатах ${\displaystyle (r,\theta ,z)}$, ${\displaystyle \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}}}\,.}$ Для аксиально-симметричной нагрузки и круглых пластин, ${\displaystyle w=w(r)}$, тогда ${\displaystyle \nabla ^{2}w\equiv {\frac {1}{r}}{\cfrac {d}{dr}}\left(r{\cfrac {dw}{dr}}\right)\,.}$

При определённых условиях нагружения плоская пластина может изгибаться, принимая форму поверхности цилиндра. Этот тип изгиба называется цилиндрическим изгибом и представляет собой особую ситуацию, когда ${\displaystyle u_{1}=u_{1}(x_{1}),u_{2}=0,w=w(x_{1})}$. В таком случае

${\displaystyle {\begin{bmatrix}N_{11}\\N_{22}\\N_{12}\end{bmatrix}}={\cfrac {2hE}{(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}u_{1,1}^{0}\\0\\0\end{bmatrix}}}$

а также

${\displaystyle {\begin{bmatrix}M_{11}\\M_{22}\\M_{12}\end{bmatrix}}=-{\cfrac {2h^{3}E}{3(1-\nu ^{2})}}~{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}w_{,11}^{0}\\0\\0\end{bmatrix}}}$

и определяющие уравнения становятся к[3]

${\displaystyle {\begin{aligned}N_{11}&=A~{\cfrac {\mathrm {d} u}{\mathrm {d} x_{1}}}\quad \implies \quad {\cfrac {\mathrm {d} ^{2}u}{\mathrm {d} x_{1}^{2}}}=0\\M_{11}&=-D~{\cfrac {\mathrm {d} ^{2}w}{\mathrm {d} x_{1}^{2}}}\quad \implies \quad {\cfrac {\mathrm {d} ^{4}w}{\mathrm {d} x_{1}^{4}}}={\cfrac {q}{D}}\\\end{aligned}}}$

Основные уравнения динамики пластины Кирхгофа — Лява:

${\displaystyle {\begin{aligned}N_{\alpha \beta ,\beta }&=J_{1}~{\ddot {u}}_{\alpha }^{0}\\M_{\alpha \beta ,\alpha \beta }+q(x,t)&=J_{1}~{\ddot {w}}^{0}-J_{3}~{\ddot {w}}_{,\alpha \alpha }^{0}\end{aligned}}}$

где для пластины с плотностью ${\displaystyle \rho =\rho (x)}$ ,

${\displaystyle J_{1}:=\int _{-h}^{h}\rho ~dx_{3}=2~\rho ~h~;~~J_{3}:=\int _{-h}^{h}x_{3}^{2}~\rho ~dx_{3}={\frac {2}{3}}~\rho ~h^{3}}$

а также

${\displaystyle {\dot {u}}_{i}={\frac {\partial u_{i}}{\partial t}}~;~~{\ddot {u}}_{i}={\frac {\partial ^{2}u_{i}}{\partial t^{2}}}~;~~u_{i,\alpha }={\frac {\partial u_{i}}{\partial x_{\alpha }}}~;~~u_{i,\alpha \beta }={\frac {\partial ^{2}u_{i}}{\partial x_{\alpha }\partial x_{\beta }}}}$

Вывод уравнений, регулирующих динамику пластин Кирхгофа — Лява

Полная кинетическая энергия пластины ${\displaystyle K=\int _{0}^{T}\int _{\Omega ^{0}}\int _{-h}^{h}{\cfrac {\rho }{2}}\left[\left({\frac {\partial u_{1}}{\partial t}}\right)^{2}+\left({\frac {\partial u_{2}}{\partial t}}\right)^{2}+\left({\frac {\partial u_{3}}{\partial t}}\right)^{2}\right]~\mathrm {d} x_{3}~\mathrm {d} A~\mathrm {d} t}$ Таким образом, вариация кинетической энергии ${\displaystyle \delta K=\int _{0}^{T}\int _{\Omega ^{0}}\int _{-h}^{h}{\cfrac {\rho }{2}}\left[2\left({\frac {\partial u_{1}}{\partial t}}\right)\left({\frac {\partial \delta u_{1}}{\partial t}}\right)+2\left({\frac {\partial u_{2}}{\partial t}}\right)\left({\frac {\partial \delta u_{2}}{\partial t}}\right)+2\left({\frac {\partial u_{3}}{\partial t}}\right)\left({\frac {\partial \delta u_{3}}{\partial t}}\right)\right]~\mathrm {d} x_{3}~\mathrm {d} A~\mathrm {d} t}$ Тут мы используем следующую нотацию ${\displaystyle {\dot {u}}_{i}={\frac {\partial u_{i}}{\partial t}}~;~~{\ddot {u}}_{i}={\frac {\partial ^{2}u_{i}}{\partial t^{2}}}~;~~u_{i,\alpha }={\frac {\partial u_{i}}{\partial x_{\alpha }}}~;~~u_{i,\alpha \beta }={\frac {\partial ^{2}u_{i}}{\partial x_{\alpha }\partial x_{\beta }}}}$ Тогда ${\displaystyle \delta K=\int _{0}^{T}\int _{\Omega ^{0}}\int _{-h}^{h}\rho \left({\dot {u}}_{\alpha }~\delta {\dot {u}}_{\alpha }+{\dot {u}}_{3}~\delta {\dot {u}}_{3}\right)~\mathrm {d} x_{3}~\mathrm {d} A~\mathrm {d} t}$ Для пластин Кирхгофа — Лява ${\displaystyle u_{\alpha }=u_{\alpha }^{0}-x_{3}~w_{,\alpha }^{0}~;~~u_{3}=w^{0}}$ Отсюда ${\displaystyle {\begin{aligned}\delta K&=\int _{0}^{T}\int _{\Omega ^{0}}\int _{-h}^{h}\rho \left[\left({\dot {u}}_{\alpha }^{0}-x_{3}~{\dot {w}}_{,\alpha }^{0}\right)~\left(\delta {\dot {u}}_{\alpha }^{0}-x_{3}~\delta {\dot {w}}_{,\alpha }^{0}\right)+{\dot {w}}^{0}~\delta {\dot {w}}^{0}\right]~\mathrm {d} x_{3}~\mathrm {d} A~\mathrm {d} t\\&=\int _{0}^{T}\int _{\Omega ^{0}}\int _{-h}^{h}\rho \left({\dot {u}}_{\alpha }^{0}~\delta {\dot {u}}_{\alpha }^{0}-x_{3}~{\dot {w}}_{,\alpha }^{0}~\delta {\dot {u}}_{\alpha }^{0}-x_{3}~{\dot {u}}_{\alpha }^{0}~\delta {\dot {w}}_{,\alpha }^{0}+x_{3}^{2}~{\dot {w}}_{,\alpha }^{0}~\delta {\dot {w}}_{,\alpha }^{0}+{\dot {w}}^{0}~\delta {\dot {w}}^{0}\right)~\mathrm {d} x_{3}~\mathrm {d} A~\mathrm {d} t\end{aligned}}}$ Определим для постоянной по толщине ${\displaystyle \rho }$ ${\displaystyle J_{1}:=\int _{-h}^{h}\rho ~dx_{3}=2~\rho ~h~;~~J_{2}:=\int _{-h}^{h}x_{3}~\rho ~dx_{3}=0~;~~J_{3}:=\int _{-h}^{h}x_{3}^{2}~\rho ~dx_{3}={\frac {2}{3}}~\rho ~h^{3}}$ Тогда ${\displaystyle \delta K=\int _{0}^{T}\int _{\Omega ^{0}}\left[J_{1}\left({\dot {u}}_{\alpha }^{0}~\delta {\dot {u}}_{\alpha }^{0}+{\dot {w}}^{0}~\delta {\dot {w}}^{0}\right)+J_{3}~{\dot {w}}_{,\alpha }^{0}~\delta {\dot {w}}_{,\alpha }^{0}\right]~\mathrm {d} A~\mathrm {d} t}$ Интегрирование по частям даёт ${\displaystyle \delta K=\int _{\Omega ^{0}}\left[\int _{0}^{T}\left\{-J_{1}\left({\ddot {u}}_{\alpha }^{0}~\delta u_{\alpha }^{0}+{\ddot {w}}^{0}~\delta w^{0}\right)-J_{3}~{\ddot {w}}_{,\alpha }^{0}~\delta w_{,\alpha }^{0}\right\}~\mathrm {d} t+\left|J_{1}\left({\dot {u}}_{\alpha }^{0}~\delta u_{\alpha }^{0}+{\dot {w}}^{0}~\delta w^{0}\right)+J_{3}~{\dot {w}}_{,\alpha }^{0}~\delta w_{,\alpha }^{0}\right|_{0}^{T}\right]~\mathrm {d} A}$ Вариации ${\displaystyle \delta u_{\alpha }^{0}}$ и ${\displaystyle \delta w^{0}}$ равны нулю при ${\displaystyle t=0}$ и ${\displaystyle t=T}$. Таким образом, после перемены последовательности интегрирования ${\displaystyle \delta K=-\int _{0}^{T}\left\{\int _{\Omega ^{0}}\left[J_{1}\left({\ddot {u}}_{\alpha }^{0}~\delta u_{\alpha }^{0}+{\ddot {w}}^{0}~\delta w^{0}\right)+J_{3}~{\ddot {w}}_{,\alpha }^{0}~\delta w_{,\alpha }^{0}\right]~\mathrm {d} A\right\}~\mathrm {d} t+\left|\int _{\Omega ^{0}}J_{3}~{\dot {w}}_{,\alpha }^{0}~\delta w_{,\alpha }^{0}\mathrm {d} A\right|_{0}^{T}}$ Интеграция по частям в срединной поверхности даёт ${\displaystyle {\begin{aligned}\delta K&=-\int _{0}^{T}\left\{\int _{\Omega ^{0}}\left[J_{1}\left({\ddot {u}}_{\alpha }^{0}~\delta u_{\alpha }^{0}+{\ddot {w}}^{0}~\delta w^{0}\right)-J_{3}~{\ddot {w}}_{,\alpha \alpha }^{0}~\delta w^{0}\right]~\mathrm {d} A+\int _{\Gamma ^{0}}J_{3}~n_{\alpha }~{\ddot {w}}_{,\alpha }^{0}~\delta w^{0}~\mathrm {d} s\right\}~\mathrm {d} t\\&\qquad -\left|\int _{\Omega ^{0}}J_{3}~{\dot {w}}_{,\alpha \alpha }^{0}~\delta w^{0}~\mathrm {d} A-\int _{\Gamma ^{0}}J_{3}~{\dot {w}}_{,\alpha }^{0}~\delta w^{0}~\mathrm {d} s\right|_{0}^{T}\end{aligned}}}$ Опять же, поскольку вариации остаются нулевыми в начале и в конце промежутка времени, то ${\displaystyle \delta K=-\int _{0}^{T}\left\{\int _{\Omega ^{0}}\left[J_{1}\left({\ddot {u}}_{\alpha }^{0}~\delta u_{\alpha }^{0}+{\ddot {w}}^{0}~\delta w^{0}\right)-J_{3}~{\ddot {w}}_{,\alpha \alpha }^{0}~\delta w^{0}\right]~\mathrm {d} A+\int _{\Gamma ^{0}}J_{3}~n_{\alpha }~{\ddot {w}}_{,\alpha }^{0}~\delta w^{0}~\mathrm {d} s\right\}~\mathrm {d} t}$ Для динамического случая вариация внутренней энергии ${\displaystyle \delta U=-\int _{0}^{T}\left\{\int _{\Omega ^{0}}\left[N_{\alpha \beta ,\alpha }~\delta u_{\beta }^{0}+M_{\alpha \beta ,\beta \alpha }~\delta w^{0}\right]~\mathrm {d} A-\int _{\Gamma ^{0}}\left[n_{\alpha }~N_{\alpha \beta }~\delta u_{\beta }^{0}+n_{\alpha }~M_{\alpha \beta ,\beta }~\delta w^{0}-n_{\beta }~M_{\alpha \beta }~\delta w_{,\alpha }^{0}\right]~\mathrm {d} s\right\}\mathrm {d} t}$ Интеграция по частям и предположение о нулевой вариации на границе срединной поверхности дает ${\displaystyle \delta U=-\int _{0}^{T}\left\{\int _{\Omega ^{0}}\left[N_{\alpha \beta ,\alpha }~\delta u_{\beta }^{0}+M_{\alpha \beta ,\beta \alpha }~\delta w^{0}\right]~\mathrm {d} A-\int _{\Gamma ^{0}}\left[n_{\alpha }~N_{\alpha \beta }~\delta u_{\beta }^{0}+n_{\alpha }~M_{\alpha \beta ,\beta }~\delta w^{0}+n_{\beta }~M_{\alpha \beta ,\alpha }~\delta w^{0}\right]~\mathrm {d} s\right\}\mathrm {d} t}$ Если имеется внешняя распределенная сила ${\displaystyle q(x,t)}$ действуя по нормали к поверхности пластины, то виртуальная внешняя работа ${\displaystyle \delta V_{\mathrm {ext} }=\int _{0}^{T}\left[\int _{\Omega ^{0}}q(x,t)~\delta w^{0}~\mathrm {d} A\right]\mathrm {d} t}$ Из принципа виртуальной работы ${\displaystyle \delta U+\delta V_{\mathrm {ext} }=\delta K}$. Таким образом, основные уравнения баланса для пластины ${\displaystyle {\begin{aligned}N_{\alpha \beta ,\beta }&=J_{1}~{\ddot {u}}_{\alpha }^{0}\\M_{\alpha \beta ,\alpha \beta }-q(x,t)&=J_{1}~{\ddot {w}}^{0}-J_{3}~{\ddot {w}}_{,\alpha \alpha }^{0}\end{aligned}}}$