油气井工程_《现代CAE技术》大作业

1、现代CAE技术作业学生姓名xxx学号xxxxxx教学院系石油与天然气工程学院专业年级油气井工程2014级完成日期2015年6月26日第一题 3D实体结构有限元力学分析结构尺寸(单位mm)见图，轴承套受压力10MPa，底板底面固定不动。材料1: E=2.0e5MPa,v=0.3，ss=325MPa；材料2: E=2.5e5MPa, v =0.28，ss=650MPa；材料3:E=1.8e5MPa, v=0.25，ss=250MPa。要求：用ANSYS的GUI建立实体模型及有限元力学模型，并分析计算结果a.交APDL文件；b.设置材料性质；c.网格模型用6面体网格(用Sweep扫描功能)；d.交结

2、果图(实体模型和网格模型显示材料性质)；e.提交Von Mises应力计算云图结果及其最大变形结果云图；f.分析三种材料的安全系数和可能的危险位置。g.模型网格节点数、单元数分别是多少？1、APDL编程3fini/clear /filnam,201421000029_Lixin_bear-housing /PREP7 ! 建立底板RECTNG,0,200,0,100, wpoff,25,25 csys,4 CYL4,0,0,12.5 CYL4,150,0,12.5 asba,1,2asba,4,3VEXT,1, , ,0,0,25, ! 建立轴套WPSTYLE,1KWPAVE,15,16 wp

3、rot,0,90 RECTNG,-50,50,0,50, ! 建立挡板CYL4,0,50,45 asba,15,16 VEXT,17, , ,0,0,75, RECTNG,50,100,0,50, LFILLT,60,59,25, , AL,62,63,64 ASBA,24,25 CYL4,75,25,12.5 ASBA,26,24 VEXT,25, , ,0,0,15, VSYMM,X,3, , , ,0,0 /triad,off ! 剖分实体WPSTYLE,1KWPAVE,33 wprot,0,0,90VSBW,1 KWPAVE,30 VSBW,6 wprot,0,0,90VSBW,1 K

4、WPAVE,66 VSBW,all vglue,all ! 定义单元类型ET,1,SOLID185 ! 定义材料属性MP,EX,1,2.0E5 MP,PRXY,1,0.3 MP,EX,2,2.5E5 MP,PRXY,2,0.28 MP,EX,3,1.8E5 MP,PRXY,3,0.25! 赋予材料属性VSEL,18VATT,3, ,1,0 VSEL,19VATT,3, ,1,0 VSEL, ,7VATT,2, ,1, 0 VSEL, ,8VATT,2, ,1, 0 VSEL, , 1VATT,1, ,1,0 VSEL, , 2VATT,1, ,1,0 VSEL, ,5VATT,1, ,1,0

5、VSEL, , ,6VATT,1, ,1,0 VSEL, , 11VATT,1, ,1,0 VSEL, , 16VATT,1, ,1,0 VSEL, , 17VATT,1, ,1,0 !* ALLSEL,ALL /REPLOT ! 划分单元LESIZE,186, , ,25, , , , ,1 LESIZE,187, , ,25, , , , ,1 LESIZE,47, , ,40, , , , ,1 LESIZE,48, , ,40, , , , ,1 LESIZE,133, , ,50, , , , ,1 LESIZE,42, , ,5, , , , ,1 LESIZE,49, , ,5,

6、 , , , ,1 LESIZE,167, , ,40, , , , ,1 LESIZE,163, , ,40, , , , ,1 LESIZE,178, , ,10, , , , ,1 LESIZE,176, , ,10, , , , ,1 Vsweep,8Vsweep,7LESIZE,52, , ,30, , , , ,1 LESIZE,4, , ,30, , , , ,1 LESIZE,56, , ,30, , , , ,1 LESIZE,138, , ,30, , , , ,1 LESIZE,96, , ,20, , , , ,1 LESIZE,69, , ,20, , , , ,1

