【小呆的力学笔记】弹塑性力学的初步认知二:应力分析(1)

文章目录

        • 1.1 一点的应力状态
        • 1.2 一点主应力状态
        • 1.3 应力偏张量、球张量、应力不变量

1.1 一点的应力状态

物体在受到外力或者自身不均匀的温度场等作用时,在其内部会产生内力,物体的内力与方向和截面都有关系。假设有一个受到外力作用的变形体,被一个平面截成A、B两个部分,B部分对A部分施加有作用力,在该截面上的dS微小面积上,作用力为dP,那么我们成dP与dS的比例极限为应力,如下式
σ = lim ⁡ Δ S → 0 d P d S (1) \boldsymbol\sigma=\lim_{\Delta S\to0} \frac{d\mathbf P}{dS}\tag{1} σ=ΔS0limdSdP(1)
上式黑体表明是方向相关量。
在这里插入图片描述
应力是有方向,我们规定垂直于截面的分量成为正应力,平行于截面的分量为剪应力。
为了分析一点的应力状态,从物体内任一点取一微小的四面体单元,如下图。其中 σ x x \sigma_{xx} σxx σ y y \sigma_{yy} σyy σ z z \sigma_{zz} σzz为正应力, τ x y \tau_{xy} τxy τ x z \tau_{xz} τxz τ y z \tau_{yz} τyz τ y x \tau_{yx} τyx τ z x \tau_{zx} τzx τ z y \tau_{zy} τzy为剪应力,在很多时候我们借鉴矩阵的应用,将这些分量放在一起来表示一点的应力状态,比如如下形式。(实际上,更方便的表示一点的应力状态就是应力张量,后面在专门学习张量的时候引入张量的概念和计算,应力张量是二阶张量,具有矩阵的显像化形式)
写成矩阵形式,如下式
[ σ x x τ x y τ x z τ y x σ y y τ y z τ z x τ z y σ z z ] \begin{bmatrix} \sigma_{xx} & \tau_{xy} & \tau_{xz}\\ \tau_{yx} & \sigma_{yy} & \tau_{yz}\\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \end{bmatrix} σxxτyxτzxτxyσyyτzyτxzτyzσzz
在这里插入图片描述

1.2 一点主应力状态

