| 查看: 2914 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
ansys焊接温度场模拟
|
|||
|
刚刚接触ansys，老师就让做焊接温度场模拟，现在有一段代码如下： /UNITS,SI /CONFIG,NRES,10000 /PREP7 ET,1,SOLID70 !单元类型，1表示此单元编号为1 MP,DENS,1,7930 !定义材料密度 MP,C,1,502 !定义材料比热 mptemp,1,20,100,500 mpdata,kxx,1,1,12.1,16.3,21.4 !kxx热导系数 /VIEW,1,1,1,1 !*************定义基板尺寸*********** lx=0.1 lz=0.1 ly=0.006 !基板大小 !*************定义焊接参数*********** layer=7 !总的焊接道数 lap=1/3 !两道之间搭接量 weld=0.01 !每道焊道宽度 wweld=layer*weld-(layer-1)*lap*weld !焊道的总宽度 wheight=0.004 !焊道的高度 remain=0.01 LSIZE=0.002 V=0.008 !速度 pi=3.1415926 U=25 I=180 Q=U*I !电源有效功率（在此，可以对Q通过UI来调整） R=0.006 !电源有效热半径 effect=0.8 Qmax=effect*Q/(pi*R*R) !中心处最大热流 tm=(lx-2*remain)/V !一道的总时间 tinc=LSIZE/V !每小段的时间间隔 !*****************建立模型，生成网格************** block,0,lx,0,ly,0,lz !创建一个块体积，可以根据基板尺寸调整坐标 wpoff,remain,ly,(lz-wweld)/2 !随着坐标系轴的变化平移工作平面 block,0,lx-2*remain,0,wheight,0,wweld vsbw,all,,delete !对block布尔运算 wpoff,,,wweld vsbw,all,,delete wprota,90,90 !转动工作面 vsbw,all,,delete wpoff,,,lx-2*remain vsbw,all,,delete save vglue,all !粘合所有的几何体 lsel,s,,,9,12,1 lsel,a,,,25,26,1 lsel,a,,,30 lsel,a,,,32 lsel,a,,,38 lsel,a,,,45 lsel,a,,,59 lsel,a,,,73 lesize,all,0.005,,,0.2,,,,1 allsel,all lsel,s,,,13,16,1 lsel,a,,,19,21,1 lsel,a,,,23,24,1 lsel,a,,,34 lsel,a,,,41,44,1 lsel,a,,,63,64,1 lsel,a,,,77,78,1 lsel,a,,,81,82,1 lesize,all,0.001,,,1,,,,1 allsel,all lsel,s,,,1 lsel,a,,,3,6,1 lsel,a,,,8 lsel,a,,,27 lsel,a,,,29 lsel,a,,,31 lsel,a,,,33 lsel,a,,,35 lsel,a,,,37 lsel,a,,,40 lsel,a,,,50 lsel,a,,,53 lsel,a,,,55 lsel,a,,,62 lsel,a,,,67 lsel,a,,,69 lsel,a,,,85 lesize,all,0.002,,,1,,,,1 allsel,all lsel,s,,,46 lsel,a,,,86 lsel,a,,,88 lesize,all,0.001,,,5,,,,1 allsel,all vsweep,all /PNUM,DEFA EPLOT fini !*******************进入后处理器********************* /solu antype,trans,new !nlgeom,on timint,0,struct timint,1,therm timint,0,mag timint,0,elect tref,25 nropt,auto autos,on kbc,0 pred,on lnsrch,on !************杀死焊缝去单元*********** nsel,s,loc,y,ly,ly+wheight esln,s,1 ekill,all allsel,all esel,s,live eplot !**********施加对流载荷************* esel,s,live nsel,s,loc,z,0 nsel,a,loc,z,lz nsel,a,loc,x,0 nsel,a,loc,x,lx sf,all,conv,10,25 !************定义数组维数********** MAX_TIME=1+tm/tinc !MAX_X=1+(lx-2*remain)/LSIZE MAX_X=1+lx/LSIZE MAX_Z=1+lz/LSIZE !*************定义table数组 ****************************************** *DIM,heatflux,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME *do,n,1,layer,1 m=mod(n,2) *if,m,eq,1,then !n=1 *DO,K,1,MAX_TIME,1 *DO,I,1,MAX_X,1 heatflux(I,0,K)=(I-1)*LSIZE *ENDDO *DO,J,1,MAX_Z,1 heatflux(0,J,K)=(J-1)*LSIZE *ENDDO *DO,I,1,MAX_X,1 *DO,J,1,MAX_Z,1 XCENTER=remain+V*(K-1)*tinc ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2 DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE- ZCENTER)**2) *IF,DISTANCE,LE,R,THEN heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2)) *ELSE heatflux(I,J,K)=0 *ENDIF *ENDDO *ENDDO heatflux(0,0,K)=(K-1)*tinc !**********激活单元********* esel,s,live nsel,s,loc,x,xcenter-0.002,xcenter+0.001 nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2 nsel,r,loc,y,ly,ly+wheight esln,s,0 ealive,all allsel,all !**********施加热流载荷********** esel,s,live nsel,s,loc,Y,ly sf,all,hflux,%heatflux% allsel,all !************求解***************** time,(n-1)*(tm+tinc)+k*lsize/v nsubst,2 SOLVE !*************删除热流载荷*********** nsel,s,loc,y,ly sfdele,all,hflux allsel,all esel,s,live eplot OUTRES,ALL,ALL, /PSF,HFLUX,,2 /REPLOT OUTRES,ALL,ALL, *ENDDO *else !n=2 *DO,K,1,MAX_TIME,1 *DO,I,1,MAX_X,1 heatflux(I,0,K)=(I-1)*LSIZE *ENDDO *DO,J,1,MAX_Z,1 heatflux(0,J,K)=(J-1)*LSIZE *ENDDO *DO,I,1,MAX_X,1 *DO,J,1,MAX_Z,1 XCENTER=lx-remain-V*(K-1)*tinc ZCENTER=2*(n-1)*weld/3+(lz-wweld)/2+weld/2 DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE- ZCENTER)**2) *IF,DISTANCE,LE,R,THEN heatflux(I,J,K)=Qmax/EXP(2*DISTANCE**2/(R**2)) *ELSE heatflux(I,J,K)=0 *ENDIF *ENDDO *ENDDO heatflux(0,0,K)=(K-1)*tinc !**********激活单元********* esel,s,live nsel,s,loc,x,xcenter-0.001,xcenter+0.002 nsel,r,loc,z,zcenter-weld/2,zcenter+weld/2 nsel,r,loc,y,ly,ly+wheight esln,s,0 ealive,all allsel,all !**********施加热流载荷********** esel,s,live nsel,s,loc,Y,ly sf,all,hflux,%heatflux% allsel,all !************求解***************** time,(n-1)*(tm+tinc)+k*lsize/v nsubst,2 SOLVE !*************删除热流载荷*********** nsel,s,loc,y,ly sfdele,all,hflux allsel,all esel,s,live eplot OUTRES,ALL,ALL, /PSF,HFLUX,,2 /REPLOT OUTRES,ALL,ALL, *ENDDO *endif *enddo FINISH 请问各位高手们，里面lsel命令是在选择线段吗？？选线段目的是什么？里面线段的标号是哪里定义的？？？？ 另外，要是手里有可用的模拟焊接温度场资料，请大家不吝施舍啊！！！ |
回复此楼» 收录本帖的淘帖专辑推荐
 焊接有限元分析 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- ansys apdl 热-应力耦合算残余应力，材料参数是哪些
