24小时热门版块排行榜    

查看: 2779  |  回复: 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的回帖

htbbzzg

铁杆木虫 (著名写手)

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

夕阳雨

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了 +11 wuchongjun 2026-08-20 15/750 2026-08-20 17:28 by wuchongjun
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见