7、LESIZE,70, , ,10, , , , ,1 LESIZE,97, , ,10, , , , ,1LESIZE,95, , ,10, , , , ,1 LESIZE,68, , ,10, , , , ,1 LESIZE,72, , ,10, , , , ,1 LESIZE,73, , ,10, , , , ,1 LESIZE,74, , ,10, , , , ,1 LESIZE,75, , ,10, , , , ,1 LESIZE,99, , ,10, , , , ,1 LESIZE,100, , ,10, , , , ,1 LESIZE,101, , ,10, , , , ,1 LE

8、SIZE,102, , ,10, , , , ,1 Vsweep,18Vsweep,19LESIZE,188, , ,40, , , , ,1 LESIZE,189, , ,40, , , , ,1 LESIZE,129, , ,25, , , , ,1 LESIZE,113, , ,25, , , , ,1 LESIZE,3, , ,50, , , , ,1 LESIZE,132, , ,10, , , , ,1 LESIZE,17, , ,10, , , , ,1 LESIZE,18, , ,10, , , , ,1 LESIZE,19, , ,10, , , , ,1 LESIZE,20

9、, , ,10, , , , ,1 LESIZE,21, , ,10, , , , ,1 LESIZE,22, , ,10, , , , ,1 LESIZE,23, , ,10, , , , ,1 LESIZE,24, , ,10, , , , ,1 LESIZE,26, , ,10, , , , ,1 LESIZE,25, , ,10, , , , ,1 LESIZE,53, , ,10, , , , ,1 LESIZE,13, , ,10, , , , ,1 LESIZE,14, , ,10, , , , ,1 Vsweep,5Vsweep,1Vsweep,2Vsweep,11Vsweep

10、,6Vsweep,16Vsweep,17VPLOT FINISH ! 施加约束条件和载荷 /SOLcsys,0nsel,s,loc,z,0cm,down,nodeD,down,ALLallsasel,s,81,1asel,a,82,1asel,a,89,1asel,a,90,1cm,zhoutao,areaSFA,zhoutao,1,PRES,10!alls!vplotsolvefini2、实体模型和网格模型结果 图1.1 显示材料性质的实体模型 图1.2 显示材料性质的的网格模型3、计算结果云图显示图1.3 Von-Mises应力计算结果云图图1.4应变计算结果云图4、安全系数计算图1.5

11、底板Von Mises应力云图及可能危险位置图1.6 轴承套Von Mises应力云图及可能危险位置图1.7 挡板Von Mises应力云图及可能危险位置表1.1 不同材料部分最大应力及安全系数材料编号123屈服应力（MPa）325650250最大应力（MPa）61.012190.439143.114安全系数5.3263.4161.747图1.5-1.7所示，红色圆圈内表示每一种材料部件最大应力位置，即最容易产生破坏的位置。由由表1.1可知：三种不同材料在题目所述的载荷下有不同的安全系数，材料一部分安全系数最大，n1=5.326，材料二部分次之，n2=3.416，材料三部分安全系数最小，n3=

12、1.747。因此相比之下，材料三部分最容易发生破坏，且破坏位置位于材料二部分和材料三部分接触位置，如图1.7圆圈内部分所示。5、模型网格节点数、单元数统计图1.5模型网格节点数、单元数统计由统计结果可知：节点数：95501 单元数：8432038第二题 井眼稳定性与井筒压力地层岩石弹性模量6.54×103MPa，泊松比：0.22，内聚力：18MPa，内摩擦角：28°。地应力s1=50MPa，s2=(1050)MPa，井筒内压p=(040)MPa，井眼直径220mm，取井眼周围约15倍井眼直径范围进行分析(即取3.0m)，用弹性和弹塑性有限元分析井眼径向位移变化、井壁周围的应

13、力变化和井眼稳定性讨论。1、APDL编程（以s1=50MPa，s2=30MPa，p=20MPa为例）fini/clear/filname,201421000029_Lixin_Wellbore-Stability/PREP7 !进入前处理ET,1,PLANE82 !定义单元类型为PLANE82KEYOPT,1,3,2 !选择平面应变分析MPTEMP, MPTEMP,1,0MPDATA,EX,1,6.54e9 !定义杨氏模量MPDATA,PRXY,1,0.22 !定义泊松比TB,DP,1, !选择Drucker-Prager准则TBMODIF,1,1,18e6 !定义内聚力TBMODIF,1,2

