计算力学课程设计课程设计弯曲段套管抗挤强度有限元分析.doc

《计算力学课程设计课程设计弯曲段套管抗挤强度有限元分析.doc》由会员分享，可在线阅读，更多相关《计算力学课程设计课程设计弯曲段套管抗挤强度有限元分析.doc（21页珍藏版）》请在三一办公上搜索。

1、东 北 石 油 大 学课 程 设 计课 程 计算力学课程设计 题 目 弯曲段套管抗挤强度有限元分析 学 院 机械科学与工程学院 专业班级 学生姓名 学生学号 指导教师 年 月 日东北石油大学课程设计任务书课程 计算力学课程设计题目 弯曲段套管抗挤强度有限元分析专业 工程力学 姓名 学号 主要内容：石油工业中，API（美国石油学会）套管强度计算公式是没有考虑任何缺陷的套管强度计算公式，井下套管的使用都是通过API套管强度标准设计的。但是在实际使用的套管都存在一定的缺陷，如出厂的不圆度、壁厚不均度等，使用过程中套管内壁被钻杆接头磨损和套管在弯曲段时的弯曲。这些因素或多或少地影响着套管抗挤强度，当套

2、管抗挤强度降低到一定程度时就会造成套管损坏，影响油气资源开发的经济效益。研究和分析弯曲段套管的抗挤强度就可以知道套管曲率对套管抗挤强度的影响关系，就可以指导现场套管的设计和选材。套管钢级为N80，屈服强度551.6MPa，泊松比0.3，弹性模量206GPa，外径177.8mm，壁厚13.72mm。套管曲率取2/100m、4/100m、6/100m、8/100m、10/100m，分析计算各种曲率条件下的套管抗挤强度。基本要求：在课程设计期间，巩固有限元理论知识，掌握边界处理方法，能够应用有限元分析软件ANSYS求解工程中的实际问题，了解力学分析软件的前后处理，掌握有限元分析流程。在3周时间内，应

3、用ANSYS软件完成课题题目的有限元分析与计算，提交所设计题目的有限元模型、结果和命令流文件，提交5000字左右论文1份（附录为分析过程命令流）。主要参考资料：1 刘巨保石油设备有限元分析M北京：石油工业出版社，19962 刘扬，刘巨保，罗敏有限元分析及应用M中国电力出版社，20083 罗敏，张强ANSYS应用基础篇M大庆石油学院自编教材，20084 祝效华，余志祥ANSYS高级工程有限元分析范例精选M电子工业出版社，2004完成期限 指导教师 专业负责人 年月日目录第1章 概述11.1 弯曲段套管抗挤强度有限元分析的研究目的和意义11.2 弯曲段套管抗挤强度有限元分析的主要研究内容1第2章

4、理论分析32.1套管抗挤强度分析32.2 SOLID45简介3第3章 偏磨套管抗挤强度有限元分析53.1 问题描述53.2 ANSYS有限元模型建立及求解5结论13第1章 概述1.1 弯曲段套管抗挤强度有限元分析的研究目的和意义套管是油井生产中重要的设施，套管损坏问题己受到国内外的普遍关注。套管柱在大井眼曲率的水平井弯曲井段将产生较大弯曲应力，使得套管抗挤强度降低，可能导致套管发生挤毁破坏。为了得到套管抗挤强度与水平井弯曲段造斜率间的定量关系。由于传统石油工程中，石油套管主要用于钻井过程中和完井后对井壁的支撑，以保证钻井过程的进行和完井后整个油井的正常运行。石油套管是维持油井运行的生命线，由于

5、地质条件不同，井下受力状态复杂，拉、压、弯、扭应力综合作用于管体，这对套管本身的质量提出了较高的要求。一旦套管本身由于某种原因而损坏，可能导致整口井的减产，甚至报废。 基于普氏自然平衡拱理论,并考虑套管与地层的相互作用,对空洞的形态进行了分析得出其油层段应尽量使用高强度、大壁厚的套管。近年来随着水平井钻井技术的发展, 高曲率井段套管柱的设计引起了人们的关注。对于直径大、壁厚相对钻柱来说较小的套管, 如果井眼曲率过高会产生另一类突出的问题, 即弯曲井段的套管会发生较大的径向变形。而下套管作业时，除了要考虑摩阻外，还要考虑套管柱的变形，尤其是弯曲井段内套管的变形有别于其他井段，因此有必要对套管柱在

