| 查看: 472 | 回复: 2 | ||
| 【悬赏金币】回答本帖问题，作者Juneok123将赠送您 5 个金币 | ||
Juneok123新虫 (初入文坛)
|
[求助]
ansys高斯热源模拟
|
|
|
/com, units: Mpa, mm - c - mW - mg -sec - mJ /title,LBW /prep7 et,1,solid70 mptemp,1,0,25,100,200,300,400 mptemp,,500,550,580,600,620,650 mptemp,,700,800,900,1000,1100,1200 mpdata,dens,1,1,1.79,1.79,1.79,1.79,1.79,1.79 mpdata,dens,1,,1.79,1.79,1.79,1.79,1.79,1.79 mpdata,dens,1,,1.79,1.79,1.79,1.79,1.79,1.79 mpdata,kxx,1,1,83.8,83.8,88,92.2,101,105 mpdata,kxx,1,,111.2,114.1,115.9,117.1,118.2,119.99 mpdata,kxx,1,,122.9,128.74,134.6,140.4,140.4,140.4 mpdata,c,1,1,1.089,1.089,1.090,1.130,1.210,1.269 mpdata,c,1,,1.269,1.269,1.269,1.269,1.269,1.269 mpdata,c,1,,1.269,1.269,1.269,1.269,1.269,1.269 mpdata,enth,1,1,0,25*1.79*1.089,100*1.79*1.090,200*1.79*1.130,300*1.79*1.210,400*1.79*1.269 mpdata,enth,1,,500*1.79*1.269,550*1.79*1.269,580*1.79*1.269,600*1.79*1.269,620*1.79*1.269+1.79*600,650*1.79*1.269+1.79*600 mpdata,enth,1,,700*1.79*1.269+1.79*600,800*1.79*1.269+1.79*600,900*1.79*1.269+1.79*600,1000*1.79*1.269+1.79*600,1100*1.79*1.269+1.79*600,1200*1.79*1.269+1.79*600 *dim,conve,table,8,1,1,temp 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 rectng,-4,0,0,10 rectng,0,4,0,10 rectng,-50,50,0,10 aovlap,all aglue,all vext,all,,,,,100 vglue,all esize,1 type,1 mat,1 vmesh,1 vmesh,2 esize,5 type,1 mat,2 vmesh,3 vmesh,4 /solu antype,4 trnopt,full pred,on solcontrol,on nropt,full,,on !cnvtol,u,,0.5 !cnvtol,f,,0.05 cutcontrol,plslimit,10, !eqslv,sparse,,2 lnsrch,on autots,on tintp,,,,1,, timint,on tref,20 kbc,0 neqit,50 outres,all *********************稳态分析处于初始条件****************** timint,off alls ic,all,temp,20 alls nsel,s,ext sf,all,conv,%conve%,20 time,1e-3 deltim,1e-3,1e-3,1e-3 kbc,1 allsel,all solve eplot ******************************热源加载***************************** lsize=1 v=20 tinc=lsize/v rb=4 pi=3.1415926 p=8000e3 K=0.7 qm=(3*k*p)/(pi*16) max_time=5 max_z=5 max_x=5 *dim,laser,table,max_z,max_x,max_time,z,x,time *do,k,1,max_time,1 *do,i,1,max_z,1 laser(i,0,k)=(i-1)*lsize *enddo *do,j,1,max_x,1 laser(0,j,k)=(j-1)*lsize *enddo *do,i,1,max_z,1 *do,j,1,max_x,1 zcenter=v*(k-1)*tinc xcenter=0 *if,k,le,6,then distance=sqrt(((i-1)*lsize-zcenter)**2+((j-1)*lsize-xcenter)**2) *if,distence,le,rb,then laser(i,j,k)=qm/exp(2*distance**2/16) *else laser(i,j,k)=0 *endif *enddo *enddo laser(0,0,k)=(k-1)*tinc *enddo ************求解*************** t=0.05 timint,on nsle,s,loc,y,10 *do,t,1,5,tinc time,t deltim,,,,1 sf,all,hflux,%laser% solve *enddo finish ************************************************************************************************* 麻烦大神帮忙看看，问题：热源加载是否正确呢？ 求解设置是否又问题啊 @beefly @月只蓝 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有77人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有121人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
Juneok123
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 50.3
- 散金: 19
- 帖子: 49
- 在线: 7.6小时
- 虫号: 22310065
- 注册: 2020-05-25
- 专业: 金属结构材料
2楼2020-07-23 23:59:42
Juneok123
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 50.3
- 散金: 19
- 帖子: 49
- 在线: 7.6小时
- 虫号: 22310065
- 注册: 2020-05-25
- 专业: 金属结构材料
3楼2020-07-24 08:00:55
