在一个面上加载高斯分布的热流密度

zcy2009

在做一个在面上加载高斯分布的热流密度，激光作用点应该是正方形的顶点处，不是中点处，可是用下面的程序计算出来以后得到的温度场却是一个类似圆环状，和预测的结果相差很大，搞了两天了都不知道原因，希望版上的高手能给予指导，谢谢！
/CLEAR,NOSTART
/FILNAM,temperature field !定义工作文件名
/TITLE,RESIDUAL STRESS ANALYSIS !定义工作标题
/UNITS,cgs!确定单位制

/PREP7 !进入前处理器
TOFFST,273
ET,1,plane55 !定义单元类型


MP,DENS,1,1.093 !输入密度
MP,C,1,3.37 !输入比热
MP,KXX,1,4.5e-3 !输入材料导热系数

/VIEW,1,1,1,1
RECTNG,0,0.2,0,0.2  !建立二维模型
SAVE
LESIZE,1,0.01 !设置线上单元的大小  
LESIZE,2,0.01
MSHKEY,1       !选择分网方式KEY=0（自由分网）1（映射分网）
               !2（自主选择）   
AATT,1,,1,0 !制定面的单元属性
AMESH,1        !划分单元
EPLOT          !元素显示

/solu
antype,trans,new!选择瞬态分析
timint,on !打开瞬态效应
autots,on !选择自动时间步进
outres,,all!输出每一步计算结果
lsize=0.01 !单元长度
t=600 !输入计算时间
tinc=1 !输入单步时间长度
n=t/tinc !需要计算的步数
p=12!输入激光的功率
w0=0.05!输入激光光束原始功率
q0=2*p/(3.14*w0*w0)

max_x=21 !x方向的节点数
max_y=21 !y方向的节点数

TUNIF,37 !初始温度
allsel,all
nsel,s,loc,y,0
sf,all,conv,6.2e-4,25
allsel,all




*DIM,flux2,TABLE,max_x,max_y,n !定义储存激光作用热流密度的数组



n=4
*DO,k,1,n,1

  *DO,i,1,max_x,1
  flux2(i,0,K)=(i-1)*lsize
    *ENDDO
   
  *DO,j,1,max_y,1
  flux2(0,j,K)=(j-1)*lsize
    *ENDDO
   
  
   *DO,i,1,max_x,1
    *DO,j,1,max_y,1
     xi=(i-1)*lsize
     yi=(j-1)*lsize
             mmm=sqrt(xi*xi+yi*yi)
flux2(i,j,k)=exp(-mmm)
*ENDDO
   *ENDDO
  
  flux2(0,0,k)=(k-1)*tinc
t=(k-1)*tinc+1e-6

time,t
autots,0
deltim,,,,1
kbc,0
allsel,all
nsel,r,loc,z,0

Sf,ALL,HFLUX,%flux2%

allsel,all
solve  
*ENDDO

zcy2009

刚刚试着用公式编辑器加载的方法得出了我想要的结果，在此发帖攒个人品 ！发现还是用公式编辑器好用，直接定义数组的方法加载真的很难控制计算机内部的运行情况。

zhibijingtian

2# zcy2009
很想知道 怎么用公式编辑器加载？？？ ANSYS运行命令跟公式编辑器有什么关系？

yanyankankan

好贴，我也在研究这方面的东西

东西南北

顶一下

无宇伦比

顶一下，

parapara

我用公式编辑器加载总是不成功