[物理学与PDEs]第5章第5节 弹性动力学方程组及其数学结构

简介: 5.5.1 线性弹性动力学方程组     1.  线性弹性动力学方程组 $$\beex \bea 0&=\rho_0\cfrac{\p{\bf v}}{\p t}-\Div_x{\bf P}-\rho_0{\bf b}\\ &=\rho_0\cfrac{\p}{\p t}\sex{\cfra...

5.5.1 线性弹性动力学方程组

 

 

1.  线性弹性动力学方程组 \beex \bea 0&=\rho_0\cfrac{\p{\bf v}}{\p t}-\Div_x{\bf P}-\rho_0{\bf b}\\ &=\rho_0\cfrac{\p}{\p t}\sex{\cfrac{\p{\bf u}}{\p t}} -\Div_x({\bf A}{\bf E})-\rho_0{\bf b}\quad\sex{{\bf u}={\bf y}-{\bf x}}\\ &=\rho_0\cfrac{\p^2{\bf u}}{\p t^2}-\Div_x({\bf A}{\bf E})-\rho_0{\bf b}, \eea \eeex

\beex \bea 0&=\rho_0\cfrac{\p{\bf v}}{\p t}-\Div_x{\bf P}-\rho_0{\bf b}\\ &=\rho_0\cfrac{\p}{\p t}\sex{\cfrac{\p{\bf u}}{\p t}} -\Div_x({\bf A}{\bf E})-\rho_0{\bf b}\quad\sex{{\bf u}={\bf y}-{\bf x}}\\ &=\rho_0\cfrac{\p^2{\bf u}}{\p t^2}-\Div_x({\bf A}{\bf E})-\rho_0{\bf b}, \eea \eeex
其分量形式为 \bee\label{5_5_1:el} \bea \rho_0\cfrac{\p ^2u}{\p t^2} &=\cfrac{1}{2}\sum_{j,k,l}\cfrac{\p}{\p x_j} \sez{a_{ijkl}\sex{\cfrac{\p u_k}{\p x_l}+\cfrac{\p u_l}{\p x_k}}} +\rho_0b_i\\ &=\cfrac{1}{2}\sum_{j,k,l}a_{ijkl}\sez{\cfrac{\p ^2u_k}{\p x_j\p x_l} +\cfrac{\p^2u_l}{\p x_j\p x_k}}+\rho_0b_i\\ &=\sum_{j,k,l}a_{ijkl}\cfrac{\p^2u_k}{\p x_j\p x_l}+\rho_0b_i.  \eea \eee
\bee\label{5_5_1:el} \bea \rho_0\cfrac{\p ^2u}{\p t^2} &=\cfrac{1}{2}\sum_{j,k,l}\cfrac{\p}{\p x_j} \sez{a_{ijkl}\sex{\cfrac{\p u_k}{\p x_l}+\cfrac{\p u_l}{\p x_k}}} +\rho_0b_i\\ &=\cfrac{1}{2}\sum_{j,k,l}a_{ijkl}\sez{\cfrac{\p ^2u_k}{\p x_j\p x_l} +\cfrac{\p^2u_l}{\p x_j\p x_k}}+\rho_0b_i\\ &=\sum_{j,k,l}a_{ijkl}\cfrac{\p^2u_k}{\p x_j\p x_l}+\rho_0b_i.  \eea \eee

 

 

2.  四阶张量 A=(aijkl) 满足强椭圆性条件, 是指 \bex α>0,\sti,j,k,laijklξiξkηjηlα|ξ|2|η|2, ξ,η\bbR3.\eex

A 满足强椭圆性条件, 则称 (???) 为二阶双曲型方程组.

 

 

3.  对各向同性材料, A 满足强椭圆性条件 \lra \bexμ>0,\lm+2μ>0.\eex

 

 

4.  Cauchy 问题、初边值问题的提法 (给定边界上的位移 u 或应力向量 (Pn)i=jklaijkl\puk\pxlnj).

 

 

5.  各向同性材料时的线性弹性动力学方程组 \bex\sedd\barl\p2u\pt2=μ\lapu+(\lm+μ)\n\Divu,u(0)=u0,\pu\pt(0)=u1.\ea\eex

 

(1)  将 u 分解为 \beeu=v+w,\rotv=0,\Divw=0.\eee

v,w 分别满足 \beex \bea \sedd{\ba{rl} \cfrac{\p^2{\bf v}}{\p t^2}=a_1^2\lap{\bf v},\\ {\bf v}(0) ={\bf u}^0_L,\quad \cfrac{\p {\bf v}}{\p t}(0) ={\bf u}^1_L; \ea},&\quad\sedd{\ba{rl} \cfrac{\p ^2{\bf w}}{\p t^2}=a_2^2\lap{\bf w},\\ {\bf w}(0) ={\bf u}^0_T,\quad\cfrac{\p {\bf w}}{\p t}(0) ={\bf u}^1_T. \ea} \eea \eeex
其中 a21=\lm+2μ, a22=μ. 由于 (???) 分解的整体依赖性 (而非点依赖性), u(t,x) 依赖于 \bex\sedy; a2t|yx|a1t.\eex

 

(2)  \sex\p2\pt2a21\lap\sex\p2\pt2a22\lapu=0.

 

 

6.  稳定性条件 \bex ˜α>0,\sti,j,k,laijkleijekl˜α|E|2,\eex

对称矩阵 E=(eij) 成立.

 

(1)  稳定性条件 \ra 强椭圆性条件 (只要取 eij=12\sexξiηj+ξjηi). 反之不然.

 

(2)  对各向同性材料, 稳定性条件 \lra \bexμ>0,κ=\lm+23μ>0.\eex

 

 

 

5.5.2 非线性弹性动力学方程组

 

 

 

 

1.  P(x)=ˆP(F(x))=detFˆT(F)FT 代入动量守恒方程有 \bexρ0\p2ui\pt2=j,k,laijkl(\nu)\puk\pxj\pxl+ρ0bi,\eex

其中 \bexaijkl(F)=\ppij\pfkl.\eex

 

 

2.  强椭圆性条件: \bexi,j,k,laijklξiξkηjηl>0, F,  ξ,ηR3\bs\sed0.\eex

 

 

 

5.5.3 非线性弹性动力学方程组的一阶守恒律形式

 

 \bee\label{5_5_3_ne} \bea \cfrac{\p f_{kl}}{\p t}-\cfrac{\p v_k}{\p x_l}&=0,\\ \rho_0\cfrac{\p v_i}{\p t}-\sum_j\cfrac{\p}{\p x_j}p_{ij}({\bf F})-\rho_0b_i&=0.  \eea \eee

 

 

1.  (???) 可化为守恒律形式的一阶拟线性方程组.

 

 

2.  若材料是超弹性的, A=(aijkl) 满足强椭圆性条件, 则 (???) 为双曲型的.

 

 

3.  在解的间断面上应满足熵不等式 \bex\p\ptη(U)+j\p\pxjqj(U)0,\eex

其中 \bexη=12|v|2+ˆW(F),qj=jpijvi.\eex

 

 

 

5.5.4 化弹性动力学方程组为一阶对称双曲组

 

 

 

 

1.  当 \lm+2μ>μ>0 时, 变形在自然状态附近的各向同性材料的非线性弹性动力学方程组可化为一阶对称双曲组; 也可通过构造一附加守恒律的方法化为具守恒律的一阶对称双曲组.

 

 

2.  对一般的非线性超弹性动力学方程组, 如果贮能函数是严格多凸的, 则也可化为具守恒律的一阶对称双曲组.

 

 

 

5.5.5 一维非线性弹性动力学方程组

 

 

 

 

1.  各向同性材料的纯轴向变形 \bexρ0\p2u1\pt2=\p\px1t11\sex\pu1\px1+ρ0b1.\eex

这是一维拟线性波动方程.

 

 

2.  各向同性材料的纯剪切变形 \bexρ0\p2u1\pt=\p\px2t12\sex\pu1\px2+ρ0b1.\eex

这也是一维拟线性波动方程.

 

目录
打赏
0
0
0
0
15
分享
相关文章
[物理学与PDEs]第2章第1节 理想流体力学方程组 1.4 一维理想流体力学方程组
1.  一维理想流体力学方程组 $$\beex \bea \cfrac{\p\rho}{\p t}+\cfrac{\p}{\p x}(\rho u)&=0,\\ \cfrac{\p}{\p t}(\rho u) +\cfrac{\p}{\p x}(\rho u^2+p)&=\rho F,\\ \cf...
830 0
[物理学与PDEs]第5章第6节 弹性静力学方程组的定解问题
5. 6 弹性静力学方程组的定解问题           5. 6. 1 线性弹性静力学方程组         1.  线性弹性静力学方程组 $$\bee\label{5_6_1_le} -\sum_{j,k,l}a_{ijkl}\cfrac{\p ^2u_k}{\p x_j\p x_l}=\rho_0b_i,\quad i=1,2,3.
750 0
[物理学与PDEs]第5章习题4 广义 Hookean 定律的张量的对称性
设材料是超弹性的, 并设参考构形为自然状态, 证明由线性化得到的张量 A=(aijkl)=\sex2\pˉpijckl 具有以下的对称性: $$\bex a_{ijkl}=a_{klij}.
578 0
[物理学与PDEs]第5章习题9 伴随矩阵的特征值
3×3A 的特征值为 \lm1,\lm2,\lm3, 证明 \cofA 的特征值为 $$\bex \lm_2\lm_3,\quad \lm_3\lm_1,\quad \lm_1\lm_2.
726 0
[物理学与PDEs]第5章习题1 矩阵的极分解
证明引理 2. 1.    证明:   (1)  先证明存在正交阵 P,Q 及对角阵 D 使得 $$\bex {\bf F}={\bf P}{\bf D}{\bf Q}.
847 0
[物理学与PDEs]第4章第2节 反应流体力学方程组 2.4 反应流体力学方程组的数学结构
1.  粘性热传导反应流体力学方程组是拟线性对称双曲 - 抛物耦合组.     2.  理想反应流体力学方程组是一阶拟线性对称双曲组 (取 u,p,S,Z 为未知函数).     3.  右端项具有间断性.
710 0
[物理学与PDEs]第3章第4节 磁流体力学方程组的数学结构
1.  在流体存在粘性、热传导及 σ 时, 磁流体力学方程组是一个拟线性对称双曲 - 抛物耦合组.     2.  在流体存在粘性、热传导但 σ= 时, 磁流体力学方程组是一个拟线性对称双曲 - 抛物耦合组.
768 0
[物理学与PDEs]第5章 弹性力学
[物理学与PDEs]第5章第1节 引言   [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.1 变形梯度张量   [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.2 Cauchy - Green 应变张量   [物理学与PDEs]第5章第2节 变形的描述, 应变张量 2.
679 0
[物理学与PDEs]第4章习题3 一维理想反应流体力学方程组的数学结构
证明: Euler 坐标系下的一维反应流体力学方程组 (3. 10)-(3. 13) 也是一个一阶拟线性双曲型方程组. 证明: 由 (3. 10), (3. 12), (3. 13) 知 $$\bex \cfrac{1}{\rho c^2}\cfrac{\p p}{\p t} +\cfrac{u}{\rho c^2}\cfrac{\p p}{\p x}+\cfrac{\p u}{\p x}=0.
861 0
[物理学与PDEs]第3章习题5 一维理想磁流体力学方程组的数学结构
试将一维理想磁流体力学方程组 (5. 10)-(5. 16) 化为一阶拟线性对称双曲组的形式. 解答: 由 (5. 12),(5. 16) 知 $$\beex \bea 0&=\cfrac{\p p}{\p \rho}\sex{\cfrac{\p \rho}{\p t}+u_1\cfrac{\p ...
798 0
AI助理

你好,我是AI助理

可以解答问题、推荐解决方案等