接触单元conta174如何设置?

浏览:108654 回答:2
大家好,我现在做热结构耦合分析,需要创建接触单元170与174。但对174的设置好像一直有问题,求解中弹出一下警告和错误:
ERROR
If one or more parts of the model are held together only by contact
verify that the contact surfaces are closed. You can check contact
status in the SOLUTION module for the converged solutions using
CNCHECK.
WARNING
The unconverged solution (identified as time 6000 substep 999999) is
output for analysis debug purposes. Results should not be used for
any other purpose.
求大家帮我看看如何解决。命令流先谢谢了。大家好,我现在做热结构耦合分析,需要创建接触单元170与174。但对174的设置好像一直有问题,求解中弹出一下警告和错误:
ERROR
If one or more parts of the model are held together only by contact
verify that the contact surfaces are closed. You can check contact
status in the SOLUTION module for the converged solutions using
CNCHECK.
WARNING
The unconverged solution (identified as time 6000 substep 999999) is
output for analysis debug purposes. Results should not be used for
any other purpose.
求大家帮我看看如何解决。先谢谢了。
命令流:
Finish
~PARAIN,'ballscrew','x_t',,SOLIDS,0,0 !导入实体模型(x_t格式)
/NOPR
/GO
/FACET,NORML
/REPLOT
/title,Thermal Analysis of Ballscrew
/prep7
et,1,solid70
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.1e11
mpdata,prxy,1,,0.3
mpdata,dens,1,,7820
mpdata,c,1,,460
mpdata,kxx,1,,36.5
uimp,1,reft,,,20
mpdata,alpx,1,,1.12e-5
esize,0.002
MSHKEY,0
MSHAPE,1,3d
vmesh,all !自由网格划分
/VIEW,1,1
/ANG,1
/REP,FAST
Asel,s,,,44 !螺母内表面
Nsla,s,1
Cm,plate1,node !生成第一个节点组件
Allsel
Asel,s,,,5 !丝杠外表面
Nsla,s,1
Cm,shaft1,node !生成第一个节点组件
Allsel
R,3
Real,3
et,2,170
et,3,174
Nsel,s,,,plate1
Type,2
Esln,s,0
Esurf,all
Allsel
Nsel,s,,,shaft1
Type,3
Esln,s,0
Esurf,all
Allsel
Keyopt,3,9,1 !对接触单元进行设置
Keyopt,3,12,4
Finish
/solu !进入求解器
Antype,trans !定义为瞬态分析
Timint,on !时间积分
Kbc,1
Nlgeom,on !包含大变形效应
Nropt,full !使用完全的NR法
Autots,on
Time,1000 !指定求解时间
Deltim,50,50,100 !定义时间子步
kwpave,44,59 !坐标系偏移(点44和59的中间位置)
wprota,0,0,180 !绕Y轴旋转180°
wpoffs,0,0,0.146 !沿转后的Z轴平移0.146
cswpla,11,1,1,1 !定义局部柱坐标系
csys,11
asel,s,area,,2,10
asel,a,area,,12,13
asel,a,area,,22,32
asel,a,area,,37,43
asel,a,area,,45,50,5
asel,a,area,,48,49
asel,a,area,,54,56
asel,a,area,,60,62
asel,a,area,,66,68
nsla,s,1
sf,all,conv,22.75,20 !在丝杠外表面上施加空气对流载荷
Allsel
asel,s,area,,1,11,10
asel,a,area,,14,21
nsla,s,1
sf,all,conv,532.7,20 ! 在丝杠内表面上施加水的对流载荷
allsel
nsel,s,loc,x,0.015
nsel,r,loc,z,-0.066,-0.036
sfdele,all, conv
allsel
nsel,s,loc,x,0.015
nsel,r,loc,z,1.232,1.262 !删去轴承处空气对流换热系数
sfdele,all,conv
allsel
nsel,s,loc,x,0.015
nsel,r,loc,z,-0.066,-0.036
sf,all,hflux,6518 ! 施加轴承处的热流密度
allsel
nsel,s,loc,x,0.015
nsel,r,loc,z, 1.232,1.262
sf,all,hflux,6518 ! 施加轴承处的热流密度
allsel
nsel,s,loc,x,0.02
nsel,r,loc,z,0.46,0.685
sfdele,all,conv !删去螺母处空气对流换热系数
allsel
nsel,s,loc,x,0.02
nsel,r,loc,z,0.46,0.685
sf,all,hflux,6590 ! 施加螺母处的热流密度
Allsel
Tunif,20 !设定初始所有节点温度
Toffst,273
Outres,,all
solve
/post1
Plnsol,temp,,0,
邀请回答 我来回答

全部回答

(2)
默认 最新
邹正刚
为楼上补充一点:首先可以增加一句:
Keyopt,3,1,1
试试,否则 174 号单元的节点自由度只是 Ux,Uy,Uz,不包含温度 Temp,热分析时不起作用。
而 solid70 单元的节点自由度只有 Temp,二者不匹配。
2010年7月11日
评论 点赞
bme
solid70是三维热单元,而你的问题是要模拟结构热耦合,应该选用solid5耦合单元;看你的做法好像是要先得到瞬态温度场,然后再施加到结构上,间接进行热结构耦合;瞬态热分析需要先设定只有一个子步的时间很小的载荷步,作稳态分析,确定初始条件,然后再删除稳态分析中定义的节点温度,最后进行瞬态求解;另外你的命令流生成接触单元也有问题,Esln,s,0这句话不需要,时间步长可能设置过大,导致接触面在变形过程中闭合。
2010年7月10日
评论 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP