求助一个简单的热分析的解答！

求助 ！光垂直于圆筒中心照在外表面上一段时间的圆筒温度场分布。一个内径0.51m，外径0.52m，长1m的空心铁圆筒，内外为空气，初始温度均为20度，材料的热传导系数45，热容为675，密度为7790，外表面的对流系数为35，不考虑内表面的对流，内表面绝热，不考虑热辐射和相变，存在（载荷，热流密度）外表面中心，也就是在柱坐标下以(z=0, Ө=0,R=0.51m)为中心的R方向有一个垂直于R分布5431exp[-2r/(0.17m*0.17m)]（0<r<0.17）的光斑，（注意光斑照在圆筒外壁上分布要发生变化，5431exp[-2r/(0.17m*0.17m)]*cos Ө）求解加热10秒的温度场分布。希望高手能给出命令流或者GUI。