14、,28 !定义摩擦角TBMODIF,1,3, !定义其它单元属性RECTNG,0,1.5,0,1.5, !建立矩形面RYL4, , ,0.11 !建立圆形面ASBA,1,2 !布尔减操作NUMCMP,ALL !压缩单元编号LESIZE,2, , ,20, , , , ,1 !设置2号线段分段数LESIZE,1, , ,20, , , , ,1 !设置1号线段分段数LESIZE,5, , ,20,0.3, , , ,1!设置5号线分段数LESIZE,4, , ,20,0.3, , , ,1!设置4号线分段数LESIZE,3, , ,40, , , , ,1 !设置4号线段分段数LCCAT,2,1

15、 !连接2号和1号线段AMESH,1 !对面进行网格划分LREFINE,3, , ,1,1,1,1 !进行网格加密FINISH !完成该处理层/SOL !进入求解处理器DL,4, ,SYMM !施加对称约束DL,5, ,SYMM !施加对称约束ALLSEL,ALL !选择所有实体SFL,2,PRES,50e6,!2号线上施加载荷50e6PaSFL,1,PRES,30e6,!1号线上施加载荷30e6PaSFL,3,PRES,20e6,!3号线上施加载荷20e6PaSOLVE !求解FINISH !完成该处理层/POST1 !进入后处理CSYS,1 !转换坐标AVPRIN,0AVRES,2,/EF

16、ACET,1 LAYER,0 FORCE,TOTAL /PLOPTS,INFO,1 /PLOPTS,MINM,0 !输图处理/TRIAD,OFF !关闭坐标系显示PLNSOL, U,X, 0,1.0 !查看径向位移云图/EXPAND,4,POLAR,HALF,90 !图形拓展PLNSOL, S,EQV, 0,1.0 !查看von Mises应力云图2、网格划分结果图2.1 网格划分（单元类型为PLANE82）3、云图输出结果（1）s1=50MPa，s2=10MPa，p=0MPa图2.2 径向位移云图 图2.3 1/4井壁处位移随弧长变化图 图2.4 von Mises应力云图 图2.5 1/4

17、井壁处应力随弧长变化图 （2）s1=50MPa，s2=10MPa，p=10 MPa图2.6 径向位移云图 图2.7 1/4井壁处位移随弧长变化图 图2.8 von Mises应力云图 图2.9 1/4井壁处应力随弧长变化图 （3）s1=50MPa，s2=10MPa，p=20 MPa图2.10 径向位移云图 图2.11 1/4井壁处位移随弧长变化图 图2.12 von Mises应力云图 图2.13 1/4井壁处应力随弧长变化图（4）s1=50MPa，s2=10MPa，p=30 MPa图2.14 径向位移云图 图2.15 1/4井壁处位移随弧长变化图 图2.16 von Mises应力云图 图2

18、.17 1/4井壁处应力随弧长变化图（5）s1=50MPa，s2=10MPa，p=40 MPa图2.18 径向位移云图 图2.19 1/4井壁处位移随弧长变化图 图2.20 von Mises应力云图 图2.21 1/4井壁处应力随弧长变化图（6）s1=50MPa，s2=20MPa，p=0 MPa图2.22 径向位移云图 2.23 1/4井壁处位移随弧长变化图 图2.24 von Mises应力云图 图2.25 1/4井壁处应力随弧长变化图（7）s1=50MPa，s2=20MPa，p=10 MPa图2.26 径向位移云图 图2.27 1/4井壁处位移随弧长变化图图2.28 von Mises应

