警告：① The MINREF field of the CNVTOL command has been assigned a default value of 1.0 for convergence on F.  Check results carefully.

       ② BEAM189 will generate a gyroscopic damping matrix contribution(stationary frame ON on the CORIOLIS command) in a dynamic analysis.The element matrix will be non-symmetric.

       ③ The gyroscopic damping matrix (stationary frame ON in the CORIOLIS command) is based on small deflection theory.  The large deflection effect (NLGEOM,ON) needs to be used with caution.                                        

APDL语言如下：

FINISH

/CLEAR

/CONFIG,NRES,10000

/PREP7

L1=0.01

L2=0.112

L3=0.0162

L4=0.034

L5=0.041

L6=0.01

L0=L3+L4+L5

D=0.01

R1=0.051

R2=0.0208

R3=0.0295

T1=0.0187

T2=0.015

T3=0.0187

MUNB1=1E-4

MUNB2=6E-5

MUNB3=2E-5

MUNBR1=0.01

MUNBR2=0.01

MUNBR3=0.01

ET,1,BEAM189

ET,2,COMBI214,,1

MP,EX,1,2.19E11

MP,PRXY,1,0.3

MP,DENS,1,7830

SECTYPE,1,BEAM,CSOLID

SECDATA,R1,16,4

SECTYPE,2,BEAM,CSOLID

SECDATA,R2,16,4

SECTYPE,3,BEAM,CSOLID

SECDATA,R3,16,4

SECTYPE,4,BEAM,CSOLID

SECDATA,D/2

R,1,2.4E8,1.8E8,3E7,3E7,2.15E5,7.5E4  

RMORE,6E4,6E4

R,2,8.2E7,7.4E7,-2E6,-2E6,2.45E5,1.50E5

RMORE,5.5E4,5.5E4

K,1

K,2,,,L1-T1/2

K,3,,,L1

K,4,,,L1+T1/2

K,5,,,L1+L2-T2/2

K,6,,,L1+L2

K,7,,,L1+L2+T2/2

K,8,,,L1+L2+L0-T3/2

K,9,,,L1+L2+L0

K,10,,,L1+L2+L0+T3/2

K,11,,,L1+L2+L0+L6

K,12,,R1

K,13,,R2

K,14,,R3

L,1,2

L,2,3

L,3,4

L,4,5

L,5,6

L,6,7

L,7,8

L,8,9

L,9,10

L,10,11

LSEL,S,,,1,10,3

LATT,1,,1,,,14,4

LSEL,S,,,2,3

LATT,1,,1,,,14,1

LSEL,S,,,5,6

LATT,1,,1,,,14,2

LSEL,S,,,8,9

LATT,1,,1,,,14,3

LSEL,ALL

ESIZE,D/2

LMESH,ALL

N1=NODE(0,0,L1+L2+L3)

N2=NODE(0,0,L1+L2+L3+L4)

N,1000,,0.1,NZ(N1)

N,2000,,0.1,NZ(N2)

TYPE,2

REAL,1

MAT,1

E,N1,1000

E,N2,2000

D,1000,ALL

D,2000,ALL

D,N1,,,UZ,,,ROTZ

D,N2,ROTZ

FINISH

PI=ACOS(-1)

SPIN=4000*PI/30

DT=5E-4

TEND=3

SPINSLO=SPIN/TEND

NPOINT1=NINT(TEND/DT)+1

NPOINT2=NPOINT1+NINT(1/DT)+1

*DIM,FXU1,TABLE,NPOINT2,,,TIME

*DIM,FYU1,TABLE,NPOINT2,,,TIME

*VFILL,FXU1(1,0),RAMP,0,DT

*VFILL,FYU1(1,0),RAMP,0,DT

*DO,I,1,NPOINT1

TI=(I-1)*DT

SPINTI=SPINSLO*TI

SPINTI2=SPINTI**2

ROTTI=SPINSLO*TI**2/2

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU1(I,1)=MUNB1*MUNBR1*(SPINTI2*COSR-SPINSLO*SINR)

FYU1(I,1)=MUNB1*MUNBR1*(SPINTI2*SINR+SPINSLO*COSR)

*ENDDO

*DO,I,NPOINT1,NPOINT2

TI=(I-1)*DT

SPIN2=SPIN**2

