(完整版)高斯热源

Introduction

在一平面上加载一个高斯分布的温度,和高斯热源的原理都是一样的,你只要修改一下就可以了.

1.define a areas

Preprocessor > Modeling > Create >Areas >Rectangle>By Dimensions

2.Define the Type of Element

Preprocessor > Element Type > Add/Edit/Delete...

3.Define Element Material Properties

Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic In the window that appears, enter the following geometric properties for steel:

KXX: 60.5

这个其实是随便给的

4.mesh

将smart size调到fine

点击mesh

然后出现选择对话框

5.定义载荷

Parameters>function>define/edit

其中{X},{Y}从TIME的那个下拉菜单中选取.编辑好公式后file>save

这只是定义了一个函数,还要把它加载

Parameters>function>read from file 然后选择你刚才存储的那个文件

会出现如下对话框

在Table parameter name 中输入表单的名字(随便起好了,只要不合ansys默认的冲突就

6.加载

Solution > Define Loads > Apply > Themal > On Area

然后选择pick all

出现如下对话框

在apply as 中选择Existing table

Ok后会出现下面的对话框

选择你刚才给表单起的那个名字

7.求解

Solution > Solve > Current LS

8.后处理

General Postproc>Plot Results>Contour Plot>Nodal Solu

点ok 可以看到刚开始的那张图

命令流如下

finish

/CLEAR,START

/prep7

ET,1,PLANE55

MPTEMP,,,,,,,,

MPTEMP,1,0

MPDATA,KXX,1,,0.07

BLC5,0,0,10,10

ESIZE,0.1

AMESH,ALL

*SET,_FNCNAME,'TEMP1'

*DIM,%_FNCNAME%,TABLE,6,13,1

!

! Begin of equation: exp(-1*({X}^2+{{Y}}^2))

*SET,%_FNCNAME%(0,0,1), 0.0, -999

*SET,%_FNCNAME%(2,0,1), 0.0

*SET,%_FNCNAME%(3,0,1), 0.0

*SET,%_FNCNAME%(4,0,1), 0.0

*SET,%_FNCNAME%(5,0,1), 0.0

*SET,%_FNCNAME%(6,0,1), 0.0

*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 0, 0, -3 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1

*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 2, 17, -1 *SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 17, -1 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -3, 1, -4 *SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,12,1), 0.0, -1, 7, 1, -3, 0, 0 *SET,%_FNCNAME%(0,13,1), 0.0, 99, 0, 1, -1, 0, 0 ! End of equation: exp(-1*({X}^2+{{Y}}^2))

!-->

!ASEL,ALL

DA,ALL,ALL, %TEMP1%

/STATUS,SOLU

/SOLU

SOLVE

FINISH

/POST1

/EFACET,1

PLNSOL, TEMP,, 0