19、力云图 图2.29 1/4井壁处应力随弧长变化图（8）s1=50MPa，s2=20MPa，p=20sMPa图2.30 径向位移云图 图2.31 1/4井壁处位移随弧长变化图 图2.32 von Mises应力云图 图2.33 1/4井壁处应力随弧长变化图（9）s1=50MPa，s2=20MPa，p=30MPa图2.34 径向位移云图 图2.35 1/4井壁处位移随弧长变化图图2.36 von Mises应力云图 图2.37 1/4井壁处应力随弧长变化图（10）s1=50MPa，s2=20MPa，p=40MPa图2.38 径向位移云图 图2.39 1/4井壁处位移随弧长变化图图2.40 von

20、Mises应力云图 图2.41 1/4井壁处应力随弧长变化图（11）s1=50MPa，s2=30MPa，p=0 MPa图2.42 径向位移云图 图2.43 1/4井壁处位移随弧长变化图图2.44 von Mises应力云图 图2.45 1/4井壁处应力随弧长变化图（12）s1=50MPa，s2=30MPa，p=10 MPa图2.46 径向位移云图 图2.47 1/4井壁处位移随弧长变化图图2.48 von Mises应力云图 图2.49 1/4井壁处应力随弧长变化图（13）s1=50MPa，s2=30MPa，p=20MPa图2.50 径向位移云图 图2.51 1/4井壁处位移随弧长变化图图2.

21、52 von Mises应力云图 图2.53 1/4井壁处应力随弧长变化图（14）s1=50MPa，s2=30MPa，p=30MPa图2.54 径向位移云图 图2.55 1/4井壁处位移随弧长变化图 图2.56 von Mises应力云图 图2.57 1/4井壁处应力随弧长变化图（15）s1=50MPa，s2=30MPa，p=40 MPa图2.58 径向位移云图 图2.59 1/4井壁处位移随弧长变化图图2.60 von Mises应力云图 图2.61 1/4井壁处应力随弧长变化图（16）s1=50MPa，s2=40MPa，p=0MPa图2.62 径向位移云图 图2.63 1/4井壁处位移随弧

22、长变化图 图2.64 von Mises应力云图 图2.65 1/4井壁处应力随弧长变化图（17）s1=50MPa，s2=40MPa，p=10MPa图2.66 径向位移云图 图2.67 1/4井壁处位移随弧长变化图 图2.68 von Mises应力云图 图2.69 1/4井壁处应力随弧长变化图（18）s1=50MPa，s2=40MPa，p=20MPa图2.70 径向位移云图 图2.71 1/4井壁处位移随弧长变化图 图2.72 von Mises应力云图 图2.73 1/4井壁处应力随弧长变化图（19）s1=50MPa，s2=40MPa，p=30MPa图2.74 径向位移云图 图2.75 1

23、/4井壁处位移随弧长变化图 图2.76 von Mises应力云图 图2.77 1/4井壁处应力随弧长变化图（20）s1=50MPa，s2=40MPa，p=40MPa图2.78 径向位移云图 图2.79 1/4井壁处位移随弧长变化图 图2.80 von Mises应力云图 图2.81 1/4井壁处应力随弧长变化图（21）s1=50MPa，s2=50MPa，p=0MPa图2.82 径向位移云图 图2.83 1/4井壁处位移随弧长变化图 图2.84 von Mises应力云图 图2.85 1/4井壁处应力随弧长变化图（22）s1=50MPa，s2=50MPa，p=10MPa图2.86 径向位移云图

24、 图2.87 1/4井壁处位移随弧长变化图 图2.88 von Mises应力云图 图2.89 1/4井壁处应力随弧长变化图（23）s1=50MPa，s2=50MPa，p=20MPa图2.90 径向位移云图 图2.91 1/4井壁处位移随弧长变化图图2.92 von Mises应力云图 图2.93 1/4井壁处应力随弧长变化图 （24）s1=50MPa，s2=50MPa，p=30MPa图2.94 径向位移云图 图2.95 1/4井壁处位移随弧长变化图 图2.96 von Mises应力云图 图2.97 1/4井壁处应力随弧长变化图（25）s1=50MPa，s2=50MPa，p=40MPa图2.

