哪位大侠帮忙加下冷却的命令流!!!

定义材料属性的APDL命令流程序刚孩口下:
/prep7
mptemp,1,20,250,500,750,1000,1500,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,c, 1,1,460,480,530,675,670,660,780,820
MPDATA,PRXY,1,1,0.28,0.29,0.31,0.35,0.4,0.45,0.48,0.5
mpdata,alpx, 1 , 1,1 .10E-05,1.22E-05,1.39E-05,1.48E-05,1.34E-05,1.33E-05,1.32E-05,1.31E-05
mpdata,ex, 1,1,2.05E+11,1.87E+11,1.50E+11,7.00E+10,2.00E+10,1.90E+10,1.80E+10,1.20E+10
!下面定义材料在各温度下的焙值
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,1533
MPTEMP,3,1590
MPTEMP,4,1670
MPDE,ENTH,1
MPDATA,ENTH,1,,0
MPDATA,ENTH,1,,7.5E+009
MPDATA,ENTH,1,,9.6E+009
MPDATA,ENTH,1,,1.1E+010
!下面在各温度点下的屈服应力和屈服后的弹性模量
TB,BISO,1,6,2,
TBTEMP,20
TBDATA,,2.20E+08,2.05E+10,,,,
TBTEMP,500
TBDATA,,1.30E+08,1.50E+10,,,,
TBTEMP,750
TBDATA,,4.00E+07,7.00E+09,,,,
TBTEMP,1000
TBDATA,,2.50E+07,2.00E+09,,,,
TBTEMP,1700
TBDATA,,1.00E+06,1.80E+09,,,,
TBTEMP,2500
TBDATA,,1.00E+05,1.20E+09
SAVE
Fini
建立几何模型及网格划分的APDL程序如下:
/prep7
/units, si
rb=0.002
length=0.12
width=0.06
height=0.006
et, 1, solid5
et, 2, solid95
block, 0, length, 0, width, 0, height
Block, 0, length, 0, 3*Rb, 0, height
Block,0, length,0,Rb*6, 0, height
Block,0, length,0, Rb*12, 0, height
VOVLAP,ALL
/view,1,1,1,1
!下面定义焊缝区的单元尺寸
Lesize, 4,rb
Lesize, 5,rb
Lesize, 9,rb
Lesize, 14,rb
Lesize, 50,rb
Lesize, 10,rb
Lesize, 52,rb
Lesize, 23,rb
Lesize, 51,rb
Lesize, 9,rb
Lesize, 49,rb
Lesize, 24,rb
!定义过渡区的单元尺寸
Lesize,31,2*rb
Lesize,26,2*rb
!定义远离焊缝区的单元尺寸
Lesize, 43,4*rb
Lesize, 38,4*rb
Lesize, 2,4*rb
Lesize, 7,4*rb
Lesize, 53,4*rb
Lesize, 54,4*rb
Lesize, 55,4*rb
Lesize, 56,4*rb
Lesize,11,1 .5*rb
Lesize,12,1 .5*rb
Lesize,47,1 .5*rb
Lesize,48,1 .5*rb
Save
!下面进行网格划分
type, 1
mshkey, 1
mshape, 0
vmesh, 5
vmesh, 6
type, 2
mshkey,0
mshape,1
vmesh,7
vmesh,8
Save
其APDL命令流的实现如下:
/SOLU
DA, 3, SYMM
/GO
DL,53,,UX,
DL,61,,UX,
DL,57,,UX,
DL,49,,UX,
DL,4,,UZ,
DL,9,,UY,
!开始求解设置
/solu
!时间步长设置
antype,trans,new
nlgeom,on
timint,1,struct
Timint,1,thern
timintm,0,mag
Timint,0, elect
Tintp,0.005,,,-1,0.5,0.2
tref,25
nropt,auto,,
autos,on
outres,,all
v=0.008
tinc=rb/v
!移动热源设置
U=25
I=180
ef=0.75
p=25*180*ef
r=0.006
qmax=p/(3.14159*r**2)
!热源移动开始
max_time=length/rb+1
max_ x= length/rb+1
max_y=width/rb+ 1
*dim,flux2,table,max_ x,max_y,max_ time,x,y,time
*do,k,1,max_ time,1
* do,i,1,max_ x,1
flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*rb
*enddo
* do,i,1,max_ x, 1
*do,j,1 ,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/EXP(3*DISTANCE**2/(R**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
nsel,s,loc,z,height
sf,all,hflux, %flux2%
allsel,all
nsel,s,loc,x,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,x,length
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,width
sf,all,conv,10,20
allsel,all
nsel,s,loc,z,0
sf,all,conv,10,20
allsel,all
time,k*rb/v
nsubst,5
solve
BFDELE,all,HGEN
*enddo
save
程序到这里都能运行正确,但进行后面的冷却后,温度越来越高了,应该是热源没有卸载,哪位大侠帮忙加下。谢谢!!!
time,20
nsubst,10
solve
time,44
nsbust,24
solve
time,644
nsbust,20
solve
time,7244
nsubst,100
solve
save
邀请回答 举报

共1个回答

默认排序
  • yanglin0616
    2012年2月26日
    收藏备用了!贴子很不错
    0评论 采纳答案 举报

  • 哪位大侠帮忙加下冷却的命令流!!!的相关案例教程

    !ANSYS命令流学习笔记6 !热应力分析 !学习重点: !1、 理解热力耦合的直接法和间接法 !间接法:先进行热分析,然后将求得的节点温度作为载荷施加到结构应力分析中。 !直接法:直接采用具有温度和位移自由度的耦合单元,同时得到热分析和结构应力分析的结果。直接法又分弱耦合和强耦合选择强耦合时,形成不对称矩阵,线性系统可以直接求解。选择弱耦合时,对称矩阵,还是把热和结构分别进行求解,并将热结果施加
    点击蓝字 关注我们 好久没有更新算例了,这期更新一个小算例,用强度折减法计算边坡的稳定性,用ANSYS软件计算,ABAQUS的算例我前面写了两篇,感兴趣的可以看看。 边坡的尺寸及参数与之前的ABAQUS的边坡算例相同。模型如下。 接下来讲解一下建模的过程。先定义一下单元的类型。这里采用平面应变单元,PLANE42单元。 !定义单元类型ET,1,PLANE82     !定义边坡围岩单元KEYOPT
    0 引言 最近刚好刷到《战狼1》的短视频,里面有一个经典的子*弹穿墙镜头,想着或许可以使用离散元的方法进行模拟。对于弹道分析、穿墙速度损失等方面应该会有一定的指导意义。案例中使用到了很多技术方法,比如子*弹的形状是用圆锥二次化后生成的,墙体的边界处理等。 1 成样 破碎方面一般只在被击中的砖周边,所以我们只建了三行三列的砖墙,其余的边界用有限元的zone模拟。砖头采用240*115*90mm尺寸,
    极限载荷法是求解结构能承受的极限载荷的一种方法。极限载荷法是一种强度分析方法,注意要与结构稳定性分析区分开。下面给出一个带有缺陷的管道结构极限载荷APDL算例,最后给出几点注意事项,供读者参考。 ! ***************环境设置************************finish/clear/filn,file/title,!基本参数h1=108!外径h2=108-2*8!内径h
    Spring comes 小伙们大家好!之前有一期讲了ABAQUS进行反应谱分析时怎么进行静动叠加,利用工况叠加的原理。这一期主要讲解了ANSYS反应谱分析时怎么进行静动叠加。其实,同样是采用了工况叠加。 我接下来以一根柱子来做这个案例。反应谱采用水工抗震规范[1]的反应谱。柱子的尺寸是1×2×5m,弹性模量为2.1E8Pa,泊松比0.2,密度2400kg/m3。 假设本案例地面最大加速度为 a=