疲劳大作业.docVIP

疲劳与断裂 大作业 姓名 学号 2011-7-41、（1）计算材料的S-N、P-S-N曲线 应力水平（MPA） 编号 寿命次数 (95%) 900 2.9542 1 142000 5.1523 5.10986 0.1100 4.92891 2 94000 4.9731 3 186000 5.2695 4 123000 5.0899 5 116000 5.0645 800 2.9031 6 206000 5.3139 5.22708 0.1430 4.99185 7 98000 4.9912 8 223000 5.3483 9 192000 5.2833 10 158000 5.1987 700 2.8451 11 987000 5.9943 5.80078 0.1466 5.5596 12 466000 5.6683 13 834000 5.9211 14 512000 5.7092 15 514000 5.7110 650 2.8129 16 992000 5.9965 5.85148 0.1182 5.6570 17 875000 5.9420 18 702000 5.8463 19 535000 5.7283 20 555000 5.7443 600 2.7782 21 960000 5.9823 5.87256 0.1441 5.6355 22 592000 5.7723 23 469000 5.6712 24 865000 5.9370 25 1000000 6.0000 表中，标准偏差，，时，。 令，，最小二乘法拟合结果为，其中： ， 通过计算，最小二乘法拟合结果如下表： A B m C 50 3.9697 -0.1994 -0.953 5.015 95 3.9658 -0.2068 -0.948 4.836 中值S-N曲线： P-S-N曲线： 表中相关系数r的求法：, , 其中为载荷水平的种数。 （2）考虑到机匣的中心对称和上下对称，在ansys中建立如图模型，先择SOLID95单元，设置单元的轴对称特性，设置材料参数。所建立的模型如下： 机匣有限元模型1/2×10o 有限元模型局部放大 划分网格后开始施加约束条件及载荷，轴向载荷换算成面上的载荷，内压，计算结果如下： 应力云图分布 最大应力处局部放大应力云图分布 由上面的结果可看出，机匣危险点出现在机匣安装边与机匣壳体的过度处，最大应力为543MPa。因此可由第一问的结果计算出该机匣的寿命，结果如下： S-N曲线：，代入最大应力可得寿命N=1559677 P-S-N曲线（存活率95%）：，代入计算出寿命N=897836 2、编程实现雨流计数法，程序功能：实现复杂载荷谱典型段的循环计数。即：顺序输入典型段各点的峰值、谷值，程序输出：（1）循环；（2）变程；（3）均值。 从上图读取数据如下： 序号 1 2 3 4 5 6 7 8 9 10 1’ 值（MPa） 500 -100 300 -400 400 -200 100 -300 0 -200 500 本程序采用c语言编制。在文件CgmaVsTime.txt中输入上表中的峰值和谷值，在result.dat中看计算结果，包括循环、变程和均值。 C语言源程序： #includeiostream.h #includefstream.h #includemath.h #define N 100 void main() { int m; float a[N],fg[N],fgn[N],fl[N]={0},Ap[N],Av[N],hmax,he; int tmax,t,nt,np,tn[N],path[N][N]={0},npath[N]={0},fpath[N][2]={0},h; ifstream is1(CgmaVsTime.txt); ofstream os1(result.txt); is1m; for(int i=0;im;i++) //读入 名义应力-时间历程 {is1a[i];} tmax=0; //找出绝对值最大的峰/谷 for(int j=0;jm;j++) { if(fabs(a[j])fabs(a[tmax])) tmax=j; } for(int k=0;km;k++) //重新安