强迫运动法 | 解决固支基础条件不能求解位移/加速度传递率的问题

FINISH/CLEAR/OUT,EMM_SUMMARY,OUT  !将输出信息保存到.out文件/TITLE, ENFORCED MOTION METHOD /UNIT,SI/FILNAME,HARMIC_EMM,1/PREP7                                                        ! ELEMENT TYPE    ! **************ET,1,SOLID185  KEYOPT,1,2,3
! MATERIALS PROPERTIES! *********************MP,EX,1,210E9     MP,PRXY,1,0.3MP,DENS,1,7850

! GEOMETRY ! *********L=10E-2  !PLATE LENGTH W=10E-2 !PLATE WIDTH T=0.1E-2 !PLATE THICKNESS 
! Acceleration inputNG = 1 ! Number of g'sG = 9.81 ! Gravity accelerationTOTG = NG*G
! PLATE MODELING! ****************BLOCK,0,L,0,W,0,T
! MESH CONTROL! *************LSEL,S,LOC,Z,0LSEL,A,LOC,Z,TLESIZE,ALL,,,60ALLSEL,ALLLSEL,S,LENGTH,,TLESIZE,ALL,,,4ALLSEL,ALL
! MESHING! *******VSEL,ALLVATT,1,,1VMESH,ALLALLSEL,ALL
! BC'S - CCCC! *****NSEL,S,LOC,X,0NSEL,A,LOC,X,LNSEL,A,LOC,Y,0NSEL,A,LOC,Y,WD,ALL,ALL,0CM,BASE_NODES,NODE ! BASE NODESALLSEL,ALL


! MODE SUPERPOSITION HARMONIC ANALYSIS SOLUTION! **********************************************! MODAL ANALYSIS! ****************/SOLUANTYPE,2NMODE=10MODOPT,LANB,NMODE   !Block Lanczos 方法EQSLV,SPAR !采用稀疏矩阵求解器MXPAND,NMODE , , ,1 LUMPM,0 PSTRES,0MODCONT,,ON ! ACTIVATE ENFORCED MOTIONCMSEL,S,BASE_NODES,NODED,ALL,UZ,1      !  !这里定义强制运动标识1,注意里1只表示标识1，不表示位移大小ALLSEL,ALLSOLVE FINISH
! HARMONIC ANALYSIS! ******************FREQ1=1      ! FIRST FREQUENCY IN RANGEFREQ2=10000      ! LAST FREQUENCY IN RANGEFREQSTEP=(FREQ2-FREQ1)*1  ! SUBSTEPS
/SOLU ANTYPE,3 HROPT,MSUP, , ,0HROUT,ONDVAL,1,ACC,TOTG  ! 在标识1基础上施加加速度激励HARFRQ,FREQ1,FREQ2,NSUBST,FREQSTEP,  KBC,1     ALPHAD,0,   BETAD,0,DMPRAT,0.005SOLVE SAVEFINISH 
! TIME HISTORY POST-PROCESSING! ****************************/POST26 file,,rfrqNSOL,2,5702,U,Z,                     !  NODE 5702 DISPLACEMENT/AXLAB,X,Frequency [Hz]  /AXLAB,Y,Absolute Response [m]  /GROPT,LOGX,OFF  /GROPT,LOGY,ON PLVAR,2