24小时热门版块排行榜

返回列表

中兵大将

新虫 (初入文坛)

应助: 2 (幼儿园)
金币: 81.5
帖子: 32
在线: 15.9小时
虫号: 3515551
注册: 2014-11-03
专业: 制造系统与自动化

[求助] 加载高斯热源的一个分析，总是出错，希望大神指点一下

/clear                !清除数据库
/UNITS,SI             ！使用国际单位
/PREP7                ！进入前处理器
ET,1,solid70
et,2,surf152

MPTEMP,1,20,250,500,750,1000,1500
MPTEMP,,1700,2500             ！定义材料1 性能的温度表
MPDATA,DENS,1,1,7820,7700,7610,7550,7490,7350
MPDATA,DENS,1,,7300,7090  ！定义各温度下的密度
MPDATA,KXX,1,1,50,47,40,27,30,35
MPDATA,KXX,1,,140,142 !定义各温度下的X上热传导系数
MPDATA,c,1,1,460,480,530,675,670,660
MPDATA,c,1,,780,820

block,0/1000,3/1000,0/1000,3/1000,,1/1000       ！建立一个矩形
block,3/1000,10/1000,0/1000,3/1000,,1/1000       ！建立另一个矩形
vglue,all                                  ！把两个面粘结起来
NUMCMP,ALL                               ！压缩所有的号码
/pnum,line,1                               !打开线编号
vplot

!常数定义
PI=3.1415
U=20                               ！电压
I=100                               ！电流
eff=0.18                            ！有效效率
q=eff*U*I                   ！最大的热流量
Ra=2.1/1000 !电弧半径
QMAX=3*Q/(pi*ra**2)
SizeE=0.05/1000 !网格大小

lesize,5,SizeE       ！定义1、3线的划分单元的尺寸
lesize,7,SizeE

lesize,8,,,40,1/2       ！定义4、2、5 的划分个数和相邻单元的比例
lesize,6,,,40,2
lesize,14,,,40,2

lesize,19,,,15,3
lesize,20,,,15,3

lesize,9,,,25,2

asel,s,loc,z,1/1000             !选中实体上表面
AATT,1, , 2,       !指定实体上表面用152号单元
MSHAPE,0,2D    !实体单元类型六面体
MSHKEY,1
amesh,all                !对上表面划分网格

allsel,all
aatt,1,,1                ！给相关的元素分配材料属性
MSHAPE,0,3D    !实体单元类型六面体
MSHKEY,1
vmesh,all
save                ！保存
eplot                ！
finish                ！从当前的处理器中退出

/solu                ！进入求解处理器
ANTYPE,TRANS          ！进入瞬态分析
tunif,25             ！定义统一温度
TOFFST,273          !定义绝对温度与所选温度的偏差

allsel
nsel,s,loc,x,0
dsym,symm,x          !施加关于x轴（yz面）的对称边界条件
nsel,s,loc,y,0
dsym,symm,y

!表面传热状态
allsel

esel,s,type,,2             !选中实体上表面的表面效应单元以方便加对流
sfe,all,,conv,2,10000
allsel

asel,s,loc,z,1/1000             !加载热流，选择高斯热源面
nsla,r,1                         !选择面上所有节点
esln,r                         !选择和这些节点有关的单元
eplot                            !所有单元可视

*do,Xlocatio,0,Ra,SizeE                ！定义do循环
*do,Ylocatio,0,Ra,SizeE                ！定义do循环
  rr=sqrt(xlocatio**2+ylocatio**2)
  *if,rr,le,ra,then
qr=qmax*exp(-3*(rr/ra)**2)          ！定义所加的热量？？？？
nsel,s,loc,y,ylocatio                ！选择Y坐标为8的节点
nsel,r,loc,x,xlocatio                   ！再选择X坐标为XLOC的节点
nsel,r,loc,z,1/1000
f,all,heat,qr                   ！在所选的节点上施加热载荷
  *endif
*enddo !逐点计算并施加热流
*enddo

allsel
NSUBST,8,16,5          ！定义载荷步的子步数量
TIME,1e-5 !时间步长(0.001,0.0005)
KBC,0 !斜坡载荷加载方式
AUTOTS,On !自动时间步长
timint,on !Turns on transient effects.
OUTRES,ALL !All solution items written to the database
OUTPR,ALL,LAST !-- Controls the solution printout
TINTP,,,,1 !Defines transient integration parameters
TRNOPT,FULL !Specifies transient analysis options
EQSLV,ICCG,,2 ！Specifies the type of equation solve save saver
LNSRCH,ON             ！打开线搜索
solve
/POST1       !进入POST1
SET,1,last
etable,tempnode,temp,,          !将总应变存入ETABLE etable,tempreture,bfe,temp
esel,s,etab,tempnode,1350,       !选择所有总应变大于或等于0.20的单元  esel,s,etab,tempreture,,
FINISH
!
/SOLU             !重新进入求解器
ANTYPE,，rest       !在分析，缺省为new
nropt,full       !打开牛顿-
estif,
EKILL,ALL       !杀死选择（超过允许值）的单元
ESEL,s,live

回复此楼