|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function [Nmatrixss,Nmatrixtt]=shapeFunction(s,z)
syms s z;
N1 = (1-s)*(1-z)*(-s-z-1)/4;
N2 = (1+s)*(1-z)*(s-z-1)/4;
N3 = (1+s)*(1+z)*(s+z-1)/4;
N4 = (1-s)*(1+z)*(-s+z-1)/4;
N5 = (1-z)*(1+s)*(1-s)/2;
N6 = (1+s)*(1+z)*(1-z)/2;
N7 = (1+z)*(1+s)*(1-s)/2;
N8 = (1-s)*(1+z)*(1-z)/2;
N=[N1 N2 N3 N4 N5 N6 N7 N8];
N1s = diff(N1,s);
N2s = diff(N2,s);
N3s = diff(N3,s);
N4s = diff(N4,s);
N5s = diff(N5,s);
N6s = diff(N6,s);
N7s = diff(N7,s);
N8s = diff(N8,s);
N1z = diff(N1,z);
N2z = diff(N2,z);
N3z = diff(N3,z);
N4z = diff(N4,z);
N5z = diff(N5,z);
N6z = diff(N6,z);
N7z = diff(N7,z);
N8z = diff(N8,z);
Nmatrixss=[N1s,N2s,N3s,N4s,N5s,N6s,N7s,N8s];
Nmatrixtt=[N1z,N2z,N3z,N4z,N5z,N6z,N7z,N8z];
return
>> [Nmatrixss,Nmatrixtt]=shapeFunction(s,z)
??? Undefined command/function 'shapeFunction'.
[ 本帖最后由 eight 于 2007-8-27 20:08 编辑 ] |
|