【转载】ANSYS完全法与模态叠加法瞬态分析实例

原文地址：http://www.caetecc.com/thread-2172-1-1.html

! 半脉冲载荷 --- 模态叠加法
fini
/clear,nostart
/PREP7
ET,1,BEAM4             
MP,EX,1,110E9   
MP,NUXY,1,0.3
MP,DENS,1,1600
R,1,0.000045,3.75e-12,8.44e-10,0.015,0.003,0    
k,1,0,0
k,2,0.2,0
LSTR,       1,       2
LATT,1,1,1, , , , 
lsel,all
LESIZE,all, , ,100, , , , ,1
lmesh,all
FINISH
/solu
antype,2
d,1,all
modopt,lanb,50
solve
fini

/solu
plusetime=0.01
pi=3.1415926
delt_time=0.0001
antype,trans
trnopt,msup,50
deltim,delt_time
d,1,all
time,0
f,101,fz,0
kbc,1
autots,off
solve
*DO,j,1,plusetime/delt_time
  TIME,delt_time*j
  KBC,0  
  autots,off
  F,101,FY,-500*sin(delt_time*j*pi*100)
  solve
*enddo
TIME,10
autots,on
F,101,FY,0 
solve

! 脉冲载荷 --- 模态叠加法
fini
/clear,nostart
/PREP7
ET,1,BEAM4             
MP,EX,1,110E9   
MP,NUXY,1,0.3
MP,DENS,1,1600
R,1,0.000045,3.75e-12,8.44e-10,0.015,0.003,0    
k,1,0,0
k,2,0.2,0
LSTR,       1,       2
LATT,1,1,1, , , , 
lsel,all
LESIZE,all, , ,100, , , , ,1
lmesh,all
FINISH
/solu
antype,2
d,1,all
modopt,lanb,100
solve
fini

/solu
plusetime=0.01
pi=3.1415926
delt_time=0.0001
antype,trans
trnopt,msup,100
deltim,delt_time
d,1,all
time,0
f,mid_node,fz,0
kbc,1
autots,off
solve
*DO,j,1,plusetime/delt_time
  TIME,delt_time*j
  KBC,0  
  autots,off
  F,101,FY,500
  solve
*enddo
TIME,1
autots,on
F,101,FY,0 
solve

! full method
fini
/clear,nostart
/PREP7
ET,1,BEAM4             
MP,EX,1,110E9   
MP,NUXY,1,0.3
MP,DENS,1,1600
R,1,0.000045,3.75e-12,8.44e-10,0.015,0.003,0    
k,1,0,0
k,2,0.2,0
LSTR,       1,       2
LATT,1,1,1, , , , 
lsel,all
LESIZE,all, , ,100, , , , ,1
lmesh,all
FINISH
/solu
ANTYPE,TRANS
d,1,all,all
outres,all,all
!alpha,10
TINTP,0
trnopt,full
deltim,1e-3
time,1e-3
kbc,1
autots,on
solve
time,1e-2
kbc,1
autots,off
F,101,FY,1
solve
time,1
fdele,all,all
kbc,1
autots,on
solve

! ***典型的用模态叠加法进行瞬态动力学分析的输入命令流：

!Build　the　Model

/FILNAM,...!　Jobname

/TITLE,...!　Title

/PREP7!　Enter　PREP7

---

---!　Generate　model

---

FINISH

!Obtain　the　Modal　Solution

/SOLU!　Enter　SOLUTION

ANTYPE,MODAL!　Modal　analysis

MODOPT,REDU!　Reduced　method

M,...!　Master　DOF

TOTAL,...

D,...!　Constraints

SF,...!　Element　loads

ACEL,...

SAVE

SOLVE

FINISH

!Obtain　the　Mode　Superposition　Transient　Solution

/SOLU!　Re-enter　SOLUTION

ANTYPE,TRANS!　Transient　analysis

TRNOPT,MSUP,...!　Mode　superposition　method

LVSCALE,...!　Scale　factor　for　element　loads

F,...!　Nodal　Loads

MDAMP,...!　Modal　damping　ratios

DELTIM,...!　Integration　time　step　sizes

LSWRITE!　Write　first　load　step　(Remember:　the　first　load　step

---!　is　solved　statically　at　time=0.)

---

---!　Loads,　etc.　for　2nd　load　step

TIME,...!　Time　at　end　of　second　load　step

KBC,...!　Ramped　or　stepped　loads

OUTRES,...!　Results-file　data　controls

---

LSWRITE!　Write　2nd　load　step　(first　transient　load　step)

SAVE

LSSOLVE!　Initiate　multiple　load　step　solution

FINISH

!Review　results　of　the　mode　superposition　solution

/POST26!　Enter　POST26

FILE,,RDSP!　Results　file　is　Jobname.RDSP

SOLU,...!　Store　solution　summary　data

NSOL,...!　Store　nodal　result　as　a　variable

PLVAR,...!　Plot　variables

PRVAR,...!　List　variables

FINISH

!Expand　the　Solution

/SOLU!　Re-enter　SOLUTION

EXPASS,ON!　Expansion　pass

NUMEXP,...!　No.　of　solutions　to　expand;　time　range

OUTRES,...!　Results-file　data　controls

SOLVE

FINISH

!Review　the　Results　of　the　Expanded　Solution

/POST1

SET,...!　Read　desired　set　of　results　into　database

PLDISP,...!　Deformed　shape

PRRSOL,...!　Reaction　loads

PLNSOL,...!　Contour　plot　of　nodal　results

PRERR!　Global　percent　error　(a　measure　of　mesh　adequacy)

---

---!　Other　postprocessing　as　desired

---

FINISH

 

 

!**下面给出的是可以概括用完全法进行瞬态动力学分析的过程的输入命令流：

!Build　the　Model

/FILNAM,...!　Jobname

/TITLE,...!　Title

/PREP7!　Enter　PREP7

---

---!　Generate　model

---

FINISH

!Apply　Loads　and　Obtain　the　Solution

/SOLU!　Enter　SOLUTION

ANTYPE,TRANS!　Transient　analysis

TRNOPT,FULL!　Full　method

D,...!　Constraints

F,...!　Loads

SF,...

ALPHAD,...!　Mass　damping

BETAD,...!　Stiffness　damping

KBC,...!　Ramped　or　stepped　loads

TIME,...!　Time　at　end　of　load　step

AUTOTS,ON!　Auto　time　stepping

DELTIM,...!　Time　step　size

OUTRES,...!　Results　file　data　options

LSWRITE!　Write　first　load　step

---

---!　Loads,　time,　etc.　for　2nd　load　step

---

LSWRITE!　Write　2nd　load　step

SAVE

LSSOLVE,1,2!　Initiate　multiple　load　step　solution

FINISH

!

!Review　the　Results

/POST26

SOLU,...!　Store　solution　summary　data

NSOL,...!　Store　nodal　result　as　a　variable

ESOL,,,,!　Store　element　result　as　a　variable

RFORCE,...!　Store　reaction　as　a　variable

PLVAR,...!　Plot　variables

PRVAR,...!　List　variables

FINISH

/POST1

SET,...!　Read　desired　set　of　results　into　database

PLDISP,...!　Deformed　shape

PRRSOL,...!　Reaction　loads

PLNSOL,...!　Contour　plot　of　nodal　results

PRERR!　Global　percent　error　(a　measure　of　mesh　adequacy)

---

---!　Other　postprocessing　as　desired

---

FINISH