CN109657334B - 一种超静定结构的受力特性分析方法 - Google Patents

一种超静定结构的受力特性分析方法 Download PDF

Info

Publication number
CN109657334B
CN109657334B CN201811536615.0A CN201811536615A CN109657334B CN 109657334 B CN109657334 B CN 109657334B CN 201811536615 A CN201811536615 A CN 201811536615A CN 109657334 B CN109657334 B CN 109657334B
Authority
CN
China
Prior art keywords
force
packer
basic
constraint
redundant
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811536615.0A
Other languages
English (en)
Other versions
CN109657334A (zh
Inventor
曹银萍
窦益华
秦彦斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Shiyou University
Original Assignee
Xian Shiyou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Shiyou University filed Critical Xian Shiyou University
Priority to CN201811536615.0A priority Critical patent/CN109657334B/zh
Publication of CN109657334A publication Critical patent/CN109657334A/zh
Application granted granted Critical
Publication of CN109657334B publication Critical patent/CN109657334B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Earth Drilling (AREA)

Abstract

一种超静定结构的受力特性分析方法，将封隔器坐封后的四封隔器管柱简化，管柱被四个固定端固定，每个固定端能提供三个反力，但是通常井中的管柱只承受轴向力，因此每个固定端只提供水平反力，从而此结构简化为3次超静定结构，如果去掉超静定结构中的三个多余约束，成为力法的基本结构；若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3，可以得到含有多余未知力的静定结构，求解后各封隔器管柱段轴向力和轴向变形的目的，本发明克服了传统分析方法中分段逐段求解各封隔器之间约束的局限性，可同时计算封隔器之间的每段管柱轴向力，所建立模型更符合多封隔器管柱坐封后的实际受力情况，计算结果也更精确。

Description

一种超静定结构的受力特性分析方法
技术领域
本发明涉及油气井开采技术领域，特别涉及一种超静定结构的受力特性分析方法。
背景技术
油气井开采过程中，为提高油气采收率，增大储层接触面积，常采用直井分层和水平井分段技术。为实现这两项技术，需要采用多封隔器(即多于三个封隔器)。封隔器的作用为定位和密封，从力学角度来说，多封隔器的存在引入了多余约束的问题，即原来的油井管结构变成了超静定结构。
从上世纪60年代开始至今，有很多国内外专家学者对油气井管柱力学进行了研究。1980年，Hammerlindl以双封隔器单一管柱为例，分析了管柱在四种基本效应下的轴向力及变形；1993年，冯建华等人建立了双封隔器复合管柱的轴向力计算模型，并介绍了封隔器受力及其承受压差的计算方法；2014年，吕占国分析了压裂过程中水平井多封隔器管柱的残余轴向力及强度安全性；2015年，练章华等人建立了多封隔器坐封等三种工况下的水平井压裂管柱力学计算模型，并对多级压裂管柱进行了安全性评价；2016年，张智等人建立了四种基本效应下的多封隔器复合管柱受力及变形计算模型；2017年，胡志强等人建立了双封隔器密闭环空压力预测模型，对不同工况下的封隔器管柱进行了安全校核。
现有的技术缺陷在于双封隔器结构中没有考虑超静定的约束条件，仅建立静力平衡方程进行求解，而多封隔器结构从变形协调条件建立超静定结构的变形方程，分段逐段计算管柱的受力和变形。
发明内容
为了克服上述现有技术的不足，本发明的目的在于提供一种超静定结构的受力特性分析方法，具有精确计算各封隔器管柱段轴向力和轴向变形的目的，特点在于充分考虑了前一个封隔器的坐封状态对后一个封隔器的影响，可同时计算封隔器间每段管柱的轴向力，而不是像传统的计算公式分段逐段进行计算管柱的受力和变形。
为了实现上述目的，本发明采用的技术方案是：
一种超静定结构的受力特性分析方法，包括以下步骤；
步骤一：模型简化；
封隔器坐封后的四封隔器管柱简化，管柱被四个固定端固定，其中从左往右第一个固定端为封隔器1对管柱的固定约束，第二个固定端为封隔器2对管柱的固定约束，第三个固定端为封隔器3对管柱的固定约束，第四个固定端为封隔器4对管柱的固定约束；
步骤二：超静定次数的判定；
每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力，由以上简化模型可知此结构为9次超静定结构，但是通常井中的管柱只承受轴向力，因此每个固定端只提供X方向的水平反力，从而此结构简化为3次超静定结构，其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束；
如果去掉超静定结构中的三个多余约束，则原结构变为静定结构，成为力法的基本结构；若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3，可以得到含有多余未知力的静定结构，此静定结构称为力法的基本体系；
步骤三：基本方程；
由于原超静定结构在任何情况下的变形量均为零，因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同，即：
Figure BDA0001907044110000031
由叠加原理可得此三次超静定结构的力法基本方程为：
Figure BDA0001907044110000032
式中：系数δi1在数值上等于基本结构在单位力X1＝1单独作用下沿Xi方向产生的位移；系数δi2在数值上等于基本结构在单位力X2＝1单独作用下沿Xi方向产生的位移；系数δi3在数值上等于基本结构在单位力X3＝1单独作用下沿Xi方向产生的位移；自由项ΔiP为基本结构在外界因素(温度、压力、外载荷等)作用下沿Xi方向产生的位移；i＝1，2，3；
方程(2)中各个系数及自由项的计算公式如下：
Figure BDA0001907044110000033
Figure BDA0001907044110000034
Figure BDA0001907044110000035
i，j＝1，2，3
式中：
Figure BDA0001907044110000036
是基本结构由于Xi＝1作用而产生的内力，N；L为管柱的长度，m；E为管柱的弹性模量，Pa；A为管柱的横截面积，m2；εi为外界因素(温度、压力、外载荷等)使管柱产生的线应变；
步骤四：基本方程的求解
在基本体系中，当X1＝1单独作用时，由平衡条件可得杆1的内力为1，杆2的内力为0；当X2＝1单独作用时，由平衡条件可得杆1和杆2的内力均为1；当X3＝1单独作用时，由平衡条件可得杆1、杆2和杆3的内力均为1，则可得方程中各个系数及自由项的大小；
将各个系数及自由项的计算所得值代入基本方程(2)，则有：
Figure BDA0001907044110000041
超静定结构中各杆的内力可以根据叠加原理用下式计算：
Figure BDA0001907044110000042
解方程组(6)，可得三段管柱的轴向力分别为：
Figure BDA0001907044110000043
由于水平段管柱温度改变量和压力改变量处处相同，即：
ΔT1＝ΔT2＝ΔT3＝ΔT4＝ΔT5＝ΔT
Δpi1＝Δpi2＝Δpi3＝Δpi4＝Δpi5＝Δpi
Δpo1＝Δpo2＝Δpo3＝Δpo4＝Δpo5＝Δpo
因此，由温度效应和鼓胀效应引起的水平段管柱线应变也处处相等，即：
ε温1＝ε温2＝ε温3＝ε温4＝ε温5＝βΔT
Figure BDA0001907044110000051
从而可知，由温度效应和鼓胀效应引起的水平段管柱轴向力处处相等，即：
F温1＝F温2＝F温3＝-EAβΔT
Figure BDA0001907044110000052
本发明的有益效果：
本发明克服了传统分析方法中分段逐段求解各封隔器之间约束的局限性，可同时计算封隔器之间的每段管柱轴向力，所建立模型更符合多封隔器管柱坐封后的实际受力情况，计算结果也更精确。
附图说明
图1为水平井四封隔器管柱简图。
图2为四封隔器管柱简化模型。
图3为基本力系示意图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
以水平井四封隔器单一尺寸管柱为研究对象(如图1所示)，从左往右分别为封隔器1、封隔器2、封隔器3、封隔器4；井口到水平段管柱的垂深为h；封隔器1、2之间的管柱长度为l1，为第一段管柱；封隔器2、3之间的管柱长度为l2，为第二段管柱；封隔器3、4之间的管柱长度为l3，为第三段管柱。
(1)模型简化
封隔器坐封后的四封隔器管柱简化模型如图2所示，由于本专利主要研究四个封隔器之间的三个管柱段所受轴向力的大小，所以封隔器1以前的管柱段以及封隔器4以后的管柱段不再考虑，略去。从图中可以看出，管柱被四个固定端固定，其中从左往右第一个固定端为封隔器1对管柱的固定约束，第二个固定端为封隔器2对管柱的固定约束，第三个固定端为封隔器3对管柱的固定约束，第四个固定端为封隔器4对管柱的固定约束。
(2)超静定次数的判定
因为每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力，而对于图2所示结构只能列出三个平衡方程，所以由以上简化模型可知此结构为9次超静定结构。但是通常井中的管柱只承受轴向力，因此每个固定端只提供X方向的水平反力，从而此结构简化为3次超静定结构，其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束。
如果去掉超静定结构中的三个多余约束，则原结构变为静定结构，成为力法的基本结构；若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3，可以得到含有多余未知力的静定结构，此静定结构称为力法的基本体系，如图3所示。
(3)基本方程
由于原超静定结构在任何情况下的变形量均为零，因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同，即：
Figure BDA0001907044110000061
由叠加原理可得此三次超静定结构的力法基本方程为：
Figure BDA0001907044110000062
式中：系数δi1在数值上等于基本结构在单位力X1＝1单独作用下沿Xi方向产生的位移；系数δi2在数值上等于基本结构在单位力X2＝1单独作用下沿Xi方向产生的位移；系数δi3在数值上等于基本结构在单位力X3＝1单独作用下沿Xi方向产生的位移；自由项ΔiP为基本结构在外界因素(温度、压力、外载荷等)作用下沿Xi方向产生的位移；i＝1，2，3。
方程(2)中各个系数及自由项的计算公式如下：
Figure BDA0001907044110000071
Figure BDA0001907044110000072
Figure BDA0001907044110000073
i，j＝1，2，3
式中：
Figure BDA0001907044110000074
是基本结构由于Xi＝1作用而产生的内力，N；L为管柱的长度，m；E为管柱的弹性模量，Pa；A为管柱的横截面积，m2；εi为外界因素(温度、压力、外载荷等)使管柱产生的线应变。
(4)基本方程的求解
由图3可知，在基本体系中，当X1＝1单独作用时，由平衡条件可得杆1的内力为1，杆2的内力为0；当X2＝1单独作用时，由平衡条件可得杆1和杆2的内力均为1；当X3＝1单独作用时，由平衡条件可得杆1、杆2和杆3的内力均为1。则可得方程中各个系数及自由项的大小如表1和表2所示。
表1系数的计算
Figure BDA0001907044110000075
Figure BDA0001907044110000081
表2自由项的计算
Figure BDA0001907044110000082
将表中各个系数及自由项的计算所得值代入基本方程(2)，则有：
Figure BDA0001907044110000091
超静定结构中各杆的内力可以根据叠加原理用下式计算：
Figure BDA0001907044110000092
解方程组(6)，可得三段管柱的轴向力分别为：
Figure BDA0001907044110000093
由于水平段管柱温度改变量和压力改变量处处相同，即：
ΔT1＝ΔT2＝ΔT3＝ΔT4＝ΔT5＝ΔT
Δpi1＝Δpi2＝Δpi3＝Δpi4＝Δpi5＝Δpi
Δpo1＝Δpo2＝Δpo3＝Δpo4＝Δpo5＝Δpo
因此，由温度效应和鼓胀效应引起的水平段管柱线应变也处处相等，即：
ε温1＝ε温2＝ε温3＝ε温4＝ε温5＝βΔT
Figure BDA0001907044110000094
从而可知，由温度效应和鼓胀效应引起的水平段管柱轴向力处处相等，即：
F温1＝F温2＝F温3＝-EAβΔT
Figure BDA0001907044110000095
由所推导的公式，可以定量求解多封隔器管柱这种超静定结构在不同工况下的温度效应和鼓胀效应，及在这两种效应下各个管柱段间的力学特性，进而分析管柱的强度安全性，为安全施工提供参数。
本专利运用结构力学中的力法对以超静定结构为特点的多封隔器管柱轴向力进行分析，可同时计算封隔器之间的每段管柱轴向力，更符合实际工况下多封隔器管柱的受力情况，计算结果也更精确。
本专利以直井多封隔器单一尺寸管柱为研究对象，考虑温度效应和鼓胀效应对管柱的影响，运用结构力学中的力法对以超静定结构为特点的多封隔器管柱轴向力进行分析，建立用于同时计算封隔器之间的每段管柱轴向力的计算模型，最后以四封隔器为例推导两种效应下的管柱轴向力计算公式，并和传统的计算方法进行比较，计算两者之间的差值。

Claims (1)

1.一种超静定结构的受力特性分析方法，其特征在于，包括以下步骤；
步骤一：模型简化；
封隔器坐封后的四封隔器管柱简化，管柱被四个固定端固定，其中从左往右第一个固定端为封隔器1对管柱的固定约束，第二个固定端为封隔器2对管柱的固定约束，第三个固定端为封隔器3对管柱的固定约束，第四个固定端为封隔器4对管柱的固定约束；
步骤二：超静定次数的判定；
每个固定端能提供水平X方向、竖直Y方向和弯矩M三个反力，由以上简化模型可知此结构为9次超静定结构，但是通常井中的管柱只承受轴向力，因此每个固定端只提供X方向的水平反力，从而此结构简化为3次超静定结构，其中三个多余约束分别为封隔器2对管柱的轴向约束、封隔器3对管柱的轴向约束以及封隔器4对管柱的轴向约束；
如果去掉超静定结构中的三个多余约束，则原结构变为静定结构，成为力法的基本结构；若去掉多余约束后在其相应位置代之以多余未知力X1、X2和X3，可以得到含有多余未知力的静定结构，此静定结构称为力法的基本体系；
步骤三：基本方程；
由于原超静定结构在任何情况下的变形量均为零，因此基本体系沿多余未知力X1、X2和X3方向的位移Δ1、Δ2和Δ3应与原结构相同，即：
Figure FDA0001907044100000011
由叠加原理可得此三次超静定结构的力法基本方程为：
Figure FDA0001907044100000021
式中：系数δi1在数值上等于基本结构在单位力X1＝1单独作用下沿Xi方向产生的位移；系数δi2在数值上等于基本结构在单位力X2＝1单独作用下沿Xi方向产生的位移；系数δi3在数值上等于基本结构在单位力X3＝1单独作用下沿Xi方向产生的位移；自由项ΔiP为基本结构在外界因素作用下沿Xi方向产生的位移；i＝1，2，3；
方程(2)中各个系数及自由项的计算公式如下：
Figure FDA0001907044100000022
Figure FDA0001907044100000023
Figure FDA0001907044100000024
i，j＝1，2，3
式中：
Figure FDA0001907044100000025
是基本结构由于Xi＝1作用而产生的内力，N；L为管柱的长度，m；E为管柱的弹性模量，Pa；A为管柱的横截面积，m2；εi为外界因素使管柱产生的线应变；
步骤四：基本方程的求解
在基本体系中，当X1＝1单独作用时，由平衡条件可得杆1的内力为1，杆2的内力为0；当X2＝1单独作用时，由平衡条件可得杆1和杆2的内力均为1；当X3＝1单独作用时，由平衡条件可得杆1、杆2和杆3的内力均为1，则可得方程中各个系数及自由项的大小。
CN201811536615.0A 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法 Active CN109657334B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811536615.0A CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811536615.0A CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Publications (2)

Publication Number Publication Date
CN109657334A CN109657334A (zh) 2019-04-19
CN109657334B true CN109657334B (zh) 2023-03-31

Family

ID=66114274

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811536615.0A Active CN109657334B (zh) 2018-12-14 2018-12-14 一种超静定结构的受力特性分析方法

Country Status (1)

Country Link
CN (1) CN109657334B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112507503A (zh) * 2019-09-16 2021-03-16 深圳市建筑设计研究总院有限公司 一种次生自反力结构的加载及预反力计算方法
CN111651922A (zh) * 2020-06-03 2020-09-11 重庆大学 力法计算的模拟手算分析方法、系统、装置及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102720185A (zh) * 2012-02-03 2012-10-10 招商局重庆交通科研设计院有限公司 h型抗滑桩桩身结构段分式计算方法
CN104166792A (zh) * 2014-08-06 2014-11-26 中国科学院工程热物理研究所 一种预应力混凝土连续刚构桥温度作用有限元分析方法
CN106548030A (zh) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 一种多封隔器复合管柱受力分析计算方法
WO2017121112A1 (zh) * 2016-01-14 2017-07-20 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种多级水力喷射压裂管柱性能评价方法
CN107563000A (zh) * 2017-08-01 2018-01-09 哈尔滨工程大学 一种内燃机曲轴系统纵扭耦合刚度计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6526819B2 (en) * 2001-02-08 2003-03-04 Weatherford/Lamb, Inc. Method for analyzing a completion system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102720185A (zh) * 2012-02-03 2012-10-10 招商局重庆交通科研设计院有限公司 h型抗滑桩桩身结构段分式计算方法
CN104166792A (zh) * 2014-08-06 2014-11-26 中国科学院工程热物理研究所 一种预应力混凝土连续刚构桥温度作用有限元分析方法
WO2017121112A1 (zh) * 2016-01-14 2017-07-20 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 一种多级水力喷射压裂管柱性能评价方法
CN106548030A (zh) * 2016-11-03 2017-03-29 中国石油天然气股份有限公司 一种多封隔器复合管柱受力分析计算方法
CN107563000A (zh) * 2017-08-01 2018-01-09 哈尔滨工程大学 一种内燃机曲轴系统纵扭耦合刚度计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"结构力学"课程中力法与位移法的对偶关系剖析;王仁华;《中国建设教育》;20170215(第01期);全文 *
超静定结构徐变效应的力法分析方法;沈蒲生等;《铁道科学与工程学报》;20060228(第01期);全文 *
超静定结构的分析与应用;丁瑞娜等;《城市地理》;20160725(第14期);全文 *

Also Published As

Publication number Publication date
CN109657334A (zh) 2019-04-19

Similar Documents

Publication Publication Date Title
CN109657334B (zh) 一种超静定结构的受力特性分析方法
WO2016192077A1 (zh) 一种致密气压裂水平井数值试井模型建立求解方法
CN101539477A (zh) 一种深水顶张式立管涡激振动与疲劳分析的方法
Liu et al. Study on failure mechanism and sealing performance optimization of compression packer
Toumi et al. Approximate Riemann solvers and flux vector splitting schemes for two-phase flow
CN108414351A (zh) 一种套管非均匀外挤能力评价方法
Li et al. The cavity profile of a diaphragm compressor for a hydrogen refueling station
CN109871507A (zh) 正交各向异性煤层裂隙绝对渗透率计算方法
Wang et al. Study on the stress and deformation of a diaphragm compressor cylinder head under extreme conditions
CN103994747A (zh) 梁形结构拉伸弯曲复合变形场的无基准分布式测量方法
Juan et al. Application of linear viscoelastic differential constitutive equation in ABAQUS
Lin et al. Numerical and experimental distribution of stress fields for double shoulder tool joint
Zhang et al. Nonlinear and non-local analytical solution for Darcy–Forchheimer flow through a deformable porous inclusion within a semi-infinite elastic medium
Xuewen et al. Finite element analysis on double-telescopic prop of hydraulic support
Guo et al. The FEA contact analysis of high pressure packers
Feng et al. Finite element analysis of the packer rubbers on sealing process
Ciarlet et al. On the existence of solutions to the generalized Marguerre-von Kármán equations
CN107063993A (zh) 基于应变的高速气流摩阻系数计算方法
CN113191102A (zh) 一种钻井液瞬态波动压力确定方法及系统
CN104808647B (zh) 无线网络环境下抽油机井的故障检测方法
胡蓝 et al. Plastic deformation analysis of thin-walled tube bending under internal pressure
Gu et al. A new method of estimating the equivalent elastic modulus of RCCD
CN209117333U (zh) 一种确定封隔器胶筒坐封力的试验装置
Yadav et al. Weight Optimization of API 6D 12”-150 Class Plug Valve Body by Finite Element Analysis and Experimental Method
CN203716956U (zh) 用于石油测井仪器的压力补偿装置和结构尺寸补偿装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant