24小时热门版块排行榜    

查看: 2238  |  回复: 8

starkyu

新虫 (初入文坛)

[求助] [求助]用ANSYS耦合平面出现的错误~设置了周期对称的节点还能再进行平面耦合约束么?

对螺栓法兰连接系统进行分析,设置的是周期对称边界条件,考虑螺栓的预紧力,其余的边界条件是想在法兰一端进行轴向约束,另一端进行一个平面耦合Y方向的约束,但是每次以计算就会出现如下错误:
*** ERROR ***                           CP =      91.885   TIME= 20:02:29
Constraint equation 940 has node 1 = 1898 degree of freedom = UY which  
is also a slave in coupled set 1.                                       

*** ERROR ***                           CP =      91.994   TIME= 20:02:30
Degree of freedom UY at node 2312 is both in constraint equation 940   
and a slave degree of freedom of coupled set 1, which may cause a      
problem.  Please replace the CP,1,UY,1896,...,2312,..., command with   
the CE,next,0.0,2312,UY,-1.0,1896,UY,1.0 command.

已经纠结了好多天,实在不知道要怎么解决,求大神帮忙~~
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夕阳雨

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
starkyu: 金币+10, ★★★很有帮助, 难得有人回答 去研究一下 谢谢了~ 2013-10-25 15:34:44
我这里有一个实例你看一下《螺栓预紧法兰盘衬垫分析》
/FILNAME,Ring Flange,0
/PREP7
H1=0.01  !法兰盘1的厚度
H2=0.002 !衬垫厚度
H3=0.01  !法兰盘2的厚度
H4=0.008 !螺帽的厚度
R1=0.05  !法兰盘内半径
R2=0.08  !法兰盘外半径
R3=(R2-R1)/2+R1 !螺栓中心位置
R4=0.5*(R2-R1)/2!螺栓半径
R5=0.7*(R2-R1)/2!螺帽半径

ET,1,186
ET,2,194

MP,EX,1,1.05E11
MP,NUXY,1,0.3

TB,GASKET,2,,13,COMP             !材料压缩曲线
TBPT,,0.508000E-04, 0.161226E+07
TBPT,,0.101600E-03, 0.520884E+07
TBPT,,0.152400E-03, 0.113134E+08
TBPT,,0.203200E-03, 0.200499E+08
TBPT,,0.254000E-03, 0.259960E+08
TBPT,,0.304800E-03, 0.290345E+08
TBPT,,0.355600E-03, 0.357453E+08
TBPT,,0.406400E-03, 0.440064E+08
TBPT,,0.457200E-03, 0.563189E+08
TBPT,,0.508000E-03, 0.748254E+08
TBPT,,0.558800E-03, 0.967287E+08
TBPT,,0.609600E-03, 0.129001E+09
TBPT,,0.683260E-03, 0.157147E+09
TB,GASKET,2,,5,LUNL              !定义衬垫的线性卸载曲线
TBPT,,0.152400E-03, 2.430000E+11
TBPT,,0.304800E-03, 3.565000E+11
TBPT,,0.406400E-03, 5.923000E+11
TBPT,,0.558800E-03, 1.088000E+12
TBPT,,0.683260E-03, 1.490000E+12

CYLIND,R1,R2,0,H1,0,360          !定义圆柱体
CYLIND,R1,R2,H1,H1+H2,0,360      !定义圆柱体
CYLIND,R1,R2,H1+H2,H1+H2+H3,0,360!定义圆柱体
NUMMRG,ALL                       !合并所有实体

CYLIND,R4,0,0,H1+H2+H3,0,360,     !定义圆柱体
VGEN,,4,,,R3,,,,,1                !移动圆柱体
CSYS,1                            !激活总体柱坐标系
!复制圆柱体
VGEN,2,4,,,,90,,,0            
VGEN,2,4,,,,180,,,0
VGEN,2,4,,,,270,,,0

!把圆柱体从法兰盘中减去
FLST,2,3,6,ORDE,2
FITEM,2,1
FITEM,2,-3
VSBV,P51X,4