假设该斜面上只有正应力分量,没有剪应力,即四面体的斜面上 σ v \boldsymbol\sigma_{v} σv为正应力,其中 σ v \boldsymbol\sigma_{v} σv的应力方向余弦为 ( l , m , n ) (l,m,n) (l,m,n),根据力的平衡,可得下式。
σ x x ⋅ S Δ B O C + τ z x ⋅ S Δ A O B + τ y x ⋅ S Δ A O C = σ v x S Δ A B C σ y y ⋅ S Δ A O C + τ z y ⋅ S Δ A O B + τ x y ⋅ S Δ B O C = σ v y S Δ A B C σ z z ⋅ S Δ A O B + τ y z ⋅ S Δ A O C + τ x z ⋅ S Δ B O C = σ v z S Δ A B C (2) \sigma_{xx}\cdot S_{\Delta BOC}+\tau_{zx}\cdot S_{\Delta AOB}+\tau_{yx}\cdot S_{\Delta AOC}=\sigma_{vx}S_{\Delta ABC}\\ \sigma_{yy}\cdot S_{\Delta AOC}+\tau_{zy}\cdot S_{\Delta AOB}+\tau_{xy}\cdot S_{\Delta BOC}=\sigma_{vy}S_{\Delta ABC}\\ \sigma_{zz}\cdot S_{\Delta AOB}+\tau_{yz}\cdot S_{\Delta AOC}+\tau_{xz}\cdot S_{\Delta BOC}=\sigma_{vz}S_{\Delta ABC}\tag{2} σxxSΔBOC+τzxSΔAOB+τyxSΔAOC=σvxSΔABCσyySΔAOC+τzySΔAOB+τxySΔBOC=σvySΔABCσzzSΔAOB+τyzSΔAOC+τxzSΔBOC=σvzSΔABC(2)
其中
S Δ A O C = 1 2 d x d z = S Δ A B C ⋅ m S Δ A O B = 1 2 d x d y = S Δ A B C ⋅ n S Δ B O C = 1 2 d y d z = S Δ A B C ⋅ l (3) S_{\Delta AOC}=\frac{1}{2}dxdz=S_{\Delta ABC}\cdot m\\ S_{\Delta AOB}=\frac{1}{2}dxdy=S_{\Delta ABC}\cdot n\\ S_{\Delta BOC}=\frac{1}{2}dydz=S_{\Delta ABC}\cdot l \tag{3} SΔAOC=21dxdz=SΔABCmSΔAOB=21dxdy=SΔABCnSΔBOC=21dydz=SΔABCl(3)
那么力平衡方程改为
σ x x ⋅ l + τ z x ⋅ n + τ y x ⋅ m = σ v x = σ v l σ y y ⋅ m + τ z y ⋅ n + τ x y ⋅ l = σ v y = σ v m σ z z ⋅ n + τ y z ⋅ m + τ x z ⋅ l = σ v z = σ v n (4) \sigma_{xx}\cdot l+\tau_{zx}\cdot n+\tau_{yx}\cdot m=\sigma_{vx}=\sigma_{v}l\\ \sigma_{yy}\cdot m+\tau_{zy}\cdot n+\tau_{xy}\cdot l=\sigma_{vy}=\sigma_{v}m\\ \sigma_{zz}\cdot n+\tau_{yz}\cdot m+\tau_{xz}\cdot l=\sigma_{vz}=\sigma_{v}n \tag{4} σxxl+τzxn+τyxm=σvx=σvlσyym+τzyn+τxyl=σvy=σvmσzzn+τyzm+τxzl=σvz=σvn(4)
合并同类相,将其改为
( σ x x − σ v ) ⋅ l + τ y x ⋅ m + τ z x ⋅ n = 0 τ x y ⋅ l + ( σ y y − σ v ) ⋅ m + τ z y ⋅ n = 0 τ x z ⋅ l + τ y z ⋅ m + ( σ z z − σ v ) ⋅ n = 0 (5) (\sigma_{xx}-\sigma_{v})\cdot l+\tau_{yx}\cdot m+\tau_{zx}\cdot n=0\\ \tau_{xy}\cdot l+(\sigma_{yy}-\sigma_{v})\cdot m+\tau_{zy}\cdot n=0\\ \tau_{xz}\cdot l+\tau_{yz}\cdot m+(\sigma_{zz}-\sigma_{v})\cdot n=0 \tag{5} (σxxσv)l+τyxm+τzxn=0τxyl+(σyyσv)m+τzyn=0τxzl+τyzm+(σzzσv)n=0(5)
写成矩阵形式,如下式
[ σ x x − σ v τ y x τ z x τ x y σ y y − σ v τ z y τ x z τ y z σ z z − σ v ] [ l m n ] = [ 0 0 0 ] (6) \begin{bmatrix} \sigma_{xx}-\sigma_{v} & \tau_{yx} & \tau_{zx}\\ \tau_{xy} & \sigma_{yy}-\sigma_{v} & \tau_{zy}\\ \tau_{xz} & \tau_{yz} & \sigma_{zz}-\sigma_{v} \end{bmatrix} \begin{bmatrix} l \\m \\n \end{bmatrix}= \begin{bmatrix} 0 \\0 \\0 \end{bmatrix} \tag{6} σxxσvτxyτxzτyxσyyσvτyzτzxτzyσzzσv lmn = 000 (6)
方向余弦存在非零解,那么系数行列式需为零,如下所示。
∣ σ x x − σ v τ y x τ z x τ x y σ y y − σ v τ z y τ x z τ y z σ z z − σ v ∣ = 0 (7) \begin{vmatrix} \sigma_{xx}-\sigma_{v} & \tau_{yx} & \tau_{zx}\\ \tau_{xy} & \sigma_{yy}-\sigma_{v} & \tau_{zy}\\ \tau_{xz} & \tau_{yz} & \sigma_{zz}-\sigma_{v} \end{vmatrix}=0 \tag{7} σxxσvτxyτxzτyxσyyσvτyzτzxτzyσzzσv =0(7)
将其展开,如下式,并利用剪应力互等关系。
∣ σ x x − σ v τ y x τ z x τ x y σ y y − σ v τ z y τ x z τ y z σ z z − σ v ∣ = ( σ x x − σ v ) ( σ y y − σ v ) ( σ z z − σ v ) + τ y x τ z y τ x z + τ z x τ x y τ y z − ( σ y y − σ v ) τ z x τ x z − ( σ z z − σ v ) τ y x τ x y − ( σ x x − σ v ) τ z y τ y z = σ x x σ y y σ z z − ( σ x x σ y y + σ x x σ z z + σ y y σ z z ) σ v + ( σ x x + σ y y + σ z z ) σ v 2 − σ v 3 + 2 τ x z τ x y τ y z + ( τ x z 2 + τ x y 2 + τ y z 2 ) σ v − σ y y τ x z 2 − σ z z τ x y 2 − σ x x τ y z 2 = − σ v 3 + ( σ x x + σ y y + σ z z ) σ v 2 + ( τ x z 2 + τ x y 2 + τ y z 2 − σ x x σ y y − σ x x σ z z − σ y y σ z z ) σ v + σ x x σ y y σ z z + 2 τ x z τ x y τ y z − σ y y τ x z 2 − σ z z τ x y 2 − σ x x τ y z 2 = − σ v 3 + I 1 σ v 2 + I 2 σ v + I 3 (8) \begin{aligned} \begin{vmatrix} \sigma_{xx}-\sigma_{v} & \tau_{yx} & \tau_{zx}\\ \tau_{xy} & \sigma_{yy}-\sigma_{v} & \tau_{zy}\\ \tau_{xz} & \tau_{yz} & \sigma_{zz}-\sigma_{v} \end{vmatrix}=&(\sigma_{xx}-\sigma_{v})(\sigma_{yy}-\sigma_{v})(\sigma_{zz}-\sigma_{v})+\tau_{yx}\tau_{zy}\tau_{xz}+\tau_{zx}\tau_{xy}\tau_{yz}\\ &-(\sigma_{yy}-\sigma_{v})\tau_{zx}\tau_{xz}-(\sigma_{zz}-\sigma_{v})\tau_{yx}\tau_{xy}-(\sigma_{xx}-\sigma_{v})\tau_{zy}\tau_{yz}\\ =&\sigma_{xx}\sigma_{yy}\sigma_{zz}-(\sigma_{xx}\sigma_{yy}+\sigma_{xx}\sigma_{zz}+\sigma_{yy}\sigma_{zz})\sigma_{v}+(\sigma_{xx}+\sigma_{yy}+\sigma_{zz})\sigma_{v}^2\\ &-\sigma_{v}^3+2\tau_{xz}\tau_{xy}\tau_{yz}+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)\sigma_{v}-\sigma_{yy}\tau_{xz}^2-\sigma_{zz}\tau_{xy}^2-\sigma_{xx}\tau_{yz}^2\\ =&-\sigma_{v}^3+(\sigma_{xx}+\sigma_{yy}+\sigma_{zz})\sigma_{v}^2+\\ &(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2-\sigma_{xx}\sigma_{yy}-\sigma_{xx}\sigma_{zz}-\sigma_{yy}\sigma_{zz})\sigma_{v}+\\ &\sigma_{xx}\sigma_{yy}\sigma_{zz}+2\tau_{xz}\tau_{xy}\tau_{yz}-\sigma_{yy}\tau_{xz}^2-\sigma_{zz}\tau_{xy}^2-\sigma_{xx}\tau_{yz}^2\\ =&-\sigma_{v}^3+I_1\sigma_{v}^2+I_2\sigma_{v}+I_3 \end{aligned}\tag{8} σxxσvτxyτxzτyxσyyσvτyzτzxτzyσzzσv ====(σxxσv)(σyyσv)(σzzσv)+τyxτzyτxz+τzxτxyτyz(σyyσv)τzxτxz(σzzσv)τyxτxy(σxxσv)τzyτyzσxxσyyσzz(σxxσyy+σxxσzz+σyyσzz)σv+(σxx+σyy+σzz)σv2σv3+2τxzτxyτyz+(τxz2+τxy2+τyz2)σvσyyτxz2σzzτxy2σxxτyz2σv3+(σxx+σyy+σzz)σv2+(τxz2+τxy2+τyz2σxxσyyσxxσzzσyyσzz)σv+σxxσyyσzz+2τxzτxyτyzσyyτxz2σzzτxy2σxxτyz2σv3+I1σv2+I2σv+I3(8)