6、弯曲井段的通过能力进行研究。发生径向变形的套管不仅会因通径的减小影响钻井工具和完井工具的下入, 而且还会因其椭圆度增加降低套管的抗挤强度。大多数油井套管的损坏是由过大的外载荷引起的，其在各种载荷作用下的强度、变形是油井工程中非常受关心的问题之一，由于套管的成本在整个油井成本中占很大的比例，所以要尽量避免油井套管所加载的外载荷超出其所承受的极限，从而保证油井的正常使用寿命。在进行油井套管的强度分析时，通常通过解析法对其进行计算校核，其中难免采用很多经验值和经验公式，因此，用该方法分析得到的结果不够准确。本文中通过大型有限元软件对套管的强度进行分析，可以得到非常精确的结果，对工程应用有非常重要的指

7、导意义。1.2 弯曲段套管抗挤强度有限元分析的主要研究内容分析该问题应该首先从没有缺陷套管的抗挤强度入手，用ANSYS对理想套管建模井求解，得到理想状态无缺陷下套管的抗挤强度。此时的抗挤强度应该和API套管强度公式计算出来的一致，主要用来验证ANSYS模型的正确性。ANSYS在本研究方向具有非常大的潜力，很多用解析法解决不了的问题都可以通过ANSYS方便的求解，其中不仅包括套管抗挤强度的问题，还有抗拉强度等，不但可以求解带有理想状态弯曲段套管的问题，还可以求解有凹槽、壁厚不均匀等初始缺陷的问题。借助ANSYS有限元软件就可以对理想和使用过程中的套管建模和求解，可以得到理想和使用过程中的缺陷对套

8、管抗挤强度的影响关系，避免井下事故的发生。在钻井及修井作业中，套管弯曲是不可避免的，尤其在深井、超深井、大位移井和水平井中更为显著。弯曲后的套管不仅抗挤强度会降低，其他使用性能也会有所下降。本研究以弯曲形式为空心柱形套管进行分析，得到这种弯曲模型的程度和套管抗挤强度之间的关系。接下来分析弯曲套管的抗挤强度，使用过程中套管内壁被钻杆接头磨损和套管在弯曲段时的弯曲。这些因素或多或少地影响着套管抗挤强度，当套管抗挤强度降低到一定程度时就会造成套管损坏，影响油气资源开发的经济效益。 建立空心圆柱形的实体模型，对套管的外表面施加1MPa的外压，接着使其一端受到全约束，另一端受到弯曲载荷的作用，让套管自然

9、就弯曲。接着改动命令流，可以得到不同曲率的套管的应力集中。这个问题涉及到弯曲的问题，属于三维问题，建立三维有限元模型，采用SOLID45单元。本文中只给出套管曲率为6的情况下的GUI实例和APDL程序，其他曲率情况套管可以通过改变APDL程序得到。第2章 理论分析2.1套管抗挤强度分析套管所能承受的最大外挤压力称为套管的抗挤强度。套管抗挤强度是国际上一个比较活跃的研究领域,多年来,API/ISO一直在组织国际合作研究。目前,API/ISO Collapse Sub-Team(美国石油学会/国际标准化组织套管挤毁工作组)正计划修订现行API5C3标准，研究的焦点集中在提高抗挤强度计算公式的精度。

10、针对美国石油协会的套管抗挤强度计算公式与实测套管挤毁压力差距较大的问题，利用统计学方法，详细分析了套管实物挤毁试验检测数据，得到了各因素影响套管抗挤强度的定量指标.研究表明，套管抗挤强度主要由径厚比决定；套管计算屈服外压和失稳外压等其他参数及其交互作用对套管的抗挤强度也有显著影响;套管外径和壁厚数据的变异系数比外径不圆度和壁厚不均度对套管抗挤强度的影响更显著.最后给出了套管抗挤强度统计计算公式,为生产厂和用户提供了预测套管抗挤强度的简单实用的计算方法。API所公布的套管强度数据是套管受相应的单一外载作用时的强度，如抗挤强度是套管仅受外挤压力作用时套管所能承受的最大外挤压力值；套管在井下主要考虑

