Rectangular Kirchhoff-Love plates
The example for a Sinusoidal load case:
clc
clear
M=50;
N=50;
q0=30;
D=1;
a=1;
b=2;
dx=a/M;
dy=b/N;
m=1;
n=1;
aaa=1/(((m/a)^2+(n/b)^2)^2);
aaa1=(pi*dx*m)/a;
aaa2=(pi*dy*n)/b;
aa=q0/((pi^4)*D);
for i=1:M;
for j=1:N;
w(i,j)=aa*aaa*sin(aaa1*i)*sin(aaa2*j);
end
end
w(1:M,1)=0;
w(1:M,N)=0;
w(1,1:N)=0;
w(M,1:N)=0;
[x,y]=meshgrid(0:dx:a,0:dy:b);
surf(x(1:M,1:N),y(1:M,1:N),w)
xlabel('x')
ylabel('y')
zlabel('')
axis equal
clear
M=50;
N=50;
q0=30;
D=1;
a=1;
b=2;
dx=a/M;
dy=b/N;
m=1;
n=1;
aaa=1/(((m/a)^2+(n/b)^2)^2);
aaa1=(pi*dx*m)/a;
aaa2=(pi*dy*n)/b;
aa=q0/((pi^4)*D);
for i=1:M;
for j=1:N;
w(i,j)=aa*aaa*sin(aaa1*i)*sin(aaa2*j);
end
end
w(1:M,1)=0;
w(1:M,N)=0;
w(1,1:N)=0;
w(M,1:N)=0;
[x,y]=meshgrid(0:dx:a,0:dy:b);
surf(x(1:M,1:N),y(1:M,1:N),w)
xlabel('x')
ylabel('y')
zlabel('')
axis equal
The equations coded are taken from the following link:
Unless otherwise noted, all content on this site is @Copyright by Ahmed Al Makky 2012-2013 - http://cfd2012.com