马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我的命令流如下ssj1是什么意思啊
/PREP7
/NUMBER,0 !0 - Color (terminal dependent) the numbered items and show numbers
/VIEW,,1,1,1
*SET,PI,3.14159
*SET,DIAO,43E-3
*SET,DIAOO,38E-3
*SET,DIAI,15E-3
*SET,HEI,4E-3
*SET,D1,5E-3
*SET,D2,8E-3
*SET,D3,19E-3
*SET,D10,D1/2+2*HEI+D2
*SET,D11,D10+4*HEI
*SET,D12,D11+D3
ET,1,SOLID227,1001 !IT IS A PIEZOELECTRIC ELEMENT
MP,DENS,1,7800 !STEEL DENSITY
MP,NUXY,1,0.3 !STEEL POISSON RATIO
MP,EX,1,205E9 !STEEL MODULUS OF ELASTICITY
MP,DENS,4,2810 !ALUMINIUM DENSITY
MP,NUXY,4,0.33 !ALUMINIUM POISSON RATIO
MP,EX,4,72E9 !ALUMINIUM MODULUS OF ELASTICITY
*SET,E31,-2.4
*SET,E33,17.3
*SET,E16,12.95
*SET,C11,15E10
*SET,C12,8.4E10
*SET,C13,8.4E10
*SET,C33,12.9E10
*SET,C66,2.67E10
*SET,C44,3.33E10
*SET,DIE11,765
*SET,DIE33,640
MP,DENS,2,7600
TB,PIEZ,2 !DEFINE PIEZO OF PZT4_POLAR_Z
TBDATA,3,E31
TBDATA,6,E31
TBDATA,9,E33
TBDATA,14,E16
TBDATA,16,E16
TB,ANEL,2 !DEFINE STRUCTURE OF PZT4_POLAR_Z
TBDATA,1,C11,C12,C13
TBDATA,7,C11,C13
TBDATA,12,C33
TBDATA,16,C44
TBDATA,19,C66
TBDATA,21,C66
MP,PERX,2,DIE11 !PERMITTIVITY(X AND Y DIRECTION ) OF PZT4_PLOAR_Z
MP,PERZ,2,DIE33 !PERMITTIVITY(Z DIRECTION) OF PZT4_POLAR_Z
MP,DENS,3,7600
TB,PIEZ,3 !DEFINE PIEZO OF PZT4_POLAR_-Z
TBDATA,3,-E31
TBDATA,6,-E31
TBDATA,9,-E33
TBDATA,14,-E16
TBDATA,16,-E16
TB,ANEL,3 !DEFINE STRUCTURE OF PZT4_POLAR_-Z
TBDATA,1,C11,C12,C13
TBDATA,7,C11,C13
TBDATA,12,C33
TBDATA,16,C44
TBDATA,19,C66
TBDATA,21,C66
MP,PERX,3,DIE11 !PERMITTIVITY(X AND Y DIRECTION ) OF PZT4_PLOAR_-Z
MP,PERZ,3,DIE33 !PERMITTIVITY(Z DIRECTION) OF PZT4_POLAR_-Z
C*** WE HAVE DEFINED ALL MATERIAL PROPERTY
C***NOW WE BEGIN TO MAKE THE GEOMETRY MODEL
CYLIN,(DIAI-3E-3)/2,,-D12,D12 !THE BOLT : 14mm LONGER THAN D11
CM,BOLT,VOLU !COMPONENT THE BOLT
VSEL,NONE
BLOCK,-DIAOO/2,DIAOO/2,-DIAO/2,DIAO/2,-D1/2,D1/2 !THE FLANGE PART
CM,FLANGE,VOLU
VSEL,NONE
SSJ1,DIAOO,DIAO,D1/2,D1/2+HEI,DIAI
SSJ1,DIAOO,DIAO,D1/2+HEI,D1/2+2*HEI,DIAI
VSEL,S,LOC,Z,D1/2,D1/2+2*HEI
VSYMM,Z,ALL
VSEL,NONE
SSJ2,-DIAOO/2,DIAO,D10,D10+HEI,DIAI
SSJ2,DIAOO/2,DIAO,D10,D10+HEI,DIAI
SSJ2,-DIAOO/2,DIAO,D10+HEI,D10+2*HEI,DIAI
SSJ2,DIAOO/2,DIAO,D10+HEI,D10+2*HEI,DIAI
SSJ2,-DIAOO/2,DIAO,D10+2*HEI,D10+3*HEI,DIAI
SSJ2,DIAOO/2,DIAO,D10+2*HEI,D10+3*HEI,DIAI
SSJ2,-DIAOO/2,DIAO,D10+3*HEI,D10+4*HEI,DIAI
SSJ2,DIAOO/2,DIAO,D10+3*HEI,D10+4*HEI,DIAI
VSEL,S,LOC,Z,D10,D11
VSYMM,Z,ALL
vsel,none
BLOCK,-DIAOO/2,DIAOO/2,-DIAO/2,DIAO/2,D11,D12 !THE BLOCK BASE OF THE HORN
CM,AAAA,VOLU
VSEL,NONE
CYLIND,(DIAI-3E-3)/2,,D11+50E-3,D11-50E-3
CM,BBBB,VOLU
CMSEL,A,AAAA
VSBV,ALL,BBBB
VSYMM,Z,ALL
VSEL,ALL
VSEL,U,VOLU,,BOLT
VGLUE,ALL
VSEL,S,LOC,Z,D11,D12
VSEL,A,LOC,Z,-D11,-D12
VSEL,A,VOLU,,BOLT
VGLUE,ALL
VSEL,S,LOC,Z,0
VADD,ALL
NUMCMP,ALL
C***NOW WE BEGIN TO ASSIGN ATTRIBUTE TO ALL THE VOLUMES
VSEL,S,LOC,Z,0
VATT,1,,1 !ASSIGN MATERIAL ATTRIBUTES TO ALL STEEL PARTS
VSEL,S,LOC,Z,D11,D12
VSEL,A,LOC,Z,-D11,-D12
VSEL,A,LOC,Z,D1/2+2*HEI,D10
VSEL,A,LOC,Z,-(D1/2+2*HEI),-D10
VATT,4,,1 !ASSIGN MATERIAL ATTRIBUTES TO ALL ALUMINUM PARTS
VSEL,S,LOC,Z,D1/2,D1/2+HEI
VATT,3,,1
VSEL,S,LOC,Z,D1/2+HEI,D1/2+2*HEI
VATT,2,,1
VSEL,S,LOC,Z,-D1/2,-D1/2-HEI
VATT,2,,1
VSEL,S,LOC,Z,-D1/2-HEI,-D1/2-2*HEI
VATT,3,,1 !ASSIGN MATERIAL ATTIBUTES TO FOUR LONGITUDINAL
VSEL,S,LOC,Z,D10,D10+HEI
VSEL,R,LOC,X,,DIAO
VATT,3,,1
VSEL,S,LOC,Z,D10,D10+HEI
VSEL,R,LOC,X,,-DIAO
VATT,2,,1
VSEL,S,LOC,Z,D10+HEI,D10+2*HEI
VSEL,R,LOC,X,,DIAO
VATT,2,,1
VSEL,S,LOC,Z,D10+HEI,D10+2*HEI
VSEL,R,LOC,X,,-DIAO
VATT,3,,1
VSEL,S,LOC,Z,D10+2*HEI,D10+3*HEI
VSEL,R,LOC,X,,DIAO
VATT,3,,1
VSEL,S,LOC,Z,D10+2*HEI,D10+3*HEI
VSEL,R,LOC,X,,-DIAO
VATT,2,,1
VSEL,S,LOC,Z,D10+3*HEI,D11
VSEL,R,LOC,X,,DIAO
VATT,2,,1
VSEL,S,LOC,Z,D10+3*HEI,D11
VSEL,R,LOC,X,,-DIAO
VATT,3,,1 !ASSIGN MATERIAL ATTRIBUTES OF LEFT SIDE BENDING
VSEL,S,LOC,Z,-D10,-D10-HEI
VSEL,R,LOC,X,,DIAO
VATT,3,,1
VSEL,S,LOC,Z,-D10,-D10-HEI
VSEL,R,LOC,X,,-DIAO
VATT,2,,1
VSEL,S,LOC,Z,-D10-HEI,-D10-2*HEI
VSEL,R,LOC,X,,DIAO
VATT,2,,1
VSEL,S,LOC,Z,-D10-HEI,-D10-2*HEI
VSEL,R,LOC,X,,-DIAO
VATT,3,,1
VSEL,S,LOC,Z,-D10-2*HEI,-D10-3*HEI
VSEL,R,LOC,X,,DIAO
VATT,3,,1
VSEL,S,LOC,Z,-D10-2*HEI,-D10-3*HEI
VSEL,R,LOC,X,,-DIAO
VATT,2,,1
VSEL,S,LOC,Z,-D10-3*HEI,-D11
VSEL,R,LOC,X,,DIAO
VATT,2,,1
VSEL,S,LOC,Z,-D10-3*HEI,-D11
VSEL,R,LOC,X,,-DIAO
VATT,3,,1 !ASSIGN MATERIAL ATTRIBUTES OF RIGHT SIDE BENDING
VSEL,ALL
ESIZE,5E-3
VMESH,ALL
C*** THE PREP/7 WORK DONE
CSYS,1
NSEL,S,LOC,Z,D1/2+HEI
NSEL,A,LOC,Z,-D1/2-HEI
NSEL,R,LOC,X,DIAI/2,DIAO
*GET,NNE1,NODE,,COUNT !NODE NUMBER OF THE ELECTRODE1
CM,ELECTRODE1,NODE !DIFINED THE ELECTRODE1
CP,1,VOLT,ALL
*GET,NELECTRODE1,NODE,0,NUM,MIN
NSEL,S,LOC,Z,D10+HEI
NSEL,A,LOC,Z,D10+3*HEI
NSEL,A,LOC,Z,-D10-HEI
NSEL,A,LOC,Z,-D10-3*HEI
NSEL,R,LOC,X,DIAI/2,DIAO
*GET,NNE2,NODE,,COUNT !NODE NUMBER OF THE ELECTRODE2
CM,ELECTRODE2,NODE !DEFINED THE ELECTRODE2
CP,2,VOLT,ALL
*GET,NELECTRODE2,NODE,,NUM,MIN
NSEL,S,LOC,Z,D1/2
NSEL,A,LOC,Z,-D1/2
NSEL,A,LOC,Z,D1/2+2*HEI
NSEL,A,LOC,Z,-D1/2-2*HEI
NSEL,A,LOC,Z,D10
NSEL,A,LOC,Z,-D10
NSEL,A,LOC,Z,D10+2*HEI
NSEL,A,LOC,Z,-D10-2*HEI
NSEL,A,LOC,Z,D11
NSEL,A,LOC,Z,-D11
NSEL,R,LOC,X,DIAI/2,DIAO
CM,ELECTRODE4,NODE !DEFINED THE ELECTRODE4
CP,4,VOLT,ALL
*GET,NELECTRODE4,NODE,,NUM,MIN
CSYS
NSEL,ALL !IMPORTANT !! NEVER FORGET!!
D,NELECTRODE2,VOLT,0
D,NELECTRODE4,VOLT,0
D,NELECTRODE1,VOLT,0
FINISH
/SOLU
ANTYP,MODAL
MODOPT,LANB,20,15000,25000
MXPAND,20
SOLVE
FINISH
/prep7
D,NELECTRODE2,VOLT,0,100
D,NELECTRODE4,VOLT,0
D,NELECTRODE1,VOLT,100
/SOLU
EQSLV,SPARSE
ANTYPE,HARMIC
OUTRES,ALL,ALL
HARFRQ,19650,19800
NSUBST,20
KBC,1
SOLVE
FINISH
/POST26
NUMVAR,20
RFORCE,2,NELECTRODE1,CHRG
RFORCE,3,NELECTRODE2,CHRG
PI2=3.14159*2
PROD,4,2,1,,AMPSL,,,PI2,
NSOL,5,NELECTRODE1,VOLT
QUOT,6,5,4,,OMUL,,,
PROD,7,3,1,,AMPSF,,,PI2,
NSOL,8,NELECTRODE2,VOLT
QUOT,9,8,7,,OMUF,,,
NSOL,10,2845,U,Z
NSOL,11,2845,U,X
NSOL,12,2845,U,Y
NSOL,13,5883,U,Z
NSOL,14,5883,U,X
NSOL,15,5883,U,Y
SAVE,1
C*** START DEFINE THE VOLTA AND VOLTB
*AFUN,RAD
*SET,PI,3.1415926
*SET,RFREQ,19720
*SET,T_P,1/RFREQ !ZHOU QI (S)
*SET,VOLTARMS,100 !ROOT MEAN SQUARE VALUE OF VOLTA
*SET,VOLTBRMS,100 !ROOT MEAN SQURAE VALUE OF VOLTB
*SET,VLTA,VOLTARMS*SQRT(2) !PEAK VALUE OF VOLTA
*SET,VLTB,VOLTBRMS*SQRT(2) !PEAK VALUE OF VOLTB
*SET,PH_D_D,90 !PHASE DIFFERENCE OF A AND B IN DEGREE
*SET,PH_D,PH_D_D/180*PI !PHASE DIFFERENCE IN RADIAS
*SET,T_S,PH_D_D/360*T_P !PHASE DIFFERENCE TIME SPAN
*SET,PRD_N,30 !PERIOD NUMBER OF SIN INPUT SIGNAL
*SET,PRT_N,64 !PARTATION NUMBER IN ONE CRICLE
*SET,T_ICR,T_P/PRT_N ! TIME INCREMENT
*SET,LST_TIME,T_P*PRD_N !SETS THE TIME FOR A LOAD STEP:TIME,LST_TIME
*SET,LSTN_1C,16 !LOAD SUBSTEP NUMBER IN ONE CIRCLE
*SET,LSTN_W,LSTN_1C*PRD_N !THE WHOLE LOAD SUBSTEP NUMBER
*DIM,V_A,TABLE,1+PRT_N*PRD_N,1,1
*SET,V_A(0,1,1),0
*SET,V_A(1,0,1),0
*SET,V_A(1,1,1),0
*DO,I,1,PRD_N
*DO,J,1,PRT_N
*SET,V_A((I-1)*PRT_N+J+1,0,1),(I-1)*T_P+J*T_ICR
*ENDDO
*DO,K,1,PRT_N
*SET,V_A((I-1)*PRT_N+K+1,1,1),VLTA*SIN((2*PI/T_P)*K*T_ICR)
*ENDDO
*ENDDO
*DIM,V_B,TABLE,1+PRT_N*PRD_N,1,1
*SET,V_B(0,1,1),0
*SET,V_B(1,0,1),0
*SET,V_B(1,1,1),VLTB*SIN(PH_D)
*DO,I,1,PRD_N
*DO,J,1,PRT_N
*SET,V_B((I-1)*PRT_N+J+1,0,1),(I-1)*T_P+J*T_ICR
*ENDDO
*DO,K,1,PRT_N
*SET,V_B((I-1)*PRT_N+K+1,1,1),VLTB*SIN((2*PI/T_P)*(K*T_ICR+T_S))
*ENDDO
*ENDDO
C***END DEFINE THE VOLTA AND VOLTB
*VPLOT,V_A(1,0,1),V_A(1,1,1)
*VPLOT,V_B(1,0,1),V_B(1,1,1)
/SOLU
ANTYPE,4
TRNOPT,FULL
D,NELECTRODE4,VOLT,0
D,NELECTRODE1,VOLT,%V_A%
D,NELECTRODE2,VOLT,%V_B%
NSUBST,LSTN_W,2*LSTN_W,LSTN_W/2
OUTRES,ALL,ALL
KBC,0
TIME,LST_TIME
SOLVE
/POST26
NUMVAR,200
NSOL,2,2847,U,Z,U_Z_L
NSOL,3,2847,U,X,U_X_L
NSOL,5,18950,U,Z,U_Z_R
NSOL,6,18950,U,X,U_X_R
NSOL,7,18950,U,Y,U_Y_R
XVAR,1
PLVAR,2,3
NSOL,8,2880,U,Z,U_Z_R_120
NSOL,9,2880,U,X,U_X_R_120
NSOL,10,2880,U,Y,U_Y_R_120
SSJ1,DIAOO,DIAO,D1/2+2*HEI,D10,DIAI
VSYMM,Z,ALL
VSEL,NONE |