11、管柱的摩阻力限制。除了要考虑这个以外，还要考虑套管的变形。如果基于刚性模型和弹性模型对套管进行套管柱在井下一般是处于复合外载作用状态（两种及两种以上外载同时作用状态）。在复合外载作用下，套管的强度要发生变化，有时套管的强度增加，有时套管的强度降低。套管抗挤强度计算:在外压力与轴向力或内压力与轴向力作用下，管柱内的应力状态为三向应力状态，三个主应力分别为轴向应力、周向应力和径向应力 ： （1）根据第四强度理论，可得套管在多向应力下的强度条件为：sc ss第四强度理论变为：2.2 SOLID45简介 单元性质：3-D维实体结构单元 单元说明：solid45单元用于构造三维实体结构.单元通过8个节点

12、来定义,每个节点有3个沿着xyz方向平移的自由度.单元具有塑性,蠕变,膨胀,应力强化,大变形和大应变能力。有关该单元的细节参看 ANSYS, 理论参考中的SOLID45部分。类似的单元有适用于各向异性材料的solid64单元。Solid45单元的更高阶单元是solid95。SOLID45假定和限制：（1）体积等于0的单元是不允许的。 单元结点编号可参照图45.1: SOLID45几何描述，面IJKL和 MNOP也可互换。 单元不能扭曲，这样单元就会有两个独立的体。这通常发生在当单元结点编号不当时。 所有单元都必须有8个结点。 可以通过定义重合的K和L、O和P来形成棱柱形单元（参见三角形、棱柱形

13、和四面体单元）。 四面体形状也是允许的。对四面体单元，额外形状被自动删除。SOLID45产品限制对于以下产品，将在上述一般假设和限制的基础上再增加一定的限制：。ANSYS Professional. 不允许有DAMP材性。 不能施加流量体荷载。 唯一允许的特殊性能是应力钢化。 KEYOPT(6)=3不可用。图1 SOLID45单元几何第3章 偏磨套管抗挤强度有限元分析3.1 问题描述套管钢级为N80，屈服强度551.6MPa，泊松比0.3，弹性模量206GPa，外径177.8mm，壁厚13.72mm。套管曲率取2/100m、4/100m、6/100m、8/100m、10/100m，分析计算各种

14、曲率条件下的套管抗挤强度。图2 套管实体示意图3.2 ANSYS有限元模型建立及求解GUI步骤：（1）进入ANSYS程序ANSYS9.0ANSYS product launcherfile managementworking directory（路径）job Name（文件名）run。图3 定义文件名（2）设置计算类型ANSYS主菜单：Preferencesselect StructuralOK（3）选择单元类型ANSYS主菜单：PreprocessorElement TypeAdd/Edit/DeletAddselectSolid Brick 8node 45Okclose（）图4 定义单元

15、类型(4)定义材料参数ANSYS主菜单：PreprocessorMaterial PropsMaterial ModelsStructuallinearElasticIsotropicinput EX：2.06e11，PRXY：0.3OK图5 定义材料属性（5）生成几何单元模型ANSYS主菜单PreprocessormodelingcreateVolumeCylinderHollow Cylinderinput：WP X（0），WPY（0），Rad-1（75.18e-3），Rad-2（88.9e-3），Depth（1）OK。如图1所示图6 建立几何模型建立未弯曲套管实体模型，加上弯曲载荷后就变

16、成为弯曲段套管抗挤压分析实体模型。(6)网格划分打开当前工作坐标，把它移动到关键点4（0 ，0.0889, 0），旋转工作平面绕Y轴旋转90度，用工作平面把实体切成两部分。ANSYS主菜单：PreprocessorModelingOperateBooleansDivideVolume by workplanepick all 把工作平面移动到关键点2（0.0889, 0, 0），旋转工作平面绕X轴旋转90度，用工作平面把实体模型切成4部分。ANSYS主菜单：PreprocessorModelingOperateBooleansDivideVolume by workplanepick all

17、图7 用工作平面切实体模型ANSYS主菜单：PreprocessorMeshingMeshtoolLines选择厚度方向上的8条线OK选择2份OK PreprocessorMeshingMeshtoolLines环向上的16条线OK选择6份OKPreprocessorMeshingMeshtoolLines长度方向上的8条线OK选择20份OK图8 实体分网控制最后如下图设置，点击Sweeppick all分网完成，保存一下分网模型。图8 对所有进行分六面体网格图9 有限元网格模型（7）模型施加约束ANSYS主菜单：SolutionDefine LoadsApplyStructuralDispl

18、acementOn Areas选择Z坐标为0的平面OkLab2：ALL DOFOk（8）模型施加载荷1.施加外载荷：ANSYS主菜单：SolutionDefine LoadsApplyStructuralPressureOn Areas选择外圆周上的四个平面7,8,15,16：Value：Load pressure value：1e6Ok2.施加弯曲载荷：弯曲载荷的施加比较复杂，首先根据套管的曲率来计算弯曲载荷的大小，再平均分配到每一个节点上。由材料力学可以知道： 式中：F是施加在断面上的弯曲载荷，N；E是为套管的材料弹性模量，Pa；I是套管截面惯性矩，是套管的转角，rad；L是套管的长度，m

19、；将参数代入计算得F为5170.137N弯曲载荷施加在断面的每个节点上，首先选择断面上的每一个节点，并且知道断面总共有几个节点。ANSYS主菜单：SelectEntitiesAreasApply选择所要选择的断面OkNodesAttached ToArea ，allOkANSYS主菜单：PlotNodesOk ListNodesO，本实例中的节点数为72，故每个节点的载荷为5170.137/72=71.807。k列出所有的节点的数量ANSYS主菜单：SolutionDefine LoadsApplyStructuralForce/MomentOnNodespick allLab：FY，VALU

20、E：-71.807OK图10 施加弯曲载荷示意图（9）分析计算ANSYS主菜单：SelecteverythingANSYS主菜单：SolutionSolveCurrent LSOk（to close the solve current Load step window）Ok图11 求解后的模型（10）结果显示在提取结果时，根据圣文南原理，取远离施加载荷和约束的位置，所以本例中应去中间位置查看结果。ANSYS主菜单：General PostprocPlot ResultContour PlotNodal SoluSelect：stress intensity SINTOk图12 结果云图显示（1

21、1）结果提取 ANSYS主菜单：General PostprocQuery ResultSubgrid SoluSelect：stress intensity SINTOk（Turn to Query Result window）SelectMaxOk 根据本实例提取最大应力值为23.755MPa，如下图13 所示。（12）退出系统ANSYS主菜单：FileExitSave everythingOk图13 提取结果图3.3 结论分析前面的有限元计算中得到套管曲率为6/100m时，外压为1MPa情况下最大应力强度为23.755MPa。根据套管的最大屈服强度是551.6MPa，所以计算得到套管曲率

22、为6/100m时的抗挤强度为29.156MPa。确定套管抗挤强度的方法和偏磨套管抗挤强度的方法是一样的。当给定外径、壁厚、弹性模量、泊松比时套管受到外压，在线性和小变形内，其最大应力强度和和外压之比是一个无量纲常数K。 （1）根据徐秉业的塑性力学书选取屈服准则，其应力强度的表达式为： （2）根据屈强度限挤毁形式套管抗挤强度的定义（使套管中某一点表达达到材料的屈服强度的外压力值为套管抗挤强度），在求得K值之后，将材料的屈服极限代入上式，即可求得套管抗挤强度： (3)将数值带入(1)和(3)式，可以得到该种套管的抗挤强度为29.156MPa同理可把套管曲率6换为2、4、8、10，即可求出弯曲套管在

23、这些曲率下的最大应力强度和套管抗挤强度。表1 弯曲套管的曲率对最大应力值和抗挤强度的影响套管曲率最大应力强度Mpa套管抗挤强度MPa弯曲载荷F/N单位节点载荷/N211.38358.8711723.37923.936415.36737.3843446.75848.872623.75529.1565170.13770.807826.00423.2696893.51695.7431035.58719.9888616.895119.679结论经过三周的时间的学习，我完成了弯曲套管抗挤强度有限元分析，同时对ANSYS有了一个新的了解。在这我要感谢我的指导老师对我耐心的指导，感谢老师给我的帮助，因为在设

24、计的过程中遇到了很多的困难，所以老师给予我极大的帮助。在设计的过程中，我通过查阅大量有关资料，与同学交流经验和自学，使自己学到了不少知识，也经历了不少艰辛，但收获颇丰。在整个设计中我懂得了许多东西，也培养了我独立工作的能力，树立了对自己工作能力的信心，相信会对今后的学习特别是毕业设计、工作、生活有非常重要的影响。虽然这个设计做的也不太好，但是在设计过程中所学到的东西是这次课程设计的最大收获和财富，使我终身受益。（1）套管曲率增大，弯曲井段套管的径向变形和应力也随之增大，套管弯曲载荷和单位节点载荷都增大，最大应力强度也增大，说明曲率与它们成正比。（2）套管曲率增大，套管的抗挤强度反而减小，说明曲

