| 查看: 965 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者lk12将赠送您 5 个金币 | ||
[求助]
ansys中使用表面单元加载对流载荷，求救如何使用
|
||
|
请教大佬表面效应单元的问题，我想在这个程序中增加一个表面效应的对流载荷，显示单元4未定义，esurf只能应用在线或者面上，请问怎么改呢？ FINISH /CLEAR,START /TITLE,Laser /units,si /VIEW,1,1,1,1 /prep7 !===定义单元 et,1,solid70 et,4,surf151 keyopt,4,4,1 keyopt,4,5,0 keyopt,4,8,4 !===定义材料 !===熔覆层为铝zl104 mp,dens,1,2630 !密度(kg/m^3) mp,hf,1,100 mptemp,,,,,,,, mptemp,1,25 mptemp,2,100 mptemp,3,200 mptemp,4,300 mptemp,5,400 mptemp,6,500 mptemp,7,600 mptemp,8,700 mptemp,9,800 mptemp,10,900 mptemp,11,1000 mptemp,12,1100 mptemp,13,1200 mptemp,14,1300 mpdata,c,1,,725 !比热容(J/(kg℃)) mpdata,c,1,,754 mpdata,c,1,,796 mpdata,c,1,,837 mpdata,c,1,,921 mpdata,c,1,,1011 mpdata,c,1,,1100 mpdata,c,1,,1130 mpdata,c,1,,1150 mpdata,c,1,,1170 mpdata,c,1,,1190 mpdata,c,1,,1210 mpdata,c,1,,1250 mpdata,kxx,1,,113 !热导率(W/(m℃)) mpdata,kxx,1,,144.9 mpdata,kxx,1,,149.1 mpdata,kxx,1,,149.1 mpdata,kxx,1,,144.9 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 mpdata,kxx,1,,145.1 !===熔覆层尺寸(m) Width_Clad=2.5e-3 !熔覆层宽度 Length_Clad=30e-3 Height_Clad=0.35e-3 Layer=5 !==基材尺寸(m) Width_Base=30e-3 Length_Base=50e-3 Height_Base=6e-3 !===建模 block,-Width_Clad/2,Width_Clad/2,0,Length_Clad,0,Height_Clad*Layer block,-Width_Clad/2,Width_Clad/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2 block,-Width_Clad/2,-Width_Base/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2 block,Width_Clad/2,Width_Base/2,-10/1000,Length_Base-10/1000,0,-Height_Base/2 vglue,all nummrg,all numcmp,all !===网格划分 type,1 mat,1 lesize,9,,,Layer lesize,10,,,Layer lesize,11,,,Layer lesize,12,,,Layer esize,2.2/1000/6 vmesh,1 type,1 mat,1 lesize,21,0.001,,,3 lesize,22,0.001,,,3 lesize,19,0.001,,,3 lesize,20,0.001,,,1/3 lesize,21,0.001,,,1/3 lesize,22,0.001,,,3 vsweep,4 type,1 mat,1 lesize,31,0.001,,,3 lesize,32,0.001,,,3 lesize,33,0.001,,,3 lesize,34,0.001,,,3 lesize,37,0.001,,,3 lesize,38,0.001,,,3 lesize,39,0.001,,,3 lesize,40,0.001,,,3 vsweep,2 vsweep,3 !=========激光参数 LaswerPower=1400 !激光功率(W) Radius1=0.85/1000 !熔覆层半径(m) Radius2=1.2/1000 !光斑半径(m) Factor=0.35 !激光利用率 TruePower=LaswerPower*Factor pi=3.14159265358979310 dt=0.00001 !小量 Init_Temp=20!初始温度 Coff_Conv=100 !对流换热系数 AirTemp=20 AirConvection=30 BaseTemp=20 BaseConvection=80 !===求解参数设置 StepDis=2.2/1000/6 !每个载荷步位移(m) V=6/1000 !扫描速度(m/s) StepTime=StepDis/V !单步时间长度(s) TotalTime=(Length_Clad-2*Radius1)/V !总时长(s) !StepNum=(Length_Clad-2*Radius1)/StepDis !载荷步数 StepNum=76 !载荷步数 !===对流 nsel,s,ext nsel,u,loc,z,Height_Clad*Layer SF,all,CONV,30,20 allsel,all !===对流 nsel,s,ext nsel,r,loc,z,Height_Clad*5 asel,s,loc,z,Height_Clad*5 esurf,all ESEL,S,TYPE,,4 SFE,all,1,CONV,1,100 SFE,all,1,CONV,2,20 allsel,all !===求解部分 /solu antype,trans nropt,full tintpr,,,,1.0 trnopt,full !指定瞬态分析选项 timint,on !时间积分效应 autots,on !指定使用自动步长跟踪 deltim,,,, nsubst,1 kbc,1 TUNIF,BaseTemp !初始温度20度 !toffst,273 !===杀死所有的熔覆层单元 VSEL,S,,,1 !选择熔覆层体 ESLV,S !选择依附于所选体上的单元 EKILL,ALL !杀死所选的单元(熔覆层) ESEL,S,LIVE !选择活着的单元(基体) EPLOT !刷新显示 LOCAL,11,0,0,0,0 !定义局部坐标系 *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNCCSYS *SET,_FNCNAME,'huanxing' *SET,_FNCCSYS,11 !将方程坐标系设定为自己定义的局部坐标系 ! /INPUT,.\function\huanxing.func,,,1 *DIM,%_FNCNAME%,TABLE,6,22,1,,,,%_FNCCSYS% ! Begin of equation: 1.1E8*EXP(-2*((sqrt({X}^2+{Y}^2)-0.00106)^2)/(0.00000025)) *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, 2, 0, 0, 3 *SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 1, 3, 17, -1 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -3, 1, -4 *SET,%_FNCNAME%(0,11,1), 0.0, -1, 16, 1, -1, 0, 0 *SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 0.00106, 0, 0, -1 *SET,%_FNCNAME%(0,13,1), 0.0, -4, 0, 1, -1, 2, -3 *SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, -4 *SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -4, 17, -1 *SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -2, 3, -3 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 0.00000025, 0, 0, -1 *SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -1, 4, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -1, 7, 1, -3, 0, 0 *SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 1.1E8, 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: 1.1E8*EXP(-2*((sqrt({X}^2+{Y}^2)-0.00106)^2)/(0.00000025)) !--> !===第一层 time=0 XC=0 YC=3*StepDis *DO,K,0,StepNum,1 LOCAL,11,0,XC,YC,0 t=0.45+StepTime*(k+1) time,t !===基于局部坐标系下复活第一层熔覆层 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,0,Height_Clad nsel,r,loc,y,-3*StepDis,3*StepDis esln,s,1 ealive,all allsel,all ESEL,S,LIVE EPLOT !===第一层加载热源 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad nsel,r,loc,y,-3*StepDis,3*StepDis SF,all,HFLUX,%huanxing% ESEL,S,LIVE /PSF,HFLUX,,2 /REPLOT EPLOT ALLSEL,ALL OUTRES,ALL,ALL SOLVE YC=YC+StepDis *ENDDO !===第二层 CSYS,0 XC=0 YC=3*StepDis *DO,k,0,StepNum,1 LOCAL,11,0,XC,YC,0 t=5.45+StepTime*(k+1) time,t !===基于局部坐标系下复活第二层熔覆层 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad,Height_Clad*2 nsel,r,loc,y,-3*StepDis,3*StepDis esln,s,1 ealive,all allsel,all ESEL,S,LIVE EPLOT !===第二层加载热源 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*2 nsel,r,loc,y,-3*StepDis,3*StepDis SF,all,HFLUX,%huanxing% ESEL,S,LIVE /PSF,HFLUX,,2 /REPLOT EPLOT ALLSEL,ALL OUTRES,ALL,ALL SOLVE YC=YC+StepDis *ENDDO !===第三层 CSYS,0 XC=0 YC=3*StepDis *DO,k,0,StepNum,1 LOCAL,11,0,XC,YC,0 t=10.45+StepTime*(k+1) time,t !===基于局部坐标系下复活第三层熔覆层 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*2,Height_Clad*3 nsel,r,loc,y,-3*StepDis,3*StepDis esln,s,1 ealive,all allsel,all ESEL,S,LIVE EPLOT !===第三层加载热源 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*3 nsel,r,loc,y,-3*StepDis,3*StepDis SF,all,HFLUX,%huanxing% ESEL,S,LIVE /PSF,HFLUX,,2 /REPLOT EPLOT ALLSEL,ALL OUTRES,ALL,ALL SOLVE YC=YC+StepDis *ENDDO !===第四层 CSYS,0 XC=0 YC=3*StepDis *DO,k,0,StepNum,1 LOCAL,11,0,XC,YC,0 t=15.45+StepTime*(k+1) time,t !===基于局部坐标系下复活第四层熔覆层 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*3,Height_Clad*4 nsel,r,loc,y,-3*StepDis,3*StepDis esln,s,1 ealive,all allsel,all ESEL,S,LIVE EPLOT !===第四层加载热源 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*4 nsel,r,loc,y,-3*StepDis,3*StepDis SF,all,HFLUX,%huanxing% ESEL,S,LIVE /PSF,HFLUX,,2 /REPLOT EPLOT ALLSEL,ALL OUTRES,ALL,ALL SOLVE YC=YC+StepDis *ENDDO !===第五层 CSYS,0 XC=0 YC=3*StepDis *DO,k,0,StepNum,1 LOCAL,11,0,XC,YC,0 t=20.45+StepTime*(k+1) time,t !===基于局部坐标系下复活第五层熔覆层 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*4,Height_Clad*5 nsel,r,loc,y,-3*StepDis,3*StepDis esln,s,1 ealive,all allsel,all ESEL,S,LIVE EPLOT !===第五层加载热源 nsel,s,loc,x,-Width_Clad/2,Width_Clad/2 nsel,r,loc,z,Height_Clad*5 nsel,r,loc,y,-3*StepDis,3*StepDis SF,all,HFLUX,%huanxing% ESEL,S,LIVE /PSF,HFLUX,,2 /REPLOT EPLOT ALLSEL,ALL OUTRES,ALL,ALL SOLVE YC=YC+StepDis *ENDDO FINISH /POST1 PLNSOL,TEMP,,0 |
回复此楼» 猜你喜欢
考研调剂
已经有7人回复
-
0856调剂
已经有4人回复
-
材料277求调剂
已经有5人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有6人回复
-
资源与环境 调剂申请(333分)
已经有8人回复
-
334分 一志愿武理 材料求调剂
已经有3人回复
-
284求调剂
已经有5人回复
-
0703化学/290求调剂/本科经历丰富/工科也可
已经有6人回复
-
北科281学硕材料求调剂
已经有18人回复
-
275求调剂
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2020-11-26 18:56:53
20