ROTTI=SPINSLO*TEND**2/2+SPIN*(TI-TEND)

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU1(I,1)=MUNB1*MUNBR1*SPIN2*COSR

FYU1(I,1)=MUNB1*MUNBR1*SPIN2*SINR

*ENDDO

*DIM,FXU2,TABLE,NPOINT2,,,TIME

*DIM,FYU2,TABLE,NPOINT2,,,TIME

*VFILL,FXU2(1,0),RAMP,0,DT

*VFILL,FYU2(1,0),RAMP,0,DT

*DO,I,1,NPOINT1

TI=(I-1)*DT

SPINTI=SPINSLO*TI

SPINTI2=SPINTI**2

ROTTI=SPINSLO*TI**2/2

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU2(I,1)=MUNB2*MUNBR2*(SPINTI2*COSR-SPINSLO*SINR)

FYU2(I,1)=MUNB2*MUNBR2*(SPINTI2*SINR+SPINSLO*COSR)

*ENDDO

*DO,I,NPOINT1,NPOINT2

TI=(I-1)*DT

SPIN2=SPIN**2

ROTTI=SPINSLO*TEND**2/2+SPIN*(TI-TEND)

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU2(I,1)=MUNB2*MUNBR2*SPIN2*COSR

FYU2(I,1)=MUNB2*MUNBR2*SPIN2*SINR

*ENDDO

*DIM,FXU3,TABLE,NPOINT2,,,TIME

*DIM,FYU3,TABLE,NPOINT2,,,TIME

*VFILL,FXU3(1,0),RAMP,0,DT

*VFILL,FYU3(1,0),RAMP,0,DT

*DO,I,1,NPOINT1

TI=(I-1)*DT

SPINTI=SPINSLO*TI

SPINTI2=SPINTI**2

ROTTI=SPINSLO*TI**2/2

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU3(I,1)=MUNB3*MUNBR3*(SPINTI2*COSR-SPINSLO*SINR)

FYU3(I,1)=MUNB3*MUNBR3*(SPINTI2*SINR+SPINSLO*COSR)

*ENDDO

*DO,I,NPOINT1,NPOINT2

TI=(I-1)*DT

SPIN2=SPIN**2

ROTTI=SPINSLO*TEND**2/2+SPIN*(TI-TEND)

SINR=SIN(ROTTI)

COSR=COS(ROTTI)

FXU3(I,1)=MUNB3*MUNBR3*SPIN2*COSR

FYU3(I,1)=MUNB3*MUNBR3*SPIN2*SINR

*ENDDO

/SOLU

ANTYPE,TRANSIENT

NLGEOM,ON

CORIOLIS,ON,,,ON

OUTRES,ALL,ALL

AUTOTS,ON

CNVTOL,F,,0.005

DMPRAT,0.005

KBC,0

TIME,1.5

DELTIM,DT,0.1*DT,10*DT

OMEGA,,,SPINSLO*1.5

N3=NODE(0,0,L1)

F,N3,FX,%FXU1%

F,N3,FY,%FYU1%

N4=NODE(0,0,L1+L2)

F,N4,FX,%FXU2%

F,N4,FY,%FYU2%

N5=NODE(0,0,L1+L2+L0)

F,N5,FX,%FXU3%

F,N5,FY,%FYU3%

SOLVE

TIME,TEND

DELTIM,DT,0.1*DT,2*DT

OMEGA,,,SPIN

SOLVE

TIME,TEND+1

SOLVE

FINISH

/POST26

NSOL,2,N3,U,X

PROD,3,2,2

NSOL,4,N3,U,Y

PROD,5,4,4

ADD,6,3,5

SQRT,7,6

NSOL,2,N4,U,X

PROD,3,2,2

NSOL,4,N4,U,Y

PROD,5,4,4

ADD,6,3,5

SQRT,8,6

NSOL,2,N5,U,X

PROD,3,2,2

NSOL,4,N5,U,Y

PROD,5,4,4

ADD,6,3,5

SQRT,9,6

PLVAR,7,8,9

XVAR,2

PLVAR,4

*GET,TMAX7,VARI,7,EXTREM,TMAX

/POST1

/VIEW,1,1,1,1

/ESHAPE,1

SET,,,,,TMAX7

PLNSOL,U,SUM

PLNSOL,S,X

@htbbzzg