西南石油大学现代CAE作业6.docVIP

现代CAE技术作业 作业一 3D实体结构有限元力学分析 结构尺寸(单位mm)见图，轴承套受压力10MPa，底板底面固定不动。材料1: E=2.0e5MPa,u=0.3，δs=325MPa；材料2:E=2.5e5MPa,u=0.28，δs=650MPa ；材料3: E=1.8e5MPa,u=0.25，δs=250MPa 。 作业二 井眼稳定性与井筒压力 地层岩石弹性模量6.54×103MPa，泊松比：0.22，内聚力：18MPa，内摩擦角：28°。地应力(1=50MPa，(2=(10～50) MPa，井筒内压p=(0～40) MPa，井眼直径220mm，取井眼周围约15倍井眼直径范围进行分析(即取3.0m)，用弹性和弹塑性有限元分析井眼径向位移变化、井壁周围的应力变化和井眼稳定性讨论。 要求： 1. 用APDL编程序。 2. 结果要求Word文件分析报告。 3. 报告中曲线用Excel生成。 APDLFinish /clear /filname,homework1 /prep7 a=1500 !mm b=1500 !mm ro=110 !mm Ex=6.54e3 !MPa ux=0.22 ch=18 !MPa angle=28 p1=50 !MPa p2=10 !MPa p0=5 !MPa ET,1,PLANE42 MP,Ex,1,ex MP,PRXY,1,ux TB,DP,1,,, TBMODIF,1,1,ch TBMODIF,1,2,angle RECTNG,0,a,0,b, PCIRC,ro, ,0,360, ASBA, 1, 2 LESIZE,10,,,10,0.2,,,,1 LESIZE,9,,,10,0.2,,,,1 LESIZE,5,,,10,,,,,1 LESIZE,2,,,10,,,,,1 LESIZE,3,,,10,,,,,1 MSHAPE,0,2D MSHKEY,0 AMESH,3 !/eof /SOL DL,10,,SYMM DL,9,,SYMM SFL,3,PRES,p1 SFL,2,PRES,p2 SFL,5,PRES,p0 solve （1） 当P1=50MPa，P2、P0赋初值分别为10、5MPa时， 1.1 Von Mises应力云图 图1.2 塑性等效应力云图 由图1.2知，在井壁附近局部区域出现塑性变形，表明在P1=50Mpa、P2=10Mpa 、P0=5Mpa应力作用下井壁是不稳定。 1.3井眼径向位移变化云图 图1.3所示上下边界位移最大，为11.258mm。 1.4 P0=5Mpa时的SEPL曲线 图1.4 P0=25Mpa时的SEPL曲线 分析结果表明当两个主应力P1、P2越接近时，井壁愈稳定，即使没有内压存在，井眼仍会保持稳定。P1为50MPa，内压P0取25MPa，当P2为20MPa时井眼是稳定的，而当P2为10MPa时会在井眼周围出现不同程度的应力集中，导致井壁失稳。 因此，井眼失稳主要是由于两个主应力的差异造成的。但也与井眼内液柱产生的内压有很大关系，内压得存在会消除部分由于地应力的各向异性带来的井眼 附近的应力集中。 . 钢板受力分析 钢板：E=2.1105MPa； v=0.3，材料屈服应力：650MPa，厚度=25mm；P=30MPa 1.求孔眼附近(x的值(用路径曲线分析)，钢板内(x、(y、( xy 以及Von Mises的应力场分布，形成Word文件分析报告。 2. P约为多少时，钢板开始屈服，钢板內最大Von Mises应力发生在什么位置？孔眼内壁附件的位移变化以及各种应力变化等值线云图，形成分析报告。 要求： 1.用APDL编程序。 2.结果要求Word文件分析报告。 3.报告中曲线用Excel生成。 APDL程序语言 Finish /clear /filnam,homework2 a=150 !mm b=150 r=13 h=25 Ex=2.1e5 !MPa ux=0.3 p=30 !MPa /PREP7 ET,1,PLANE42 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,ex MPDATA,PRXY,1,,ux RECTNG,0,a,0,b, CYL4, , ,r ASBA,1, 2 LESIZE,10, , ,20,0.2, , , ,1 LESIZE,5, , ,10, , , , ,1 LESIZE,9, , ,20,0.2, , , ,