应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 关于使用ansys做胀接接触分析不收敛的问题的求出
81/1返回列表
查看: 1597  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

marsltp

银虫 (初入文坛)

[求助] 关于使用ansys做胀接接触分析不收敛的问题的求出 已有1人参与

FINISH
/CLEAR
/FILENAME,3D
/TITLE.,zhangjie
/PREP7
ET,1,solid185
ET,2,TARGE169
ET,3,CONTA172

MP,EX,1,1.9E5  !换热管性能参数
MP,PRXY,1,0.3
TB,MISO,1
TBPT,,0.001,190
TBPT,,0.008,275
TBPT,,0.027,370
TBPT,,0.054,430
TBPT,,0.080,475
TBPT,,0.107,510
TBPT,,0.134,530
TBPT,,0.157,545
MP,MU,1,0.1

MP,EX,2,2.1E5   !母材性能参数
MP,PRXY,2,0.3
TB,MISO,2
TBPT,,0.0023,483
TBPT,,0.013,490
TBPT,,0.033,570
TBPT,,0.053,630
TBPT,,0.076,670
TBPT,,0.116,715
TBPT,,0.155,730
MP,MU,2,0.1

MP,EX,3,1.9E5   !堆焊层性能参数
MP,PRXY,3,0.3
TB,MISO,3
TBPT,,0.001,190
TBPT,,0.008,275
TBPT,,0.027,370
TBPT,,0.054,430
TBPT,,0.080,475
TBPT,,0.107,510
TBPT,,0.134,530
TBPT,,0.157,545
MP,MU,3,0.1

CYLIND,12,8.04,0,7,330,360,
CYLIND,33,12,0,7,330,360,
CYLIND,12,8.04,7,352,330,360,
CYLIND,33,12,7,352,330,360,
CYLIND,7.94,6.59,0,352,330,360,

VGLU,1,2,3,4,
NUMCMP,ALL

LESIZE,5,,,10,,,,,1
LESIZE,7,,,10,,,,,1
LESIZE,6,,,4,,,,,1
LESIZE,8,,,4,,,,,1
LESIZE,10,,,7,,,,,1
MAT,3
VSWEEP,1

LESIZE,26,,,10,,,,,1
LESIZE,28,,,10,,,,,1
LESIZE,27,,,2,,,,,1
LESIZE,29,,,2,,,,,1
LESIZE,30,,,704,,,,,1
MAT,1
VSWEEP,2

LESIZE,17,,,10,,,,,1
LESIZE,19,,,10,,,,,1
LESIZE,18,,,4,,,,,1
LESIZE,20,,,4,,,,,1
LESIZE,36,,,345,,,,,1
MAT,3
VSWEEP,3

LESIZE,5,,,10,,,,,1
LESIZE,14,,,10,,,,,1
LESIZE,41,,,5,,,,,1
LESIZE,40,,,5,,,,,1
LESIZE,15,,,7,,,,,1
MAT,3
VSWEEP,4

LESIZE,17,,,10,,,,,1
LESIZE,21,,,10,,,,,1
LESIZE,44,,,5,,,,,1
LESIZE,45,,,5,,,,,1
LESIZE,43,,,345,,,,,1
MAT,2
VSWEEP,5


WPSTYLE,,,,,,,,0
NUMMRG,ALL
NUMCMP,ALL


ASEL,S,,,11,
NSLA,R,1
CM,CONTACT_node,NODE
ALLSEL,ALL

ASEL,S,,,4
ASEL,A,,,17
NSLA,R,1
CM,TARGET_node,NODE
ALLSEL,ALL

MAT,2
R,3
REAL,3
ET,2,TARGE170
ET,3,CONTA174
KEYOPT,1,3,1  
KEYOPT,3,5,3 !CONTACT设置icont adjustment为gap/penetration
KEYOPT,3,9,0  !包含初始化穿透
KEYOPT,3,10,1 !
R,3,,,0.1,0.1,-0.1,
RMORE,  
RMORE,,0
RMORE,0
NSEL,S,,,TARGET_node
TYPE,2
ESLN,S,0
ESURF
NSEL,S,,,CONTACT_node
TYPE,3
ESLN,S,0
ESURF
ALLSEL,ALL


NSEL,S,LOC,Z,15,342
ASEL,S,,,12
NSLA,R,1
NSLA,R,1
CM,P_node,NODE
ALLSEL,ALL

/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
WPCSYS,-1,0
WPSTYL,,,,,,,,0
FINISH
/SOLU
ASEL,S,LOC,Y,O
DA,ALL,SYMM
CSYS,1
ASEL,S,LOC,Y,330
DA,ALL,SYMM
CSYS,O
DA,1,UZ
DA,8,UZ
DA,19,UZ
DA,24,UZ

