|
我刚好在算曲率模态,这是一个matllab的子程序,你看看吧
function =qulv(FI)
% 模态曲率差值法,需要比较密的测点,以便使用中心差分。
global pi noe use_femorder use_emaorder rou E l
h=l/noe; % !!!!!!!!!!!!!!!
row=length(FI(:,1));
col=length(FI(1,:));
U=zeros(row,col);
for j=1:col
for i=1
U(i,j)=(FI(i+2,j)+FI(i,j)-2*FI(i+1,j))/h^2;
end
for i=2:row-1
U(i,j)=(FI(i-1,j)+FI(i+1,j)-2*FI(i,j))/h^2;
end
for i=row
U(i,j)=(FI(i-2,j)+FI(i,j)-2*FI(i-1,j))/h^2;
end
end
end
只需要损伤前后的模态,就可以利用这个子程序分别算出损伤前后的模态曲率。
要是有不明白的,可以联系我:lwming2001@sjtu.edu.cn
大家互相交流一下。 |
|