24小时热门版块排行榜

返回列表

shuoqi

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 109
帖子: 107
在线: 25.5小时
虫号: 780986
注册: 2009-05-27
性别: MM
专业: 零件成形制造

[求助] 200金币求助，为什么我的生死单元没起作用？已有1人参与

目的：

类似焊接，喷涂多层材料到一基材表面，每喷涂一层后，停2分钟左右，基材表面的热量向空气辐射；再向表面喷涂第二层，再停2分钟向空气辐射，如此直到喷涂4层为止，观察基材的温度场分布及热应力分布，但我用的生死单元不知为什么，没起作用，请哪位大侠帮忙看下，解决问题的话，奉上多半家产，200金币。

Keyw,thermal,1 !过滤图形用户界面
Keyw,structural,1
/prep7
!基材参数
ET, 1, PLANE55 !T选择热分析单元
!ET,2,LINK34 !选择用于2D的对流线单元
MP,KXX,1, 0.558 !定义热传导系数 5 W/mC
!注意此处是一个“，”设置，如果是Mp,c,1,,200，则错误，是指随温度变化的性能情况，最终结果错误，只能出现一个，无法出现瞬态各个时间段的结果。
Mp,c, 1,2550
!注意此处是一个“，”设置，如果是Mp,dens,1,,5000，则错误，最终结果只能出现一个，无法出现瞬态各个时间段的结果。
Mp, dens, 1, 900
Mp,alpx,1,6.43e-6
Mp,ex,1,1.57e7
Mp,prxy,1,0.3
!涂层性能参数
!ET,2,PLANE55
Mp,kxx,2,131
Mp,c,2,475
Mp,dens,2,6610
Mp,alpx,2,2.858e-5 !热膨胀系数
Mp,ex,2,1.7e10
Mp,prxy,2,0.3

!Et,2,surf152    !152适用于三维体单元分析时的空间点
Et,2,surf151    !151适用于二维面单元分析时的空间点

Keyopt,2,4,1 !定义单元关键字,此处说明surf152单元没有中间节点.2表示单元2，4表示
!单元第4个选项，1表示单元没有中间节点。
Keyopt,2,5,1 !surf152单元有额外脱离基体的节点
Keyopt,2,9,1 !采用辐射形状因子实常数辐射
r,2,1,5.67e-8 !定义单元实常数辐射形状因子

!两种材料空气吸收辐射率
mp,emis,3,0.8 !输入辐射率

!建立分析模型
K,1,0,0
K,2,200e-6,0
K,3,200e-6,200e-6
K,4,0,200e-6
K,5,200e-6,240e-6
K,6,0,240e-6
!K,7,200e-6,280e-6
!K,8,0,280e-6

A,1,2,3,4
A,3,5,6,4
!A,5,7,8,6
Esize,2e-6

Mat,1
Asel,s,,,1
Amesh,all

Allsel
Mat,2
!Asel,s,,,2,3
Asel,s,,,2
Amesh,all

!生成表面单元
Type,2
Real,2
Mat,3 !指定材料类型
N,60000,100e-6,250e-6,0    !生成节点N, NODE, X, Y, Z, THXY, THYZ, THZX
!若是处节点设置坐标不同时，对最终结果有较大影响，应该是没有将生成体的面单元清除干净。
!Space,60000 !在此处加这一个命令对最终结果影响不大。
Esurf,60000 !生成表面单元，当额外节点时，可以以esurf定义

/solu
Antype,trans
TRNOPT,FULL
Nlgeom,on  !激活大变形选项
Nropt,full  !设定牛顿-拉夫森选项
Outres,all,all

ic,all,temp,25
!Asel,s,,,1
!Nsla,s,1
!D,all,temp,25

allsel
!asel,s,,,2,3
asel,s,,,2
nsla,s,1
d,all,temp,453

!*
Allsel
asel,s,,,3
nsla,s,1
d,all,temp,20

Asel,s,,,3
Esla,s
Ekill,all
!*

TIME,0.01
!Time,0.4e-3
Autots,on
Deltim,0.02,0.01,0.05
TIMINT,on
allsel
!Esel,s,live
!Nsle,s
!NSEL,INVE                         !选择所有非活动结点（不与活动单元相连的结点）
!D,ALL,temp,0                      !约束所有不活动的结点自由度（可选）
!ncnv, 2  ! 终止分析选项
Solve

allsel
asel,s,,,2
nsla,s,1
Ddele,all,temp
allsel
Toffst,273    !定义温度偏移量
D,60000,temp,25  !设定周围介质温度
!Time,60
TIME,60

Autots,on
!Deltim,2,1,5
Deltim,2,1,5
TIMINT,on
Allsel
Solve

Allsel
Asel,s,,,1,2
Nsla,s,1
Ddele,all,temp
Allsel
Lsel,s,,,6
Nsll,s,1
Sf,all,conv,8,25 ! Specifies surface loads on nodes

allsel
TIME,10
Autots,on
Deltim,3,1,5
TIMINT,on
Allsel
Solve

Allsel
Lsel,s,,,6
Nsll,s,1
Sfdele,all,conv
allsel
Asel,s,,,3
Esla,s
Ealive,all
Allsel
Asel,s,,,3
Nsla,s,1
Ddele,all,temp
D,all,temp,100
Allsel

Time,6.01
Autots,on
!Deltim,0.2,0.1,0.5
Timint,on
Solve

allsel
Asel,s,,,3
Nsla,s,1
Ddele,all,temp
Lsel,s,,,9
Nsll,s,1
Sf,all,conv,8,20
allsel
Time,12
Autots,on
Deltim,0.2,0.1,0.5
Timint,on

ncnv, 2  ! Terminate the analysis, but not the program execution, if the solution fails to converge
Solve

回复此楼