Error in my FEM Matlab code. Please help

3 次查看（过去 30 天）

Virajan Verma 2018-10-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/422108-error-in-my-fem-matlab-code-please-help

评论： Torsten 2018-10-4

for i=1:2
    fepsn(i)=i;
end
%%Loop for epsilon values
for icounter=1:2
    epsn=fepsn(icounter)
end
%%Shape function for axial effect
N1(epsn)=(-1/2)*epsn*(1-epsn)
N2(epsn)=(1+epsn)*(1-epsn)
N3(epsn)=(1/2)*epsn*(1+epsn)
syms epsn
%%First Derivatives of Shape Function
N1=diff(N1(epsn),epsn)   (Error here)
N2=diff(N2(epsn),epsn)
N3=diff(N3(epsn),epsn)
%%%%Second Derivatives of Shape Function
N11=diff(N1(epsn));
N22=diff(N2(epsn));
N33=diff(N3(epsn));

8 个评论
显示 6更早的评论隐藏 6更早的评论

Virajan Verma 2018-10-4

I have removed it now. But error is there as stated in my previous comment

Torsten 2018-10-4

在 MATLAB Online 中打开

syms epsn N1(epsn) N2(epsn) N3(epsn)
N1(epsn) = (-1/2)*epsn*(1-epsn);
N2(epsn) = (1+epsn)*(1-epsn);
N3(epsn) = (1/2)*epsn*(1+epsn);
DN1 = diff(N1,epsn);
DN2 = diff(N2,epsn);
DN3 = diff(N3,epsn);
D2N1 = diff(DN1,epsn);
D2N2 = diff(DN2,epsn);
D2N3 = diff(DN3,epsn);