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轴为频率）