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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

本帖最后由 zouxinfree 于 2012-6-6 11:23 编辑

QQ截图20120606111544.jpg

这个式方程[img=621,86]file:///d:/Documents/Tencent%20Files/125937112/Image/I3I[%7BZ5W1%60[75(V[ZIT8SV8.jpg[/img]
其中
QQ截图201206061116161.jpg

求高人指点，。。。
这是我编的程序。。但是积分老出错。。。
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

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

不给力啊。。。

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 求这个二重积分的解、急。急、、、

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

本帖被以下淘专辑推荐: