zouxinfree 发表于 2012-6-6 11:10

求这个二重积分的解、急。急、、、

本帖最后由 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)


zouxinfree 发表于 2012-6-6 11:21

各位大神。。。。

zouxinfree 发表于 2012-6-6 13:31

来位大神指导下啊。。。。

zouxinfree 发表于 2012-6-10 20:06

{:{17}:}

zouxinfree 发表于 2012-6-15 12:48

不给力啊。。。
页: [1]
查看完整版本: 求这个二重积分的解、急。急、、、