| 查看: 714 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者Su peng将赠送您 25 个金币 | ||
[求助]
ANSYS 高斯函数型热源+生死单元求助已有1人参与
|
||
 求助论坛大神，小弟刚入门apdl来模拟激光熔覆温度场，想用高斯移动热源+生死单元模拟同步送粉，根据看的一些资料 利用apdl里面加载函数的形式，加载了高斯热源函数 ，函数是这样的2.866242e7*exp(-2*(({X})^2+({Z}+0.003*{TIME})^2)/0.004^2)，其中0.003是热源移动速度，但是编写出来进行加载，移动热源只在我激活的第一个单元加载，其他单元激活不了，而且热源也不移动，不知道程序出了什么问题。自己实在找不到问题在哪，恳请论坛大神出手相助，帮忙看一下问题出在哪，小弟不胜感激。APDL代码如下: /CLEAR /units,SI /PREP7 mptemp,1,0,400,600,800,1000,1500 mpdata,kxx,1,1,18.9,28.35,31.08,32.76,31.5,28.2 mpdata,c,1,1,580,756,840,882,756,797 mpdata,hf,1,1,9.66,59.43,84.32,109.2,180.6,359 mpdata,dens,1,1,7800,7550,7550,7550,7550,7550 *dim,conve,table,8,1,1,temp !定义了一个表格 8行1列一面的网格 对流和辐射总换热系数，按照对流加载 conve(0,1,1)=1 conve(1,0,1)=0 conve(1,1,1)=12e-3 conve(2,0,1)=200 conve(2,1,1)=18.23e-3 conve(3,0,1)=400 conve(3,1,1)=40.38e-3 conve(4,0,1)=600 conve(4,1,1)=72.69e-3 conve(5,0,1)=800 conve(5,1,1)=125.00e-3 conve(6,0,1)=1000 conve(6,1,1)=150.38e-3 conve(7,0,1)=1500 conve(7,1,1)=200.00e-3 conve(8,0,1)=2000 conve(8,1,1)=250.65e-3 mpcopy,,1,2 !################################定义数据 et,1,solid70 et,2,solid70 r=0.0025 P=1300 v=0.003 pi=3.1415926 efficient=0.5 Qm=3*P*efficient/(pi*r*2) kk=0.85 Qmax=Qm*kk lx=0.015 Lz=0.06 ly=0.008 cladding=0.002 lzsize=0.003 !z方向上单元为3mm lxsize=cladding/6 !x方向上单元长度 将x方向6等分 tz=lz/v !z方向上热源走完需要的总时间 tl=lzsize/v !z方向上 走一个单元所需要的时间 单元时间 也就是3/3=1s step_time=tz/tl !时间步长 需要走的时间步数 step_z=lz/lzsize!z方向上的总步长数（总长度除一个单元长度） step_x=cladding/lxsize !x方向的步长 熔覆层半圆形总共分了6等分 !##########################################建立模型 BLOCK,0,cladding,0,-ly,0,lz BLOCK,cladding,lx,0,-ly,0,lz CYL4,0,0,0.002,90,0,0,lz VGLUE,1,2,3 NUMCMP,ALL lplot !##########################################################################################网格划分 lesize,11,,,20 lesize,17,,,20 lesize,12,,,20 lesize,2,,,6 lesize,7,,,6 lesize,18,,,6 lesize,21,,,6 mshape,0,3d mshkey,1 type,1 mat,1 VMESH,2 lplot lesize,1,,,4,2 lesize,8,,,4,2 lesize,9,,,20 type,1 mat,2 vmesh,1 lplot lesize,25,,,13,2 lesize,24,,,13,2 lesize,22,,,13,2 lesize,23,,,13,2 type,1 mat,2 vmesh,3 !##########################################################################################选择生死单元 alls vsel,s,volu,,2 ESLV,s,2 eplot ekill,all alls esel,s,live eplot !########################定义变量 zmin=-1e-5 zmax=lzsize zd=lzsize *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNCCSYS *SET,_FNCNAME,'SU' *SET,_FNCCSYS,0 ! /INPUT,.\Desktop\1.func,,,1 *DIM,%_FNCNAME%,TABLE,6,22,1,,,,%_FNCCSYS% ! ! Begin of equation: 2.866242e7*exp(-2*(({X})^2+({Z}+0.003*{TIME})^2)/0.004^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, 2, 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, 0.003, 0, 0, 1 *SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, -1, 3, 1 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, 4, 1, -4 *SET,%_FNCNAME%(0,11,1), 0.0, -4, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, -1, 17, -4 *SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 1, -3, 1, -5 *SET,%_FNCNAME%(0,14,1), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 0.004, 0, 0, 0 *SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,17,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 1, -3, 4, -4 *SET,%_FNCNAME%(0,19,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 2.866242e7, 0, 0, -1 *SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, -2, 3, -1 *SET,%_FNCNAME%(0,22,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: 2.866242e7*exp(-2*(({X})^2+({Z}+0.003*{TIME})^2)/0.004^2) !--> !####################################非求解选项 /solu antype,4 ! ANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action, --, PRELP（Antype:0静态分析/1瞬态分析） trnopt,full !完全瞬态积分法 !TRNOPT, Method, MAXMODE, --, MINMODE, MCout, TINTOPT, VAout, DMPSFreq 定义瞬态分析选项（Method,1完全的方法） pred,on !打开预测校正 !PRED, Sskey, --, Lskey，在非线性分析中激活一个预测器（On,接下来分析都使用预测器） solcontrol,on !自动控制 !autots, key 是否使用自动时间步长，key  n: 当solcontrol为on时缺省为on，当solcontrol为off时缺省为offnropt,full,,on !完全牛顿-拉普森方法 !NROPT, Option1, Option2, Optva ，在静态分析和动态分析中，定义牛顿-拉普拉斯选项，Option1：full(采用完全牛顿拉普莱斯)，Option1为full时Option2无效，Optva ：on(使用自适应系统) timint,on !瞬态效果 tref,25 !参考温度 !TREF, tref(定义热应变计算参考温度) kbc,0 !连续性载荷 !KBC, KEY（指定加载步骤内的倾斜加载或步进加载），KEY：默认为线性加载，值为0 neqit,25 !规定每个子步中最大迭代次数,默认为25 ！NEQIT, NEQIT, FORCEkey（指定非线性分析中的最大迭代步数） outres,all !OUTRES, Item, Freq, Cname, -- , NSVAR, DSUBres ：控制求解数据写到数据库中，all是默认值 !###########################稳态分析 timint,off alls ic,all,temp,25 alls asel,s,area,,1,2 asel,a,area,,12,13 asel,a,area,,10 asel,a,area,,15 asel,a,area,,7,9 aplot nsla,s,1 sf,all,conv,%conve%,25 time,1e-5 deltim,1e-5,1e-5,1e-5 kbc,1 allsel,all solve eplot !############################瞬态分析 m=1 timint,on *do,c,1,20 /solu esel,s,mat,,1 nsel,r,loc,z,zmin,zmax nsel,r,loc,x,zmin,zmax nsel,r,loc,y,zmin,zmax esln,s,all ealive,all sf,all,hflux,%SU% time,m deltim,0.5,0.5,0.5 alls solve esel,s,mat,,1 nsel,r,loc,z,zmin,zmax nsel,r,loc,x,zmin,zmax nsel,r,loc,y,zmin,zmax esln,s,all SFDELE,all,hflux m=m+tl zmin=zmin+zd zmax=zmax+zd *enddo alls eplot |
