马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
最近需要从ansys模型中提取结构的刚度和质量矩阵,利用论坛中的方法进行了尝试,并与理论值进行了比较验证。利用弹簧单元combin14和质量单元mass21组成简单的振动模型。验证了该方法的正确性。
finish
/clear
/prep7
Et,1,combin14
Et,2,Mass21
keyopt,2,3,4
k1=100
k2=200
Mass=100
R,1,k1
R,2,k2
R,3,Mass
N,1
N,2,0.1
N,3,0.2
Type,1
real,1
E,1,2
real,2
E,2,3
Type,2
real,3
E,1
E,2
E,3
/solu
d,1,all
nsel,s,,,2,3
d,all,uy,,,,,uz
F,3,fx,50
allsel
antype,substr
!seopt,file,1,1
seopt,matrixname,2
m,all,all
/output,matrix
solve
/output
selist,matrixname,3
但是,当用beam3建立模型提取某一单元的质量矩阵时与理论值不符,但刚度矩阵正确。
finish
/clear
/CONFIG,NRES,2e4
/prep7
/title, Static, Modal, Harmonic response of a MEMS structure
/com
L=0.01 ! beam length
b=0.02 ! beam width
h=0.002 ! beam height
I=b*h**3/12 ! beam moment of inertia
E=70e9 ! modulus
dens=2700 ! density
et,1,3 ! 2-D beam element
r,1,b*h,I,h ! beam properties
mp,ex,1,E
mp,dens,1,dens
mp,prxy,1,0.345
K,1,0,
K,2,L
LSTR,1,2
LEsize,1,,,10
Lmesh,1
nplot
/solu
d,1,ux,,,,,uy
d,2,ux,,,,,uy
allsel
antype,substr
!seopt,file,1,1
seopt,matrixname,2
esel,s,,,5
nsle,s,
m,all,UY,,,Rotz
/output,matrix
solve
/output
selist,matrixname,3
下面是根据根据beam单元理论公式,用matlab 编写的beam单元的刚度和质量矩阵数值
clear all;
%close all;
format short e
b=0.02;h=0.002;l=0.001;
dens=2700;E=70e9;
I=b*h^3/12;
A=b*h;
K=(E*I/l^3)*[12 6*l -12 6*l;
6*l 4*l^2 -6*l 2*l^2;
-12 -6*l 12 -6*l;
6*l 2*l^2 -6*l 4*l^2];
M=(dens*A*l/420)*[156 22*l 54 -13*l;
22*l 4*l^2 13*l -3*l^2;
54 13*l 156 -22*l;
-13*l -3*l^2 -22*l 4*l^2]
不知道哪里出现错误,希望大家给与帮助!
[ 本帖最后由 xinglong-liu 于 2008-6-8 21:33 编辑 ] |