25、98 径向位移云图 图2.99 1/4井壁处位移随弧长变化图图2.100 von Mises应力云图 图2.101 1/4井壁处应力随弧长变化图4、井壁稳定性分析 图2.102 不同情况井壁处最大应力变化线 图2.103 不同情况井壁处最大位移变化曲线由上述的不同情况下井壁岩石的有限元分析，结合不同情况下井壁处最大应力变化曲线（图101）以及不同情况下井壁处最大位移变化曲线（图102）可以得出：（1）在s1=50MPa以及井筒内的有效液柱压力一定时，随着水平主应力s2不断变大，井壁处岩石所受到的最大von Mises应力整体呈减小趋势，同时井壁处的最大径向位移由正值变为负值，即井壁处的岩石的受

26、力状态逐渐从受拉变为受压。同时可以看出最大和最小水平主应力之间的差值越大，井壁岩石所受的应力也就越大，井壁就越容易失稳。（2）在s1=50MPa以及水平主应力s2一定时，随着有效液柱压力p的不断增大，井壁处岩石所受到的最大von Mises应力整体呈减小趋势，由此可以看出，钻井过程中的有效液柱压力对于维持井壁稳定具有重要作用，当有效液柱压力处于较低范围值时（钻井过程中小于地层压力从而处于欠平衡状态），单从von Mises应力分析可以得出井壁越容易失稳。但是随着有效液柱压力的增大，井壁围岩则易发生张性破裂，从而产生井漏。从上述分析可以看出，最大水平主地应力、最小水平主地应力以及有效液柱压力的之

27、间差值较大时，井壁在一定程度上就越容易失稳。因此在钻井过程中，需要分析所钻地层地应力状态，保持有效液柱压力在合理的范围之内。第三题 钢板受力分析钢板：E=2.1*105MPa，v=0.3，材料屈服应力：650MPa，厚度=25mm；P=30MPa1.求孔眼附近sx的值(用路径曲线分析)，钢板内sx、sy、txy以及Von Mises的应力场分布，形成Word文件分析报告。2.P约为多少时，钢板开始屈服，钢板內最大Von Mises应力发生在什么位置？孔眼内壁附件的位移变化以及各种应力变化等值线云图，形成分析报告。1、APDL编程fini/clear/filname,201421000029_L

28、ixin_Gangban/PREP7 !进入前处理ET,1,PLANE42 !定义单元类型为PLANE42KEYOPT,1,3,3 !选择平面应力分析R,1,25, !设置实常数 厚度25mmMP,EX,1,2.1e5 !定义杨氏模量MP,PRXY,1,0.3 !设置泊松比RECTNG,0,150,0,150, !建立矩形面CYL4, , ,13 !建立圆形面ASBA,1,2!布尔减操作NUMCMP,ALL !压缩编号LESIZE,5, , ,20,0.1, , , ,1!设置5号线分段数LESIZE,4, , ,20,0.1, , , ,1!设置4号线分段数LESIZE,1, , ,20,

29、, , , ,1!设置1号线分段数LESIZE,2, , ,20, , , , ,1!设置2号线分段数LESIZE,3, , ,40, , , , ,1!设置3号线分段数LCCAT,2,1 !连接2号和1号线段AMESH,1 !对面进行网格划分LREFINE,3, , ,1,1,1,1 !进行网格加密FINISH !完成该处理层/SOL !进入求解处理器DL,5, ,SYMM !施加对称约束DL,4, ,SYMM !施加对称约束SFL,1,PRES,213.62, !1号线上施加载荷SOLVE !求解FINISH !完成该处理层/PLOPTS,INFO,1/PLOPTS,MINM,0 !输图处

30、理/TRIAD,OFF !转换坐标/EXPAND,4,POLAR,HALF,90 !图形拓展/POST1 !进入后处理PLNSOL, S,X, 0,1.0 !查看x方向应力云图PLNSOL, S,Y, 0,1.0 !查看y方向应力云图PLNSOL, S,XY, 0,1.0!查看切应力云图PLNSOL, S,EQV, 0,1.0!查看von Mises云图2、网格划分结果图3.1 网格划分（单元类型为PLANE42）3、云图输出结果（1）P=30MPa钢板受力分析 图3.2 sx应力云图 图3.3 1/4孔壁处sx随弧长变化图 图3.4 sy应力云图 图3.5 1/4孔壁处sy随弧长变化图 图3

