大家好,欢迎来到IT知识分享网。
1.热应力概述
物体的变形不仅仅由外力作用引起,温度的变化也能够引起变形,称为热变形:如热胀冷缩。
但是单有温度的变化,不一定就在物体内部产生应力,只有当温度变化所引起的膨胀或收缩受到约束时,才会在物体内产生应力。这种无外力作用而是由于温度变化引起的热变形受到约束而产生的应力,称为热应力或温度应力。
- 外部变形约束:物体受到外部条件影响,不能自由伸缩,如两端固定等;
- 内部变形约束:物体内部温度分布不均匀,如物体内外、上下存在温差,导致零件翘曲;
- 相互变形约束:不同零件组合在一起,热膨胀系数不同,相互影响;
2.广义胡克定律
对于各向同性物体,在自由膨胀情况下,假设温度变化为 t,热膨胀率为 α \alpha α,则:
ε x = ε y = ε z = α t γ x y = γ y z = γ z x = 0 (1) \begin{aligned} &\varepsilon_x=\varepsilon_y=\varepsilon_z=\alpha t\\ &\gamma_{xy}=\gamma_{yz}=\gamma_{zx}=0 \end{aligned}\tag1 εx=εy=εz=αtγxy=γyz=γzx=0(1)
已知体积应力 Θ = σ x + σ y + σ z (3) \Theta =\sigma_x+\sigma_y+\sigma_z\tag3 Θ=σx+σy+σz(3),
剪切弹性模量 G = E 2 ( 1 + μ ) (4) G=\frac{E}{2(1+\mu)}\tag4 G=2(1+μ)E(4)
式(2)中正应变可整理为:
{ ε x = 1 2 G ( σ x − μ 1 + μ Θ ) + α t ε y = 1 2 G ( σ y − μ 1 + μ Θ ) + α t ε z = 1 2 G ( σ z − μ 1 + μ Θ ) + α t (5) \begin{cases} \varepsilon_x=\frac{1}{2G}(\sigma_x-\frac\mu{1+\mu}\Theta)+\alpha t\\[1.5ex] \varepsilon_y=\frac{1}{2G}(\sigma_y-\frac\mu{1+\mu}\Theta)+\alpha t\\[1.5ex] \varepsilon_z=\frac{1}{2G}(\sigma_z-\frac\mu{1+\mu}\Theta)+\alpha t \end{cases}\tag5 ⎩⎪⎪⎪⎨⎪⎪⎪⎧εx=2G1(σx−1+μμΘ)+αtεy=2G1(σy−1+μμΘ)+αtεz=2G1(σz−1+μμΘ)+αt(5)
由式(5),对于体积应变 e e e,
e = ε x + ε y + ε z = 1 − 2 μ 1 + μ Θ 2 G + 3 α t (7) e=\varepsilon_x+\varepsilon_y+\varepsilon_z=\frac{1-2\mu}{1+\mu}\frac{\Theta}{2G}+3\alpha t\tag7 e=εx+εy+εz=1+μ1−2μ2GΘ+3αt(7)
根据式(6)、式 (7),将正应力公式整理为以下形式:
{ σ x = 2 G ε x + λ e − β t σ y = 2 G ε y + λ e − β t σ z = 2 G ε y + λ e − β t (9) \begin{cases}\sigma_x=2G\varepsilon_x+\lambda e -\beta t\\[1.5ex] \sigma_y=2G\varepsilon_y+\lambda e -\beta t\\[1.5ex] \sigma_z=2G\varepsilon_y+\lambda e -\beta t \end{cases}\tag9 ⎩⎪⎪⎪⎨⎪⎪⎪⎧σx=2Gεx+λe−βtσy=2Gεy+λe−βtσz=2Gεy+λe−βt(9)
其中, λ \lambda λ 为拉梅常数:
λ = E μ ( 1 + μ ) ( 1 − 2 μ ) (10) \lambda=\frac{E\mu}{(1+\mu)(1-2\mu)}\tag{10} λ=(1+μ)(1−2μ)Eμ(10)
β \beta β 为热应力系数:
β = α ( 3 λ + 2 G ) (11) \beta=\alpha(3\lambda+2G)\tag {11} β=α(3λ+2G)(11)
3.位移平衡方程
我们假设热应力为 p p p,则:
∇ p = [ ∂ σ x ∂ x + ∂ τ y x ∂ y + ∂ τ z x ∂ z ∂ τ x y ∂ x + σ τ y ∂ y + ∂ τ z y ∂ z ∂ τ x z ∂ x + ∂ τ y z ∂ y + ∂ σ z ∂ z ] (12) \nabla p = \left[\begin{matrix} \frac{\partial\sigma_{x}}{\partial x}+\frac{\partial\tau_{yx}}{\partial y}+\frac{\partial\tau_{zx}}{\partial z} \\[1.5ex] \frac{\partial\tau_{xy}}{\partial x}+\frac{\sigma\tau_{y}}{\partial y}+\frac{\partial\tau_{zy}}{\partial z} \\[1.5ex] \frac{\partial\tau_{xz}}{\partial x}+\frac{\partial\tau_{yz}}{\partial y}+\frac{\partial\sigma_{z}}{\partial z} \end{matrix}\right] \tag{12} ∇p=⎣⎢⎢⎢⎡∂x∂σx+∂y∂τyx+∂z∂τzx∂x∂τxy+∂yστy+∂z∂τzy∂x∂τxz+∂y∂τyz+∂z∂σz⎦⎥⎥⎥⎤(12)
仅受应力影响的应变微分方程位移形式表述如下:
ε = [ ε x ε y ε z γ x y γ y z γ z x ] = [ ∂ u ∂ x ∂ v ∂ y ∂ w ∂ z ∂ v ∂ x + ∂ u ∂ y ∂ w ∂ y + ∂ v ∂ z ∂ u ∂ z + ∂ w ∂ x ] (13) \varepsilon=\begin{bmatrix}\varepsilon_x \\\varepsilon_y \\\varepsilon_z\\ \gamma_{xy}\\\gamma_{yz}\\\gamma_{zx} \end{bmatrix}=\begin{bmatrix} \frac{\partial u}{\partial x} \\\frac{\partial v}{\partial y}\\ \frac{\partial w}{\partial z} \\ \frac{\partial v}{\partial x}+\frac{\partial u}{\partial y} \\\frac{\partial w}{\partial y}+\frac{\partial v}{\partial z} \\\frac{\partial u}{\partial z}+\frac{\partial w}{\partial x} \end{bmatrix} \tag{13} ε=⎣⎢⎢⎢⎢⎢⎢⎡εxεyεzγxyγyzγzx⎦⎥⎥⎥⎥⎥⎥⎤=⎣⎢⎢⎢⎢⎢⎢⎢⎡∂x∂u∂y∂v∂z∂w∂x∂v+∂y∂u∂y∂w+∂z∂v∂z∂u+∂x∂w⎦⎥⎥⎥⎥⎥⎥⎥⎤(13)
结合式(6)、(9)、(12)、(13)可得:
∇ p = [ ( λ + G ) ∂ e ∂ x + G ∇ 2 u − β ∂ t ∂ x ( λ + G ) ∂ e ∂ y + G ∇ 2 v − β ∂ t ∂ y ( λ + G ) ∂ e ∂ z + G ∇ 2 w − β ∂ t ∂ z ] (14) \nabla p = \left[\begin{matrix} (\lambda+G)\frac{\partial e}{\partial x}+G\nabla^2u-\beta\frac{\partial t}{\partial x} \\[1.5ex] (\lambda+G)\frac{\partial e}{\partial y}+G\nabla^2v-\beta\frac{\partial t}{\partial y} \\[1.5ex] (\lambda+G)\frac{\partial e}{\partial z}+G\nabla^2w-\beta\frac{\partial t}{\partial z} \end{matrix}\right] \tag{14} ∇p=⎣⎢⎢⎡(λ+G)∂x∂e+G∇2u−β∂x∂t(λ+G)∂y∂e+G∇2v−β∂y∂t(λ+G)∂z∂e+G∇2w−β∂z∂t⎦⎥⎥⎤(14)
除了必须满足微分平衡方程外,也需要满足边界条件。当忽略表面力(即表面力为零)时:
{ σ x l + τ y x m + τ z x n = 0 σ y m + τ y x n + τ z x l = 0 σ z n + τ y x l + τ z x m = 0 (15) \begin{cases} \sigma_{x}l+\tau_{yx}m+\tau_{zx}n=0\\[1.5ex] \sigma_{y}m+\tau_{yx}n+\tau_{zx}l =0\\[1.5ex] \sigma_{z}n+\tau_{yx}l+\tau_{zx}m=0 \end{cases} \tag{15} ⎩⎪⎪⎪⎨⎪⎪⎪⎧σxl+τyxm+τzxn=0σym+τyxn+τzxl=0σzn+τyxl+τzxm=0(15)
其中,l、m、n 分别为边界表面在x、y、z方向的法向余弦。
免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://haidsoft.com/135628.html