FLST,2,3,6,ORDE,2
FITEM,2,8
FITEM,2,-10
VSBV,P51X,5


FLST,2,3,6,ORDE,2
FITEM,2,1
FITEM,2,-3
VSBV,P51X,6


FLST,2,3,6,ORDE,3
FITEM,2,4
FITEM,2,-5
FITEM,2,8
VSBV,P51X,7
!定义圆柱体
CYLIND,R4*0.95,0,0,H1+H2+H3,0,360
CYLIND,R5,0,H1+H2+H3,H1+H2+H3+H4,0,360
CYLIND,R5,0,-H4,0,0,360
!把体4,5和6进行加运算
VADD,4,5,6

CSYS,0  !激活总体笛卡尔坐标系
VGEN,,7,,,R3,,,,,1
CSYS,1
VGEN,2,7,,,,90,,,0
VGEN,2,7,,,,180,,,0
VGEN,2,7,,,,270,,,0

!设置总体单元尺寸
ESIZE,0.005

TYPE,1
MAT,1
VSWEEP,2
!划分衬垫单元
TYPE,2
MAT,2
IMESH,AREA,19,25, ,0,0,H2,0.0005

TYPE,1
MAT,1
VSWEEP,3


MSHKEY,0
MSHAPE,1,3D
VMESH,4
VMESH,5
VMESH,6
VMESH,7

!定义预拉伸单元
PSMESH, ,P1, ,VOLU,4,   0,Z,(H1+H2+H3)/2,,,,,
PSMESH, ,P2, ,VOLU,5,   0,Z,(H1+H2+H3)/2, , , , ,  
PSMESH, ,P3, ,VOLU,6,   0,Z,(H1+H2+H3)/2, , , , ,  
PSMESH, ,P4, ,VOLU,7,   0,Z,(H1+H2+H3)/2, , , , ,  


!定义接触单元节点组1
ASEL,S, , ,48
NSLA,S,1
CM,CN1,NODE
ALLSEL,ALL   
!定义接触单元节点组2
ASEL,S, , ,28
NSLA,S,1
CM,CN2,NODE
ALLSEL,ALL   
!定义目标单元节点组1
FLST,5,4,5,ORDE,4   
FITEM,5,50  
FITEM,5,54  
FITEM,5,64  
FITEM,5,74  
ASEL,S, , ,P51X
NSEL,ALL
NSLA,S,1
CM,TN1,NODE
ALLSEL,ALL   
!定义目标单元节点组1
FLST,5,4,5,ORDE,4   
FITEM,5,70  
FITEM,5,80  
FITEM,5,53  
FITEM,5,60  
ASEL,S, , ,P51X
NSEL,ALL
NSLA,S,1
CM,TN2,NODE
ALLSEL,ALL  
!通过接触向导建立接触对
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.25
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3  
ET,4,170
ET,5,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,5,4,0
KEYOPT,5,5,0
NROPT,UNSYM
KEYOPT,5,7,0
KEYOPT,5,8,0
KEYOPT,5,9,0
KEYOPT,5,10,2   
KEYOPT,5,11,0   
KEYOPT,5,12,0   
KEYOPT,5,2,0
KEYOPT,4,5,0
! Generate the target surface   
NSEL,S,,,TN1
CM,_TARGET,NODE
TYPE,4  
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,S,,,CN1
CM,_CONTACT,NODE
TYPE,5  
ESLN,S,0
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  
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.25
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,4
REAL,4  
ET,6,170
ET,7,174
R,4,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,7,4,0
KEYOPT,7,5,0
NROPT,UNSYM
KEYOPT,7,7,0
KEYOPT,7,8,0
KEYOPT,7,9,0
KEYOPT,7,10,2   
KEYOPT,7,11,0   
KEYOPT,7,12,0   
KEYOPT,7,2,0
KEYOPT,6,5,0
! Generate the target surface   
NSEL,S,,,TN2
CM,_TARGET,NODE
TYPE,6  
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
NSEL,S,,,CN2
CM,_CONTACT,NODE
TYPE,7  
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,4  
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  




