原文地址：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