25、率和套管抗挤强度成反比。（3）套管的最大应力强度和抗挤强度成反比，通过以上研究，对套管曲率下弯曲井段套管柱强度设计提供了理论依据。参考文献1 刘巨保石油设备有限元分析M北京：石油工业出版社，19962 刘扬，刘巨保，罗敏有限元分析及应用M中国电力出版社，20083 罗敏，张强ANSYS应用基础篇M大庆石油学院自编教材，20084 祝效华，余志祥ANSYS高级工程有限元分析范例精选M电子工业出版社，20045 刘鸿文材料力学I第四版高等教育出版社，20086 朱伯芳有限单元法原理与应用M第二版中国水利水电出版社，19977 祝效华，高原水平井弯曲段套管下入可行性分析，20118 李黔，施太和水平

26、井弯曲井段套管径向变形计算，19989 蔡恩宏弯曲井段中套管通过能力分析研究 ，20101011附录分析过程APDL程序方法/PREP7 ET,1,SOLID45MPTEMP, MPTEMP,1,0 MPDATA,EX,1,2.06e11MPDATA,PRXY,1,0.3 CYL4, , ,75.18e-3, ,88.9e-3, ,1 KWPAVE, 4 wprot,0,0,90VSBW, 1KWPAVE, 2 wprot,0,0,-90 wprot,0,90,0FLST,2,2,6,ORDE,2 FITEM,2,2 FITEM,2,-3 VSBW,P51X LPLOT FLST,5,4,4,

27、ORDE,4 FITEM,5,25 FITEM,5,27 FITEM,5,29 FITEM,5,31 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y LESIZE,_Y1, , ,2, , , , ,1 FLST,5,4,4,ORDE,4 FITEM,5,26 FITEM,5,28 FITEM,5,30 FITEM,5,32 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y !* LESIZE,_Y1, , ,2, , , , ,1 FLST,5,8,4,ORDE,2 FITEM,5,9 FITEM,

28、5,-16 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y !* LESIZE,_Y1, , ,6, , , , ,1 !*FLST,5,8,4,ORDE,2 FITEM,5,1 FITEM,5,-8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y !* LESIZE,_Y1, , ,6, , , , ,1FLST,5,8,4,ORDE,2 FITEM,5,17 FITEM,5,-24 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,_Y !* LESIZE

29、,_Y1, , ,20, , , , ,1 FLST,5,4,6,ORDE,3 FITEM,5,1 FITEM,5,4 FITEM,5,-6 CM,_Y,VOLU VSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 FINISH /SOLFLST,2,4,5,ORDE,4 FITEM,2,1 FITEM,2,4 FITEM,2,6 FITEM,2,23 DA,P51X,ALL,FINISH /PREP7 FINISH /SOLFLST,2,4,5,ORDE,4

30、 FITEM,2,7 FITEM,2,-8 FITEM,2,15 FITEM,2,-16 /GO !* SFA,P51X,1,PRES,1e6 FLST,5,4,5,ORDE,4 FITEM,5,2 FITEM,5,5 FITEM,5,21 FITEM,5,24 ASEL,S, , ,P51X NSLA,S,1NPLOT FLST,2,72,1,ORDE,20 FITEM,2,476 FITEM,2,-480FITEM,2,505 FITEM,2,526 FITEM,2,623 FITEM,2,-627FITEM,2,652 FITEM,2,749 FITEM,2,-753FITEM,2,77

31、8 FITEM,2,875 FITEM,2,-879FITEM,2,899 FITEM,2,-912FITEM,2,963 FITEM,2,-974FITEM,2,1025FITEM,2,-1036 FITEM,2,1123FITEM,2,-1132 /GO F,P51X,FY,-71.807ALLSEL,ALL /STATUS,SOLUSOLVE/POST1 PLNSOL, S,INT, 0,1.0东北石油大学课程设计成绩评价表课程名称计算力学课程设计题目名称弯曲段套管抗挤强度有限元分析学生姓名黄保富学号080403240127指导教师姓名张强职称序号评价项目指 标满分评分1工作量、工作态度和出勤率按期圆满的完成了规定的任务，难易程度和工作量符合教学要求，工作努力，遵守纪律，出勤率高，工作作风严谨，善于与他人合作。202课程设计质量课程设计选题合理，计算过程简练准确，分析问题思路清晰，结构严谨，文理通顺，撰写规范，图表完备正确。453创新工作中有创新意识，对前人工作有一些改进或有一定应用价值。54答辩能正确回答指导教师所提出的问题。30总分评语：指导教师： 年 月 日