/SOL
CSYS,1           !激活总体柱坐标系
NSEL,S,LOC,X,R1  !选择法兰盘内空上的所有节点
/PREP7
NROTAT,ALL       !把选中的节点转到总体柱坐标系中
/SOL
D,ALL,,,,,,UX,,,,,  !约束法兰盘内径的径向位移
D,ALL,,,,,,UY,,,,,  !约束法兰盘内径的周向位移
ALLSEL,ALL
CSYS,0
!施加预拉伸力
*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,1   
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 1000,1,2
*ENDDO  

*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,2   
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 1000,1,2
*ENDDO  

*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,3   
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 1000,1,2
*ENDDO

*DEL,_z
*DIM,_z,ARRAY,1,1,1
*SET,_z(1) ,4   
*DO,_z1,1,1
SLOAD,_z(_z1),9,LOCK,FORC, 1000,1,2
*ENDDO  

SOLVE
学无止境
2楼2013-10-25 10:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

错误 1:
    一个节点只能在一个耦合方程中做为从节点。因为节点耦合的意思是使从节点的位移自由度与主节点的自由度相对,如果一个节点同时存在于两个耦合中,那么它到底和哪一个主节点的自由度相等呢?
错误 2:
    一个节点的 Uy 自由度既出现在某个约束方程中,又是一个节点耦合中的从节点,处理时会造成混乱。软件建议将节点耦合也修改为约束方程。
3楼2013-10-26 08:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

starkyu

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by htbbzzg at 2013-10-26 08:24:55
错误 1:
    一个节点只能在一个耦合方程中做为从节点。因为节点耦合的意思是使从节点的位移自由度与主节点的自由度相对,如果一个节点同时存在于两个耦合中,那么它到底和哪一个主节点的自由度相等呢?
错误 2: ...

*** WARNING ***   SUPPRESSED MESSAGE    CP =    2159.818   TIME= 17:16:19
CE 8194 has undefined node 144434                                       
  Equation not generated.      

*** WARNING ***   SUPPRESSED MESSAGE    CP =    2160.130   TIME= 17:16:19
Constraint equation 8194 uses only undefined nodes, it has been         
deleted.

这种又是什么原因呢~
4楼2013-10-27 09:16:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

估计是你在求解前做过实体选择,因而求解时没有包括所有实体,即上述信息中列出的节点不在当前选择集中。
在求解前执行一下命令:
    allsel,all
试试。
5楼2013-10-28 17:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

starkyu

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by htbbzzg at 2013-10-28 17:38:08
估计是你在求解前做过实体选择,因而求解时没有包括所有实体,即上述信息中列出的节点不在当前选择集中。
在求解前执行一下命令:
    allsel,all
试试。

试过了还是不行
就是我的模型里设置了四对接触对 然后计算收敛失败之后我再检查的时候发现有八对接触了 就是重复了但是我也没有重新定义啊 不知道是什么情况。。
6楼2013-10-29 12:46:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

程序不会自动给你添加接触的,应该还是你自己重复定义了。
可以把多出来的接触对删除掉,然后再试试。
可能你的模型中还有其它问题,需要仔细检查一下。
7楼2013-10-29 17:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

starkyu

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by htbbzzg at 2013-10-29 17:50:33
程序不会自动给你添加接触的,应该还是你自己重复定义了。
可以把多出来的接触对删除掉,然后再试试。
可能你的模型中还有其它问题,需要仔细检查一下。

我用GUI和命令流都单独的定义过 计算热分析的时候挺正常的 进行结构分析的时候就会出现这种情况
8楼2013-10-29 22:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

最后一个办法,把你的模型发给我看看。只是我很忙,可能回答要慢一点。163 邮箱,用户: htbbzzg。
9楼2013-10-30 07:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 starkyu 的主题更新
信息提示
请填处理意见