回复此楼» 猜你喜欢
普通院校药学硕士，做合成的，感觉找不到工作
已经有17人回复
-
如何分离原料与产品
已经有13人回复
-
LWT投
已经有3人回复
-
真的好想读博！
已经有22人回复
-
博士白读了
已经有49人回复
-
基金开始函评了吗？
已经有16人回复
-
没学上
已经有12人回复
-
研二光催化6月底四篇二区什么水平
已经有15人回复
-
环肽的合成
已经有5人回复
-
学硕专硕
已经有5人回复
2楼2022-11-17 20:28:20
【答案】应助回帖
感谢参与，应助指数 +1
|
简单看了一下你的命令流，有一些问题： 首先，你定义的对流+辐射的换热系数表数组，换热系数的值太小。 更重要的是，你使用函数进行热流的加载，但是函数定义的部分放错了为止。按照 ANSYS 的规定，函数的定义部分，必须放在紧靠使用该函数的命令前面，原因是，你是在瞬态分析的循环中使用的该函数，而函数中使用的参数值 x, y, z 是变化的。你把函数定义放在循环的前面，那么其中的参数值都是取得循环之前的值，且始终不变，而不是循环中，热源位置变化之后的值，所有计算是错误的。 再有，在循环命令 *do,c,1,20 下面，多了一个 /solu，也需要去掉。这是因为 ANSYS 中，/solu 表示出现进入求解模块，将开始一个新的求解过程，即不在从上一个求解过程的结果继续进行分析。这也和你的医院不一致吧？所以需要去掉该 /solu 命令，以便计算是连续进行的。 可能还会有其他问题，这里不是详细讨论的地方，先说这些吧。 |
» 本帖已获得的红花（最新10朵）
Su peng3楼2022-11-20 09:01:26
Cfdstudio8
禁虫 (初入文坛)
|
本帖内容被屏蔽 |
4楼2022-11-20 11:38:27
5楼2022-11-21 16:22:40
6楼2022-11-22 08:26:19
5