ansys（bottom to up）子结构

ansys（bottom to up）子结构

fini
/cle
/filn,full
/prep7
et,1,185
mp,ex,1,2e11
mp,prxy,1,0.3
k,1,,0.5
k,2,,0.5,0.25
k,3,,0.4,0.25
k,4,,0.4,0.05
k,5,,-0.4,0.05
k,6,,-0.4,0.25
k,7,,-0.5,0.25
k,8,,-0.5
a,1,2,3,4,5,6,7,8
vext,all,,,3
k,,2,0.25
k,,2,-0.25
k,,1,-0.25
k,,1,0.25
l,17,18
l,18,19
l,19,20
l,20,17
lfillt,28,25,0.1
lfillt,25,26,0.1
lfillt,26,27,0.1
lfillt,27,28,0.1
al,25,30,26,31,27,32,28,29
vext,11,,,,,0.1
vsbv,1,2
wpoffs,,,0.05
vsbw,all
wpoffs,,0.4
wprota,,90
vsbw,all
wpoffs,,,0.8
vsbw,all
wpoffs,1.5
wprota,,,90
vsbw,all
vsymm,z,all
vglue,all
esize,0.05
vsweep,all
save
fini

!生成超单元1
/filn,gen1
resume,full,db
/sol
antype,substr
seopt,SE1
nsel,s,loc,x,0
nsel,a,loc,x,3
m,all,all
nsel,all
save
solve
fini

!生成超单元2
/filn,gen2
resume,full,db
/prep7
vgen,,all,,,3,,,,,1
fini
/sol
antype,substr
seopt,SE2
nsel,s,loc,x,3
nsel,a,loc,x,6
m,all,all
nsel,s,loc,x,4.5
nsel,r,loc,y,0.5
m,all,uy
nsel,all
save
solve
fini


!生成超单元3
/filn,gen3
resume,full,db
/prep7
vgen,,all,,,6,,,,,1
fini
/sol
antype,substr
seopt,SE3
nsel,s,loc,x,6
nsel,a,loc,x,9
m,all,all
nsel,all
save
solve
fini

!使用部分
/cle
/filn,full
resume
fini
/filn,use
/prep7
et,2,50
type,2
setran,SE1,,100000,SE4!由于已经存在非超级单元。因此不能直接使用se,se1
se,SE4
setran,SE2,,200000,SE5
se,SE5
setran,SE3,,300000,SE6
se,SE6
esel,s,type,,2！仅仅选择超级单元部分节点耦合
nsle,s,all
nsel,r,loc,x,3
nsel,a,loc,x,6
cpintf,all,0.03
type,1
vgen,3,all,,,3！生成比较实体模型
vgen,,all,,,,-6,,,,1
nslv,s,1
nummrg,all
numcmp,all
alls
save
fini
/sol
nsel,s,loc,x,0
nsel,a,loc,x,9
d,all,all
alls
nsel,s,loc,x,4.5
nsel,r,loc,y,0.5
f,all,fy,-1000/ndinqr(0,13)!超级单元部分加载
nsel,s,loc,x,4.5
nsel,r,loc,y,-5.5
f,all,fy,-1000/ndinqr(0,13)!比较模型部分加载
alls
esel,s,type,,2
*get,emax,elem,0,num,max
*get,emin,elem,0,num,min
alls
sfe,emin,1,selv,,1
sfe,emin+1,2,selv,,1
sfe,emax,3,selv,,1
alls
solve
/post1
plnsol,s,eqv
esel,s,type,,1
nsort,u,sum,0,0,all 
!获取比较模型最大位移
*get,max_u,sort,0,max
nsort,s,eqv,0,0,all 
!获取比较模型最大应力
*get,max_s,sort,0,max
save
fini

!扩展超单元1
/cle
/filn,gen1
resume
/sol
expass,on
seexp,SE4,use
expsol,1
solve
save
fini
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all 
!获取超级单元1最大位移
*get,max_u1,sort,0,max
nsort,s,eqv,0,0,all 
!获取超级单元1最大应力
*get,max_s1,sort,0,max
save
fini


/cle
/filn,gen2
resume
/sol
expass,1
seexp,SE5,use
expsol,1
solve
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all 
!获取超级单元2最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all 
!获取超级单元2最大应力
*get,max_s1,sort,0,max
fini
save

/cle
/filn,gen3
resume
/sol
expass,1
seexp,SE6,use
expsol,1
solve
/post1
plnsol,s,eqv
nsort,u,sum,0,0,all 
!获取超级单元3最大位移
*get,max_u2,sort,0,max
nsort,s,eqv,0,0,all 
!获取超级单元3最大应力
*get,max_s1,sort,0,max
save
fini

基础单元网格：

第一个/三个超级单元：


第二个超级单元：

超级单元与整体比较模型：

中部应力结果：

比较模型位移结果：

比较模型最大应力结果：

超级单元扩展最大应力结果：


列表显示各应力与位移：