求这个二重积分的解、急。急、、、
本帖最后由 zouxinfree 于 2012-6-6 11:23 编辑这个式方程file:///d:/Documents/Tencent%20Files/125937112/Image/I3I[%7BZ5W1%60
其中
求高人指点,。。。
这是我编的程序。。但是积分老出错。。。
syms xi;
R=0.015;L_1=0.057;
a=0.9*R;
b=(R^2-(R-a)^2)^(1/2); %裂纹宽度
eta=a-R+(R^2-xi^2)^(1/2); %局部裂纹深度
h=2*(R^2-xi^2)^(1/2);
v=0.3;
E=2.1*10^11;
L_R=L_1/(R*3*(1-v^2));
F_2=((2*h/pi*eta*tan(pi*eta/2*h))^(1/2))*(0.923+0.199*(1-sin(pi*eta)/(2*h))^4)/cos(pi*eta/(2*h))
F_1=((2*h/pi*eta*tan(pi*eta/2*h))^(1/2))*(0.752+2.02*(eta/h)+0.37*(1-sin(pi*eta/(2*h)))^3)/cos(pi*eta/(2*h))
y1=quadl(@(xi)arrayfun(@(xi)quadl(@(eta)(32*pi*R^3*E/(1-v^2))*(1-(xi/R)^2)*(eta/R)*F_2^2,0,eta/R),xi),b/R,-b/R)
各位大神。。。。 来位大神指导下啊。。。。 {:{17}:} 不给力啊。。。
页:
[1]