下载文档原格式
热采井套管应力的弹塑性分析我国稠油资源是油气资源的重要组成部分,在总资源中的比例达到20%-25%。
稠油有多种开采方式,其中以蒸汽驱、蒸汽吞吐和火烧油层三种方法为主,我国一般以蒸汽吞吐方法为主。
稠油开采中套管损坏现象时常发生,影响稠油开采进程,对油田经济造成严重影响。
在稠油的开采过程中注入高温高压蒸汽,套管产生热应力,多个生产周期循环作用,套管冷热交替还可能产生残余应力,加快套管的破坏程度。
为了更好的优化稠油井的井筒结构、高效安全的进行稠油活动,准确预测稠油开采过程中井筒内温度场和热载荷的分布,分析每个生产周期内套管的应力分布是具有重要意义的。
本文针对热采井套管每个生产周期内的工况,对热采井套管进行了以下几方面分析:在分析热采井采油方式和其井筒结构基本特点的基础上,以热采井井筒-地层危险段为研究对象,分析了一个生产周期内的耦合温度场,并建立温度场的非稳态数学模型;得到套管、水泥环及地层三者的温升函数,分析了稠油热采井内套管、水泥环及地层在整个生产周期内的温度变化。
在热弹性力学和均匀内外压作用的厚壁筒理论的基础上,建立了井筒-地层耦合应力场的力学模型,在弹性条件下分析了套管、水泥环及地层的应力情况,并计算出三者的应力和位移表达式;根据某油田热采井的常规参数,分析出单个生产周期内三种不同钢级套管的应力分布及屈服情况。
在屈雷斯卡屈服条件的基础上,建立了井筒-地层耦合应力场弹塑性力学模型,分别推导了套管、水泥环弹性区及塑性区应力及位移的分布公式;根据套管、水泥环、地层三者紧密接触,满足位移协调方程,得到套管-水泥环-地层系统的接触压力及屈服半径,并求得一个生产周期后套管、水泥环发生的不可恢复塑性变形,重复计算得到套管的安全生产周期。
运用VB编程软件编写了套管-水泥环-地层系统弹性及弹塑性应力分析软件,并结合某油田热采常规参数计算了三种钢级套管的生产周期。
稠油热采井套管HYP110H的开发【摘要】本文利用有限元模拟的方法，借鉴国内外研制热采井套管所采用的先进技术，从材料耐热稳定性方面出发，考虑材料和扣型的抗热变形性能，选择高温下密封性能良好且低应力的特殊扣型，开发了优质热采井用套管。
【关键词】稠油开采热采井套管热载荷套管强度目前，国内外稠油开采应用最广泛且效益较高的方法是注蒸汽单井吞吐的开采工艺，这是一种高温高压下开采的工艺，注蒸汽的平均温度在320℃左右，有的高达375℃；压力在10MPa～15MPa之间。
套管在这种井中工作环境非常恶劣，套管柱受力较普通井套管柱更为复杂。
在向稠油热采井注入蒸汽过程中，井中套管温度剧烈变化，发生热胀冷缩，但套管受到地锚和井壁的约束，不能自由伸长或缩短，因此套管内部产生很大的热应力。
而在停注采油的过程中，由于金属材料的松弛和蠕变作用，套管又受到较大的拉应力作用。
根据理论计算，这种拉-压应力非常大，严重时超过700MPa[1]。
套管周期性承受如此高的拉-压应力作用，是热采井套管损坏的根本原因。
另外，普通的API偏梯形螺纹在抵抗热应力方面存在很大缺陷，且偏梯形螺纹的密封性能很差，即使选用优质的耐高温螺纹脂，其密封性能也不可靠[2]。
据有关资料统计，稠油热采井的套损率平均为30%以上，局部区块达到70%，列油田套管损坏率之首[3]。
本文借鉴国内外研制热采井套管所采用的先进技术，从材料耐热稳定性方面出发，考虑材料和扣型的抗热变形性能，选择高温下密封性能良好且低应力的特殊扣型，开发了优质热采井用套管。
1 研制内容1.1 耐热钢种研制热采井中注入蒸汽过程中，套管本身随温度上升而强度降低，而温度升高将导致套管膨胀引起的压缩应力和弯曲应力使套管强度降低，且封闭环空流体膨胀导致套管内压增加而引起套管强度降低；文献研究表明，温度对金属材料线膨胀系数、钢材弹性模量等性能参数都有很大的影响。
随温度升高，钢材弹性模量呈直线减小，线膨胀系数则先升高，到达一定温度后降低。
收稿日期：2005-12-02；修回日期：2006-06-12基金项目：该项目受油气藏地质及开发工程国家重点实验室基金项目（P L N 0141）资助。
作者简介：曾玉强（1979-），在读博士研究生，油气田开发。
地址：（610500）四川成都新都区西南石油大学研究生院博士2004级，电话：028-*******，83030404，E -m a i l ：z e n g3407@163.c o m 注蒸汽开发稠油油藏中的井筒热损失分析曾玉强1，李晓平1，陈礼2，鲁小会3，王琴1（1西南石油大学2辽河油田分公司欢喜岭采油厂工艺研究所3新疆油田分公司重油开发公司）曾玉强等.注蒸汽开发稠油油藏中的井筒热损失分析.钻采工艺，2006，29（4）：44-46摘要：注蒸汽开发稠油油藏过程中，为了预测沿井深和随时间变化的蒸汽温度分布、干度分布和蒸汽压力分布、套管和地层温度分布，以及焖井、开井生产过程中温度、压力的变化，首先必须建立注入过程中的井筒-地层温度分布模型。
而这一模型建立的基础是井筒热损失的计算分析。
文中应用热传递基本理论，通过井筒内能量守恒定理建立起注蒸汽井井筒热损失计算预测模型，重点分析了井筒总传热系数U t o 和注汽速度对井筒热损失和井底蒸汽干度的影响，对注蒸汽开发稠油油藏有一定的指导作用。
关键词：稠油油藏；注蒸汽；井筒；热损失；温度中图分类号：T E345文献标识码：A文章编号：1006-768X （2006）04-0044-03稠油热采时，井筒热损失分析是一项很重要的内容。
分析和预测蒸汽温度分布、干度分布和蒸汽压力分布，都以此为基础。
蒸汽注入阶段中，湿蒸汽沿井筒向下流动为气液两相流，在建立注入过程中的温度分布模型中，有以下的基本假设条件：（1）湿蒸汽沿油管向下稳定流动，不对外做功；（2）油管、隔热管和套管同心；（3）使用耐热封隔器，蒸汽不窜入油套环空；（4）不考虑接箍、扶正器等的影响；（5）初始地层温度按地温梯度分布；（6）井筒为一维径向稳定传热，地层为一维稳定传热。
双管热采井口的抗冲蚀磨损性能裴峻峰;王兵;殷舜时;秦志坚;邵金海;徐延海【期刊名称】《中国石油大学学报（自然科学版）》【年(卷),期】2018(042)003【摘要】针对蒸汽辅助重力泄油(SAGD)双管热采井口装置内冲蚀磨损严重问题,研究冲蚀磨损的影响因素.采用正交试验设计对SAGD井口装置进行流体动力学模拟和优化设计,应用ANSYS软件中FSI(fluid-structure interaction)单向流固耦合模块,对正交优化设计试验进行轴向速度、径向速度、全域最大磨损率及壁面剪切应力模拟,判定井口装置中冲蚀破坏的危险区域.结果表明:入口端速度为0.1 m/s,颗粒粒度为0.125 mm,过渡圆角半径为8 mm,管壁壁厚为21.8 mm,管材种类取35CrMo基体与WC-12Co涂层.优化方案使多相流体减轻了对井口管壁的冲蚀磨损作用.【总页数】9页(P122-130)【作者】裴峻峰;王兵;殷舜时;秦志坚;邵金海;徐延海【作者单位】常州大学机械工程学院,江苏常州213016;常州大学机械工程学院,江苏常州213016;江苏金石机械集团有限公司,江苏金湖211600;江苏金石机械集团有限公司,江苏金湖211600;江苏金石机械集团有限公司,江苏金湖211600;江苏金石机械集团有限公司,江苏金湖211600【正文语种】中文【中图分类】TE832【相关文献】1.双管热采井口的有限元应力分析及结构优化 [J], 胡承军;蒙永立;黄晓东;向河;刘清友2.含凹坑缺陷稠油热采井口用四通管的冲蚀数值模拟研究 [J], 邱福寿;王国涛;彭辉;孙亮;傅登伟;刘岑凡3.新疆石油管理局采油工艺研究院双管热采井口销往海外油田 [J], 裘新农4.基于试验和CFD模拟的稠油热采井口四通管冲蚀规律分析 [J], 王国涛;朱丽云;刘岑凡;王振波;李安俊;石景元5.SAGD双管热采井口装置研制与应用 [J], 彭辉因版权原因，仅展示原文概要，查看原文内容请购买。
热采井口装置全尺寸热-流-固强度计算及安全评定张喜庆【摘要】采用数值模拟的方法,对KR14-337型热采井口装置进行全尺寸热-流-固强度计算,分析其流场、温度场以及应力场.基于弹性分析评价准则,对其主要的承压部件(上法兰、阀体、油管头四通、小四通、螺栓)进行了安全评定.结果表明:该热采井口装置中压力波动并不大,可以看作为均匀内压.阀体和小四通位置存在较大的温度梯度,强度校核时必须考虑不均匀温度场所产生的热应力.主要承压部件中,上法兰和螺栓处的安全裕量较小,在以后的在役检测评价过程中应给予特别关注,以保证设备长周期安全稳定运行.【期刊名称】《石油矿场机械》【年(卷),期】2019(048)003【总页数】7页(P12-18)【关键词】热采;井口装置;应力分析;安全评定【作者】张喜庆【作者单位】胜利油田胜机石油装备有限公司,山东东营257067【正文语种】中文【中图分类】TE931.102采油井口装置主要由采油树、油管头、套管头及相关配件组成，其性能的优劣直接决定了油气井能否长期安全稳定地高效生产[1]。
随着石油工业的不断发展，钻井深度日益加深，采油井口装置的工作环境日益恶劣。
特别是热采井口装置，除了本身结构复杂，易出现结构不连续而引起的应力集中之外，服役过程中高温、高压水蒸汽所带来的流体冲击作用和热载荷作用，使得其强度评定工作难度显著增加。
龚凯[2]、华剑[3]和李斌[4]等分别对采油井口装置中闸阀、阀体和大四通等承压部件进行了有限元应力分析和结构优化。
邱福寿等[5]对热采井口装置六通进行了三维应力分析，但并未考虑热载荷和流场冲击的作用。
胡成军等[6-7]采用热力顺次耦合的方法研究了热采井口装置大四通结构的热应力。
综合现有研究，一方面，流体冲击作用和不均匀热载荷作用对装置强度安全性能的影响尚不明确；另一方面，目前的研究多局限于采油井口装置的局部部件，并未从整个系统的角度进行深入的分析。
本文基于专业流场分析和结构分析软件，采用数值计算的方法，对全尺寸热采井口装置进行流场、温度场以及应力场分析，并基于弹性分析评价准则对几个关键承压部件进行强度评定。
稠油热采特种井口工艺技术研究稠油油藏开采进入中后期，为了进一步提高原油采收率，进行了多种新工艺、新技术的现场推广应用，为了适应新形势的需求，辽河油田科研人员针对当前开发中应用的蒸汽驱、SAGD、水平井等新技术自主研制了一系列特种井口，从注汽、生产、测试几个方面，形成了具有辽河油田特色的工艺技术，达到了国内领先水平。
该系列特种井口的研制成功为蒸汽驱、SAGD、水平井技术的顺利进行提供了技术保障。
标签：蒸汽驱；特种井口；堵塞装置针对普通井口在新工艺，新技术应用中所存在的问题，我们在特种井口研制过程中，首先进行理论研究计算，其次对井口材料进行了优选，第三和生产实际相结合，对井口结构进行了优化设计。
使之既能够满足新工艺新技术的特殊要求，又能够满足现场常规作业要求，做到操作简单，安全实用。
1 蒸汽驱注汽井口研制1.1 蒸汽驱注汽井口结构蒸汽驱注汽井口主要由大四通、异径法兰、隔热管悬挂器及热采闸阀四大部分组成。
热采闸阀的阀体材料采用30CrMo，密封圈用螺纹连接在阀体上，采用楔式钢性单闸板结构设计，更适用于高温状态，便于补偿材料热胀冷缩，使阀门性能可靠、寿命长。
1.2 蒸汽驱注汽井口特点蒸汽驱注汽井口采用油管悬挂器和异径法兰采用分体式设计，可实现不压井作业。
油管悬挂器与大四通之间采用多级密封。
油管悬挂器内堵塞器设计，可直接从测试阀门下入堵塞器堵头堵住油管悬挂器，实现蒸汽驱注汽过程中采油树的带压更换。
2 双管注汽井口研制2.1 双管注汽井口结构该双管注汽井口主要由大四通、双管四通、隔热管悬挂器及热采闸阀四大部分组成。
内外管双流道互不干扰，在井口处不需采用任何密封件就使两流道有效分割开。
2.2 双管注汽井口特点双管注汽井口装置是水平井注汽开采地面控制设备的重要组成部分，是实现水平井段的脚跟和脚尖的均匀注汽的关键技术设备。
水平井双管注汽井口在高温、高压注汽过程中具有可靠的长效密封与耐冲击性能，其结构设计新颖，体积小，操作方便，闸门开关灵活。
稠油热采井口装置的热应力分析胡承军①彭辉蒙永利新疆油田公司采油工艺研究院，834000新疆克拉玛依摘要论述了新疆油田浅层稠油注蒸汽开采的抽油井采油树主要承压部件在高温状态下力学分析及有限元计算过程。
利用Ansys软件对实体建立三维模型，应用有限元的方法分析、计算应力分布，掌握其应力应变的变化规律，为今后持续改进设计研发提供科学的理论计算依据和方法。
主题词井口装置热应力稠油热采建模有限元分析1前言目前，国内油田老区浅层稠油注蒸汽生产抽油井、注汽井等普遍使用的采油树在高温状态下，温度对金属材料的机械性能有很大影响。
本文针对油田注蒸汽开采的采油树的使用工况，对采油树主要承压部件进行了有限元分析计算，掌握其在高温状态下的应力应变的变化规律，找到在不同工况条件下的应力分布情况，并在此基础上提出了一些建议，对采油树今后的改进设计和生产提供了较为科学的理论依据和方法。
2热—结构耦合分析的基本原理和方法步骤浅层稠油抽油井采油树在高温状态下，材料的塑性指标随温度升高而增大，强度指标随温度升高而下降。
其中，强度指标影响材料的许用应力的大小，因此需要结合传热学的知识计算采油树的温度分布，确定采油树在该温度分布下的等效应力分布，使其计算更符合实际工况，从而得到更准确的计算结果；其次，当物体各点的温度分布不同时，物体每一部分都将由于温度的升高或降低而趋于膨胀或收缩，但是由于弹性体所受的外在约束及各个部分之间的相互约束，这种膨胀或收缩并不能自由地发生，于是产生了由于各部分变形大小不同而引起的温度应力，因此也有必要计算出采油树的温度分布，从而计算采油树的温度应力。
首先要对采油树主要承压部件大四通进行热分析。
热分析是热—结构耦合分析的基础，目的是用于计算部件的温度分布及其他热物理参数，原理是基于能量守恒原理的平衡方程、用有限元法计算各点的温度。
它包括稳态传热和瞬态传热两种过程，其中稳态传热用于研究稳定的热载荷对系统或者部件的影响，通过稳态分析可以了解由于稳定的热载荷引起的温度、热梯度等参数。
瞬态分析用于计算一个系统的随时间变化的温度场及其他热参数。
热分析过程表述为以下3个步骤。
2.1建模利用Pro/E建立三维实体几何模型，由于Ansys 软件直接调用Pro/E建立的三维实体几何模型难度较大，采用的方式是将Pro/E建立的三维实体几何模型转化为Igs通用三维几何模型通用格式，从而实现在Ansys软件中的调用。
完成三维实体几何模型的建立后，还需要指定材料的物理属性和单元类型。
在稳态传热分析中，仅需要指定导热系数；在瞬态传热工程中，还需要指定比热和密度等参数。
然后通过制定的单元类型划分网格，完成建模。
2.2施加载荷施加载荷包括分析有关的边界温度分布，空气和水蒸汽的对流系数，也就是指定分析的边界条件。
2.3求解及结果后处理Ansys将求解结果放入*.rth文件中，进行热结构耦合分析时，可将其代入进行计算，其主要包括节点温度、单元热流密度、节点及单元的热梯度、单元热流率、节点的反作用热流率等，其结果可通过Ansys 软件的POST1进行结果的显示和处理，可分别通过①作者简介：高级工程师，1994-07毕业于天津大学机械系彩色云图、矢量图和列表对结果进行显示。
热—结构耦合分析可以采用直接法和间接法两种方法。
直接法是使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应力分析的结果。
这种方法一般适用于热与结构耦合是双向的情况下，即热分析结果影响结构分析结果，同时，结构分析结果引起的大变形、接触等问题又影响热分析结果。
间接法首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中。
这种方法可以使用热分析的功能和结构分析的功能，考虑到了结构由于在加热和冷却的情况下发生膨胀和收缩，由于各部分膨胀和收缩程度不同或者受到了限制而引起热应力的影响。
其中对于瞬时变化的温度梯度，只需要找到温度梯度最大的时间点即可。
由于本文研究的采油树主要承压部件热与结构的耦合是单向的，即热分析结果虽然影响结构分析的结果，但是由于没有出现大变形、接触等问题，结构变形不会影响热分析的结果。
因此，采用间接法进行处理。
通过对工作状态的分析，大四通的热传导过程是一个缓慢达到热平衡的过程。
由于其过程的缓慢性，可将其动态的热传导过程简化为一个初始阶段和最后平衡阶段的综合。
初始阶段由初始通入高温水蒸气阶段开始，大四通热传导刚刚开始还没有达到平衡状态，其主要特点是温度梯度较大，由弹塑性力学知识可知，其温度应力最大。
此时，温度应力是热量对大四通应力分布产生影响的主要因素。
由于大四通材料具有较好的导热性，在平衡状态温度梯度很小，此时温度分布基本处处相等。
在这个状态下，温度对材料性能的影响是大四通应力分布在热作用下变化的主要因素。
由于初始阶段属于一个瞬态过程，其应力分布仅具有参考价值，因此，本文的分析主要以平衡状态下应力分布为主。
以下采用Ansys软件进行大四通的热-结构耦合分析。
3计算分析过程及结果3．1实体三维模型的建立过程在用Pro/E软件进行零件三维建模时，必需先把表征这个零件的二维图形设计出来，然后将这个二维图形按照某种方式生成符合要求的三维形体。
本文利用双管壳体和大四通的原始图纸，首先在草绘器（Sketcher）中完成三维形体的二维截面，然后将这个截面按一定的方式运动而生成实体。
Pro/E提供了4种基本运动方式将面生成实体：拉伸（Extrude）、旋转（Revolve）、混成（Blend）和扫描（Sweep）。
建立的三维实体模型见图1。
3．2大四通的力学分析及有限元计算采油树大四通的工作参数和性能参数如下：设计压力：0～21MPa；工作温度：-40～350℃；材料：ZG30CrMo（σs=420MPa）。
采油树大四通结构为前后、左右对称的结构。
依据模型简化理论，取结构整体的1/8作为计算对象。
本文在计算中采用精度较高的空间四面体十节点单元，根据有限元网格划分原则，应用Pro/E分析软件，自动把估计应力梯度大的部位，如孔壁、倒角、相贯线等处的单元尽量划分得细一些。
大四通边界条件处理：由对称性规定XOY平面上不允许有Z方向位移，YOZ平面上不允许有X方向位移，XOZ平面上不允许有Y方向位移。
大四通的设计额定工作压力为21MPa，作用在1/4大四通上的5个载荷分别为：21MPa内压，M30螺栓拉力，钢圈顶紧力，M24螺栓拉力，上法兰压力。
下面介绍有关载荷的计算。
安装时，由《阀门设计计算手册》知，钢圈在预紧状态下需要的压紧力按下式计算：F G=3.14DG by；式中D G—钢圈压紧力作用中心圆直径，mm；b—钢圈有效密封宽度，b=ω，ω为椭圆形钢圈的环厚，查标准，取为11.1mm；y—钢圈比压，查手册取为150.3MPa.通过对工作过程的分析，大四通的热传导过程是一个缓慢达到热平衡的过程。
由于其过程的缓慢性，可将其动态的热传导过程简化为一个初始阶段和最后平衡阶段的综合。
初始阶段由初始通入高温水蒸汽阶段开始，大四通热传导刚刚开始还没有达到平衡状态，其主要特点是温度梯度较大，由弹塑性力学知识可知，其温度应力最大，此时，温度应力是热对大四通图1大四通三维实体模型图图5大四通应力分布图NODAL SOLUTION STRP=1SUB=1TIME=1SRQV (AVG)DMX=.663669SMN=14.077SMX=727.45∧NSYS JUL 10200422∶07∶3414.077172.604331.131489.659648.18693.34251.868410.395568.922727.45图6大四通稳定温度场分布NODAL SOLUTION STEP=1SUB=1TIME=1TRMP (AVG)RSYS=0SMN=263.007SMX=269.281∧NSYS JUL 10200422∶56∶00263.007264.401265.795267.19268.584263.704265.098266.493267.887269.281应力分布产生影响的主要因素。
由于大四通材料较好的导热性，在平衡状态，温度梯度很小，此时温度分布基本处处相等，在这个状态下，温度对材料性能的影响是大四通应力分布在热作用下变化的主要因素。
3．3大四通的Ansys 有限元分析结果将Pro/E 建立的大四通模型导出为Igs 文件，在Ansys 中调入，制定单元类型为solid97热分析单元，采用自适应网格划分，网格划分结果如图2所示。
3.3.1大四通热传导初始状态温度分布在初始阶段，给定内壁与外壁的壁面温度，从而计算出大四通温度梯度最大时的应力分布情况。
图3是按照第一边界条件施加在大四通上的温度分布示意图，其边界条件为外壁温度T w =300℃，内壁温度T n =20℃。
由于温度梯度对温度应力有很大影响，因此，图4显示了温度梯度的分布图。
其最大温度梯度集中在内壁和外壁倒角附近。
由此，在原有载荷不变的前提下，引入温度载荷，其应力分布如图5所示，由于温度梯度较大，因此产生较大的温度应力，出现在内外壁上。
由图5可知，在瞬时温度场下，大四通的应力高达631.864MPa ，高于材料的屈服强度，但由于此过程很短，不足于引起材料的破坏，因此可以忽略。
3.3.2大四通热传导稳态温度分布在稳定阶段，给定内壁与外壁的流体温度和对流换热系数，从而计算出大四通稳态温度分布情况。
内壁流体温度T b w=350℃，对流换热系数βw =6000W/m 3K ；外壁流体温度T bn =20℃，对流换热系数βw =800W/m 2K ；由此得到大四通稳定温度场分布，如图6所示。
从图上可见，温度场稳定以后，当介质温度为350℃时，大四通的温度处处相同，其温度为289.33℃。
得到温度场分布后，就可以在原有载荷不变的前提下，引入温度载荷，计算大四通在高温下的应力，其应力分布如图7所示。
由图可以看出，其应力分布与常温下大四通的应力分布趋势相同，这是由于不存在温度应力的影响。
稳定温度场分布下的最大应力比不考虑温度场的应力大5MPa ，这是由于所施加的稳定图3大四通温度场分布图NODAL SOLUTION STEP=1SUB=1TIME=1TEMP （AVG ）RSYS=0SMN=20SMX=350∧NSYS JUL9200412∶03∶3256.667130203.333276.6673502093.333166.667240313.333图4大四通温度梯度分布图NODAL SOLUTION SUB=1TOME=1TGSUM (AVG)RSYS=0SMN=.622E-03SMX=28.766SUM=.002151∧NSYS JUL 10200414∶47∶12.622E -03 6.39312.78519.17725.573.1979.58915.98122.37328.766TGSUM=28.72MN图7大四通施加稳定温度载荷的应力分布NODAL SOLUTION STEP=2SUB=1TIME=2SEQV (AVG)DMX=1.184SMN=263.007SMX=251.627∧NSYS JUL 10200422∶56∶00.88970356.609112.328168.048223.76728.749140.188195.907251.627图2Ansys 划分大四通的网格图∧NSYS JUL10200413∶51∶58ELEMENTS表1大四通有限元分析结果分析软件名称Ansys 有限元分析软件727.45MPa 249.712MPa瞬时温度场无温度最大应力最大应力263.007～269.281℃251.629MPa 稳定温度场温度范围稳定场最大应力温度在250～300℃左右，ZG30GrMo 的σb 略有下降，因此计算出的等效应力比常温下应力分布略微增加。