其中 I 1 I_1 I1 I 2 I_2 I2 I 3 I_3 I3成为应力不变量,如下所示。
I 1 = σ x x + σ y y + σ z z = t r [ σ ] I 2 = τ x z 2 + τ x y 2 + τ y z 2 − σ x x σ y y − σ x x σ z z − σ y y σ z z I 3 = σ x x σ y y σ z z + 2 τ x z τ x y τ y z − σ y y τ x z 2 − σ z z τ x y 2 − σ x x τ y z 2 = ∣ σ ∣ (9) I_1=\sigma_{xx}+\sigma_{yy}+\sigma_{zz}=tr[\boldsymbol \sigma]\\ I_2=\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2-\sigma_{xx}\sigma_{yy}-\sigma_{xx}\sigma_{zz}-\sigma_{yy}\sigma_{zz}\\ I_3=\sigma_{xx}\sigma_{yy}\sigma_{zz}+2\tau_{xz}\tau_{xy}\tau_{yz}-\sigma_{yy}\tau_{xz}^2-\sigma_{zz}\tau_{xy}^2-\sigma_{xx}\tau_{yz}^2=|\boldsymbol \sigma| \tag{9} I1=σxx+σyy+σzz=tr[σ]I2=τxz2+τxy2+τyz2σxxσyyσxxσzzσyyσzzI3=σxxσyyσzz+2τxzτxyτyzσyyτxz2σzzτxy2σxxτyz2=σ(9)

难么,公式(7)就变为
∣ σ x x − σ v τ y x τ z x τ x y σ y y − σ v τ z y τ x z τ y z σ z z − σ v ∣ = − σ v 3 + I 1 σ v 2 + I 2 σ v + I 3 = 0 (10) \begin{aligned} \begin{vmatrix} \sigma_{xx}-\sigma_{v} & \tau_{yx} & \tau_{zx}\\ \tau_{xy} & \sigma_{yy}-\sigma_{v} & \tau_{zy}\\ \tau_{xz} & \tau_{yz} & \sigma_{zz}-\sigma_{v} \end{vmatrix}=-\sigma_{v}^3+I_1\sigma_{v}^2+I_2\sigma_{v}+I_3=0 \end{aligned}\tag{10} σxxσvτxyτxzτyxσyyσvτyzτzxτzyσzzσv =σv3+I1σv2+I2σv+I3=0(10)

那么 σ v \boldsymbol\sigma_{v} σv可以通过公式(10)来求解,回代公式(6),可以解的 ( l , m , n ) (l,m,n) (l,m,n)
根据三次方程的韦达定理,有三个方程的根(也就是主应力)的和等于下式。
σ 1 + σ 2 + σ 3 = − I 1 − 1 = I 1 (11) \sigma_1+\sigma_2+\sigma_3=-\frac{I_1}{-1}=I_1\tag{11} σ1+σ2+σ3=1I1=I1(11)

其实观察上面的计算,不难发现,正应力 σ v \boldsymbol\sigma_{v} σv和方向余弦 ( l , m , n ) (l,m,n) (l,m,n)为应力矩阵的特征值和特征向量。

1.3 应力偏张量、球张量、应力不变量

下图为能反映一点的应力状态的六面体,众多的金属实验表明(在常见的应力范围内)当六面体各个面上只有相等正应力无切应力时,物体不发生塑性变形和形状变化。由此,定义了这么一种应力状态,即
[ σ m ] = [ σ m 0 0 0 σ m 0 0 0 σ m ] (12) [\boldsymbol\sigma_m]=\begin{bmatrix} \sigma_{m} & 0 & 0\\ 0 & \sigma_{m} & 0\\ 0 & 0 & \sigma_{m} \end{bmatrix}\tag{12} [σm]= σm000σm000σm (12)
σ m = σ 1 = σ 2 = σ 3 = 1 3 ( σ 1 + σ 2 + σ 3 ) = 1 3 I 1 (13) \sigma_m=\sigma_1=\sigma_2=\sigma_3=\frac{1}{3}(\sigma_1+\sigma_2+\sigma_3)=\frac{1}{3}I_1\tag{13} σm=σ1=σ2=σ3=31(σ1+σ2+σ3)=31I1(13)
在这里插入图片描述
那么将应力张量减去应力应力球张量,可得应力偏张量(同时可以确定的是[s]也是对称矩阵,由于剪应力互等)
[ s ] = [ σ x x − σ m τ x y τ x z τ y x σ y y − σ m τ y z τ z x τ z y σ z z − σ m ] = [ s x x s x y s x z s y x s y y s y z s z x s z y s z z ] (14) [\boldsymbol s]=\begin{bmatrix} \sigma_{xx}-\sigma_{m} & \tau_{xy} & \tau_{xz}\\ \tau_{yx} & \sigma_{yy}-\sigma_{m} & \tau_{yz}\\ \tau_{zx} & \tau_{zy} & \sigma_{zz}-\sigma_{m} \end{bmatrix} =\begin{bmatrix} s_{xx} & s_{xy} & s_{xz}\\ s_{yx} & s_{yy} & s_{yz}\\ s_{zx} & s_{zy} & s_{zz} \end{bmatrix}\tag{14} [s]= σxxσmτyxτzxτxyσyyσmτzyτxzτyzσzzσm = sxxsyxszxsxysyyszysxzsyzszz (14)

同样,按照1.2的过程,可以得到应力偏张量的主应力的公式如下所示。
∣ s x x − s v s y x s z x s x y s y y − s v s z y s x z s y z s z z − s v ∣ = ∣ s x x − s v s x y s x z s y x s y y − s v s y z s z x s z y s z z − s v ∣ = − s v 3 + I ^ 1 s v 2 + I ^ 2 s v + I ^ 3 = 0 (15) \begin{aligned} \begin{vmatrix} s_{xx}-s_{v} & s_{yx} & s_{zx}\\ s_{xy} & s_{yy}-s_{v} & s_{zy}\\ s_{xz} & s_{yz} & s_{zz}-s_{v} \end{vmatrix}&=\begin{vmatrix} s_{xx}-s_{v} & s_{xy} & s_{xz}\\ s_{yx} & s_{yy}-s_{v} & s_{yz}\\ s_{zx} & s_{zy} & s_{zz}-s_{v} \end{vmatrix}\\ &=-s_{v}^3+\hat I_1s_{v}^2+\hat I_2s_{v}+\hat I_3=0 \end{aligned}\tag{15} sxxsvsxysxzsyxsyysvsyzszxszyszzsv = sxxsvsyxszxsxysyysvszysxzsyzszzsv =sv3+I^1sv2+I^2sv+I^3=0(15)
其中,有应力偏张量的不变量如下所示。
I ^ 1 = s x x + s y y + s z z = σ x x − σ m + σ y y − σ m + σ z z − σ m = 0 I ^ 2 = s x z 2 + s x y 2 + s y z 2 − s x x s y y − s x x s z z − s y y s z z I ^ 3 = s x x s y y s z z + 2 s x z s x y s y z − s y y s x z 2 − s z z s x y 2 − s x x s y z 2 (16) \hat I_1=s_{xx}+s_{yy}+s_{zz}=\sigma_{xx}-\sigma_{m} + \sigma_{yy}-\sigma_{m}+ \sigma_{zz}-\sigma_{m}=0\\ \hat I_2=s_{xz}^2+s_{xy}^2+s_{yz}^2-s_{xx}s_{yy}-s_{xx}s_{zz}-s_{yy}s_{zz}\\ \hat I_3=s_{xx}s_{yy}s_{zz}+2s_{xz}s_{xy}s_{yz}-s_{yy}s_{xz}^2-s_{zz}s_{xy}^2-s_{xx}s_{yz}^2 \tag{16} I^1=sxx+syy+szz=σxxσm+σyyσm+σzzσm=0I^2=sxz2+sxy2+syz2sxxsyysxxszzsyyszzI^3=sxxsyyszz+2sxzsxysyzsyysxz2szzsxy2sxxsyz2(16)
对公式(15)进行展开合并同类项等过程,如下所示,
− s v 3 + I ^ 1 s v 2 + I ^ 2 s v + I ^ 3 = − s v 3 + I ^ 2 s v + I ^ 3 = − s v 3 + ( s x z 2 + s x y 2 + s y z 2 − s x x s y y − s x x s z z − s y y s z z ) s v + s x x s y y s z z + 2 s x z s x y s y z − s y y s x z 2 − s z z s x y 2 − s x x s y z 2 = − s v 3 + ( τ x z 2 + τ x y 2 + τ y z 2 ) s v − [ ( σ x x − σ m ) ( σ y y − σ m ) + ( σ x x − σ m ) ( σ z z − σ m ) + ( σ y y − σ m ) ( σ z z − σ m ) ] s v + ( σ x x − σ m ) ( σ y y − σ m ) ( σ z z − σ m ) + 2 τ x z τ x y τ y z − ( σ y y − σ m ) τ x z 2 − ( σ z z − σ m ) τ x y 2 − ( σ x x − σ m ) τ y z 2 = − s v 3 + ( τ x z 2 + τ x y 2 + τ y z 2 ) s v − ( σ x x σ y y + σ x x σ z z + σ y y σ z z ) s v + ( σ x x + σ y y + σ x x + σ z z + σ y y + σ z z ) σ m s v − 3 σ m 2 s v + σ x x σ y y σ z z − ( σ x x σ y y + σ x x σ z z + σ y y σ z z ) σ m + ( σ x x + σ y y + σ z z ) σ m 2 − σ m 3 + 2 τ x z τ x y τ y z − ( σ y y τ x z 2 + σ z z τ x y 2 + σ x x τ y z 2 ) + ( τ x z 2 + τ x y 2 + τ y z 2 ) σ m = − s v 3 + ( τ x z 2 + τ x y 2 + τ y z 2 ) ( s v + σ m ) − ( σ x x σ y y + σ x x σ z z + σ y y σ z z ) ( s v + σ m ) + 6 σ m 2 s v − 3 σ m 2 s v + σ x x σ y y σ z z + 3 σ m 3 − σ m 3 + 2 τ x z τ x y τ y z − ( σ y y τ x z 2 + σ z z τ x y 2 + σ x x τ y z 2 ) = − s v 3 + 3 σ m 2 s v + 2 σ m 3 + ( τ x z 2 + τ x y 2 + τ y z 2 − σ x x σ y y − σ x x σ z z − σ y y σ z z ) ( s v + σ m ) + σ x x σ y y σ z z + 2 τ x z τ x y τ y z − ( σ y y τ x z 2 + σ z z τ x y 2 + σ x x τ y z 2 ) = − s v 3 − σ m 3 + 3 σ m 2 s v + 3 σ m 3 + I 2 ( s v + σ m ) + I 3 = − ( s v + σ m ) ( s v 2 − s v σ m + σ m 2 − 3 σ m 2 ) + I 2 ( s v + σ m ) + I 3 = − ( s v + σ m ) ( s v + σ m ) ( s v − 2 σ m ) + I 2 ( s v + σ m ) + I 3 = − ( s v + σ m ) 2 ( s v + σ m − 3 σ m ) + I 2 ( s v + σ m ) + I 3 = − ( s v + σ m ) 3 + 3 σ m ( s v + σ m ) 2 + I 2 ( s v + σ m ) + I 3 = − ( s v + σ m ) 3 + I 1 ( s v + σ m ) 2 + I 2 ( s v + σ m ) + I 3 = 0 (17) \begin{aligned} -s_{v}^3+\hat I_1s_{v}^2+\hat I_2s_{v}+\hat I_3&=-s_{v}^3+\hat I_2s_{v}+\hat I_3\\ &=-s_{v}^3+(s_{xz}^2+s_{xy}^2+s_{yz}^2-s_{xx}s_{yy}-s_{xx}s_{zz}-s_{yy}s_{zz})s_{v}\\ &\quad+s_{xx}s_{yy}s_{zz}+2s_{xz}s_{xy}s_{yz}-s_{yy}s_{xz}^2-s_{zz}s_{xy}^2-s_{xx}s_{yz}^2\\ &=-s_{v}^3+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)s_{v} \\ &\quad- [(\sigma_{xx}-\sigma_{m})(\sigma_{yy}-\sigma_{m})+(\sigma_{xx}-\sigma_{m})(\sigma_{zz}-\sigma_{m})+(\sigma_{yy}-\sigma_{m})( \sigma_{zz}-\sigma_{m})]s_{v}\\ &\quad+(\sigma_{xx}-\sigma_{m})(\sigma_{yy}-\sigma_{m})(\sigma_{zz}-\sigma_{m})+2\tau_{xz}\tau_{xy}\tau_{yz}\\ &\quad-(\sigma_{yy}-\sigma_{m})\tau_{xz}^2-(\sigma_{zz}-\sigma_{m})\tau_{xy}^2-(\sigma_{xx}-\sigma_{m})\tau_{yz}^2\\ &=-s_{v}^3+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)s_{v}-(\sigma_{xx}\sigma_{yy}+\sigma_{xx}\sigma_{zz}+\sigma_{yy}\sigma_{zz})s_{v}\\ &\quad+(\sigma_{xx}+\sigma_{yy}+\sigma_{xx}+\sigma_{zz}+\sigma_{yy}+\sigma_{zz})\sigma_{m}s_{v}-3\sigma_{m}^2s_{v}+\sigma_{xx}\sigma_{yy}\sigma_{zz}\\ &\quad-(\sigma_{xx}\sigma_{yy}+\sigma_{xx}\sigma_{zz}+\sigma_{yy}\sigma_{zz})\sigma_{m}+(\sigma_{xx}+\sigma_{yy}+\sigma_{zz})\sigma_{m}^2-\sigma_{m}^3+2\tau_{xz}\tau_{xy}\tau_{yz}\\ &\quad-(\sigma_{yy}\tau_{xz}^2+\sigma_{zz}\tau_{xy}^2+\sigma_{xx}\tau_{yz}^2)+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)\sigma_{m}\\ &=-s_{v}^3+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)(s_{v}+\sigma_{m})-(\sigma_{xx}\sigma_{yy}+\sigma_{xx}\sigma_{zz}+\sigma_{yy}\sigma_{zz})(s_{v}+\sigma_{m})\\ &\quad+6\sigma_{m}^2s_{v}-3\sigma_{m}^2s_{v}+\sigma_{xx}\sigma_{yy}\sigma_{zz}+3\sigma_{m}^3-\sigma_{m}^3+2\tau_{xz}\tau_{xy}\tau_{yz}\\ &\quad-(\sigma_{yy}\tau_{xz}^2+\sigma_{zz}\tau_{xy}^2+\sigma_{xx}\tau_{yz}^2)\\ &=-s_{v}^3+3\sigma_{m}^2s_{v}+2\sigma_{m}^3+(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2-\sigma_{xx}\sigma_{yy}-\sigma_{xx}\sigma_{zz}-\sigma_{yy}\sigma_{zz})(s_{v}+\sigma_{m})\\ &\quad+\sigma_{xx}\sigma_{yy}\sigma_{zz}+2\tau_{xz}\tau_{xy}\tau_{yz}-(\sigma_{yy}\tau_{xz}^2+\sigma_{zz}\tau_{xy}^2+\sigma_{xx}\tau_{yz}^2)\\ &=-s_{v}^3-\sigma_{m}^3+3\sigma_{m}^2s_{v}+3\sigma_{m}^3+I_2(s_{v}+\sigma_{m})+I_3\\ &=-(s_{v}+\sigma_{m})(s_{v}^2-s_{v}\sigma_{m}+\sigma_{m}^2-3\sigma_{m}^2)+I_2(s_{v}+\sigma_{m})+I_3\\ &=-(s_{v}+\sigma_{m})(s_{v}+\sigma_{m})(s_{v}-2\sigma_{m})+I_2(s_{v}+\sigma_{m})+I_3\\ &=-(s_{v}+\sigma_{m})^2(s_{v}+\sigma_{m}-3\sigma_{m})+I_2(s_{v}+\sigma_{m})+I_3\\ &=-(s_{v}+\sigma_{m})^3+3\sigma_{m}(s_{v}+\sigma_{m})^2+I_2(s_{v}+\sigma_{m})+I_3\\ &=-(s_{v}+\sigma_{m})^3+I_1(s_{v}+\sigma_{m})^2+I_2(s_{v}+\sigma_{m})+I_3=0 \end{aligned}\tag{17} sv3+I^1sv2+I^2sv+I^3=sv3+I^2sv+I^3=sv3+(sxz2+sxy2+syz2sxxsyysxxszzsyyszz)sv+sxxsyyszz+2sxzsxysyzsyysxz2szzsxy2sxxsyz2=sv3+(τxz2+τxy2+τyz2)sv[(σxxσm)(σyyσm)+(σxxσm)(σzzσm)+(σyyσm)(σzzσm)]sv+(σxxσm)(σyyσm)(σzzσm)+2τxzτxyτyz(σyyσm)τxz2(σzzσm)τxy2(σxxσm)τyz2=sv3+(τxz2+τxy2+τyz2)sv(σxxσyy+σxxσzz+σyyσzz)sv+(σxx+σyy+σxx+σzz+σyy+σzz)σmsv3σm2sv+σxxσyyσzz(σxxσyy+σxxσzz+σyyσzz)σm+(σxx+σyy+σzz)σm2σm3+2τxzτxyτyz(σyyτxz2+σzzτxy2+σxxτyz2)+(τxz2+τxy2+τyz2)σm=sv3+(τxz2+τxy2+τyz2)(sv+σm)(σxxσyy+σxxσzz+σyyσzz)(sv+σm)+6σm2sv3σm2sv+σxxσyyσzz+3σm3σm3+2τxzτxyτyz(σyyτxz2+σzzτxy2+σxxτyz2)=sv3+3σm2sv+2σm3+(τxz2+τxy2+τyz2σxxσyyσxxσzzσyyσzz)(sv+σm)+σxxσyyσzz+2τxzτxyτyz(σyyτxz2+σzzτxy2+σxxτyz2)=sv3σm3+3σm2sv+3σm3+I2(sv+σm)+I3=(sv+σm)(sv2svσm+σm23σm2)+I2(sv+σm)+I3=(sv+σm)(sv+σm)(sv2σm)+I2(sv+σm)+I3=(sv+σm)2(sv+σm3σm)+I2(sv+σm)+I3=(sv+σm)3+3σm(sv+σm)2+I2(sv+σm)+I3=(sv+σm)3+I1(sv+σm)2+I2(sv+σm)+I3=0(17)
对比公式(17)和公式(10),如下所示。
− ( s v + σ m ) 3 + I 1 ( s v + σ m ) 2 + I 2 ( s v + σ m ) + I 3 = 0 − σ v 3 + I 1 σ v 2 + I 2 σ v + I 3 = 0 -(s_{v}+\sigma_{m})^3+I_1(s_{v}+\sigma_{m})^2+I_2(s_{v}+\sigma_{m})+I_3=0\\ -\sigma_{v}^3+I_1\sigma_{v}^2+I_2\sigma_{v}+I_3=0 (sv+σm)3+I1(sv+σm)2+I2(sv+σm)+I3=0σv3+I1σv2+I2σv+I3=0
不难发现, s v + σ m = σ v s_{v}+\sigma_{m}=\sigma_{v} sv+σm=σv,即应力偏张量主值和应力张量主值存在以上转换关系。
同时,从公式(16)的 I ^ 1 \hat I_1 I^1
I ^ 1 = s x x + s y y + s z z = 0 \hat I_1=s_{xx}+s_{yy}+s_{zz}=0 I^1=sxx+syy+szz=0
那么可以得到
( s x x + s y y + s z z ) 2 = s x x 2 + s y y 2 + s z z 2 + 2 ( s x x s y y + s x x s z z + s y y s z z ) = 0 (18) (s_{xx}+s_{yy}+s_{zz})^2=s_{xx}^2+s_{yy}^2+s_{zz}^2+2(s_{xx}s_{yy}+s_{xx}s_{zz}+s_{yy}s_{zz})=0\tag{18} (sxx+syy+szz)2=sxx2+syy2+szz2+2(sxxsyy+sxxszz+syyszz)=0(18)
于是
− 6 ( s x x s y y + s x x s z z + s y y s z z ) = 2 s x x 2 + 2 s y y 2 + 2 s z z 2 − 2 ( s x x s y y + s x x s z z + s y y s z z ) = ( s x x − s y y ) 2 + ( s x x − s z z ) 2 + ( s y y − s z z ) 2 = ( σ x x − σ y y ) 2 + ( σ x x − σ z z ) 2 + ( σ y y − σ z z ) 2 (19) \begin{aligned} -6(s_{xx}s_{yy}+s_{xx}s_{zz}+s_{yy}s_{zz})&=2s_{xx}^2+2s_{yy}^2+2s_{zz}^2-2(s_{xx}s_{yy}+s_{xx}s_{zz}+s_{yy}s_{zz})\\ &=(s_{xx}-s_{yy})^2+(s_{xx}-s_{zz})^2+(s_{yy}-s_{zz})^2\\ &=(\sigma_{xx}-\sigma_{yy})^2+(\sigma_{xx}-\sigma_{zz})^2+(\sigma_{yy}-\sigma_{zz})^2 \end{aligned}\tag{19} 6(sxxsyy+sxxszz+syyszz)=2sxx2+2syy2+2szz22(sxxsyy+sxxszz+syyszz)=(sxxsyy)2+(sxxszz)2+(syyszz)2=(σxxσyy)2+(σxxσzz)2+(σyyσzz)2(19)
把公式(16)的 I ^ 2 \hat I_2 I^2,那么
− 6 ( s x x s y y + s x x s z z + s y y s z z ) = 6 [ I ^ 2 − ( τ x z 2 + τ x y 2 + τ y z 2 ) ] = ( σ x x − σ y y ) 2 + ( σ x x − σ z z ) 2 + ( σ y y − σ z z ) 2 (20) \begin{aligned} -6(s_{xx}s_{yy}+s_{xx}s_{zz}+s_{yy}s_{zz})&=6[\hat I_2-(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)]\\ &=(\sigma_{xx}-\sigma_{yy})^2+(\sigma_{xx}-\sigma_{zz})^2+(\sigma_{yy}-\sigma_{zz})^2 \end{aligned}\tag{20} 6(sxxsyy+sxxszz+syyszz)=6[I^2(τxz2+τxy2+τyz2)]=(σxxσyy)2+(σxxσzz)2+(σyyσzz)2(20)

