文章目录
- 2. 简单应力状态下的应力应变关系
- 2.1 简单拉伸的应力应变关系
- 2.2 真实应力应变关系
- 2.3 应力-应变关系简化模型
2. 简单应力状态下的应力应变关系
我们在高中就学过,弹簧拉伸力和变形量成比例,对于一般的金属材料,在一定载荷以内这种结论也是成立的,这种情况称之为弹性。在下面的材料力学单向拉伸试验的结果中,我们可以看到材料先发生弹性变形,超过一定极限后产生塑性变形。
2.1 简单拉伸的应力应变关系
材料拉伸试样如下图1所示,在试验时在试样两端夹持,施加载荷,试验在载荷作用下会伸长,记录相应时刻的载荷-位移数据,并以此绘制曲线如下图2所示。
图
1
拉伸试样
图1 \quad拉伸试样
图1拉伸试样
图
2
载荷
−
变形曲线
图2 \quad载荷-变形曲线
图2载荷−变形曲线
但一般会进行以下变换
σ
=
F
A
0
,
ε
=
Δ
l
0
(1)
\sigma=\frac{F}{A_0}, \varepsilon=\frac{\Delta}{l_0}\tag{1}
σ=A0F,ε=l0Δ(1)
那么可以将载荷-变形曲线变换成
σ
−
ε
\sigma-\varepsilon
σ−ε应力应变曲线,如下图3。
图
3
σ
−
ε
曲线
图3\quad \sigma-\varepsilon曲线
图3σ−ε曲线
其中A点称为上屈服点,B为下屈服点,材料在两者之间呈现流动状态,应力不发生显著变化只增加应变,一般上屈服点和下屈服点区别不大(此屈服流动的现象一般在低碳钢中存在,合金钢等往往没有明显的特征),用
σ
s
\sigma_s
σs来代替,称材料的屈服强度。材料在
σ
s
\sigma_s
σs以下呈现弹性变形,也就是载荷去除之后变形能够完全恢复,且应力和应变成比例。同时,金属材料压缩的应力应变曲线基本上与拉伸的应力应变曲线接近。
有一些金属可能没有上下屈服流动阶段,如下图4。
在应力超过弹性极限后,会产生塑性应变 ε p \varepsilon^p εp,缓慢卸去载荷,变形也不能完全恢复(见图4中 ε p \varepsilon^p εp),这种现象成为屈服。同时卸载曲线也是线性的,并且斜率和刚开始的弹性段一样,直到反向屈服。
图4中A’为材料的压缩屈服点,假设A与A’对称(即压缩屈服等于拉伸屈服),左图为多晶材料,多晶材料反向屈服点M’一般绝对值小于A’,称为包晶格效应(Bauschingerx effect),即在拉伸方向的强化导致压缩方向的弱化,这种效应在后文中还会应用(就是随动强化模型)。
右图为单晶材料,材料反向屈服点M’一般绝对值大于A’,即在拉伸方向的强化导致压缩方向的同样的强化,这种效应在后文中还会应用(就是等向强化模型)。
图
4
加载
/
卸载应力应变曲线
图4\quad 加载/卸载应力应变曲线
图4加载/卸载应力应变曲线
2.2 真实应力应变关系
在图3中,在应力达到最高点C以前,应力和应变一同增加,到达C点后应变增加,应力却下降了。事实上,在C点前由于泊松比的存在,变形前试样的截面随着载荷增加,会慢慢减小,但是在随着变形的继续,某一时刻横截面会较迅速的减小,这种现象称为颈缩(也称塑性失稳),由于截面的迅速缩小,试样的承载能力随之下降,相应的名义应力也下降。因此实际上,名义应力在变形量小的时候跟试样的真实应力差别不大,但是在颈缩时,名义应力和真实应力较大差别。
按照定义,定义真实应力
σ
~
\widetilde{\sigma}
σ
如下
σ
~
=
P
A
(2)
\widetilde{\sigma}=\frac{P}{A}\tag{2}
σ
=AP(2)
其中
A
A
A为试样瞬时截面,
P
P
P为试样瞬时载荷。
瞬时的应变增量
d
ε
~
d\widetilde{\varepsilon}
dε
如下所示
d
ε
~
=
d
l
′
l
′
(3)
d\widetilde{\varepsilon}=\frac{dl'}{l'}\tag{3}
dε
=l′dl′(3)
其中
l
′
l'
l′为试样瞬时长度,
d
l
′
dl'
dl′为试样瞬时伸长量。
那么真实应变
ε
~
\widetilde{\varepsilon}
ε
为
ε
~
=
∫
l
0
l
d
l
′
l
′
=
ln
(
l
′
)
∣
l
0
l
=
ln
(
l
l
0
)
=
ln
(
l
−
l
0
+
l
0
l
0
)
=
ln
(
1
+
ε
)
(4)
\widetilde{\varepsilon}=\int_{l_0}^{l}\frac{dl'}{l'}=\ln(l')|_{l_0}^{l}=\ln(\frac{l}{l_0})=\ln(\frac{l-l_0+l_0}{l_0})=\ln(1+\varepsilon)\tag{4}
ε
=∫l0ll′dl′=ln(l′)∣l0l=ln(l0l)=ln(l0l−l0+l0)=ln(1+ε)(4)
在材料进入塑性阶段,材料表现出塑性流动的特征,这里需要引用材料几乎不可压缩的假设,因此有
A 0 l 0 = A l (5) A_0l_0=Al\tag{5} A0l0=Al(5)
那么真实应力变
σ
~
\widetilde{\sigma}
σ
为
σ
~
=
P
A
=
P
A
0
⋅
A
0
A
=
P
A
0
⋅
l
l
0
=
σ
e
ε
~
=
σ
(
1
+
ε
)
(6)
\widetilde{\sigma}=\frac{P}{A}=\frac{P}{A_0}\cdot\frac{A_0}{A}=\frac{P}{A_0}\cdot\frac{l}{l_0}=\sigma e^{\widetilde{\varepsilon}}=\sigma (1+\varepsilon)\tag{6}
σ
=AP=A0P⋅AA0=A0P⋅l0l=σeε
=σ(1+ε)(6)
在图5右侧曲线中C为名义应力达到最大值,在此时,有
d
σ
d
ε
=
0
(7)
\frac{d\sigma}{d\varepsilon}=0\tag{7}
dεdσ=0(7)
在图5右侧曲线中C为名义应力达到最大值对应的真实应力点C’’,其中真实应力有式(6),那么有
d
σ
~
d
ε
=
σ
(8)
\frac{d\widetilde{\sigma}}{d\varepsilon}=\sigma\tag{8}
dεdσ
=σ(8)
在图5左侧曲线中C为名义应力达到最大值对应的真实应力点C’,真实应力和真实应变应满足的条件如下
d
σ
~
d
ε
~
=
(
d
σ
d
ε
⋅
d
ε
d
ε
~
)
e
ε
~
+
σ
e
ε
~
=
σ
e
ε
~
=
σ
~
(9)
\frac{d\widetilde{\sigma}}{d\widetilde{\varepsilon}}=(\frac{d\sigma}{d\varepsilon}\cdot\frac{d\varepsilon}{d\widetilde{\varepsilon}})e^{\widetilde{\varepsilon}}+\sigma e^{\widetilde{\varepsilon}}=\sigma e^{\widetilde{\varepsilon}}=\widetilde{\sigma}\tag{9}
dε
dσ
=(dεdσ⋅dε
dε)eε
+σeε
=σeε
=σ
(9)
图
5
应力应变曲线
图5\quad 应力应变曲线
图5应力应变曲线
2.3 应力-应变关系简化模型
在理论分析中,常常采用简化的应力应变模型来分析具体问题(实际工程中也多有应用)。
简化模型一:理想弹塑性模型
如下图6所示,那么应力应变关系可以写为
σ
=
{
E
ε
,
ε
≤
ε
s
σ
s
s
i
g
n
ε
,
ε
>
ε
s
(10)
\sigma=\begin{cases}E\varepsilon &\quad,\quad \varepsilon\le\varepsilon_s\\ \sigma_s sign \varepsilon &\quad,\quad \varepsilon\gt\varepsilon_s \end{cases}\tag{10}
σ={Eεσssignε,ε≤εs,ε>εs(10)
图
6
理想弹塑性模型
图6\quad 理想弹塑性模型
图6理想弹塑性模型
简化模型二:线性强化弹塑性模型
如下图7所示,那么应力应变关系可以写为
σ = { E ε , ε ≤ ε s σ s + E ’ ( ε − ε s ) , ε > ε s (11) \sigma=\begin{cases}E\varepsilon &\quad,\quad \varepsilon\le\varepsilon_s\\ \sigma_s +E’( \varepsilon-\varepsilon_s) &\quad,\quad \varepsilon\gt\varepsilon_s \end{cases}\tag{11} σ={Eεσs+E’(ε−εs),ε≤εs,ε>εs(11)
图
7
线性强化弹塑性模型
图7\quad 线性强化弹塑性模型
图7线性强化弹塑性模型
当然,上式当进入塑性后,还可以写成另外一种形式,如下所示
ε
=
ε
e
+
ε
p
(12)
\varepsilon=\varepsilon^e+\varepsilon^p\tag{12}
ε=εe+εp(12)
σ
=
σ
s
+
h
ε
p
(13)
\sigma=\sigma_s+h\varepsilon^p\tag{13}
σ=σs+hεp(13)
可以通过下图来确定h,下图为应力-塑性应变图,由图7可知
E
′
=
d
σ
d
ε
(14)
E'=\frac{d\sigma}{d\varepsilon}\tag{14}
E′=dεdσ(14)
同时由图8可知
h
=
d
σ
d
ε
p
=
d
σ
d
ε
−
d
ε
e
=
1
d
ε
d
σ
−
d
ε
e
d
σ
=
1
1
E
′
−
1
E
(15)
h=\frac{d\sigma}{d\varepsilon^p}=\frac{d\sigma}{d\varepsilon-d\varepsilon^e}=\frac{1}{\frac{d\varepsilon}{d\sigma}-\frac{d\varepsilon^e}{d\sigma}}=\frac{1}{\frac{1}{E'}-\frac{1}{E}}\tag{15}
h=dεpdσ=dε−dεedσ=dσdε−dσdεe1=E′1−E11(15)
代入(13),那么有
ε
p
=
σ
−
σ
s
h
=
(
σ
−
σ
s
)
(
1
E
′
−
1
E
)
\varepsilon^p=\frac{\sigma-\sigma_s}{h}=(\sigma-\sigma_s)(\frac{1}{E'}-\frac{1}{E})
εp=hσ−σs=(σ−σs)(E′1−E1)
那么相应的式(11)可以改写为
ε
=
ε
e
+
ε
p
=
σ
E
+
(
σ
−
σ
s
)
(
1
E
′
−
1
E
)
(16)
\varepsilon=\varepsilon^e+\varepsilon^p=\frac{\sigma}{E} +( \sigma-\sigma_s)(\frac{1}{E'}-\frac{1}{E})\tag{16}
ε=εe+εp=Eσ+(σ−σs)(E′1−E1)(16)
同时,有
0
<
1
h
=
1
E
′
−
1
E
<
1
E
′
(17)
0\lt\frac{1}{h}=\frac{1}{E'}-\frac{1}{E}\lt\frac{1}{E'}\tag{17}
0<h1=E′1−E1<E′1(17)
那么有以下结论,
h
>
E
′
(18)
h\gt E'\tag{18}
h>E′(18)
而其物理意义如下图
图
8
h
的物理意义(在
σ
−
ε
p
图中)
图8\quad h的物理意义(在\sigma-\varepsilon^p 图中)
图8h的物理意义(在σ−εp图中)