| 查看: 1513 | 回复: 4 | ||
[求助]
求助高斯热源中对流边界的设置 已有1人参与
|
| 高斯热源半径是是时变的 那么半径之外的对流条件怎么施加？利用表面效应单元设置对流系数 怎么设置对流边界 ？有没有方法把对流转换成热流直接放在高斯热源的公式里 ？还是把半径之外的选在一个组里？本人新手 求大神给个思路 可能还得加在曲面上 |
回复此楼» 猜你喜欢
一志愿北交大材料工程总分358求调剂
已经有6人回复
-
312求调剂
已经有7人回复
-
材料334求调剂
已经有19人回复
-
生物与医药273求调剂
已经有6人回复
-
301求调剂
已经有10人回复
-
专硕0854初试考材科基，求调剂
已经有5人回复
-
0703化学
已经有15人回复
-
0703调剂，一志愿天津大学319分
已经有9人回复
-
085600材料与化工301分求调剂院校
已经有10人回复
-
求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 加载一个高斯热源项之后，材料是空气的时候，可以看到迭代，材料换成不锈钢，迭代不了
已经有7人回复
- 施加高斯热源可以用temperature添加吗 还是就得是heat flux函数
已经有11人回复
- 在施加旋转高斯热源时总提示错误，但感觉程序没问题，不知错在哪？
已经有5人回复
- 2种焊接移动高斯热源的程序，有什么区别吗
已经有7人回复
- 【求助】关于高斯热源峰值功率随时间变化的问题
已经有7人回复
shiranjiang
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 2856.9
- 散金: 50
- 红花: 6
- 帖子: 686
- 在线: 86.6小时
- 虫号: 3041741
- 注册: 2014-03-12
- 性别: MM
- 专业: 金属结构材料
2楼2014-06-27 08:44:57
3楼2014-06-27 11:50:41
|
有没有大神救救命 建立表面效应单元 加载一个随时间变化的对流系数 但是加不上去啊 不能这样加还是我的命令流是错误的 卡在这里好几天了 finish /clear /title,transient analysis(gauss heat flux) /filnam,transientgauss /TRIAD,RBOT /UNITS,SI /prep7 et,1,plane55 et,2,surf151 keyopt,1,3,1 keyopt,2,8,2 dividenum=60 plane_x=6e-5 plane_y=3e-5 mptemp,1,298,500,700,900,1028,1100,1500,1809,3135 mpdata,c,1,1,460,529,615,822,1064,820,652,822,821 mpdata,kxx,1,1,48,46,41,35,25,26,29,29,29 mpdata,dens,1,1,7890,7700,7490,7350,7090,7840,7840,7840,7840 mpdata,enth,1,1,0,1.21e9,1.68e9,2.76e9,4.02e9,4.29e9,6.50e9,10.3e9,18.9e9 rectng,0,plane_x,0,plane_y lesize,3,,,dividenum,0.2 lesize,1,,,dividenum,5 lesize,2,,,120,0.5 lesize,4,,,120,2 amesh,1 p=0.1 u=25 i=20 pi=3.1415926 q=p*u*i*3 maxr=60 maxenum=dividenum+1 lisannum=8 ton=30 *dim,R,array,lisannum,1,1 *dim,flux2,table,maxenum,2,lisannum,x,y,k *dim,heatcf,table,maxenum,2,lisannum,x,y,k *DO,k,1,lisannum,1 t=k*Ton/lisannum R(k)=2.4*20**0.53*t**0.38 *IF,R(k),GE,maxR,then R(k)=maxR *ENDIF *do,i,1,maxEnum,1 flux2(i,0,k)=(i-1)*plane_x/dividenum heatcf(i,0,k)=(i-1)*plane_x/dividenum *enddo *do,i,1,maxEnum,1 xcenter=0 distance=((i-1)*plane_x/dividenum-xcenter)**2 *if,distance,le,(R(k)*(1E-6))**2,then flux2(i,1,k)=Q/(pi*R(k)*R(k)*1E-12*exp(3*distance/(R(k)*R(k)*1E-12)))!-Q/3/(pi*R(k)*R(k)*1E-12*exp(3*distance/(R(k)*R(k)*1E-12))) heatcf(i,1,k)=0 *else flux2(i,1,k)=0 heatcf(i,1,k)=10000 *endif *enddo *enddo nsel,s,loc,y,3e-5 sf,all,hflux,%flux2% nsel,s,loc,y,3e-5 nsel,r,loc,x,0,6e-5 type,2 esurf nsel,all esel,s,type,,2 sfe,all,,conv,,%heatcf% sfe,all,,conv,2,298 allsel /solu antype,4 trnopt,full pred,on solcontrol,on nropt,full,,on cutcontrol,plslimit,10, lnsrch,on autots,on tintp,,,,1,, timint,on tref,298 kbc,1 neqit,50 outres,all ic,all,temp,298 k=30e-6 time,k deltim,5E-6,5E-6,5E-6 alls OUTRES,ALL,ALL, solve |
4楼2014-07-06 11:17:12
5楼2014-10-13 14:11:00