I ^ 2 = 1 6 [ ( σ x x − σ y y ) 2 + ( σ x x − σ z z ) 2 + ( σ y y − σ z z ) 2 + 6 ( τ x z 2 + τ x y 2 + τ y z 2 ) ] (21) \hat I_2=\frac{1}{6}[(\sigma_{xx}-\sigma_{yy})^2+(\sigma_{xx}-\sigma_{zz})^2+(\sigma_{yy}-\sigma_{zz})^2+6(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)]\tag{21} I^2=61[(σxxσyy)2+(σxxσzz)2+(σyyσzz)2+6(τxz2+τxy2+τyz2)](21)

应力偏张量的第二不变量更常用的还有一个形式如下。
I ^ 2 = 1 6 [ ( σ 1 − σ 2 ) 2 + ( σ 1 − σ 3 ) 2 + ( σ 2 − σ 3 ) 2 ] (22) \hat I_2=\frac{1}{6}[(\sigma_{1}-\sigma_{2})^2+(\sigma_{1}-\sigma_{3})^2+(\sigma_{2}-\sigma_{3})^2]\tag{22} I^2=61[(σ1σ2)2+(σ1σ3)2+(σ2σ3)2](22)
看公式(21)(22)是不是特别的眼熟,其实他就是等效应力的来源,定义应力强度 σ e \sigma_e σe
σ e = 3 I 2 = 1 2 [ ( σ 1 − σ 2 ) 2 + ( σ 1 − σ 3 ) 2 + ( σ 2 − σ 3 ) 2 ] = 1 2 [ ( σ x x − σ y y ) 2 + ( σ x x − σ z z ) 2 + ( σ y y − σ z z ) 2 + 6 ( τ x z 2 + τ x y 2 + τ y z 2 ) ] (22) \sigma_e=\sqrt{3I_2}=\sqrt{\frac{1}{2}[(\sigma_{1}-\sigma_{2})^2+(\sigma_{1}-\sigma_{3})^2+(\sigma_{2}-\sigma_{3})^2]}\\ =\sqrt{\frac{1}{2}[(\sigma_{xx}-\sigma_{yy})^2+(\sigma_{xx}-\sigma_{zz})^2+(\sigma_{yy}-\sigma_{zz})^2+6(\tau_{xz}^2+\tau_{xy}^2+\tau_{yz}^2)]} \tag{22} σe=3I2 =21[(σ1σ2)2+(σ1σ3)2+(σ2σ3)2] =21[(σxxσyy)2+(σxxσzz)2+(σyyσzz)2+6(τxz2+τxy2+τyz2)] (22)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/234821.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