31、.6 txy应力云图 图3.7 1/4孔壁处txy随弧长变化图 图3.8 von Mises应力云图 图3.9 1/4孔壁处von Mises应力随弧长变化图（2）钢板开始屈服受力分析图3.10 sx应力云图 图3.11 1/4孔壁处sx随弧长变化图 图3.12 sy应力云图 图3.13 1/4孔壁处sy随弧长变化图 图3.14 txy应力云图 图3.15 1/4孔壁处txy随弧长变化图 图3.16 von Mises应力云图 图3.17 1/4孔壁处von Mises应力随弧长变化图图3.18 x方向位移变化图 图3.19 1/4孔壁处x方向位移随弧长变化图图3.20 y方向位移变化图 图3

32、.21 1/4孔壁处y方向位移随弧长变化图 图3.22 合位移变化图 图3.23 1/4孔壁处合位移随弧长变化图 从图3.16可以看出，当p=213.623MPa时，孔壁处所受到的von Mises应力为649.988MPa，而材料的屈服应力为650MPa，因此当p=213.623MPa左右时，钢板开始屈服，屈服位置在孔眼的上部和下部。第四题 桁架受力分析结构均为圆形钢材桁架。直径5mm，材料特性：E=2.0*1011N/m2，v=0.3。求反力、杆的内力、节点位移，哪些杆受压、哪些杆受拉？20N10N3.6m1.8m3.6m3.6m2.4m1、APDL编程fini/clear/filname

33、,201421000029_Lixin_Hangjia/PREP7 !进入前处理ET,1,LINK1 !选择单元LINK1R,1,1.96e-5 , !定义实常数MP,EX,1,2.0e11 !定义杨氏模量MP,PRXY,1,0.3 !定义泊松比N,1,0, !定义节点1坐标为(0,0)N,2,1, !定义节点2坐标为(1,0)CSYS,1 !切换到柱坐标系NGEN,4,1,2,60,1 !旋转并复制节点CSYS,0 !切换到直角坐标系E,1,2 !把1、2号节点相连构成单元E,1,3 !把1、3号节点相连构成单元E,1,4 !把1、4号节点相连构成单元E,1,5 !把1、5号节点相连构成单元

34、E,2,3 !把2、3号节点相连构成单元E,3,4 !把3、4号节点相连构成单元E,4,5 !把4、5号节点相连构成单元FINISH !完成该处理层/SOL !进入求解处理器D,5, , , , , ,ALL, , , , , !对5号节点施加约束D,2, , , , , ,UY, , , , , !对2号节点施加约束F,1,FY,-10000 !在1节点加集中外力向下10000NSOLVE !求解FINISH !完成该处理层/PLOPTS,INFO,1 /PLOPTS,MINM,0 !输图处理/TRIAD,OFF !关闭坐标系显示/POST1 !进入后处理PLNSOL, U,X, 0,1.0

35、!查看x方向位移云图PLNSOL, U,Y, 0,1.0!查看y方向位移云图PLNSOL, U,SUM,0,1.0!查看合位移云图PRRSOL,F !输出节点支反力ETABLE,axial-f,SMISC,1!定义轴向力单元表ETABLE,axial-p,LS, 1!定义轴向应力单元表PLETAB,AXIAL-F,NOAV!显示轴向力图/PNUM,NODE,1/PNUM,SVAL,1!显示节点号和单元应力数值PLETAB,AXIAL-P,NOAV !显示轴向应力图/PNUM,NODE,1/PNUM,SVAL,1!显示节点号和单元应力数值PRETAB,AXIAL-F,AXIAL-P!显示单元轴向