ANTYPE,0
NLGEOM,1
AUTOT，ON
OUTRES,ERASE
OUTRES,ALL,ALL

NSUBST,10,10000,10
PRED,ON,,ON
TIME,1
SFA,16,1,PRES,300
ALLSEL,ALL
LSWRITE,1

TIME,2
!NSUBST,10
AUTOTS,1
SFA,16,1,PRES,120
ALLSEL,ALL
LSWRITE,2

TIME,2
SFA,16,1,PRES,0
ALLSEL,ALL
LSWRITE,3


ALLSEL,ALL
LSSOLVE,1,1,1

SAVE,'ZHANGJIE1','db','F:\ANSYS DATA\'
LSSOLVE,2,1,1
SAVE,'ZHANGJIE2','db','F:\ANSYS DATA\'
LSSOLVE,3,1,1
SAVE,'ZHANGJIE3','db','F:\ANSYS DATA\'
FINISH

!/POST1  
!plnsol,cont,pres,0,1
!esel,s,ename,,172   
!plnsol,cont,pres,0,1
!allsel,below,elem   
!prnsol,cont

以上是我用ansys做管板与换热管胀接的接触分析命令流，调整了好多次接触参数，计算始终都不收敛，请问各位大神该如何解决啊？
非常着急，谢谢了！！
回复此楼

» 猜你喜欢
1楼 2017-02-11 20:41:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
你的第一个 lssolve 求解收敛吗？
  你的第二个和第三个 lssolve 命令是错误的，无法求解。应该分别改为：
    LSSOLVE,1,2,1  和 LSSOLVE,1,3,1
  另外，你的第三个载荷步的 time 应该改为 3
赞一下 回复此楼
2楼2017-02-12 08:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marsltp

银虫 (初入文坛)
第一个都不收敛的，都计算不到后面！

发自小木虫IOS客户端
赞一下 回复此楼
3楼2017-02-12 09:08:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

试了你的命令流，确实第一步就不收敛。对你的命令流提了一点建议，见下面命令流中用 !*** 标注的地方。不过可能还需要更多的改进。
  另外，你的原始模型是轴对称的吗？中为什么在侧面而不是内表面施加压力，你的胀接是如何产生的？对你的问题还不是很清楚，最好对你的具体问题做一个详细的介绍。
  对你的命令流提了一点建议，见下面命令流中用 !*** 标注的地方。不过在了解详细情况后，可能还需要更多的改进。

!  [求助] 关于使用ansys做胀接接触分析不收敛的问题的求出  
!  初始金币 30

FINISH
/CLEAR
/FILENAME,3D, 1
/TITLE.,zhangjie
/PREP7
ET,1,solid185
ET,2,TARGE169
ET,3,CONTA172

MP,EX,1,1.9E5  !换热管性能参数
MP,PRXY,1,0.3
TB,MISO,1
TBPT,,0.001,190
TBPT,,0.008,275
TBPT,,0.027,370
TBPT,,0.054,430
TBPT,,0.080,475
TBPT,,0.107,510
TBPT,,0.134,530
TBPT,,0.157,545
MP,MU,1,0.1

MP,EX,2,2.1E5   !母材性能参数
MP,PRXY,2,0.3
TB,MISO,2
TBPT,,0.0023,483
TBPT,,0.013,490
TBPT,,0.033,570
TBPT,,0.053,630
TBPT,,0.076,670
TBPT,,0.116,715
TBPT,,0.155,730
MP,MU,2,0.1

MP,EX,3,1.9E5   !堆焊层性能参数
MP,PRXY,3,0.3
TB,MISO,3
TBPT,,0.001,190
TBPT,,0.008,275
TBPT,,0.027,370
TBPT,,0.054,430
TBPT,,0.080,475
TBPT,,0.107,510
TBPT,,0.134,530
TBPT,,0.157,545
MP,MU,3,0.1

CYLIND,12,8.04,0,7,330,360,
CYLIND,33,12,0,7,330,360,
CYLIND,12,8.04,7,352,330,360,
CYLIND,33,12,7,352,330,360,
CYLIND,7.94,6.59,0,352,330,360,

VGLU,1,2,3,4,
NUMCMP,ALL

LESIZE,5,,,10,,,,,1
LESIZE,7,,,10,,,,,1
LESIZE,6,,,4,,,,,1
LESIZE,8,,,4,,,,,1
LESIZE,10,,,7,,,,,1
MAT,3
VSWEEP,1

LESIZE,26,,,10,,,,,1
LESIZE,28,,,10,,,,,1
LESIZE,27,,,2,,,,,1
LESIZE,29,,,2,,,,,1
LESIZE,30,,,704,,,,,1
MAT,1
VSWEEP,2