ros2 学习04 工作空间说明及示例

ros2 术语说明: 大家在学习其他的开发语言之前的学习和开发中,应该有接触过某些集成开发环境,比如Visual Studio、Eclipse、Qt Creator等,当我们想要编写程序之前,都会在这些开发环境的工具栏中,点击一个“…

c语言力扣题目:消失的数字(有关时间复杂度O(N²)O(N))以及对异或操作符的更深入的理解(如何用人脑的十进制去考量二进制)

目录 Way One :暴力求解,时间复杂度为 O(N) 代码1 Way Two : 时间复杂度限制到 O(N) 代码及其详解 如题 Way One :暴力求解,时间复杂度为 O(N) 大体思路:比如这里我们需要处理的整型数组是"3,0,1",我们可以用冒泡排序或者 qsort函数将他从大到小进行排序成"…

SpringBlade export-user SQL 注入漏洞复现

0x01 产品简介 SpringBlade 是一个由商业级项目升级优化而来的 SpringCloud 分布式微服务架构、SpringBoot 单体式微服务架构并存的综合型项目。 0x02 漏洞概述 SpringBlade v3.2.0 及之前版本框架后台 export-user 路径存在安全漏洞,攻击者利用该漏洞可通过组件customSqlS…

手把手教你在 windows 上安装 Docker

