test0.m

薄壁结构有限元计算

% 输入转动惯量（x向、y向、二向混合、翘曲、x向翘曲、y向翘曲）
Ix=input('Please input the x-direction of inertia-moment of the section: ');
Iy=input('Please input the y-direction of inertia-moment of the section: ');
Ixy=input('Please input the xy-direction of inertia-moment of the section: ');
If=input('Please input the warp of inertia-moment of the section: ');
Ixf=input('Please input the x-direction warp of inertia-moment of the section: ');
Iyf=input('Please input the y-direction warp of inertia-moment of the section: ');

% 输入面积（截面、x向剪切、y向剪切、混合剪切、x向翘曲剪切、y向翘曲剪切、翘曲剪切）
A=input('Please input the area of the section: ');
Ax=input('Please input the x-direction of shear-area of the section: ');
Ay=input('Please input the y-direction of shear-area of the section: ');
Axy=input('Please input the mixed shear-area of the section: ');
Axr=input('Please input the x-direction warp of shear-moment of the section: ');
Ayr=input('Please input the y-direction warp of shear-moment of the section: ');
Ar=input('Please input the warp of shear-moment of the section: ');

% 输入x向静矩、y向静矩、弹性模量、剪切模量、扭转常数、密度、固有频率
Sx=input('Please input the x-direction of shear-moment of the section: ');
Sy=input('Please input the y-direction of shear-moment of the section: ');
E=input('Please input the elastic modulus of the material: ');
G=input('Please input the shear modulus of the material: ');
J=input('Please input the torsion-constant of the section: ');
rho=input('Please input the density of the material: ');
omega=input('Please input the inherence frequency of the construction: ');

% 计算各参数与AA矩阵
b2=rho*omega*omega*A;
b3=rho*omega*omega*Sy;
b4=rho*omega*omega*Sx;
b5=G*Ax;
b6=G*Axy;
b7=G*Axr;
b8=-G*Ax+rho*omega*omega*Iy;
b9=G*Axy-rho*omega*omega*Ixy;
b10=G*Axr-rho*omega*omega*Iyf;
b11=G*Ay;
b12=G*Ayr;
b13=-G*Ay+rho*omega*omega*Ixy;
b14=-G*Ayr+rho*omega*omega*Ixf;
b15=rho*omega*omega*If;
b16=G*Ar;
b17=-G*Ar+rho*omega*omega*If;
AA=[0 1 0 0 0 0 0 0 0 0 0 0 0 0;
   b2 0 0 0 -b3 0 0 0 b4 0 0 0 0 0;
   0 0 0 1 0 0 0 0 0 0 0 0 0 0;
   0 0 b2 0 0 -b5 0 0 0 b6 -b4 0 0 b7;
   0 0 0 0 0 1 0 0 0 0 0 0 0 0;
   -b3 0 0 b5 b8 0 0 b6 b9 0 0 b7 b10 0;
   0 0 0 0 0 0 0 1 0 0 0 0 0 0;
   0 0 0 0 0 -b6 b2 0 0 b11 b3 0 0 b12;
   0 0 0 0 0 0 0 0 0 1 0 0 0 0;
   b4 0 0 -b6 b9 0 0 -b11 b13 0 0 -b12 b14 0;
   0 0 0 0 0 0 0 0 0 0 0 1 0 0;
   0 0 -b4 0 0 -b7 b3 0 0 b12 b15 0 0 b16;
   0 0 0 0 0 0 0 0 0 0 0 0 0 1;
   0 0 0 -b7 b10 0 0 -b12 b14 0 0 -b16 b17 0];

% 计算各参数与BB矩阵
k2=-E*A;
k3=E*Sy;
k4=-E*Sx;
k5=-G*Ax;
k6=-G*Axy;
k7=-G*Axr;
k8=-E*Iy;
k9=E*Ixy;
k10=E*Iyf;
k11=-G*Ay;
k12=-G*Ayr;
k13=-E*Ix;
k14=-E*Ixf;
k15=-G*J-G*Ar;
k16=-E*If;
BB=[1 0 0 0 0 0 0 0 0 0 0 0 0 0;
   0 k2 0 0 0 k3 0 0 0 k4 0 0 0 0;
   0 0 1 0 0 0 0 0 0 0 0 0 0 0;
   0 0 0 k5 0 0 0 k6 0 0 0 k7 0 0;
   0 0 0 0 1 0 0 0 0 0 0 0 0 0;
   0 k3 0 0 0 k8 0 0 0 k9 0 0 0 k10;
   0 0 0 0 0 0 1 0 0 0 0 0 0 0;
   0 0 0 k6 0 0 0 k11 0 0 0 k12 0 0;
   0 0 0 0 0 0 0 0 1 0 0 0 0 0;
   0 k4 0 0 0 k9 0 0 0 k13 0 0 0 k14;
   0 0 0 0 0 0 0 0 0 0 1 0 0 0;
   0 0 0 k7 0 0 0 k12 0 0 0 k15 0 0;
   0 0 0 0 0 0 0 0 0 0 0 0 1 0;
   0 0 0 0 0 k10 0 0 0 k14 0 0 0 k16];

% 计算特征值与特征向量
[Z,lambda]=eig(AA,BB,'qz');