코드에 대해서

Question

chae sungwon 2021-10-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1569963-

回答： Angelo Yeo 2023-12-27

syms a b c d e;
f1=-C1*c*b+C2*a*b^2+2*C3*a*b-B11*c+D11*a-C4*c*b+2*C5*a*b^2-C6*d*b+D12*b-C7*d*b+C8*a*b^2+C9*b^2+(Ly.^2.*Nx.*b)./48+Mx+(Lx.^2.*Ny.*b)./48;
f2=-C1*a*c+C2*a^2*b+C3*a^2-C4*a*c+2*C5*a^2*b+D12*a-C6*d*a-C7*d*a+C8*a^2*b+2*C9*a*b-B22*d+D22*b+(Ly.^2.*Nx.*a)./48+(Ly.^2.*Ny.*a)./48+My;
f3=A11*c-C1*a*b-B11*a+A12*d-C4*a*b-Nx;
f4=A12*c-C6*a*b-B22*b+A22*d-C7*a*b-Ny;
S=solve(f1,f2,f3,f4,a,b,c,d);
a=double(S.a);
b=double(S.b);
c=double(S.c);
d=double(S.d);
for j=1:length(a)
    if isreal(a(j))==1
        m=0;
        n=0;
        for x=0:Lx
            m=m+1;
            for y=0:Ly
                n=n+1;
                u(m,n,j)=double(c(j)*x-((a(j)^2*x^3)/6)-((a(j)*b(j)*x*y^2)/4));
                v(m,n,j)=double(d(j)*y-((b(j)^2*y^3)/6)-((a(j)*b(j)*x^2*y)/4));
                w(m,n,j)=double(0.5*((a(j)*x^2)+(b(j)*y^2))); 
                if y==Ly
                    n=0;
                end
            end
        end
    end
end
e=w;
Df=diff(e,x);
Aa=e(:,:);
plot3(u(:,:,1),v(:,:,1),w(:,:,1),'r')
xlabel('x (mm)')
ylabel('y (mm)')
zlabel('z (mm)')
title('Laminate')