36、力与轴向应力表2、约束结果图4.1 约束结果3、云图输出结果图4.2 变形前后x轴方向位移对比图4.3 变形前后y轴方向位移对比图4.4 变形前后合位移对比4、计算结果（1）节点支反力表4-1 节点支反力情况NODEFx（N）FY（N）205000.05-0.45475E-125000.0通过ANSYS计算可以得出，节点2处的支反力Fy2=5000N，方向沿y轴正向；节点5处的支反力Fx5=-0.45475E-12 N，方向沿x轴负向，Fy2=5000N，方向沿y轴正向。（2）节点位移节点1节点5的位移计算结果如表1所示。表1 节点位移计算结果节点x轴方向位移（m）y轴方向位移（m)合位移（m

37、）10.736E-03-0.0046770.00473420.00147300.00147330.434E-18-0.0025510.0025514-0.001473-0.0025510.0029465000（3）杆的轴向力与轴向应力图4.5 杆单元轴向力图图4.6 杆单元轴向应力图图4.7 各杆单元应力图表2 各杆单元轴向力和轴向应力计算结果单元编号轴向力（AXIAL-F）/N轴向应力（AXIAL-F）/Pa杆受力状态12886.80.14728E+09拉2-5773.5-0.29457E+09压3-5773.5-0.29457E+09压4-5773.5-0.29457E+09压52886.

38、80.14728E+09拉65773.50.29457E+09拉75773.50.29457E+09拉第五题 管道层流和紊流层流模拟：1.2D管道，模拟雷诺数为90的层流问题,速度为1 in/s 。2.获得解后，增加速度到50 in/s，研究其对流场的影响获得新的解。3.增加管道长度30in，研究流场变化。尺寸和特性入口长度4 in入口高度1 in过渡区长度2 in出口高度2.5 in初始出口长度4 in扩展出口长度30 in空气密度1.21x10-7 lbf-s2/in4空气粘度2.642x10-9 lbf-s/in2入口速度1 in/sec*出口压力0 psi*计算50步后重起动1、APD

39、L编程Finish/clear/filname, 201421000029_Lixin_Fluid_CFD/PREP7ET,1,FLUID141 !定义单元类型FLDATA2,ITER,EXEC,50, !设置迭代步数50FLDATA12,PROP,DENS,0!设置空气流体FLDATA13,VARY,DENS,0设置空气流体FLDATA12,PROP,VISC,0FLDATA13,VARY,VISC,0FLDATA12,PROP,COND,0FLDATA13,VARY,COND,0FLDATA12,PROP,SPHT,0FLDATA13,VARY,SPHT,0FLDATA7,PROT,DEN

40、S,CONSTANTFLDATA8,NOMI,DENS,1.21e-7,FLDATA9,COF1,DENS,0 FLDATA10,COF2,DENS,0FLDATA11,COF3,DENS,0FLDATA7,PROT,VISC,CONSTANTFLDATA8,NOMI,VISC,2.642e-9,FLDATA9,COF1,VISC,0 FLDATA10,COF2,VISC,0FLDATA11,COF3,VISC,0FLDATA12,PROP,IVIS FLDATA7,PROT,COND,CONSTANT FLDATA8,NOMI,COND,-1, FLDATA9,COF1,COND,0 FLD

41、ATA10,COF2,COND,0FLDATA11,COF3,COND,0FLDATA7,PROT,SPHT,CONSTANT FLDATA8,NOMI,SPHT,-1, FLDATA9,COF1,SPHT,0 FLDATA10,COF2,SPHT,0FLDATA11,COF3,SPHT,0RECTNG,0,4,0,1, !创建流入区域RECTNG,6,10,0,2.5, !创建流出区域L2TAN,-3,-7A,2,5,8,3 !创建过渡区域！RECTNG,10,40,0,2.5,!-第3问执行，（增加管道长度）！NUMMRG,ALL, , , ,LOW!-第3问执行LESIZE,3, , ,15,-2, , , ,1!创建网格模式LESIZE,1, , ,15,-2, , , ,1LESIZE,9, , ,12,1, , , ,1 LES