LESIZE,17,,,10,,,,,1
LESIZE,19,,,10,,,,,1
LESIZE,18,,,4,,,,,1
LESIZE,20,,,4,,,,,1
LESIZE,36,,,345,,,,,1
MAT,3
VSWEEP,3

LESIZE,5,,,10,,,,,1
LESIZE,14,,,10,,,,,1
LESIZE,41,,,5,,,,,1
LESIZE,40,,,5,,,,,1
LESIZE,15,,,7,,,,,1
MAT,3
VSWEEP,4

LESIZE,17,,,10,,,,,1
LESIZE,21,,,10,,,,,1
LESIZE,44,,,5,,,,,1
LESIZE,45,,,5,,,,,1
LESIZE,43,,,345,,,,,1
MAT,2
VSWEEP,5


WPSTYLE,,,,,,,,0
NUMMRG,ALL
NUMCMP,ALL


ASEL,S,,,11,
NSLA,R,1
CM,CONTACT_node,NODE
ALLSEL,ALL

ASEL,S,,,4
ASEL,A,,,17
NSLA,R,1
CM,TARGET_node,NODE
ALLSEL,ALL

MAT,2
R,3
REAL,3
ET,2,TARGE170
ET,3,CONTA174
KEYOPT,1,3,1  
KEYOPT,3,5,3 !CONTACT设置icont adjustment为gap/penetration
KEYOPT,3,9,0  !包含初始化穿透
KEYOPT,3,10,1 !
R,3,,,0.1,0.1,-0.1,      !***  可以把 第一个0.1 改为 0.01；第二个 0.1 改为 0.2；去掉 -0.1 试试
RMORE,  
RMORE,,0
RMORE,0
NSEL,S,,,TARGET_node
TYPE,2
ESLN,S,0
ESURF
NSEL,S,,,CONTACT_node
TYPE,3
ESLN,S,0
ESURF
ALLSEL,ALL


NSEL,S,LOC,Z,15,342
ASEL,S,,,12
NSLA,R,1
NSLA,R,1
CM,P_node,NODE
ALLSEL,ALL

/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
WPCSYS,-1,0
WPSTYL,,,,,,,,0
FINISH
!*
/SOLU
ASEL,S,LOC,Y,O   !***  O 应该改为 0
DA,ALL,SYMM
!*
CSYS,1
ASEL,S,LOC,Y,330    !***  在 csys,1 中，y 是圆周角度，这个命令可能不起作用
DA,ALL,SYMM
cm, a-y0, area      !***  增加此命令，可以检查选择的 areas 是否正确
!*
CSYS,O        !***  O 应该改为 0？
DA,1,UZ
DA,8,UZ
DA,19,UZ
DA,24,UZ

ANTYPE,0
NLGEOM,1
AUTOT，ON
OUTRES,ERASE
OUTRES,ALL,ALL

NSUBST,10,10000,10    !***   子步数太少，需要增加，比如   NSUBST,10000,1E8,10
PRED,ON,,ON
TIME,1
SFA,16,1,PRES,300
ALLSEL,ALL
!  LSWRITE,1
SOLVE

TIME,2
!NSUBST,10
AUTOTS,1
SFA,16,1,PRES,120
ALLSEL,ALL
!  LSWRITE,2
SOLVE

TIME,3
SFA,16,1,PRES,0
ALLSEL,ALL
!  LSWRITE,3
SOLVE
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

marsltp
4楼2017-02-13 08:39:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marsltp

银虫 (初入文坛)
送红花一朵
引用回帖:
4楼: Originally posted by shgao20 at 2017-02-13 08:39:12
试了你的命令流，确实第一步就不收敛。对你的命令流提了一点建议，见下面命令流中用 !*** 标注的地方。不过可能还需要更多的改进。
  另外，你的原始模型是轴对称的吗？中为什么在侧面而不是内表面施加压力，你的胀 ...

非常感激！我再试一试！

发自小木虫IOS客户端
赞一下 回复此楼
5楼2017-02-13 16:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

好的，如果仍有问题，可以发邮件给我，因为还有一些问题需要请你仔细介绍和商量。
        htbbzzg@163.com
赞一下(2人) 回复此楼
6楼2017-02-14 08:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marsltp

银虫 (初入文坛)
引用回帖:
6楼: Originally posted by shgao20 at 2017-02-14 08:21:30
好的，如果仍有问题，可以发邮件给我，因为还有一些问题需要请你仔细介绍和商量。
        htbbzzg@163.com

请看邮件！

发自小木虫IOS客户端
回复此楼
7楼2017-02-18 20:55:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

已经回复了你的邮件，边界条件中有一点问题。
赞一下 回复此楼
8楼2017-02-19 08:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 marsltp 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