已经有19人回复
- 在ansys中做齿轮瞬态温度场分析，遇到问题，求大神指点
已经有9人回复
- ANSYS如何把已做出结构分析的结果倒入温度场分析中？
已经有12人回复
- ANSYS热模拟高温拉伸，结构场计算出现问题，麻烦各位帮忙看看，附详细命令流
已经有14人回复
- 求助！Ansys脉冲激光辐照材料的温度场模拟！
已经有12人回复
- 请教各位大侠，关于ABAQUS单位量纲的问题~
已经有13人回复
- simufact多道焊模拟（生死单元法）
已经有9人回复
- 岩土工程ANSYS经典72例（含命令代码）
已经有285人回复
- 请教下各位达人，学习ansys或者abaqus，对焊接温度场进行模拟，来得及吗？
已经有7人回复
- ANSYS10.0在流体计算中遇到的问题
已经有6人回复
- 本人研究生做材料加工ANSYS温度场，应力场数值模拟，考上交材料加工什么方向适合？
已经有9人回复
- ansys可以模拟温度场么
已经有9人回复
- ansys 焊接温度场模拟命令流
已经有4人回复
5楼2013-07-26 11:31:53
2楼2013-07-25 17:01:10
3楼2013-07-25 17:35:13
4楼2013-07-25 19:09:14
30