前言 大家好,我是潇潇雨声,今天为大家带来一篇关于在 Windows 环境下使用 Docker 的教程。对于 Docker,我们可以简单地将其看作一种方便的软件安装方式,而无需深入涉及其复杂的概念。选择使用 Docker 主要是为了省事,比…

动态规划系列 | 一文搞定区间DP

文章目录 特点石子合并题目描述问题分析程序代码复杂度分析 环形石子合并题目描述问题分析程序代码复杂度分析 能量项链题目描述问题分析程序代码复杂度分析 加分二叉树题目描述问题分析程序代码复杂度分析 凸多边形的划分题目描述问题分析程序代码复杂度分析 棋盘分割题目描述…

STM32内部是怎么工作的

STM32是怎么工作的 1 从孩子他妈说起2 早期计算机的组成2.1 五大元件(1)第一个出场的是电容元件(2)第二个出场的是二极管(3)第三个出场的是电阻元件(4)第四个出场的是电感&#xff0…

现代雷达车载应用——第3章 MIMO雷达技术 3.1节 基于MIMO雷达的虚拟阵列合成

经典著作,值得一读,英文原版下载链接【免费】ModernRadarforAutomotiveApplications资源-CSDN文库。 3 MIMO雷达技术 自20世纪90年代末以来,带有少量天线的汽车雷达已被用于高级驾驶员辅助系统(ADAS)的目的。这些早期的汽车雷达主要提供目标…

