ansys 子结构法《附命令流》
欢迎与大家一起探讨有关ansys动态子结构主题，我最近用子结构法和整体建模法分别对一个结构做了谐响应分析，位移频率相差甚远，看了help数遍还是解决不了问题，先把命令流及结果附上，希望能和大家一起共勉，进步，也可以加我qq 1460864015 一起探讨……子结构法-finish
/clear
/filname,full
/prep7
!1.定义变量
l=14
w=24
r=sqrt(l*l/4+w*w/4)
an=1
at=0.5
h=0.5
e=2.6e8 !介质弹性模量
v=0.42 !介质泊松系数
g=E/2*(1+v) !介质剪切模量
d=1900 !介质密度
a=an/at
v1=(a-2)/2*(a-1) !等效一致粘弹性边界单元等效泊松系数
lame=(1-v1)/((1+v1)*(1-2*v1))
g1=at*h*g/r !等效一致粘弹性边界单元等效剪切模量
e1=an*h*g*lame/r !等效一致粘弹性边界单元等效弹性模量
lame1=(1-v)/((1+v)*(1-2*v))
Cp=sqrt(lame1/d) !压缩波速
Cs=sqrt(g/d) !剪切波速
kt=at*g/r !切向粘弹性人工边界等效物理系统的弹簧系数Kt
ct=d*cs !切向粘弹性人工边界等效物理系统的阻尼系数Ct
kn=an*g/r !法向粘弹性人工边界等效物理系统的弹簧系数Kn
cn=d*cp !法相粘弹性人工边界等效物理系统的阻尼系数Cn
b=d*r*(cs/at+cp/an)/2*g !等效一致粘弹性边界单元阻尼比例系数b
!2.建立模型
et,1,plane42
mp,ex,1,2.6e8
mp,prxy,1,0.42
mp,dens,1,1900
r,1
mp,ex,2,e1
mp,gxy,2,g1
mp,damp,2,b
mp,prxy,2,v1
rectng,-0.5,0.5,0,16
rectng,0.5,7,0,16
rectng,0.5,7,16,24
rectng,-0.5,-7,0,16
rectng,-0.5,-7,16,24
aglue,all
esize,0.5
mshape,0,2d
mshkey,1
amesh,all
!3.施加粘弹性边界
lsel,s,loc,y,7,9$lsel,a,loc,y,17,21
lsel,u,loc,x,-0.5,0.5
nsll,s,1
esln,s
mpchg,2,all
lsel,s,loc,y !选择水平线
nsll,s,1 !选择水平线上的节点
esln,s !选择水平线上与节点相关的单元
mpchg,2,all !改变所选单元的材料属性，即施加粘弹性边界
allsel,all
save
finish
!生成超单元
/filname,gen
resume,full,db
/sol
antype,substr
seopt,cd,3,1
alphad,0.5
betad,0.0002
lsel,s,loc,y,17,21$lsel,r,loc,x,-0.5,0.5
nsll,s,1
m,all,all
lsel,s,loc,y,15,17$lsel,r,loc,x,-0.5,0.5
nsll,s,1
m,all,all
lsel,s,loc,y,7,9$lsel,r,loc,x,-7
nsll,s,1
d,all,all
lsel,s,loc,y,7,9$lsel,r,loc,x,7
nsll,s,1
d,all,all
lsel,s,loc,y,17,21$lsel,r,loc,x,-7
nsll,s,1
d,all,all
lsel,s,loc,y,17,21$lsel,r,loc,x,7
nsll,s,1
d,all,all
lsel,s,loc,y !选择水平线
nsll,s,1 !选择水平线上的节点
d,all,all
allsel,all
save
solve
finish
!使用部分
/clear
/filname,use
/prep7
et,1,matrix50
et,2,plane42
mp,ex,1,2.6e10$mp,prxy,1,0.167$mp,dens,1,2500
!mp,ex,2,2.6e8
!mp,prxy,2,0.42
!mp,dens,2,1900
r,1
r,2
rectng,-0.5,0.5,16,24
esize,0.5
mshape,0,2d
mshkey,1
amesh,all
*GET,bb,NODE,,NUM,MAX !获得最大节点号从而确定最小偏移量
*stat,bb
TYPE,1
SETRAN,cd,0,51,cd1
se,cd1
nsel,all
CPINTF,ALL
eplot
selist,cd1
allsel
save
FINISH
/solu$antype,3$hropt,full !定义谐响应分析
alphad,0.5
betad,0.0002 !定义质量阻尼系数
f,21,fy,10e4
harfrq,0,100 !定义谐响应计算强制频率范围
nsubst,300 !定义谐响应解的数目，数目大可以使频率-响应曲线光滑
elist,all,,,1 !查看超单元号
sfe,33,1,selv,,1 !施加单元荷载向量
allsel
solve
finish
/post1
set,list
finish
!扩展超单元
/clear
/filname,gen
resume
/solu
outres,all,all
outpr,all,all
expass,on
seexp,cd1,use,off
!expsol,1
numexp,all !扩展当前载荷步内所有频率范围内的子步解
solve
finish
!时程后处理
/post26
nsol,2,980,u,y
nsol,3,980,u,x
/grid,1
/gropt,view,1
plcplx,0
plva,2,3
!进行模态分析
/solu
antype,modal
modopt,lanb,30
mxpand,30,,,yes
solve
finish
/post1
set,list
整体建模法-
finish
/clear
/filname,full
/prep7
!1.定义变量
l=14
w=24
r=sqrt(l*l/4+w*w/4)
an=1
at=0.5
h=0.5
e=2.6e8 !介质弹性模量
v=0.42 !介质泊松系数
g=E/2*(1+v) !介质剪切模量
d=1900 !介质密度
a=an/at
v1=(a-2)/2*(a-1) !等效一致粘弹性边界单元等效泊松系数
lame=(1-v1)/((1+v1)*(1-2*v1))
g1=at*h*g/r !等效一致粘弹性边界单元等效剪切模量
e1=an*h*g*lame/r !等效一致粘弹性边界单元等效弹性模量
lame1=(1-v)/((1+v)*(1-2*v))
Cp=sqrt(lame1/d) !压缩波速
Cs=sqrt(g/d) !剪切波速
kt=at*g/r !切向粘弹性人工边界等效物理系统的弹簧系数Kt
ct=d*cs !切向粘弹性人工边界等效物理系统的阻尼系数Ct
kn=an*g/r !法向粘弹性人工边界等效物理系统的弹簧系数Kn
cn=d*cp !法相粘弹性人工边界等效物理系统的阻尼系数Cn
b=d*r*(cs/at+cp/an)/2*g !等效一致粘弹性边界单元阻尼比例系数b
!2.建立模型
et,1,plane42
mp,ex,1,2.6e8
mp,prxy,1,0.42
mp,dens,1,1900
r,1
mp,ex,2,2.6e10$mp,prxy,2,0.167$mp,dens,2,2500
r,2
mp,ex,3,e1
mp,gxy,3,g1
mp,damp,3,b
mp,prxy,3,v1
rectng,-7,7,0,24
wpoffs,0.5,16
wprota,,90
asbw,all
wprota,,,90
asbw,all
wpoffs,8,,-1
asbw,all
wpcsys,-1
numcmp,all
asel,s,area,,6
aatt,2,2,1
allsel,all
esize,0.5
mshape,0,2d
mshkey,1
amesh,all
!施加粘弹性边界
lsel,s,loc,y !选择水平线
nsll,s,1 !选择水平线上的节点
esln,s !选择水平线上与节点相关的单元
mpchg,3,all !改变所选单元的材料属性，即施加粘弹性边界
lsel,s,loc,y,7,9$lsel,a,loc,y,17,21
lsel,u,loc,x,-0.5,0.5
nsll,s,1
esln,s
mpchg,3,all
allsel
!定义约束
lsel,s,loc,y !选择水平线
nsll,s,1 !选择水平线上的节点
d,all,all
lsel,s,loc,y,7,9$lsel,a,loc,y,17,21
lsel,u,loc,x,-0.5,0.5
nsll,s,1
d,all,all
allsel
!进行模态分析
/solu
antype,modal
modopt,lanb,30
mxpand,30,,,yes
solve
finish
/post1
set,list
!谐响应分析
/solu$antype,3$hropt,full
alphad,0.5
betad,0.0002 !定义质量阻尼系数
f,1406,fy,10e4
harfrq,0,100 !定义谐响应计算强制频率范围
nsubst,100 !定义谐响应解的数目，数目大可以使频率-响应曲线光滑
allsel
solve
finish
!时程后处理
/post26$nsol,2,1171,u,y !位移点（4,24,0）
nsol,3,1171,u,x
/grid,1
/gropt,view,1
plcplx,0
plva,2,3 !绘制响应-频率曲线（x轴为频率）
页:
[1]