lchylinchy 发表于 2012-4-6 17:32:28

求助脉冲激光热源加载问题

因课题需要,我是最近才接触ANSYS的,我要做脉冲激光加热物体,温度\热应力场的变化情况的模拟.要将脉冲激光热源加载到一个面上,脉冲激光的的加载方式为:
1、激光持续时间为1s
2、激光关闭时间为10s
3、激光再次打开
4、依次重复上述过程
不知该怎么在ANSYS中实现这一过程,有哪位可以帮我的?
真的好急啊!!!!
程序如下:请各位大侠看看有什么错误没有啊?
fini
/clear
/filename,danliang
/title,thermal analysis of beam
/prep7
et,1,solid226,11
et,2,solid226,11
et,3,surf152         !表面效应单元
KEYOPT,3,3,0         ! OMEG 用于全局笛卡儿坐标系X轴
KEYOPT,3,4,0         !有中间节点
KEYOPT,3,5,0         !关于额外节点的辐射或对流计算: 无额外节点
!KEYOPT3,6,0          !使用体积温度,额外节点温度作为体积温度
KEYOPT,3,7,0         !对于膜层散热系数不增加经验项
KEYOPT,3,8,3         !根据表面温度计算对流系数
KEYOPT,3,9,0         !不包括热辐射

block,0,6e-3,-0.75e-3,0.75e-3,-0.015e-3,0.015e-3
block,0,6e-3,-0.75e-3,0.75e-3,0.015e-3,0.016e-3

mp,kxx,1,140         !硅片的导热系数
mp,dens,1,2230         !硅的密度
mp,ex,1,1.3e+11      !硅的杨氏模量
mp,nuxy,1,0.28         !硅的泊松比
mp,alpx,1,2.8e-6       !硅的热膨胀系数

mp,kxx,2,237         !AL的导热系数
mp,dens,2,2690         !AL的密度
mp,ex,2,7.17e+11       !AL的杨氏模量
mp,nuxy,2,0.33         !AL的泊松比
mp,alpx,2,2.5e-5       !AL的热膨胀系数
mp,c,2,880             !AL的比热
VGLUE,ALL            !将AL镀在硅上
numcmp,all             !压缩线、面、体等   
!划分梁体,映射划分
lesize,2,,,30
lesize,3,,,15
lesize,11,,,2
mshape,0,3d
mshkey,1
vsel,s,,,1
vatt,1,,1
vmesh,all

lesize,19,,,1
vsel,s,,,2
mshape,0,3d
mshkey,1   !指定自由划分或映射划分(0为自由,1为映射,2先选自由,否则有映!射)
type,2
mat,2
vmesh,all
!生成表面效应单元
asel,s,,,2
nsla,s,1
type,3
esurf,all
nsel,s,loc,x,0
d,all,all,0                  !固定左端位置
            
fini
/solu
antype,trans                  !设置为瞬态求解
tref,25                     !定义初始温度
!施加对流载荷
定义固定侧面对流15000
nsel,s,loc,x,0
nsel,r,loc,y,-0.75e-3,0.75e-3
nsel,r,loc,z,-0.01e-3,0.016e-3
sf,all,conv,15000,25         

!定义右端端面对流150
nsel,s,loc,x,6e-3
nsel,r,loc,y,-0.75e-3,0.75e-3
nsel,r,loc,z,-0.015e-3,0.016e-3   
sf,all,conv,150,25      
!定义水平面对流150
nsel,s,loc,x,0,6e-3
nsel,r,loc,y,0.75e-3
nsel,r,loc,z,-0.015e-3,0.016e-3   
sf,all,conv,150,25         
nsel,s,loc,x,0,6e-3
nsel,r,loc,y,-0.75e-3
nsel,r,loc,z,-0.015e-3,0.016e-3   
sf,all,conv,150,25

!定义前后端面对流50
nsel,s,loc,x,0,6e-3
nsel,r,loc,y,-0.75e-3,0.75e-3
nsel,r,loc,z,-0.015e-3   
sf,all,conv,50,25         
nsel,s,loc,x,0,6e-3
nsel,r,loc,y,-0.75e-3,0.75e-3
nsel,r,loc,z,0.016e-3   
sf,all,conv,50,25
alls

kbc,1                         !定义为阶跃载荷
time,1                  !激光脉冲宽度(时间从0到1s)
nsubst,10                      !指定计算子步数

!施加一个直径为1.2mm的激光均匀热流密度
esel,s,type,,1
nsle,s
nsel,s,loc,x,2.4e-3,3.6e-3
nsel,r,loc,y,-0.6e-3,0.6e-3
nsel,r,loc,z,-0.015e-3
sf,all,hflux, 75            !施加均匀热流密度
!outres,,1                     !存储求解结果选项设置
lswrite,1                     !写入到载荷步文件1

time,10                  !无激光宽度(时间从1s到10s)
nsubst,10                     !指定计算子步数
!deltim,1e-6,0.1e-6,2e-6,1    !定义载荷步
!施加一个直径为150μm的激光均匀热流密度
esel,s,type,,1
nsle,s
nsel,s,loc,x,2.4e-3,3.6e-3
nsel,r,loc,y,-0.6e-3,0.6e-3
nsel,r,loc,z,-0.015e-3
sf,all,hflux, 0               !施加均匀热流密度
!outres,,1                     !存储求解结果选项设置
lswrite,2                     !写入到载荷步文件2

time,11                  !激光脉冲宽度(时间从10s到11s)
nsubst,10                      !指定计算子步数
!deltim,1e-6,0.1e-6,2e-6,1    !定义载荷步
!施加一个直径为150μm的激光均匀热流密度
esel,s,type,,1
nsle,s
nsel,s,loc,x,2.4e-3,3.6e-3
nsel,r,loc,y,-0.6e-3,0.6e-3
nsel,r,loc,z,-0.015e-3
sf,all,hflux, 75             !施加均匀热流密度
!outres,,1                     !存储求解结果选项设置
lswrite,3                     !写入到载荷步文件3

time,21                  !无激光宽度(时间从10s到21s)
nsubst,10                      !指定计算子步数
esel,s,type,,1
nsle,s
nsel,s,loc,x,2.4e-3,3.6e-3
nsel,r,loc,y,-0.6e-3,0.6e-3
nsel,r,loc,z,-0.015e-3
sf,all,hflux, 0               !施加均匀热流密度
!outres,,1                     !存储求解结果选项设置
lswrite,4                     !写入到载荷步文件4
alls
lssolve,1,4,1                !按载荷步文件进行求解

fini

/post26
cntr_c=node(3e-3,0,-0.015e-3)!获取位置为(5e-3,0,-0.015e-3)的节点编号,并赋给cntr_c
nsol,2,cntr_c,temp             !将cntr节点的温度变化值设为2号变量
nsol,3,cntr_c,u,x            !将cntr节点的位移变化值设为3号变量
plvar,2                        !绘制温度变化曲线
plvar,3                        !绘制位移变化曲线
prvar,2                        !列表显示变量
prvar,3


/image,save,wendu,bmp
/image,save,weiyi,bmp
fini

页: [1]
查看完整版本: 求助脉冲激光热源加载问题