断裂力学作业.docVIP

PAGE 3 PAGE 17 理论与应用断裂力学作业 1、解析法求裂尖应力强度因子 如图所示, 无限大板, 中心穿透裂纹, 位于裂纹表面处作用集中法向力P和集中切向力Q, 已知： 用极限公式求解裂纹尖端应力强度因子。 解：已知： 将题目中的带入上式得： 所以，得： 2、Green函数法求裂尖应力强度因子 如图, 已知 （平面应力） （平面应变） P为作用在裂纹中心处y向对称集中力，L为集中力作用点到裂纹中点的y向距离。 当a＝5 mm、15 mm、25mm，L＝5 mm、15 mm、25mm 解：取P＝1000N，泊松比为0.3。分别计算平面应力和平面应变两种情况下的应力强度因子。 由 得 采用复化辛普森数值积分公式。 分析：复化辛普森数值积分公式需要用到被积函数在积分上界的函数值，由K1的表达式可以看出被积函数在积分上界处的函数值是无穷大。但是通过被积函数的导数得知当x→a时导数为无穷大，说明被积函数存在垂直渐近线，且渐近线是y＝a。根据函数的特性（根号下的式子互为倒数），使得在x→a的很小的左临域（a-δ，a）（δ0）内函数值急剧趋于无穷大。也就导致了被积函数在左临域（a-δ，a）内函数值趋于无穷大。因此，可以用f(a-δ) 代替f(a),所得的积分误差应该不会太大。 C＋＋程序如下： #include stdio.h #include math.h #define A0 0.0 //积分初始点 #define P 1000 //外载 #define N 5000 //积分区间等分数 ！！！N 必须为偶数。 #define PSB 0.3 double A[3]={0.005,0.015,0.025}; double L[3]={0.005,0.015,0.025}; int main(int argc, char* argv[]) { double sun=0.0,y,x,a1,a2,a3,a4,t,buchang,terta; int i,j=1,k,m; FILE *fp1,*fp2; if((fp1=fopen(平面应力状态.txt,w))==NULL) printf(无法打开这个文件!\n); if((fp2=fopen(平面应变状态.txt,w))==NULL) printf(无法打开这个文件!\n); loop1: printf(请输入1或2：1表示平面应力状态，2表示平面应变状态\n); scanf(%d,i); if(!(i==1||i==2)) { printf(您的输入有误，请重新输入\n); goto loop1; } else { if(i==1) t=(3.0-PSB)/(1.0+PSB); if(i==2) t=3.0-4.0*PSB; for(k=0;k3;k++) { for(m=0;m3;m++) { buchang=(A[k]-A0)/N; x=A0+buchang; terta=0.07*buchang; //临域半径 sun=0.0; j=1; while(1) { if(j=N-1) //进行N－1次循环 { if(j%2) //奇数项求和 { a1=sqrt(double(A[k]+x)/double(A[k]-x))+sqrt(double(A[k]-x)/double(A[k]+x)); a2=(t-1)/(t+1)+(4.0/(t+1.0))*(L[m]*L[m]/double(x*x+L[m]*L[m])); a3=3(x*x+L[m]*L[m]); a4=sqrt(3A[k]); y=(a1*(double(P))*L[m]*a2/a3)/double(a4); //被积函数表达式 sun+=4*y; x+=buchang; } if(!(j%2)) //偶数项求和 {