三参数 S-N 曲线拟合及MATLAB程序

更多阅读：sppy.site

S-N 曲线：应力 S 与疲劳寿命 N 之间关系的曲线。

三参数 S-N 曲线方程
$$(S-S_{\mathrm{f}}{)}^m\times N=C\text{\hspace{1em}(1)}$$

线性拟合方法

将式(1)取对数，可得

$$m\lg (S-S_{\mathrm{f}})+\lg N=\lg C\text{\hspace{1em}(2)}$$
令 $x=\lg (S-S_{\mathrm{f}})$ 、$y=\lg N$ 、$a=\lg C$ 、$b=-m$ ，可得
$$y=a+bx\text{\hspace{1em}(3)}$$
对于一系列的应力值 $S_i$ 、疲劳寿命值 $N_i$ ($i=1,2,3,\cdots ,n$)，由最小二乘法可得
$$\{\begin{array}{l}b=L_{xy}/L_{xx}\\ a=\bar{y}-\bar{x}b\end{array}\text{\hspace{1em}(4a)}$$
线性相关系数 $R$ 的平方为
$$R^2=\frac{L_{xy}^2}{L_{xx}{L}_{yy}}\text{\hspace{1em}(4b)}$$
其中
x i = lg ⁡ ( S i − S f ) y i = lg ⁡ N i x ˉ = 1 n ∑ i = 1 n x i y ˉ = 1 n ∑ i = 1 n y i L x x = ∑ i = 1 n x i 2 − 1 n ( ∑ i = 1 n x i ) 2 L y y = ∑ i = 1 n y i 2 − 1 n ( ∑ i = 1 n y i ) 2 L x y = ∑ i = 1 n x i y i − 1 n ( ∑ i = 1 n x i ) ( ∑ i = 1 n y i ) x_i=\lg(S_i-S_\mathrm{f}) \qquad y_i=\lg N_i\\[8pt] \bar{x}=\frac{1}{n}\sum_{i=1}^nx_i \qquad \bar{y}=\frac{1}{n}\sum_{i=1}^ny_i\\[8pt] L_{xx}=\sum_{i=1}^nx_i^2-\frac{1}{n}\Big(\sum_{i=1}^nx_i\Big)^2\\[8pt] L_{yy}=\sum_{i=1}^ny_i^2-\frac{1}{n}\Big(\sum_{i=1}^ny_i\Big)^2\\[8pt] L_{xy}=\sum_{i=1}^nx_iy_i-\frac{1}{n}\Big(\sum_{i=1}^nx_i\Big)\Big(\sum_{i=1}^ny_i\Big) xi​=lg(Si​−Sf​)y