ANSYS APDL 进行不平衡激励启动过程的瞬态分析时,出现下列警告是哪里出问题了?
警告:① 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