| 查看: 1663 | 回复: 3 | |||
| 【悬赏金币】回答本帖问题，作者yunshan216将赠送您 5 个金币 | |||
yunshan216新虫 (正式写手)
|
[求助]
边界设置问题已有1人参与
|
||
|
我在一个模型上施加一个随机热源，理论上温度场应该是左右对称的（圆形分布），但是由于靠近边界，导致不对称了，需要怎样设置边界条件呀？ 我只加了载荷，其他的面都没有设置。 命令流如下： /NOPR KEYW,PR_SET,1 KEYW,PR_STRUC,0 KEYW,PR_THERM,1 /GO /COM, Thermal /PREP7 ET,1,SOLID70 ET,2,SURF152 KEYOPT,2,4,1 KEYOPT,2,8,1 /units,si TOFFST,273 MP,KXX,1,15 MP,DENS,1,7930 MP,C,1,510 MPTEMP,1,0,1400,1445,2849 !突变太大不容易收敛 MPDATA,ENTH,1,,0,5.75906e9,7.92e9,1.3675e10 *set,rj,200e-6 *set,u,20 *set,id,0.5 *set,td,30 *set,t_sub1,1 !脉宽间子步长度、微秒 *set,tran_r,1 !脉宽步长 *set,F,0.4 !极间能量分配系数 *set,t_off,970e-6 !脉冲间隙 *set,t_sub2,10 !脉冲间隙子步长度 *set,t_fre,1 !脉冲周期次数 cylind,rj,,0,-2*rj,-90,90 wprota,,90 vsbw,all wprota,,-90 /pnum,line,1 lplot mshape,0,3d !指定划分单元的形状 mshkey,1 !指定映射网格划分 lsel,s,loc,z,0 lesize,all,,,30 lsel,s,length,,2*rj lesize,all,,,40 vmesh,all asel,s,,,8 asel,a,,,12 nsla,s,1 /VIEW,1,1,2,3 eplot ! /ANG,1 ! /REP,FAST TYPE, 2 MAT, 1 REAL, ESYS, 0 SECNUM, esurf,all alls FINISH /SOLU ANTYPE,4 ALLSEL,ALL nsel,s,,,all IC,all,TEMP,25, !环境温度25度 TIME,1e-7 !设置时间 AUTOTS,1 KBC,1 TIMINT,0 DELTIM,0.1e-6, , ,1 OUTRES,ALL,ALL, ALLSEL,ALL SOLVE *do,m,1,t_fre,1 !每步计算循环注意 *set,rand_z,rand(-2*rj,0) !设置随机数 *do,n,tran_r,td,tran_r esel,s,type,,2 *get,nemax_2,elem,,num,max !获取单元最大序号 *get,nemin_2,elem,,num,min !获取单元最小序号 *do,i,nemin_2,nemax_2,1 !每一个单元循环 i为运算时单元号 *if,esel(i),eq,1,then !如果单元被选择则执行 获取中心点坐标 *SET,xsy,centrx(i) *SET,ysy,centry(i) *SET,zsy,centrz(i) *SET,rr,abs(sqrt((xsy/1)*(xsy/1)+(zsy-rand_z)*(zsy-rand_z))) *set,rt,14.823*n**0.2648*id**0.3733*1e-6 *if,ysy,gt,0,then *if,rr,le,rt,then *SET,qr,4.55*F*u*id/3.1415926/rt/rt*exp(-4.5*rr*rr/rt/rt) sfe,i,1,hflux,,qr !在i单元表面施加面载荷 *endif *endif *endif *enddo alls TIME,(m-1)*(t_off+td*1e-6)+n*1e-6 !设置时间 AUTOTS,1 DELTIM,T_sub1*1e-6, , ,1 KBC,1 TSRES,ERASE TIMINT,1 TINTP,0.005, , ,1,0.5,0.2, !*IF,N,LT,td,THEN !OUTRES,all,NONE, !*ELSEIF,N,EQ,td OUTRES,BASIC,ALL !*ENDIF ALLSEL,ALL LNSRCH,1 ALLSEL,ALL SOLVE !脉冲结束t=td *ENDDO esel,s,type,,2 sfedele,all,all,hflux !sfa,3,1,conv,10,25 TIME,m*(t_off+td*1e-6) AUTOTS,1 DELTIM,t_sub2*1e-6, , ,1 KBC,1 TSRES,ERASE TIMINT,1 TINTP,0.005, , ,1,0.5,0.2, OUTRES,BASIC,all, ALLSEL,ALL LNSRCH,1 ALLSEL,ALL SOLVE *enddo  QQ图片20180522130254.jpg |
回复此楼» 猜你喜欢
博士白读了
已经有50人回复
-
with editor 两个月了，什么原因？
已经有16人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有17人回复
-
如何分离原料与产品
已经有13人回复
-
LWT投
已经有3人回复
-
真的好想读博！
已经有22人回复
-
基金开始函评了吗？
已经有16人回复
-
没学上
已经有12人回复
-
研二光催化6月底四篇二区什么水平
已经有15人回复
-
环肽的合成
已经有5人回复
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13377.8
- 红花: 265
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
2楼2018-05-22 17:34:19
yunshan216
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1285
- 红花: 4
- 沙发: 1
- 帖子: 510
- 在线: 34.6小时
- 虫号: 8588912
- 注册: 2018-04-22
- 专业: 机械摩擦学与表面技术
3楼2018-09-15 22:12:13
shgao20
专家顾问 (职业作家)
-
专家经验: +809 - 仿真EPI: 1
- 应助: 1311 (讲师)
- 金币: 13377.8
- 红花: 265
- 帖子: 3531
- 在线: 473.4小时
- 虫号: 3644138
- 注册: 2015-01-14
- 专业: 动力学与控制
- 管辖: 仿真模拟
4楼2018-09-16 07:10:55
5