| 查看: 809 | 回复: 0 | |||
[交流]
求助，利用table加载高斯热源模型求解没用
|
|
模型大致是，在一个基板上有一个0.3mm熔覆层，想在熔覆层上表面加载移动的高斯热源模型。在设置时，我用到了设定了一个局部坐标系12，不知道是否有影响 命令流如下： Finish /filename,TA15 Ti2AlNb SINGLE /prep7 et,1,solid70 et,2,surf152 keyopt,2,4,1 !solid 70单元没有中间节点 keyopt,2,5,0 !表面效应单元加载对流，0和1均可? keyopt,2,8,4 !根据流体温度计算对流系数 !根据Jmat Pro 模拟数据输入材料参数 !2号材料为薄壁 MPREAD,'TA15','dat','C:\Users\laser\Desktop\' MPREAD,'Ti-22Al-25Nb','dat','C:\Users\laser\Desktop\' block,0,0.15,0,0.006,0,0.1 wpave,0.05,0.006,0.05 block,0,0.048,0,0.0003,0,0.005 WPAVE,0,0,0 CSYS,0 vglue,all wpave,0.05,0.006,0.05 vsbw,all,,delete wpave,0.05,0.006,0.055 vsbw,all,,delete wprot,0,0,90 vsbw,all,,delete wpave,0.098,0.006,0.055 vsbw,all,,delete nummrg,all numcmp,all vsel,s,,,2 vatt,2,,1,0 vsel,s,,,1 vsel,a,,,3,10,1 vatt,1,,1,0 allsel vsel,s,mat,,2 esize,0.001 MSHAPE,0,3d MSHKEY,1 vmesh,all vsel,s,mat,,1 esize,0.003 MSHAPE,0,3d MSHKEY,1 vmesh,all !建模 划分网格 /solu antype,trans !瞬态分析 OUTRES,,ALL !输出设置 toffst,273 !温度偏移量设置 nropt,full timint,1,therm !打开瞬态热分析 tunif,25 !初始温度 autos,on !自动时间步长 solcontrol,on kbc,0 !渐变载荷 pred,on !打开时间预测器 lnsrch,on !瞬态热分析设置 CNVTOL,HEAT, ,0.01,1,0.000001, !控制收敛准则 LUMPM,0 !不使用集中质量矩阵 asel,s,loc,x,0 sfa,all,,conv,30,25 asel,s,loc,x,0.15 sfa,all,,conv,30,25 asel,s,loc,z,0 sfa,all,,conv,30,25 asel,s,loc,z,0.1 sfa,all,,conv,30,25 asel,s,loc,y,0 sfa,all,,conv,30,25 asel,s,loc,y,0.006 asel,u,,,9 sfa,all,,conv,30,25 local,12,0,0.05,0.0063,0.05 csys,12 asel,s,loc,y,0 lsize=0.001 v=0.003 tinc=lsize/v r=0.002 qmax=1.31e8 max_time=49 max_x=49 max_z=6 *dim,flux2,table,max_x,max_z,max_time,x,z,time *do,k,1,max_time,1 *do,i,1,max_x,1 flux2(i,0,k)=(i-1)*lsize *enddo *do,j,1,max_z,1 flux2(0,j,k)=(j-1)*lsize *enddo *do,i,1,max_x,1 *do,j,1,max_z,1 xcenter=v*(k-1)*tinc zcenter=0.0025 distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-zcenter)**2) *if,distance,le,r,then flux2(i,j,k)=qmax/exp(3*(distance**2)/(r**2)) *else flux2(i,j,k)=0 *endif *enddo *enddo flux2(0,0,k)=(k-1)*tinc *enddo *do,t,tinc,16,tinc asel,s,loc,y,0 TSRES,ERASE sfa,all,1,hflux,%flux2% allsel,all deltim,0.1,0.05,1/3 time,t solve *enddo 求指点，哪出错了  模型.jpg |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有10人回复
-
基金系统什么内容也没有
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
30
