wangliu 发表于 2012-8-7 09:55

求助ansys命令流解释,SSJ1一句是什么意思啊


我的命令流如下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

少少凌 发表于 2013-1-29 13:53

看不懂啊楼主。。。。

阳阳的飞翔 发表于 2013-3-27 00:06

顶一个,支持发问
页: [1]
查看完整版本: 求助ansys命令流解释,SSJ1一句是什么意思啊