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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » ANSYS & ABAQUS » 求教一个ansys滑动摩擦的问题，采用弧长法，限制Y向位移。
51/1返回列表
查看: 1207  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lolo123

木虫 (小有名气)

[求助] 求教一个ansys滑动摩擦的问题，采用弧长法，限制Y向位移。已有1人参与

总是计算到1.33e-9秒就终止，根本达不到设置的限定Y向位移0.5mm

finish
/clear
/PREP7
ET,1,82,,,3
R,1,1,  
!***********************单位制(N,mm,Mpa)

PCIRC,800,800+120,0,-180,
wpoff,,-756
csys,4
wprot,,90
asbw,all
wpoff,,,16
asbw,all
kwpave,12
wprot,,,90
asbw,all
kwpave,8
asbw,all
kwpave,6
asbw,all
kwpave,10
asbw,all
adele,1,,,1
adele,4,,,1
adele,6,,,1
wpoff,-3
k,
kwpave,12
wpoff,-3
k,
l,19,20
k,,13
kwpave,19
k,,13
l,21,20
l,19,22
kwpave,21
k,,,,-50
kwpave,22
k,,,,50
l,22,24
l,21,23
kwpave,23
k,,100
kwpave,24
k,,100
l,24,26
l,23,25
l,25,26
al,3,15,18,19,28,30,31,32
wpcsys,-1,0
kwpave,14
k,600,1.856358392
kwpave,23
k,700
kwpave,21
k,800
kwpave,12
k,900

a,600,700,800,900
wpcsys,-1,0
ARSYM,X,4, , , ,0,0  

kwpave,14
wpoff,,-97
wprot,,90
ASBW,ALL
adele,12,,,1
adele,13,,,1
adele,17,,,1
adele,19,,,1
kwpave,600
wprot,,,90
lsbw,20
kwpave,27
lsbw,16


/prep7
mp,ex,1,1.78e5
mp,nuxy,1,0.29
MP,ALPX,1,13.58e-6
MP,DENS,1,7800/1000000000

mp,ex,2,9500e5
mp,nuxy,2,0.25
MP,ALPX,1,5.5e-6
MP,DENS,2,1550/1000000000
!!!!!!!!!!!!mesh
alls
asel,,,,4
asel,a,,,6
AATT,2, 1, 1, 0,    !(MAT,REAL,TYPE,ESYS,SECN)
!!!!defining mat
alls
asel,u,,,4
asel,u,,,6
AATT,1, 1, 1, 0,    !(MAT,REAL,TYPE,ESYS,SECN)
alls
AESIZE,all,8,  
amesh,all

!********
MP,MU,1,0.1
MAT,1   
MP,EMIS,1,7.88860905221e-031
MP,MU,2,0.1
MAT,2   
MP,EMIS,2,7.88860905221e-031
R,3
REAL,3  
ET,2,169
ET,3,172
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
RMORE,10.0  
KEYOPT,3,3,0
KEYOPT,3,4,2
KEYOPT,3,5,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,2   
KEYOPT,3,11,0   
KEYOPT,3,12,0   
KEYOPT,3,2,4
! Generate the target surface   
LSEL,S,,,12
lsel,a,,,20
CM,_TARGET,LINE
TYPE,2  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
LSEL,S,,,36
lsel,a,,,40
CM,_CONTACT,LINE
TYPE,3  
NSLL,S,1
ESLN,S,0
ESURF   
ALLSEL  


!********
MP,MU,1,0.1
MAT,1  
MP,EMIS,1,7.88860905221e-031
MP,MU,2,0.1
MAT,2   
MP,EMIS,2,7.88860905221e-031
R,4
REAL,4  
ET,4,169
ET,5,172
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
RMORE,2.0   
KEYOPT,5,3,0
KEYOPT,5,4,2
KEYOPT,5,5,0
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,4
! Generate the target surface   
LSEL,S,,,19
lsel,a,,,28
CM,_TARGET,LINE
TYPE,4  
NSLL,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
LSEL,S,,,34
lsel,a,,,38
CM,_CONTACT,LINE
TYPE,5  
NSLL,S,1
ESLN,S,0
ESURF   
ALLSEL  
FINISH  


/SOL
alls
LSCLEAR,ALL
OUTRES,ALL,ALL,

SFL,              6         ,PRES,21,  
SFL,              9         ,PRES,21,  
SFL,             11         ,PRES,21,  
SFL,             13         ,PRES,21,  
SFL,             17         ,PRES,21,  
SFL,             22         ,PRES,21,  
SFL,             30         ,PRES,21,  
SFL,             31         ,PRES,21,  
SFL,             32         ,PRES,21,  
SFL,             33         ,PRES,21,  
SFL,             37         ,PRES,21,  

