声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1831|回复: 0

[结构分析] 转子梁和实体连接及mpc问题

[复制链接]
发表于 2012-9-24 23:32 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
问题有几个:
(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

image1.png

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-18 06:30 , Processed in 0.060896 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表