镜片热分析问题请教

cengjing5566

镜片热变形分析 请教

et,1,5
mp,ex,1,6.86e10
mp,prxy,1,0.37
mp,alpx,1,7.7e-6
mp,kxx,1,13

cyl4,,,0.011,,,,0.006
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
vglue,all
numcmp,all
lsel,s,,,13,14
lsel,a,,,18,20,2
lesize,all,,,16
lsel,s,,,5,8,1
lesize,all,,,16
lsel,s,,,9,10
lesize,all,,,20
vmesh,all

csys,1
!use function editor define    function=3*exp(-{x}^2)
!the range of x : 0——0.01
!将定义的函数保存为 I
nsel,s,loc,x,,0.01
bf,all,hgen,%I%
alls
nsel,s,loc,x,0.011
d,all,ux,,,,,uy,uz
alls

/solu
antype,4
trnopt,full

timint,1,struct                 !打开结构分析时间积分选项
timint,1,therm                  !打开热  分析时间积分选项
timint,0,mag                    !关闭磁场分析时间积分选项
timint,0,elect                  !关闭电场分析时间积分选项

neqit,100                       !最大平衡迭代次数
tref,300
kbc,1
time,2
deltim,0.05,0.01,0.1
autots,on
pred,on                         ! 打开预测矫正
lnsrch,on                       ！打开线性搜索
outres,basic,all
solve

说明：一束半径大小为0.01m的平行光从圆柱形镜片（半径0.011m，长0.006m）的一端面射入，从另一端射出，光线透射过程中能量密度保持不变，镜片的最外圆周固定，光的功率密度为函数形式=30*exp(-{x}^2)   （W/m^2)，镜片材料吸收系数设为0.1/m，  施加热载荷时我以生热率Hgen=30*0.1 exp(-{x}^2)  =3exp(-{x}^2) （W/m^3)  的形式加载，但求解不收敛。        请高手赐教！！！

iomega

原帖由 cengjing5566 于 2007-10-29 16:48 发表
镜片热变形分析 请教

et,1,5
mp,ex,1,6.86e10
mp,prxy,1,0.37
mp,alpx,1,7.7e-6
mp,kxx,1,13

cyl4,,,0.011,,,,0.006
wprota,,90
vsbw,all
wprota,,,90
vsbw,all
vglue,all
numcmp,all
lsel,s, ...

At least your heat generation function format is not correct.  It should has exp decay as a function of z, as well as x (radius)

Check this link:

http://www.simwe.com/forum/viewthread.php?tid=586443&highlight=%2Biomega