转子梁和实体连接及mpc问题
问题有几个：(1) 无论用combin14或combi214单元，均提示节点 i and j not concident，在网上查到说是因为坐标系不一致，这个问题有解决的么？
(2) mpc约束问题，命令流很大一块MPC,复制粘贴log中的，本例中做梁和solid45的连接，点和面连接，试了很多约束方法，包括all dof 约束，仅ux，uy约束，均求不出结果，这个有研究过的么，给指点一二
不胜感激！！！！！！！！！！！！！！！！
命令流如下
/CLEAR
/BATCH,list
/PREP7
/UNITS,SI ! Define units used: meter
nbstep = 46
spin = 4500
dspin = spin/(nbstep-1)
*dim,spins,,nbstep
*vfill,spins,ramp,0.,dspin
spins(1) = 0.1 !! non zero to easy the
nb=6
INNER_RAD = 0.05
DIAM = 2*INNER_RAD
L = 1.5
N_n=7 ! number of beam nodes
ET,100,BEAM188,,,2
KEYOPT,100,3,2
SECTYPE,100,BEAM,CSOLID
SECDATA,INNER_RAD,
MP,EX,100,2.11e+11 ! Young's modulus
MP,DENS,100,7810 ! Density
MP,NUXY,100,0.3 ! Possion ratio
!DEFINE PARAMETERS OF DISK
OUTER_RAD = 0.28
DISK_THICK = 0.07
DENS = 7800
N,1,0,0,0
N,N_n,0,0,L
FILL,1,N_n
TYPE,100
MAT,100
real,100
*DO,I,1,N_n-1
E,I,I+1
*END DO
et,200,45
cylind,INNER_RAD,0.28,0.5-0.035,0.5+0.035
type,200
mat,100
real,200
SMRT,8
MSHAPE,1,3D
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
cm,Spool,elem
ET,300,14
KEYOPT,300,2,1
KEYOPT,300,3,0
R,300,1E+6
N,10000,0.3,0
N,20000,0.3,0,1.5
type,300
real,300
E,1,10000
E,7,20000
D,all,uz
d,all,rotz
d,10000,all
d,20000,all
numcmp,all
allsel
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
! Define surface-based constraint type of pair
MAT,1
R,3
REAL,3
ET,4,170
ET,5,175
KEYOPT,5,12,5
KEYOPT,5,4,0
KEYOPT,5,2,2
KEYOPT,4,2,0
KEYOPT,4,4,111111
TYPE,4
! Create a pilot node
TSHAP,PILO
E,3
! Generate the contact surface
ASEL,S,,,5
ASEL,A,,,6
CM,_CONTACT,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
/MREP,EPLOT
allsel
/solu
antype,static
coriolis,on,,,on
CMOMEGA,Spool,,,,spins(5)
pstr,on
campbell,on,nbstep !! prestress Campbell analysis
solve
fini
/solu
antype,modal
modopt,qrdamp,nb,,,on
mxpand,nb
CMOMEGA,Spool,,,,spins(5)
pstr,on
solve
fini
页:
[1]