ANSYS经典例子---（4）
31.带缺口梁的屈曲分析
一带缺口的梁,一端受压确定其临界载荷.梁尺寸10&10&300，缺口尺寸10&4&10，受力Fz=5000N.
得位移----载荷如下:得临界载荷Fcr=3606.74N
*set,length,10
*set,width,10
*set,height,300
*set,crac_len,4
*set,crac_wid,10
*set,crac_hei,10
*set,force,5000
ET,1,solid45
MP,EX,1,200e3
MP,PRXY,1,0.3
TBDATA,1,400
TBDATA,2,20E3
BLOCK,0,length,0,width,0,height
K,9,0,0,(height/2)-(crac_hei/2)
BLOCK,0,crac_len,0,crac_wid,0,crac_hei
KWPAVE, 13
KWPAVE, 15
wpro,,,90.000000
FLST,2,2,6,ORDE,2
LESIZE,9, , ,20, , , , ,1
LESIZE,12, , ,20, , , , ,1
*do,i,25,28
LESIZE,i, , ,20, , , , ,1
*do,i,35,36
LESIZE,i, , ,20, , , , ,1
LESIZE,i+3, , ,20, , , , ,1
LESIZE,43, , ,20, , , , ,1
LESIZE,44, , ,20, , , , ,1
LESIZE,1, , ,6, , , , ,1
LESIZE,3, , ,6, , , , ,1
LESIZE,6, , ,6, , , , ,1
LESIZE,8, , ,6, , , , ,1
LESIZE,13, , ,6, , , , ,1
LESIZE,15, , ,6, , , , ,1
LESIZE,18, , ,6, , , , ,1
*do,i,20,21
LESIZE,i, , ,6, , , , ,1
LESIZE,32, , ,6, , , , ,1
LESIZE,37, , ,6, , , , ,1
LESIZE,42, , ,6, , , , ,1
LESIZE,14, , ,2, , , , ,1
LESIZE,16, , ,2, , , , ,1
LESIZE,17, , ,2, , , , ,1
LESIZE,19, , ,2, , , , ,1
*do,i,45,48
LESIZE,i, , ,2, , , , ,1
*do,i,10,11
LESIZE,i, , ,4, , , , ,1
LESIZE,24, , ,4, , , , ,1
LESIZE,29, , ,4, , , , ,1
*do,i,30,31
LESIZE,i, , ,4, , , , ,1
LESIZE,i+10, , ,4, , , , ,1
!Define constrain and load
FLST,2,2,5,ORDE,2
FITEM,2,10
FITEM,2,27
DA,P51X,ALL,0
F,1072,FZ,-5000
ANTYPE,static
NSUBST,12,0,0
ARCLEN,1,0,0
AUTOTS,-1.0
OUTRES,ALL,all
/STATUS,SOLU
.幕墙企业玻璃简化计算
/encrypt,yuhui4,Magic Glass,mac,d:\
/CLEAR,NOSTART
multipro,'start',10
*cset,1,3,H,'The height of glass(unit:m)',2
*cset,4,6,W,'The width of glass(unit:m)',2
*cset,7,9,a1,'Angle left side to down(unit:C)',90
*cset,10,12,a2,'Angle right side to down(unit:C)',90
*cset,13,15,hole,'the number of hole',4
*cset,16,18,type,'Layered Glass Type(!SeePeak!)',1
*cset,19,21,thi,'Thickness of inner Glass(unit:m)',0.01
*cset,22,24,tho,'Thickness of outer Glass(unit:m)',0.008
*cset,25,27,Wk,'Designed wind pressure(unit:N/m2)',1000
*cset,61,62,'Please Fill in related blanks to',' build FEA model
And apply loads'
*cset,28,30,Seismic,'Designed seismic accelerate(m/s2)',0.08
*cset,63,64,'FOR LAYERED GLASS
TYPE:1--[dan','pian],2--[jiajiao],3--[ganghua]'
multipro,'end'
*if,hole,eq,0,then
&&&&multipro,'start',8
&&&&*cset,1,3,doftop,'constrain
type on the top side',2
&&&&*cset,4,6,dofdown,'constrain
type on the down side',1
&&&&*cset,7,9,dofleft,'constrain
type on the left side',2
&&&&*cset,10,12,dofright,'constrain
type on the right side',1
&&&&*cset,13,15,Wdiv,'divide
number on the width side',40
&&&&*cset,16,18,Hdiv,'divide
number on the heigth side',40
&&&&*cset,19,21,gravity,'The
gravity Acceleration(unit:m/s2)',0
&&&&*cset,22,24,switch,'Choose
solution type',1
&&&&*cset,61,62,'Please
Fill in related blanks to',' Get FEA solution'
&&&&*cset,63,64,'SWITCH:0--[By
Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
&&&&multipro,'end'
*elseif,hole,eq,2,then
&&&&multipro,'start',9
&&&&*cset,1,3,kr,'Radium
of glass hole(unit:m)',0.0175
&&&&*cset,4,6,kbw,'Distance
to width side(unit:m)',0.106
&&&&*cset,7,9,kbh,'Distance
to height side(unit:m)',0.106
&&&&*cset,10,12,dofdown,'constrain
type on the down side',1
&&&&*cset,13,15,Wdiv,'divide
number on the width side',40
&&&&*cset,16,18,Hdiv,'divide
number on the heigth side',40
&&&&*cset,19,21,holediv,'divide
number on the hole side',32
&&&&*cset,22,24,gravity,'The
gravity Acceleration(unit:m/s2)',0
&&&&*cset,25,27,switch,'Choose
solution type',1
&&&&*cset,61,62,'Please
Fill in related blanks to',' Get FEA solution'
&&&&*cset,63,64,'SWITCH:0--[By
Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
&&&&multipro,'end'
*elseif,hole,eq,4,then
&&&&multipro,'start',8
&&&&*cset,1,3,kr,'Radium
of glass hole(unit:m)',0.0175
&&&&*cset,4,6,kbw,'Distance
to width side(unit:m)',0.106
&&&&*cset,7,9,kbh,'Distance
to height side(unit:m)',0.106
&&&&*cset,10,12,Wdiv,'divide
number on the width side',40
&&&&*cset,13,15,Hdiv,'divide
number on the heigth side',40
&&&&*cset,16,18,holediv,'divide
number on the hole side',32
&&&&*cset,19,21,gravity,'The
gravity Acceleration(unit:m/s2)',0
&&&&*cset,22,24,switch,'Choose
solution type',1
&&&&*cset,61,62,'Please
Fill in related blanks to',' Get FEA solution'
&&&&*cset,63,64,'SWITCH:0--[By
Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
&&&&multipro,'end'
*elseif,hole,eq,6,then
&&&&multipro,'start',9
&&&&*cset,1,3,kr,'Radium
of glass hole(unit:m)',0.0175
&&&&*cset,4,6,kbw,'Distance
to width side(unit:m)',0.106
&&&&*cset,7,9,kbh,'Distance
to height side(unit:m)',0.106
&&&&*cset,10,12,Dbottom,'hole
central distance to bottom',1
&&&&*cset,13,15,Wdiv,'divide
number on the width side',40
&&&&*cset,16,18,Hdiv,'divide
number on the heigth side',40
&&&&*cset,19,21,holediv,'divide
number on the hole side',32
&&&&*cset,22,24,gravity,'The
gravity Acceleration(unit:m/s2)',0
&&&&*cset,25,27,switch,'Choose
solution type',1
&&&&*cset,61,62,'Please
Fill in related blanks to',' Get FEA solution'
&&&&*cset,63,64,'SWITCH:0--[By
Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
&&&&multipro,'end'
*if,hole,eq,0,then!!!!!!!!!!平板建模
&&&k,1,1,1
&&&k,2,1+w,1
&&&*afun,deg
&&&k,3,1+H*cos(a1)/sin(a1),1+H
&&&k,4,1+w-H*cos(a2)/sin(a2),1+H
&&&a,1,2,4,3
&&&et,1,shell63
*if,type,eq,1,then
&&thforcal=1.1*thi
*elseif,type,eq,2,then
&&thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
&&thforcal=1.1*1.20*thi
&&R,1,thforcal,
&&MP,EX,1,7.2e10
&&MP,PRXY,1,0.21
&&mp,dens,1,2560
&&lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
*if,doftop,eq,0,then
&&&&*elseif,doftop,eq,1,then
&&&&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uX
&&&&&allsel
&&&&*elseif,doftop,eq,2,then
&&&&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uy
&&&&&allsel
*if,dofdown,eq,0,then
&&&&&*elseif,dofdown,eq,1,then
&&&&&lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uX
&&&&&allsel
&&&&&*elseif,dofdown,eq,2,then
&&&&&lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uy
&&&&&allsel
*if,dofleft,eq,0,then
&&&&*elseif,dofleft,eq,1,then
&&&&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uY
&&&&&allsel
&&&&*elseif,dofleft,eq,2,then
&&&&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,ux
&&&&&allsel
*if,dofright,eq,0,then
&&&&&*elseif,dofright,eq,1,then
&&&&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uY
&&&&&allsel
&&&&&*elseif,dofright,eq,2,then
&&&&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,ux
&&&&&allsel
*elseif,hole,eq,2,then!!!!!!!!!!!!两点建模
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
&&thforcal=1.1*thi
*elseif,type,eq,2,then
&&thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
&&thforcal=1.1*1.20*thi
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
&&et,1,shell63
! :do loop1
*do,i,2,kld,1
e,1,kld+1,2
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,6,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
CYL4,kx(5),ky(5),kr
CYL4,kx ,ky ,kr
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
LSEL,S,LENGTH,,2*kr,2*(h+w)
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx ,ky
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
nsel,s,loc,x,kx -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +0.001
d,all,uy,0,,,,uz
CPINTF,ALL,0.0001,
*if,dofdown,eq,0,then
&&&&&*elseif,dofdown,eq,1,then
&&&&&lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&dl,all,,uz
&&&&&allsel
&&&&&*elseif,dofdown,eq,2,then
&&&&&lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&dl,all,,uz
&&&&&dl,all,,uy
&&&&&allsel
*elseif,hole,eq,4,then!!!!!!!!!!!!!!!!!四点建模
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
&&thforcal=1.1*thi
*elseif,type,eq,2,then
&&thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
&&thforcal=1.1*1.20*thi
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
&&et,2,shell63
! :do loop1
*do,i,2,kld,1
e,1,kld+1,2
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1)
k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2)
k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
CYL4,kx(5),ky(5),kr
CYL4,kx ,ky ,kr
CYL4,kx(7),ky(7),kr
CYL4,kx ,ky ,kr
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
LSEL,S,LENGTH,,2*kr,l+w
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx ,ky
*get,nnn3,node,0,count
egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7)
*get,nnn4,node,0,count
egen,2,nnn4+kld,1,kld,1,,,,,,kx ,ky
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
nsel,s,loc,x,kx -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +0.001
d,all,uy,0,,,,uz
nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001
nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001
d,all,ux,0,,,,uz
nsel,s,loc,x,kx -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +0.001
d,all,uz,0,,,,
CPINTF,ALL,0.0001,
*elseif,hole,eq,6,then!!!!!!!!!!!!!!!!!六点建模
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
&&thforcal=1.1*thi
*elseif,type,eq,2,then
&&thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
&&thforcal=1.1*1.20*thi
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
&&et,2,shell63
! :do loop1
*do,i,2,kld,1
e,1,kld+1,2
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1)
k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2)
k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
hp=ky(7)-ky(5)
xx9=(hp-dbottom)*(kx(7)-kx(5))/hp
yy9=ky(5)+dbottom
xx10=(hp-dbottom)*(kx -kx )/hp
yy10=ky +dbottom
k,9,kx(7)-xx9,yy9
k,10,kx +xx10,yy10
CYL4,kx(5),ky(5),kr
CYL4,kx ,ky ,kr
CYL4,kx(7),ky(7),kr
CYL4,kx ,ky ,kr
CYL4,kx(9),ky(9),kr
CYL4,kx(10),ky(10),kr
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
&&&&&lesize,all,,,wdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
&&lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
&&&&&lesize,all,,,hdiv
&&&&&lsel,all
LSEL,S,LENGTH,,2*kr,l+w
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx ,ky
*get,nnn3,node,0,count
egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7)
*get,nnn4,node,0,count
egen,2,nnn4+kld,1,kld,1,,,,,,kx ,ky
*get,nnn5,node,0,count
egen,2,nnn5+kld,1,kld,1,,,,,,kx(9),ky(9)
*get,nnn6,node,0,count
egen,2,nnn6+kld,1,kld,1,,,,,,kx(10),ky(10)
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
nsel,s,loc,x,kx -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +0.001
d,all,uy,0,,,,uz
nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001
nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001
d,all,ux,0,,,,uz
nsel,s,loc,x,kx -0.001,kx +0.001
nsel,r,loc,y,ky -0.001,ky +0.001
d,all,uz,0,,,,
nsel,s,loc,x,kx(9)-0.001,kx(9)+0.001
nsel,r,loc,y,ky(9)-0.001,ky(9)+0.001
d,all,ux,0,,,,uz
nsel,s,loc,x,kx(10)-0.001,kx(10)+0.001
nsel,r,loc,y,ky(10)-0.001,ky(10)+0.001
d,all,uz,0,,,,
CPINTF,ALL,0.0001,
/solu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解部分
Ek=3*seismic*(thi+tho)*26.0*
q1=0.6*ek+wk
q2=w+0.6*e
/UIS,MSGPOP,3
esel,s,real,,1,1,1
sfe,all,1,pres,,q1
acel,,gravity,
*if,switch,eq,1,then
nlgeom,off
*elseif,switch,eq,2,then
esel,s,real,,1,1,1
plnsol,u,sum,0,1
*get,numall,node,0,count
*get,numstart,node,0,num,min
&&&inuse=numstart
*do,i,1,numall,1
&&&*get,usum,node,inuse,u,sum
&&&*if,umax,ge,usum,then
&&&&&&&&&&&&&&&umax=umax
&&&*elseif,umax,lt,usum,then
&&&&&&&&&&&&&&&umax=usum
&&&&*endif
inuse=ndnext(inuse)
/UIS,MSGPOP,3
esel,s,real,,1,1,1
sfe,all,1,pres,,q2
*if,switch,eq,1,then
nlgeom,off
esel,s,real,,1,1,1
plnsol,s,eqv,0,1
esel,s,real,,1,1,1
*get,numall,node,0,count
*get,numstart,node,0,num,min
&&&inuse=numstart
*do,i,1,numall,1
&&&*get,ssum,node,inuse,s,eqv
&&&*if,smax,ge,ssum,then
&&&&&&&&&&&&&&&smax=smax
&&&*elseif,smax,lt,ssum,then
&&&&&&&&&&&&&&&smax=ssum
&&&&*endif
inuse=ndnext(inuse)
*if,thforcal,le,0.012,then
&&&&&&&dmq=84
&&&&&&&kbq=58.8
*elseif,thforcal,gt,0.012,then
&&&&&&&dmq=59
&&&&&&&kbq=41.3
*if,h,lt,w,then
*elseif,h,ge,w,then
str=smax/1e6
esel,s,real,,1,1,1
/UIS,MSGPOP,1
*msg,note,str,umax*1000,dmq,10*duanb
Ultimate Limit States Glass STRESS=%gMa,%/&
Serviceabilitylimits DISPLACEMENT=%gmm,%/&
[Toughened Glass Allowble Stress]=%gMa,%/&
[Control Displacement]=%gmm,
/UIS,MSGPOP,3
*elseif,switch,eq,2,then
esel,s,real,,1,1,1
plnsol,s,eqv,0,1
esel,s,real,,1,1,1
*get,numall,node,0,count
*get,numstart,node,0,num,min
&&&inuse=numstart
*do,i,1,numall,1
&&&*get,ssum,node,inuse,s,eqv
&&&*if,smax,ge,ssum,then
&&&&&&&&&&&&&&&smax=smax
&&&*elseif,smax,lt,ssum,then
&&&&&&&&&&&&&&&smax=ssum
&&&&*endif
inuse=ndnext(inuse)
*if,thforcal,le,0.012,then
&&&&&&&dmq=84
&&&&&&&kbq=58.8
*elseif,thforcal,gt,0.012,then
&&&&&&&dmq=59
&&&&&&&kbq=41.3
*if,h,lt,w,then
*elseif,h,ge,w,then
str=smax/1e6
esel,s,real,,1,1,1
/UIS,MSGPOP,1
*msg,note,str,umax*1000,dmq,10*duanb
Ultimate Limit States Glass STRESS=%gMa,%/&
Serviceabilitylimits DISPLACEMENT=%gmm,%/&
[Toughened Glass Allowble Stress]=%gMa,%/&
[Control Displacement]=%gmm,
/UIS,MSGPOP,3
*elseif,swtich,eq,0
.等截面杆单元生死应用实例
等截面杆两端固定，承受均匀的温度载荷时将其中间1/3段移去。过程是将其应变自由化并移去均匀温度。分析其热应力和应变情况。
&&/TITLE, ELEMENT BIRTH/DEATH IN
A FIXED BAR WITH THERMAL LOADING
&&ET,1,LINK1 !二维杆单元
&&MP,EX,1,30E6 !材料特性
&&MP,ALPX,1,.00005
&&MP,EX,2,30E6
&&MP,ALPX,2,.00005
!重新‘出生’单元的特性
&&MP,REFT,2,100 !单元出生的参考温度
&&EGEN,3,1,-1 !生成三个单元
&&ANTYPE,STATIC
&&D,1,ALL,,,4,3 !固定杆的两端
&&TREF,0 !参考温度0度
&&TUNIF,100 !均匀温度载荷
&&NROPT,FULL
&&OUTPR,BASIC,ALL
&&EKILL,2 !‘杀死’中间的单元
&&EALIVE,2 !重新激活中间单元
&&MPCHG,2,2
!将材料特性改为2以生成自由应变单元
&&TUNIF,0 !删除温度载荷
.梁板建模联系
/input,beamandshell.inp
/title,beamandshell
et,1,shell63
et,2,beam188
MP,EX,1,3.0e10
MP,PRXY,1,0.16667
MP,DENS,1,2.551e3
rectng,,4,,4
sectype,1,beam,rect,fucksec
secdata,1,0.5,10,5
secoffset,user,0,0.5
wprota,,90
wprota,,,90
aatt,1,1,1
lsel,s,,,4,12,8
lsel,a,,,10,14,4
k,100,2,2,-1
latt,1,2,2,,100,,1
lesize,all,0.1
/eshape,1!查看梁板的外形
lesize,all,0.1
/eshape,1!查看梁板的外形
allsel,all
lsel,s,,,5,6
lsel,a,,,8,9
lsel,a,,,2,11,9
lsel,a,,,13,15,2
dl,all,,all
allsel,all
plnsol,s,eqv,0,1
35.圆盘转动例子
RADIUS_INNER = 10
RADIUS_OUTER = 30
LOAD_OMEGA = 100
/triad,lbot
/pnum,type,1
mp,ex ,1,10e6
mp,nuxy,1,0.3
mp,dens,1,0.1/386.1
keyopt,2,1,1
pcirc,RADIUS_OUTER,RADIUS_INNER,0,90
pcirc,RADIUS_OUTER,RADIUS_INNER,90,180
pcirc,RADIUS_OUTER,RADIUS_INNER,180,270
pcirc,RADIUS_OUTER,RADIUS_INNER,270,360
aatt,1,1,1
! create rigid beams inside
lsel,s,loc,x,RADIUS_INNER
type,2$real,2$mat,2
*get,NCOUNT,node,,count
NEXTNODE=0
*do,ICOUNT,1,NCOUNT-1
&&NEXTNODE=ndnext(NEXTNODE)
&&e,node(0,0,0),NEXTNODE
! create dummy element
n,,RADIUS_OUTER+5,0,0
e,node(RADIUS_OUTER,0,0),node(RADIUS_OUTER+5,0,0)
d,node(0,0,0),all
ddele,node(0,0,0),rotz
allsel,all
antype,trans
rescon,define,none
outres,all,all
deltim,6e-4,1e-5,0.1
nsel,u,node,,node(0,0,0)
*get,NCOUNT,node,,count
NEXTNODE=0
*do,ICOUNT,1,NCOUNT
&&NEXTNODE=ndnext(NEXTNODE)
&&ic,NEXTNODE,uy,,LOAD_OMEGA*nx(NEXTNODE)
allsel,all
plns,u,sum
ANDATA,0.1,,2,0,0,1,0,1
简单的例子－如何对结构的振动控制分析
第一步：首先建立模型，如图，一个很简单的角形结构梁
其中上面的号码是节点号，连接1，10的就是调频器，节点2与1均有质量。
/clear,start
et,1,beam44
R,1,78,680,,
RMODIF,1,21,,
RMODIF,1,25,,0,0,7.,
UIMP,1,EX, , ,210000,
UIMP,1,DENS, , ,1.1261e-5,
et,2,mass21
第二步：进行模态分析
antype,modal
modopt,subsp,10
自己进入后处理器，察看结果
&&&&&&&&&&&&&&&&&&TIME
再进入扩展模态，找出对应想抑制震动的主自由度
在上面就以二阶和三阶为例吧，他们的振型是在y，z平面，这也好理解，这个方向刚度很低，所以低频率就出现在这个振型上。所以上面加的tmd的方向也很讲究，也就是节点10的位置是垂直于y，z平面，不然，在其他方向就没有好的效果了。
et,3,combin14
r,3,100,10
r,4,100,100,100
n,10,0,0,-1000
然后进入阻尼模态分析，由于关心的频率范围是前两阶，因此，设置了频率范围在0.1到0.2，
在此强调一下，在damp中频率范围一定要设置，不然感兴趣的结果有很多遗漏。
antype,modal
modopt,damp,10
MXPAND,10, , ,0
MODOPT,DAMP,10,0.1,0.2, ,OFF
然后进入后处理器
有兴趣者可以对tmd命令流中的r,3,100,10中的阻尼进行变化
就可发现阻尼的作用了
对于低频作用不大，但对高频有很好的抑制效果
弹簧刚度的作用相反。
这个简单的例子是对其方法的一个简单描述，可能有不少错误和不详细的地方，希望大家共同探讨指正，谢谢
.模态分析结果的输出实例
在gui里看到的既不能修改，也不能copy
现用命令流导出
考虑斑竹以积分鼓励大家，以一分为限
例子如下：
/clear,start
mass=83500
et,1,beam44
et,2,mass21
R,1,78,680,,
RMODIF,1,21,,
RMODIF,1,25,,0,0,7.,
R,2,mass-j,mass-j,mass-j,0,0,0,
UIMP,1,EX, , ,210000,
UIMP,1,DENS, , ,1.1261e-5,
UIMP,2,EX, , ,210000,
UIMP,2,DENS, , ,0.78e-5,
latt,1,1,1
lesize,1,,,5,1
katt,2,2,2
d,1,,,,,,ux,uy,uz,rotx,roty,rotz
antype,modal
modopt,subsp,5
FILE !确定数据源
FILE,'file','rst','.'
/UI,COLL,1 !打开对话框
NUMVAR,200 !post26中最大允许存100个
SOLU,191,NCMIT !确定每时间解的存储
STORE,MERGE ! 存数据
FILLDATA,191,,,,1,1 !以斜坡函数来填充数据
REALVAR,191,191 !只使用实部
! Save time history variables to file b.csv
*DEL,_P26_EXPORT !删除参数
*DIM,_P26_EXPORT,TABLE,5,1 ! 定义数组大小
VGET,_P26_EXPORT(1,0),1 !移动数据
*if,i,eq,1,then
/OUTPUT,'a1','csv','.' ! 生成文件
*elseif,i,eq,2,then
/OUTPUT,'a2','csv','.'
*elseif,i,eq,3,then
/OUTPUT,'a3','csv','.'
*elseif,i,eq,4,then
/OUTPUT,'a4','csv','.'
*elseif,i,eq,5,then
/OUTPUT,'a5','csv','.'
*VWRITE,'TIME' ! 按一定次序写入数据
*VWRITE,_P26_EXPORT(1,0)
/OUTPUT,TERM
! End of time history save
.火车过桥动态加载实例（部分）
! 用对话框提示用户输入基本数据
!multipro,'start',3
! *cset,1,3,y_bot,'enter the value of y_bot',0
! *cset,4,6,y_top,'enter the value of y_bot',0
! *cset,7,9,ttype,'enter the value of train type',0
!输入火车活载的类型
!multipro,'end'
y_bot=13.9
y_top=14.1
*if,ttype,eq,0, then
&&&q1 = 147000
&&&L2 = 30
*elseif,ttype,eq,1,then
&&&*msg,ui
&&&error value
of ttype(1 or 2),please enter again
! 识别桥面过列车的单元，并将单元按节点坐标从小到大排序
lsel,s,loc,y,y_bot,y_top
esll,s ! 选中桥面单元
elmax = elmiqr(0,14) ! 最大单元号
ndmax = ndinqr(0,14) ! 最大节点号
*dim,eldk,,elmax ! 存放桥面过车单元
*dim,ndx,,elmax ! 存放过车单元的节点x坐标（两节点坐标中的较小值）
*dim,nodes,,ndmax ! 存放桥面过车的节点
*do,i,1,elmax
&&*if,esel ,gt,0, then
&&&&eldk(j)
*enddo ! 记录桥面单元的标号，保存于数组eldk中
eldknum = j ! 总桥面单元数
! 对eldk数组排序，按节点坐标从小到大排列
*do,i,1,eldknum
nelem(ii,1)
nelem(ii,2)
&&&*if,i2,lt,i1,then
*do,i,1,eldknum-1
&&*do,j,i+1,eldknum
&&&&*if,ndx(j),lt,ndx
&&&&&&temp
&&&&&&ndx(j)
&&&&&&temp
&&&&&&eldk
&&&&&&eldk(j)=
&&&&*endif
nsle,s ! 选中桥面单元的节点
*do,i,1,ndmax
&&*if,nsel ,gt,0,then
&&&&&j=j+1
&&&&&nodes(j)=i
ndnum = j ! 总桥面节点数
xmin = nx(nelem(eldk(1),1))
temp = nx(nelem(eldk(1),2))
*if,temp,lt,xmin,then
&&xmin = temp
xmax = nx(nelem(eldk(eldknum),1))
temp = nx(nelem(eldk(eldknum),2))
*if,temp,gt,xmax,then
&&xmax = temp
qth = xmax -xmin ! 桥面过车长度
lth1 = 2*qth + L1+L2 +182.5
lth = (nint(lth1/5)+8)*5 ! 总长度
dim = lth/5 +1
*dim,disy,table,dim,ndnum ! 桥面节点的y方向的位移值
*dim,disy1,table,dim,ndnum ! 桥面节点的y方向的位移值(从大到小)
*dim,disx,table,dim,ndnum ! 桥面节点的x方向的位移值
*dim,disa,table,dim,ndnum ! 桥面节点的转角
*dim,disymax,,ndnum
*dim,diy,table,dim
*dim,ffx,table,elmax,dim ! 单元轴力
*dim,fmz,table,elmax,dim ! 单元弯矩
*dim,detaf,,elmax ! 每个单元轴力变化值
*dim,detaf1,,elmax ! 每个单元轴力变化值（从小到大）
*dim,elnum,,elmax ! 轴力变化最大的单元号
*dim,detaf2,,20 ! 前20个轴力变化最大值
*dim,elnum2,,20 ! 前20个轴力变化最大值的单元
*dim,fmax,,elmax ! 每个单元的轴力最大值
*dim,fmin,,elmax ! 每个单元的轴力最小值
*dim,fabs,,elmax ! 每个单元轴力最大值（绝对值）
*dim,elnum3,,elmax ! 轴力最大的单元号
*dim,elnum4,,20 ! 前20个轴力最大的单元号
*dim,nmax,,20
*dim,dfmaxel,table,dim,20 ! 前20各轴力变化最大单元的轴力
*dim,dmmaxel,table,dim,20 ! 前20各轴力变化最大单元的弯矩
*dim,fmaxel,table,dim,20 ! 前20各轴力最大单元的轴力
*dim,mmaxel,table,dim,20 ! 前20各轴力最大单元的弯矩
eplot ! plot element
*do,kk,0,lth,5,
time = time+1
!/uis,msgpop,3 ! 不显示数组被重新定义的警告
!/uis,msgpop,2
/solu ! 进入求解器
antype,static
sfedele,all,1,pres, !先删除上一次加上的所有单元荷载
! 计算列车各部分在桥上的单元号
tpos1 = kk ! 火车头所在的位置
*if,tpos1,lt,(qth+L1+L2), then !火车进桥情况
&&&xlth1 =
&&&*do,i,1,eldknum
= nelem(ii,1)
= nelem(ii,2)
&&&&&elth1
= abs(nx(ndi)-nx(ndj))
&&&&&xlth1
= xlth1 + elth1
&&&&&*if,xlth1,ge,tpos1,then
&&&&&&&&*exit
&&&&&*endif
&&&tele1 = i !
火车头所在的单元
&&&*if,tpos1,gt,L1,then
&&&&&&tpos2
=tpos1 -L1 !煤水车头所在的位置
&&&&&&xlth2
&&&&&&*do,i,1,eldknum
&&&&&&&&ii
&&&&&&&&ndi
= nelem(ii,1)
&&&&&&&&ndj
= nelem(ii,2)
&&&&&&&&elth2
= abs(nx(ndi)-nx(ndj))
&&&&&&&&xlth2
= xlth2 + elth2
&&&&&&&&*if,xlth2,ge,tpos2,then
&&&&&&&&&&*exit
&&&&&&&&*endif
&&&&&&*enddo
&&&&&&tele2
= i !旅客列车头所在的单元
&&&&&&*if,tpos1,gt,(L1+L2),then
&&&&&&&&tpos3
=tpos1 -(L1+L2) !旅客车头所在的位置
&&&&&&&&xlth3
&&&&&&&&*do,i,1,eldknum
&&&&&&&&&&ii
&&&&&&&&&&ndi
= nelem(ii,1)
&&&&&&&&&&ndj
= nelem(ii,2)
&&&&&&&&&&elth3
= abs(nx(ndi)-nx(ndj))
&&&&&&&&&&xlth3
= xlth3 + elth3
&&&&&&&&&&*if,xlth3,ge,tpos3,then
&&&&&&&&&&&&&*exit
&&&&&&&&&&*endif
&&&&&&&&*enddo
&&&&&&&&tele3
= i !旅客列车头所在的单元
&&&&&&*endif
*if,tpos1,gt,(qth+L1+L2+182.5),then ! 火车过桥情况(182.5为任意取的一个常数)
&&*if,tpos1,lt,lth1, then
= tpos1-(qth+L1+L2+182.5) !旅客列车尾部所在的位置
&&&&*do,i,1,eldknum
= nelem(ii,1)
= nelem(ii,2)
&&&&&&elth4
= abs(nx(ndi)-nx(ndj))
&&&&&&xlth4
= xlth4 + elth4
&&&&&&*if,xlth4,ge,tpos4,then
&&&&&&&&*exit
&&&&&&*endif
&&&&*enddo
= i !旅客列车尾部在的单元
! 加荷载于相应的单元
! SFBEAM, ELEM, LKEY, Lab, VALI, VALJ, VAL2I, VAL2J, IOFFST,
*if,tpos1,le,L1, then ! 火车头进桥
&&c1 = xlth1 -tpos1
&&c2 = elth1 - c1
&&*if,tele1,eq,1,then
= q1*c2/elth1
&&&&&sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
&&&&&*do,j,1,tele1-1
&&&&&&&&sfbeam,eldk(j),1,pres,q1,q1,,,,,
&&&&&*enddo
= q1*c2/elth1
&&&&&sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*elseif,tpos1,le,(L1+L2),then ! 煤水车进桥
&&c1 = xlth1 - tpos1
&&*if,tele1,eq,tele2,then
!火车头在同一个单元上
= xlth2 - tpos2
= elth2 - c2
= (q1*L1+q2*c3)/elth1
&&&&sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
&&&&*do,j,1,tele1-1
&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&*enddo
= tele1 - tele2
&&&&*if,dtele,eq,1,
= xlth2 - tpos2
= elth2 - c2
&&&&&&&qq0
= q1*(elth1-c1)/elth1
= (q1*c2+q2*c3)/elth2
&&&&&&&sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
&&&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
&&&&&&&*do,j,1,tele2-1
&&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&&&*enddo
= xlth2 - tpos2
= elth2 - c2
&&&&&&&qq0
= q1*(elth1-c1)/elth1
&&&&&&&sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
&&&&&&&*do,j,tele2+1,tele1-1
&&&&&&&&&&sfbeam,eldk(j),1,pres,q1,q1,,,,,
&&&&&&&*enddo
= (q1*c2+q2*c3)/elth2
&&&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
&&&&&&&*do,j,1,tele2-1
&&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&&&*enddo
&&&&&*endif
*elseif,tpos1,le,qth,then ! 旅客列车进桥
&&c1 = xlth1 - tpos1
&&*if,tele1,eq,tele2, then
= xlth2 - tpos2
= elth2 - c2
= (q1*L1+q2*c3)/elth1
&&&&sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
= xlth3 - tpos3
= elth3 - c4
&&&&dtele1
= tele2 - tele3
&&&&*if,dtele1,eq,1,then
&&&&&&&&qq1
= (q2*c4+q3*c5)/elth3
&&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&&*enddo
&&&&&*else
&&&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&&&*enddo
&&&&&&&qq1
= (q2*c4+q3*c5)/(elth3)
&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&*enddo
&&&&&*endif
&&&&dtele1
= tele1 -tele2
&&&&*if,dtele1,eq,1,then
= xlth2 - tpos2
= elth2 -c2
&&&&&&&qq0
= q1*(elth1-c1)/elth1
&&&&&&&sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
= xlth3- tpos3
= elth3 - c4
&&&&&&&dtele2
= tele2 - tele3
&&&&&&&*if,dtele2,eq,0,then
&&&&&&&&&&qq0
= (q1*c2+q2*L2+q3*c3)/elth2
&&&&&&&&&&sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
&&&&&&&&&&*do,j,1,tele2-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&&&&*enddo
&&&&&&&*elseif,dtele2,eq,1,then
&&&&&&&&&&qq0
= (q1*c2+q2*c3)/elth2
&&&&&&&&&&sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
&&&&&&&&&&qq1
= (q2*c4+c3*c5)/elth3
&&&&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&*else
&&&&&&&&&&qq0
= (q1*c2+q2*c3)/elth2
&&&&&&&&&&sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
&&&&&&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&&&&&&*enddo
&&&&&&&&&&qq1
= (q2*c4+q3*c5)/(elth3)
&&&&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&&&&*enddo
&&&&&&&*endif
= xlth2 - tele2
= elth2 - c2
&&&&&&&qq0
= q1*(elth1-c1)/elth1
&&&&&&&sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
&&&&&&&*do,j,tele2+1,tele1-1,
&&&&&&&&&&sfbeam,eldk(j),1,pres,q1,q1,,,,,
&&&&&&&*enddo
= (q1*c2+q2*c3)/(elth2)
&&&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
&&&&&&&dtele2
= tele2 - tele3
&&&&&&&*if,dtele2,eq,1,then
&&&&&&&&&&qq1
= (q2*c4+q3*c5)/(elth3)
&&&&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&&&&*enddo
&&&&&&&*else
&&&&&&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&&&&&&*enddo
&&&&&&&&&&qq1
= (q2*c4+q3*c5)/(elth3)
&&&&&&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&&&&&&*do,j,1,tele3-1
&&&&&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&&&&&&*enddo
&&&&&&&*endif
&&&&*endif
*elseif,tpos1,le,(qth+L1),then !火车头出桥
&&&c1 = tpos1 -
qth !火车头出桥长度
&&&dtele1 =
eldknum - tele2
&&&*if,dtele1,eq,0,then
= xlth2 - tpos2
= elth2 - c2
= (q1*c2+q2*c3)/(elth2)
&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
&&&&&dtele2
= tele2 - tele3
= xlth3 -tpos3
= elth3 - c4
&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&*enddo
= (q2*c4+q3*c5)/elth3
&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&*do,j,1,tele3-1
&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&*enddo
&&&*elseif,dtele1,eq,1,then
&&&&&sfbeam,eldk(eldknum),1,pres,q1,q1,,,,,
= (q1*c2+q2*c3)/(elth2)
&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
&&&&&dtele2
= tele2- tele3
= xlth3 - tpos3
= elth3 - c4
&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&*enddo
= (q2*c4+q2*c5)/elth3
&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&*do,j,1,tele3-1
&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&*enddo
&&&*elseif,dtele1,gt,1,then
&&&&&*do,j,tele2+1,eldknum
&&&&&&&&sfbeam,eldk(j),1,pres,q1,q1,,,,,
&&&&&*enddo
= (q1*c2+q2*c3)/(elth2)
&&&&&sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
= xlth3 - tpos3
= elth3 - c4
&&&&&*do,j,tele3+1,tele2-1
&&&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&*enddo
= (q2*c4+q2*c5)/elth3
&&&&&sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
&&&&&*do,j,1,tele3-1
&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&*enddo
*elseif,tpos1,le,(qth+L1+L2),then !煤水车头出桥
&&&c4 = xlth3 -
&&&c5 = elth3 -
&&&dtele1 =
eldknum - tele3
&&&*if,dtele1,eq,0,then
= (q2*c4+q3*c5)/elth3
&&&&&sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
&&&&&*do,j,1,tele3-1
&&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&*enddo
&&&&&*do,j,tele3+1,eldknum
&&&&&&&sfbeam,eldk(j),1,pres,q2,q2,,,,,
&&&&&*enddo
= (q2*c4+q3*c5)/elth3
&&&&&sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
&&&&&*do,j,1,tele3-1
&&&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
&&&&&*enddo
*elseif,tpos1,le,(qth+L1+L2+182.5),then ! 全桥均布列车荷载(182.5为任意取的
&&&*do,j,1,eldknum
&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
*elseif,tpos1,le,lth1,then ! 均布列车荷载出桥
&&&c1 = xlth4 -
&&&c2 = elth4 -
&&&*do,j,tele4+1,eldknum
&&&&&&sfbeam,eldk(j),1,pres,q3,q3,,,,,
q3*c1/elth4
&&&sfbeam,eldk(tele4),1,pres,qq,qq,,,,,
/psf,pres,norm,1
diy(time,0) = tpos1
diy(time) = uy(91) ! 得到中间节点的y坐标
*vplot,diy(1,0),diy(1,1) ! 动态显示中间节点位移结果
disx(time,0) = tpos1
disy(time,0) = tpos1
disa(time,0) = tpos1
*do,j,1,ndnum
&&disx(0,j) = nodes(j)
&&disy(0,j) = nodes(j)
&&disa(0,j) = nodes(j)
&&disx(time,j) =
ux(nodes(j))
&&disy(time,j) =
uy(nodes(j))
&&disa(time,j) =
rotz(nodes(j))
etable,forcex,smisc,1 ! 定义单元轴力表
etable,momentz,smisc,6 ! 定义单元弯矩表
*vget,ffx(1,time),elem,1,etab,forcex,,,,
*vget,fmz(1,time),elem,1,etab,momentz,,,,
ffx(0,time) = tpos1
fmz(0,time) = tpos1
*if,time,eq,1,then
&&*do,j,1,elmax
&&&&&ffx(j,0)
&&&&&fmz(j,0)
save,,,,all
! 结果处理
*do,i,1,ndnum
&&nd = nodes
&&disymax = abs(disy(1,nd))
&&*do,j,2,dim
&&&&*if,abs(disy(j,nd)),gt,disymax
&&&&&&disymax
= abs(disy(j,nd))
&&&&*endif
*do,i,1,ndnum-1
&&&*do,j,i+1,ndnum
&&&&&*if,disymax(j),gt,disymax
&&&&&&&&temp
&&&&&&&&disymax
= disymax(j)
&&&&&&&&disymax(j)
&&&&&&&&temp
&&&&&&&&nodes
= nodes(j)
&&&&&&&&nodes(j)
&&&&&*endif
*do,i,1,ndnum ! 前5个位移最大的节点
&&&numnd =
&&&disy1(0,i) =
&&&*do,j,1,dim
&&&&&&disy1(j,0)
&&&&&&disy1(j,i)
= disy(kk,numnd)
*do,i,1,elmax
&&&*do,j,2,dim
&&&&&*if,ffx(i,j),gt,fmax
&&&&&&&&fmax
= ffx(i,j)
&&&&&*endif
&&&&&*if,ffx(i,j),lt,fmin
&&&&&&&&fmin
= ffx(i,j)
&&&&&*endif
detaf = fmax - fmin ! 每个单元的内力变化值
detaf1 = detaf
*do,i,1,elmax
&&&elnum =
*do,i,1,elmax-1
&&&*do,j,i+1,elmax
&&&&&*if,detaf1(j),gt,detaf1
&&&&&&&temp
&&&&&&&detaf1
= detaf1(j)
&&&&&&&detaf1(j)
&&&&&&&temp
&&&&&&&elnum
= elnum(j)
&&&&&&&elnum(j)=temp
&&&&&*endif
*do,i,1,20
&&detaf2 = detaf1
!前20个内力变化最大的单元
&&elnum2 = elnum
*do,i,1,20 ! 前20个内力变化最大的单元
&&&numel =
&&&dfmaxel(0,i)
&&&dmmaxel(0,i)
&&&*do,j,1,dim
&&&&&&dfmaxel(j,0)
&&&&&&dfmaxel(j,i)
= ffx(numel,kk) ! 轴力变化最大单元的轴力
&&&&&&dmmaxel(j,0)
&&&&&&dmmaxel(j,i)
= fmz(numel,kk) ! 相应的弯矩
! 前20个内力最大的单元
*do,i,1,elmax
=abs(fmax )
&&&*if,fabs
,lt,abs(fmin ),then
&&&&&&&fabs
=abs(fmin )
*do,i,1,elmax
*do,i,1,elmax-1
&&&*do,j,i+1,elmax
&&&&&*if,fabs(j),gt,fabs
&&&&&&&&&temp
&&&&&&&&&fabs
&&&&&&&&&fabs(j)
&&&&&&&&&temp
&&&&&&&&&elnum3
= elnum3(j)
&&&&&&&&&elnum3(j)
&&&&&*endif
&&&&*enddo
*do,i,1,20
&&elnum4 = elnum3
&&nmax = fabs
*do,i,1,20 ! 前20个内力最大的单元
&&&numel =
&&&fmaxel(0,i) =
&&&mmaxel(0,i) =
&&&*do,j,1,dim
&&&&&&fmaxel(j,0)
&&&&&&fmaxel(j,i)
= ffx(numel,kk) ! 轴力最大单元的轴力
&&&&&&mmaxel(j,0)
&&&&&&mmaxel(j,i)
= fmz(numel,kk) ! 相应的弯矩
save,,,,all
!/title,the 91 node displament vs train position ! 显示图形的名称
!*vplot,diy(1,0),diy(1,1)
!/ui,copy,save,bmp,graph
*do,i,1,30
cc = nodes
/title, y displacement of node �% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
/sys,rename njq*.bmp aaa*.bmp
*do,i,31,ndnum
cc = nodes
/title, y displacement of node �% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
/sys,rename njq*.bmp aa*.bmp
!*do,i,1,20
!cc = elnum2
!/title,axis of element �% VS. train position
! *vplot,dfmaxel(1,0),dfmaxel(1,i)
!/ui,copy,save,bmp,graph
!/title,moment of element �% VS. train position
! *vplot,dmmaxel(1,0),dmmaxel(1,i)
!/ui,copy,save,bmp,graph
*do,i,1,20
cc = elnum4
/title,axis of element �% vs train position
*vplot,fmaxel(1,0),fmaxel(1,i)
/ui,copy,save,bmp,graph
/title,moment of element �% vs train position
*vplot,mmaxel(1,0),mmaxel(1,i)
/ui,copy,save,bmp,graph
*cfopen,result,txt,,,
('**********************南京桥计算结果*****************************')
(' ****节点最大位移**** ')
ndnum1 = 5*nint(ndnum/5)
*do,i,1,ndnum1,5
b=nodes(i+1)
c=nodes(i+2)
d=nodes(i+3)
e=nodes(i+4)
*vwrite,a,b,c,d,e
('荷载位置',5(f8.0,5x))
*vwrite,disy1(1,0),disy1(1,i),disy1(1,i+1),disy1(1,i+2),disy1(1,i+3),disy1(
(f5.0,,5(3x,f10.5))
(' ****前20个轴力变化最大的单元和相应的轴力值**** ')
*vwrite,elnum2(1),detaf2(1)
(f15.0,' ', f15.1)
*do,i,1,18,3
b= elnum2(i+1)
c= elnum2(i+2)
*vwrite,a,b,c
('荷载位置',9x,'单元',f4.0,14x,'单元',f4.0,14x,'单元',f4.0)
*vwrite,dfmaxel(1,0),dfmaxel(1,i),dmmaxel(1,i),dfmaxel(1,i+1),dmmaxel(1,i+1)
,dfmaxel(1,i+2),dmmaxel(1,i+2)
(f5.0,3(f11.1,' ',f11.1))
a= elnum2(19)
b= elnum2(20)
c= elnum2(20)
*vwrite,a,b,c
('荷载位置',9x,'单元',f4.0,14x,'单元',f4.0,14x,'单元',f4.0)
*vwrite,dfmaxel(1,0),dfmaxel(1,19),dmmaxel(1,19),dfmaxel(1,20),dmmaxel(1,20)
,dfmaxel(1,20),dmmaxel(1,20)
(f5.0,3(f11.1,' ',f11.1))
(' ***前20个轴力最大的单元和相应的弯矩值*** ')
*vwrite,elnum4(1),detaf2(1)
(f15.0,' ', f15.1)
*do,i,1,18,3
b= elnum4(i+1)
c= elnum4(i+2)
*vwrite,a,b,c
('荷载位置',9x,'单元',f4.0,14x,'单元',f4.0,14x,'单元',f4.0)
*vwrite,fmaxel(1,0),fmaxel(1,i),mmaxel(1,i),fmaxel(1,i+1),mmaxel(1,i+1),fmax
el(1,i+2),mmaxel(1,i+2)
(f5.0,3(f11.1,' ',f11.1))
a= elnum4(19)
b= elnum4(20)
c= elnum4(20)
*vwrite,a,b,c
('荷载位置',9x,'单元',f4.0,14x,'单元',f4.0,14x,'单元',f4.0)
*vwrite,fmaxel(1,0),fmaxel(1,19),mmaxel(1,19),fmaxel(1,20),mmaxel(1,20),fmax
el(1,20),mmaxel(1,20)
(f5.0,3(f11.1,' ',f11.1))
.悬索结构的找形和计算的例题
取自&悬索结构设计&沈世钊，中国建筑工业出版设，1997.8
p149例题3-1
进行了找形计算和荷载作用下的静力计算
找形时取较小的弹性模量，静力计算时采用弹性模量的真值
这个例题说明，ANSYS还是可以进行悬索结构的分析的
ET,1,LINK10
!找形时采用很小的弹性模量
MP,EX,1,0.60551e9
MP,NUXY,1,0.3
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
n,k,(j-i)*9.15,36.6-(i-1)*9.15,0
*do,j,1,2*i
*do,j,1,2*i
k=i*i+j+25
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
*do,i,18,24
!施加位移约束
d,1,uz,3.66
d,2,uz,1.83
d,4,uz,1.83
d,10,uz,-1.83
d,16,uz,-1.83
d,17,uz,-3.66
d,25,uz,-3.66
d,35,uz,-1.83
d,41,uz,-1.83
d,27,uz,1.83
d,29,uz,1.83
d,26,uz,3.66
!设定时间步
NSUBST,20, , ,1
!输出每个子步的结果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!设定大变形，应力刚化
/STATUS,SOLU
PRNSOL,DOF,
PRESOL,SMISC,1
PLESOL,SMISC,1
PLNSOL,U,X,0,1
!更新几何形状
!a.rst为计算结果文件名，最后一个为目录
!这两个参数应根据你的计算情况定
UPGEOM,1,LAST,LAST,a,rst,E:\JZD\1\
!弹性模量恢复为真值
R,1,0..0027248,
MP,EX,1,2.0e11
!重新施加位移约束
/STATUS,SOLU
!观察施加荷载之前结构的位移
!可见位移基本为零,预张力基本不变
PLNSOL,U,Z,0,1
PLESOL,SMISC,1
!写荷载工况文件01
LSWRITE,01,
!施加节点荷载
*do,i,1,41
f,i,fz,-167445
!设置第2荷载步
NSUBST,20, , ,1
LSWRITE,02,
LSSOLVE,1,2,1,
!可以看到，该结果和书中结果一致
NSOL,2,21,U,Z,
PLVAR,2, , , , , , , , , ,
ESOL,3,33,1,F,Y,
PLVAR,3, , , , , , , , , ,
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。