无法连接弹簧单元
最近在做一个主轴的模态分析，将轴承等效为弹簧单元进行处理，但是在最后连接轴承单元的过程中出现问题。想请教一下大家程序源代码如下：
前辈您好，现在有一个关于COMBIN14单元的问题想咨询您一下。我现在想添加多个方向的弹簧刚度，但调用到E,NODE1,NODE2时出现问题，有的弹簧单元单元可以添加，有的无法添加。想请您指点一二，多谢！
附源代码：
! ===================== Local constants. ==========================
C_ANALYSIS_TYPE_MODAL = 1 ! Modal analysis
C_ANALYSIS_TYPE_UNB_RES = 2 ! Unbalanced response analysis
! ===================== User defined parameters ===================
P_RO = 8860 ! density of the rotor material, kg.m^-3
P_PEX = 1.95e11 ! elastic modulus of the rotor material, Pa
P_POISSON = 0.31 ! Poisson's ratio of the rotor material
P_BEARING_STIFFNESS = 8e7! Stiffness coefficients of the roller bearing.
! P_ANALYSIS_TYPE = C_ANALYSIS_TYPE_MODAL
P_ANALYSIS_TYPE = C_ANALYSIS_TYPE_UNB_RES
! ================== End of user defined parameters ===============
! ========================= 1. Pre-process ===========================
/PREP7
! == 1.2 == rotor material configuration
P_PGXY = P_PEX/(2*(1+P_POISSON))
MP, EX, 1, P_PEX
MP, DENS, 1, P_RO
MP,GXY, 1, P_PGXY
! MP, DMPR, 1, 0.005
! == 1.3 == prepare data for BEAM188 elements
! Set element type
ET, 1, BEAM188, , , 2
P_NELEMENTS = 10 ! total number of rotor elements
P_NRADII = P_NELEMENTS
! outer radii array
*DIM, P_ORADIUS, ARRAY, P_NRADII
P_ORADIUS(1)= 0.058
P_ORADIUS(2)= 0.058
P_ORADIUS(3)= 0.058
P_ORADIUS(4)= 0.062
P_ORADIUS(5)= 0.070
P_ORADIUS(6)= 0.076
P_ORADIUS(7)= 0.066
P_ORADIUS(8)= 0.066
P_ORADIUS(9)= 0.066
P_ORADIUS(10) = 0.066
! inner radii array
*DIM, P_IRADIUS, ARRAY, P_NRADII
P_IRADIUS(1)= 0.032
P_IRADIUS(2)= 0.032
P_IRADIUS(3)= 0.032
P_IRADIUS(4)= 0.032
P_IRADIUS(5)= 0.032
P_IRADIUS(6)= 0.032
P_IRADIUS(7)= 0.032
P_IRADIUS(8)= 0.040
P_IRADIUS(9)= 0.048
P_IRADIUS(10) = 0.054
! set data for BEAM188 element
*DO, I, 1, P_NRADII
*IF, P_IRADIUS(I), EQ, 0, THEN
SECTYPE, I, BEAM, CSOLID
SECDATA, P_ORADIUS(I)/2, 8, 8
*ELSE
SECTYPE, I, BEAM, CTUBE
SECDATA, P_IRADIUS(I)/2, P_ORADIUS(I)/2
*ENDIF
*ENDDO
! == 1.4 == prepare data for MASS21 elements
! set element type
ET, 2, MASS21
KEYOPT, 2, 1, 0
KEYOPT, 2, 2, 0
KEYOPT, 2, 3, 0
! prepare element nodes and real numbers
P_NMASSES = 1
P_MASS_REAL= 1
! real constants for MASS21 elements
! Mxx, Myy, Mzz, Ixx(polar moment of inertia), Iyy(diameter moment of inertia), Izz(polar moment of inertia)
! impeller & nut
R, P_MASS_REAL, 4, 4, 4, 1.3e-2, 8.32e-3, 8.32e-3
P_MASS_NODE= 3
! == 1.5 == prepare data for COMBI14 elements
! set element type for bearings/seals
ET, 11, COMBIN14
KX = 1.e+5 !Example stiffness value
KEYOPT,11,,1,0 ! X direction
R,21,KX
R,22,KX
ET, 12, COMBIN14
KY = 1.e+5 !Example stiffness value
KEYOPT,12,,2,0 ! X direction
R,23,KY
R,24,KY
ET, 13, COMBIN14
KZ = 1.e+5 !Example stiffness value
KEYOPT,13,,3,0 ! X direction
R,25,KZ
R,26,KZ
ET, 14, COMBIN14
KRY = 1.e+5 !Example stiffness value
KEYOPT,14,,5,0 ! X direction
R,27,KRY
R,28,KRY
ET, 15, COMBIN14
KRZ = 1.e+5 !Example stiffness value
KEYOPT,15,,6,0 ! X direction
R,29,KRZ
R,30,KRZ
! == 1.5 == geometry
! node list for the rotor
N, 1,0.00,0,0
N, 2,0.04,0,0
N, 3,0.08,0,0
N, 4,0.11,0,0
N, 5,0.185,0,0
N, 6,0.251,0,0
N, 7,0.351,0,0
N, 8,0.389,0,0
N, 9,0.415,0,0
N, 10, 0.441,0,0
N, 11, 0.467,0,0
! node list for the bearings/seals
N, 1000, 0.04,0,0
N, 1001, 0.04,0,0
N, 1002, 0.04,0,0
N, 1003, 0.04,0,0
N, 1004, 0.04,0,0
N, 9000,0.441,0,0
N, 9001,0.441,0,0
N, 9002,0.441,0,0
N, 9003,0.441,0,0
N, 9004,0.441,0,0
! == 1.6 == elements creation
! rotor elements
TYPE, 1
MAT,1
*DO, I, 1, P_NRADII
SECNUM, I
E, I, I+1
*ENDDO
! mass elements
TYPE, 2
REAL, P_MASS_REAL
E, P_MASS_NODE
! bearing elements
TYPE, 11
REAL, 21
E, 1000,2
REAL, 22
E, 9000,10
TYPE, 12
REAL, 23
E, 1001,2
REAL, 24
E, 9001,10
TYPE, 13
REAL, 25
E, 1002,2
REAL, 26
E, 9002,10
TYPE, 14
REAL, 27
E, 1003,2
REAL, 28
E, 9003,10
TYPE, 15
REAL, 29
E, 1004,2
REAL, 30
E, 9004,10
! element component
CM, PC_ROTOR, ELEM
ALLSEL
FINI
! ====================== 2. Solve =====================
/SOLU
应该是对单元类型 combin14 的定义有问题，看看以下资料吧：
页:
[1]