Flask学习四:补充

插件 flask-caching 简介 Flask-Caching 是一个 Flask 扩展,旨在为 Flask 应用程序添加缓存功能。缓存是一种提高应用性能的技术,通过将常用数据暂时存储在一个快速访问的位置(如内存或磁盘),从而减少对较慢资源&…

HarmonyOS应用开发实战—开箱即用的活动创建页面【ArkTS】【鸿蒙专栏-33】

一.HarmonyOS应用开发实战—开箱即用的个人主页页面【ArkTS】【鸿蒙专栏-32】 1.1 项目背景 HarmonyOS(鸿蒙操作系统)是华为公司推出的一种分布式操作系统。它被设计为一种全场景、全连接的操作系统,旨在实现在各种设备之间的无缝协同和共享,包括智能手机、平板电脑、智能…

Android排队预约系统(Java+SqLite+ZXing)

自己写的排队预约系统,可改写,添加功能,如管理用户,查询排队人数等功能。(由于是选修课课设,所以写的比较粗糙) 使用方法: 1.使用Android studio导入项目。 2.使用gradle加载build.gradle.kts中的依赖。…

【MySQL工具】pt-online-schema-change 详细介绍

下面文章基于 pt-online-schema-change 3.5.5 版本 简介 pt-online-schema-change - 无锁表修改表结构工具,这里无锁表也不是绝对的,在交互原表与中间表 表名的时候也会有元数据锁,只不过事件很短 语法 pt-online-schema-change [OPTIONS…

计算BMI指数-第11届蓝桥杯选拔赛Python真题精选

[导读]:超平老师的Scratch蓝桥杯真题解读系列在推出之后,受到了广大老师和家长的好评,非常感谢各位的认可和厚爱。作为回馈,超平老师计划推出《Python蓝桥杯真题解析100讲》,这是解读系列的第19讲。 计算BMI指数&…

智能优化算法应用:基于混沌博弈算法3D无线传感器网络(WSN)覆盖优化 - 附代码

智能优化算法应用:基于混沌博弈算法3D无线传感器网络(WSN)覆盖优化 - 附代码 文章目录 智能优化算法应用:基于混沌博弈算法3D无线传感器网络(WSN)覆盖优化 - 附代码1.无线传感网络节点模型2.覆盖数学模型及分析3.混沌博弈算法4.实验参数设定5.算法结果6.…

SQLturning:定位连续值范围起点和终点

在上一篇blog说到,如何去优化查询连续值范围,没看过的朋友,上篇blog链接[在此]。(https://blog.csdn.net/weixin_42575078/article/details/135067645?spm1001.2014.3001.5501) 那么今天来说说怎么将连续的数据合并,然后返回合并…

SpringSecurity入门

前言 Spring Security是一个用于在Java应用程序中提供身份验证和授权功能的强大框架。它构建在Spring框架之上,为开发人员提供了一套灵活且全面的安全性服务,本篇将为大家带来Spring Security的详细介绍及入门 一.安全框架 在学习了解Spring Security之…

MeterSphere files 任意文件读取漏洞复现 (CVE-2023-25573)

0x01 产品简介 MeterSphere 是一站式开源持续测试平台, 涵盖测试跟踪、接口测试、UI 测试和性能测试等功能,全面兼容 JMeter、Selenium 等主流开源标准。 0x02 漏洞概述 MeterSphere /api/jmeter/download/files 路径文件存在文件读取漏洞,攻击者可通过该漏洞读取系统重要…

[总线仲裁]

目录 一. 集中仲裁方式1.1 链式查询方式1.2 计数器查询方式1.3 独立请求方式 二. 分布式仲裁方式 总线仲裁是为了解决多个设备争用总线这个问题 \quad 一. 集中仲裁方式 \quad 集中仲裁方式: 就像是霸道总裁来决定谁先获得总线控制权 分布仲裁方式: 商量着谁先获得总线控制权 …

国产670亿参数的DeepSeek:超越Llama2,全面开源

模型概述 DeepSeek,一款国产大型语言模型(LLM),凭借其670亿参数的规模,正引领着人工智能领域的新浪潮。这款模型不仅在多项中英文公开评测榜单上超越了700亿参数的Llama 2,而且在推理、数学和编程能力方面…

JDK各个版本特性讲解-JDK16特性

JDK各个版本特性讲解-JDK16特性 一、JAVA16概述二、语法层面变化1.JEP 397:密封类(第二次预览)2.JEP 394:instanceof 的模式匹配3.JEP 395:记录4_JEP 390:基于值的类的警告 三、API层面变化1.JEP 338&#…

(备战2024)三天吃透Java面试八股文,面试通过率高达90%

什么样的求职者能够获得面试官的青睐?求职者需要准备哪些内容来面对形形色色的面试官?这两份资料是我在几十场面试中被面试官问到的问题,比其他复制粘贴的面试题强一百倍,堪称全网最强(我不太喜欢“全网最强”这样的字…