DL,2, ,ALL,  
DL,4, ,ALL,  
DL,33, ,UX,   
DL,37, ,UX,   
D,217, , , , , ,UX, , , , ,   

antype,0                  
nlgeom,0               
outres,all,all           
arclen,1,0               
arctrm,u,0.5,56,uy     
TIME,1  
AUTOTS,1
nsubst,200,,,1
OUTRES,ALL,ALL,
alls
SOLVE

求教一个ansys滑动摩擦的问题，采用弧长法，限制Y向位移。
assss.jpg
回复此楼

» 猜你喜欢
1楼2017-11-29 15:20:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lolo123(月只蓝代发): 金币+20, 感谢应助！ 2017-12-11 13:14:50
看来你的材料性能有问题，包括具体数据和单位两个方面。其中：
  材料 1 似乎是钢材
    mp,ex,1,1.78e5
    mp,nuxy,1,0.29
    MP,ALPX,1,13.58e-6
    MP,DENS,1,7800/1000000000    !  7.8e-6  什么单位？  
  材料 2 不知是什么材料？
    mp,ex,2,9500e5               !  此弹性模量是 材料 1 的 5000 多倍，什么材料有这么大的弹性模量？
    mp,nuxy,2,0.25
    MP,ALPX,1,5.5e-6
    MP,DENS,2,1550/1000000000    !  1.55e-6  什么单位？
  你说你的单位制是 N,mm,MPa，这似乎不是一个协调的单位制。一般将单位制，可以是：
      Kg, m, s 单位制 - 长度单位 m。
    如果长度单位只要 mm，相应的常用单位制是：
      T, mm, s    (T - 吨)
    这样，密度单位应该是   T/mm^3，你的密度值应该分别改为：
      材料 1 - 7.8e-9；    材料 2 - 1.55e-9。

  另外，你的问题是一个接触分析问题，创建了接触对，但是没有设置合适的接触属性。而且，接触问题不应该采用弧长法进行分析，使用一般的静力计算方法就可以。但是你的子步数设置太小了，建议改为：
    NSUBST,200,  1, 10
  先修改材料参数，其它问题慢慢解决吧。
赞一下(1人) 回复此楼

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

lolo123
2楼2017-11-30 08:39:46
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lolo123

木虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by shgao20 at 2017-11-30 08:39:46
看来你的材料性能有问题，包括具体数据和单位两个方面。其中：
  材料 1 似乎是钢材
    mp,ex,1,1.78e5
    mp,nuxy,1,0.29
    MP,ALPX,1,13.58e-6
    MP,DENS,1,7800/1000000000    !  7.8e-6  什么单位？ ...

材料2确实是很刚，那是个石墨密封，刚度很大。至于说没设置好合适的接触属性，我觉得应该没有问题，因为之前我用这个设置计算过一个滑移密封的问题，结果是正确的。单位制问题，可能是我行业的问题，我们比较习惯这种单位制。至于弧长法计算接触问题不容易收敛，这个也是没有办法才试着使用的，可能是限于我所知，在ansys里还找不到能做到滑移摩擦条件终止的计算方法。
其实归根到底还是找不到合适的条件终止计算的方法。
不知说的对不对，还请这位大虫指正。
赞一下 回复此楼
3楼2017-11-30 14:39:12
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lolo123(月只蓝代发): 金币+10, 感谢应助！ 2017-12-11 13:15:02
石墨的弹性模量，已知的最大值是 15000 MPa；即使是金刚石，其弹性模量 (可能是目前弹性模量最大的材料) 也只有  11e5 MPa，都远小于你的材料 2 的数据。
  至于接触属性，你的命令流中根本没有设置过，而是采用 ANSYS 的默认设置。不过，默认的接触设置对某些问题是有效的，但不是对任何问题都有效。
  弧长法用于接触分析是不合适的。对接触问题，一般可以按照静力分析处理，设置合适的接触属性和求解参数即可。当然设置也是有一些技巧的，需要多次调试以找到对你的问题最合适的设置。
赞一下 回复此楼
4楼2017-12-01 08:05:06
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lolo123

木虫 (小有名气)
引用回帖:
4楼: Originally posted by shgao20 at 2017-12-01 08:05:06
石墨的弹性模量，已知的最大值是 15000 MPa；即使是金刚石，其弹性模量 (可能是目前弹性模量最大的材料) 也只有  11e5 MPa，都远小于你的材料 2 的数据。
  至于接触属性，你的命令流中根本没有设置过，而是采用 A ...

Y向位移作为一个限制条件，应该咋设置
赞一下 回复此楼
5楼2017-12-04 09:04:17
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lolo123 的主题更新
51/1返回列表
信息提示
关闭
请填处理意见
关闭