CN107063648A - 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法 - Google Patents

一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法 Download PDF

Info

Publication number
CN107063648A
CN107063648A CN201710141266.1A CN201710141266A CN107063648A CN 107063648 A CN107063648 A CN 107063648A CN 201710141266 A CN201710141266 A CN 201710141266A CN 107063648 A CN107063648 A CN 107063648A
Authority
CN
China
Prior art keywords
mrow
mtd
msubsup
mtr
msub
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.)
Granted
Application number
CN201710141266.1A
Other languages
English (en)
Other versions
CN107063648B (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.)
CETC North semiconductor technology (Shandong) Co.,Ltd.
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201710141266.1A priority Critical patent/CN107063648B/zh
Publication of CN107063648A publication Critical patent/CN107063648A/zh
Application granted granted Critical
Publication of CN107063648B publication Critical patent/CN107063648B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/08—Testing mechanical properties
    • G01M11/083—Testing mechanical properties by using an optical fiber in contact with the device under test [DUT]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • G01L1/246—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre using integrated gratings, e.g. Bragg gratings

Abstract

一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，步骤如下：一:智能光纤光栅法兰螺栓加载监测装置构建；二：利用安置好的加载监测装置计算相关系数；三：计算载荷；四：使法兰螺栓结构进入预装配状态，此时螺栓上不施加载荷；五：执行第一轮次；六：第二轮次螺栓的目标载荷为七：根据装置显示的法兰螺栓加载最优路径，第三轮次采用反向顺序拧紧的加载顺序；八：通过步骤一构建的智能光纤光栅法兰螺栓加载监测装置监测法兰盘的应力分布，保证均匀性；通过以上步骤，达到了用最小轮数实现法兰盘螺栓加载均匀的要求，解决螺栓装配过程中无规定流程，操作不确定性大，无法保证螺栓加载均匀性的问题。

Description

一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法
技术领域
本发明涉及一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，它可以根据工况的紧密度要求计算螺栓载荷并通过装置以较少的加载轮数控制螺栓加载顺序实现法兰盘各螺栓载荷的均匀性，保证法兰盘的密封性和拧紧效率。适用于法兰螺栓加载等技术领域。
背景技术
获得均匀的预紧载荷是螺栓法兰密封系统实现长期密封的前提。然而，目前关于法兰螺栓加载方式仍无统一的规范可循，工人多根据经验交叉预紧螺栓。由于工业上无法实现同时对所有螺栓直接施加轴向载荷，故普遍采用控制扭矩的方法对螺栓逐个加载。螺栓载荷的均匀性主要由两个因素决定：扭矩转化为螺栓轴向力的精度、依次加载在过程中后加载螺栓对已加载螺栓载荷的影响。前者主要取决于加载工具以及密封组件之间的摩擦力，后者主要取决于螺栓的弹性交互作用。
螺栓间弹性交互作用是影响螺栓载荷分布的重要原因，但是进一步研究发现，合理的螺栓加载方式及加载轮次能够有效地降低螺栓间载荷的离散程度。研究表明，弹性交互作用有两种效果：弹性交互正作用和弹性交互负作用，这两种作用使最终螺栓载荷低于目标载荷，且非均匀分布。交叉加载时每轮载荷分布呈”W”形，顺次加载时每轮载荷分布呈“N”形。
根据基于紧密度的螺栓目标载荷计算以及最优路径寻优方式提出了简化的法兰螺栓拧紧方法。发明人提供了下文所述的加载监测装置的加载监测方法。
发明内容
(1)发明目的
在实际的生产制造过程中，法兰盘的密封性要求很高，而法兰螺栓的加载方法不完善且并不统一。基于气密性要求，本发明提供了一种智能光纤光栅法兰螺栓加载监测装置及方法，它利用智能光纤光栅法兰螺栓监测装置计算螺栓加载过程中同侧负影响参数并完成对螺栓加载顺序的控制，利用路径寻优算法给出最优加载顺序。
(2)技术方案
本发明是一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法。利用装置计算同侧负影响系数，根据“同侧负影响”“对侧正影响”方法建立法兰螺栓加载顺序理论并通过装置进行路径规划，从而确定螺栓法兰盘扭紧方法，达到使装配后螺栓预紧力满足目标值并形成均匀分布的目标。针对弹性交互作用的对侧正影响和同侧负影响采用的优化方法分别为：1.在对螺栓进行对称拧紧后通过一次正向顺序加载和一次反向顺序加载弱化对侧正影响；2.通过计算得出修正系数，抵消螺栓载荷引起的同侧负影响。我们仅以八螺栓法兰盘为例，具体步骤如下:
本说明采用2片尺寸相同的8螺栓法兰结构为例作为被连接件进行试验，法兰结构材料为304钢，螺栓材料为1Cr11Ni2W2MoV。法兰结构剖面见图1。
本发明一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，具体步骤如下：
步骤一:智能光纤光栅法兰螺栓加载监测装置构建：
A、将螺栓打孔见图1，把光纤光栅灌入到螺栓内进行螺栓预紧力辅助监测；在法兰盘内侧两孔之间或法兰盘表面两孔之间贴装光纤光栅应变传感器，共选取Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ、Ⅷ八个位置如图2图3；将光纤光栅栅区用胶粘紧在法兰盘内壁选取的定位点上；
B、光纤光栅应变传感器接口连接光纤光栅解调仪，同时连接安装测量应变配套软件的计算机，实时监控装配过程中各位置应变情况，分析连接件结合面的应力分布情况；
C、根据光纤光栅所获取的波长和材料杨氏模量E＝194GPa计算检测点的应变和应力，计算的结果用于监测应力分布是否均匀。
步骤二：利用安置好的加载监测装置计算相关系数；
A、计算螺栓间的弹性模量系数；本发明将法兰螺栓间弹性相互作用系数矩阵[A]n×n定义如下：
{Fvf}n×1-{Fv0}n×1＝[A]n×n({Ftf}n×1-{Ft0}n×1) (1)
弹性交互作用系数矩阵[A]n×n中，系数aij表示受影响螺栓的载荷变化量与加载螺栓加载量的比值，定义为同侧负影响系数；aij定义如下：
B、利用路径寻优及多次螺栓加载试验，在法兰螺栓加载一轮后算出最优路径，即给出下一轮的加载顺序，并将法兰螺栓加载轮数与紧密度建立关系；
步骤三：计算载荷，设定每一轮次的目标载荷值，将计算的第一轮的目标载荷设定为步骤四的目标载荷；螺栓的目标载荷值得取值范围如图4所示，即在确定安装螺栓载荷后，就可以确定工作垫片应力也就可以预测其是否能到达相应的紧密度等级保留；并根据紧密度要求选择轮数，现以三轮次为例，按照目标载荷的百分比设定三轮次每一轮次的目标载荷；其中Gb，Gs，a为垫片长度，M为设计常数：
WP＝PD(A0m0+AH)＝A0max[Sm1,Sm2]+AHPD (3)
其中：
步骤四：使法兰螺栓结构进入预装配状态，此时螺栓上不施加载荷；在法兰螺栓结构达到预装配状态后，执行步骤五到步骤七，分轮次对法兰螺栓进行拧紧，以使螺栓上预紧力达到目标值；
步骤五：执行第一轮次，首次均采用对称拧紧顺序施加载荷；第一轮轮次初始载荷为0，引入步骤二计算的同侧负影响系数，即：
其中，{F1 vf}8×1为第一轮加载的目标载荷，弹性交互作用系数矩阵[A1]8×8为：
代入计算，求得：
{F1 tf}8×1＝{(1-2a)F1 vf,F1 vf,(1-2a)F1 vf,F1 vf,(1-2a)F1 vf,F1 vf,(1-2a)F1 vf,F1 vf}1×8 (7)
步骤六：第二轮次螺栓的目标载荷为初始载荷为第一轮次加载的目标载荷，即：
根据装置显示的法兰螺栓加载最优路径，第二轮次采用正向顺序拧紧的加载顺序，其弹性交互作用系数矩阵[A2]8×8为：
在第二轮次正向顺序加载过程中，螺栓一在加载前不受到其他螺栓影响，螺栓二到七位置加载前分别受到一次同侧负影响，螺栓八受到两次同侧负影响；因此，有：
代入计算，求得：
步骤七：根据装置显示的法兰螺栓加载最优路径，第三轮次采用反向顺序拧紧的加载顺序，弹性交互作用系数矩阵[A3]8×8＝[A2]8×8，各螺栓受同侧负影响作用顺序与第二轮次相反，因此有：
代入式计算，求得：
步骤八：通过步骤一构建的智能光纤光栅法兰螺栓加载监测装置监测法兰盘的应力分布，保证均匀性。
其中，在步骤二中所述的“路径寻优算法”，是利用动态A*算法进行路径寻优，其作法如下：将八螺栓法兰盘构建成一个8*8的二维矩阵，设一号螺栓为路径的起始点，要求需要遍历所有螺栓并回到一号螺栓，求得完成要求并保证最大载荷最小载荷之差较小的最优路径；其中，每一螺栓加载后载荷变化最大的螺栓设定为A*算法中的障碍物，根据弹性系数构建障碍物的动态方程；
所述的动态A*算法即：
从初始状态经由状态n到目标状态的代价估计＝在状态空间中从初始状态到状态n的实际代价+从状态n到目标状态的最佳路径的估计代价；
其中，在步骤三中所述的“轮次”，是指需要加载到目标载荷所需要进行加载的次数，每一轮中每一螺栓仅加载一次；其中轮数的确定由所设定的紧密度等级确定，对应关系见表2。
通过以上步骤，达到了用最小轮数实现法兰盘螺栓加载均匀的要求，解决螺栓装配过程中无规定流程，操作不确定性大，无法保证螺栓加载均匀性的问题。
(3)本发明的优点
i.本发明是针对制造过程中常见又难监测的法兰螺栓加载过程，提出了一种基于紧密度的法兰螺栓加载方法，该方法可根据工况的紧密度需求计算载荷。
ii.本发明提出的基于紧密度的法兰螺栓加载方法，加载顺序轮次少且加载均匀，对于多螺栓的加载可以提高效率。
附图说明
图1光纤光栅灌注入螺栓剖面图。
图2法兰盘结构剖面图。
图3本发明法兰盘螺栓序号及光纤光栅粘贴位置图。
图4本发明目标螺栓载荷取值范围图。
图5本发明优化装配序列工艺路线图。
图6本发明优化装配序列各轮次预紧力分布图。
图7本发明优化装配序列与传统装配序列装配结果对比。
图8本发明所述方法流程图。
图中序号、符号、代号说明如下：
D:法兰盘外径；
B:法兰盘内径；
C:法兰盘厚度；
K:法兰盘孔距；
d:密封面；
f:密封面；
z:螺栓孔数；
L:螺栓孔径；
y:安全裕量；
a:同侧负影响系数；
Gb,a,Gs：垫片常数；
PD:设计压力；
P*:大气压力，0.1013MPa；
参比质量泄漏率,1mg/s；
Lr:150mm外径垫片的质量泄漏率，mg/s；
M：设计常数；
Tp:紧密性参数；
TPmin:最小紧密性参数；
AH:流体压力作用面积；
A0:垫片面积；
Sm1:设计垫片应力1；
Sm2:设计垫片应力2；
n:螺栓数目；
{Fv0}n×1:螺栓初始载荷列矩阵；
{Fvf}n×1:螺栓目标载荷列矩阵；
{Ft0}n×1:加载螺栓的初始载荷列矩阵；
{Ftf}n×1:螺栓加载载荷列矩阵；
[A]n×n:弹性交互作用系数矩阵；
fi v0:受影响螺栓的初始载荷；
fi vf:受影响螺栓的最终载荷；
加载螺栓的初始载荷；
加载螺栓的最终载荷。
具体实施方式
本发明不限于以下说明的实施方式，还可以应用于其他规格法兰盘，如六螺栓、十二螺栓等。采用2片尺寸相同的8螺栓法兰结构为例作为被连接件进行试验，法兰结构材料为304钢，螺栓材料为1Cr11Ni2W2MoV并如图1所示灌入光纤光栅进行辅助测量。法兰盘的尺寸见表1，光纤光栅贴合位置及法兰盘剖面图见图2图3。
表1法兰结构尺寸表
本发明一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，流程见图8所示，其具体步骤如下：
步骤一：连接法兰盘及相应接口，构件装置。
步骤二：根据表2紧密度等级求得螺栓间的同侧负影响系数a＝-0.15。
表2法兰设计用垫片紧密度分级
步骤三：根据紧密度要求及图4所示的载荷范围计算目标载荷并按照百分比设定3个轮次目标载荷分别为10kN、15kN、20kN。
步骤四：进一步计算修正后各轮次载荷列矩阵{F1 tf}8×1、可得：
{F1tf}8×1＝[13 10 13 10...10]TkN (14)
步骤五～七：根据A*算法计算最优路径，优化装配序列装配工艺路线见图5。根据装置给出的最优路径对法兰螺栓进行装配，各轮次螺栓的预紧力分布情况如表3。
表3优化装配序列各轮次螺栓预紧力分布结果
(kN)
步骤八：以螺栓序号为X轴螺栓载荷为Y轴建立坐标生成图6图7，如图所示采用优化装配序列对法兰螺栓进行装配，各轮次螺栓预紧力分布较为均匀，螺栓间预紧力波动较小且不再具有N型预紧力分布不平衡现象。且计算法兰盘八个螺栓的最大载荷与最小载荷，保证差值在小数点后三位，即验证该方法可以保证法兰盘螺栓加载均匀。装配完成后，各螺栓处预紧力值均达到预期目标值，保证法兰螺栓结构具有足够的预紧力。此外，与顺序拧紧和对称拧紧装配序列相比，各螺栓处预紧力的方差明显小于顺序拧紧和对称拧紧后各螺栓处预紧力的方差值。
以上结论说明本发明基于法兰螺栓装配工艺弹性交互作用提出的智能光纤光栅法兰螺栓加载监测装置及方法具备合理性与可行性。

Claims (2)

1.一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，其特征在于：具体步骤如下：
步骤一:智能光纤光栅法兰螺栓加载监测装置构建：
A、将螺栓打孔，把光纤光栅灌入到螺栓内进行螺栓预紧力辅助监测；在法兰盘内侧两孔之间及法兰盘表面两孔之间贴装光纤光栅应变传感器，共选取Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ、Ⅶ、Ⅷ八个位置；将光纤光栅栅区用胶粘紧在法兰盘内壁选取的定位点上；
B、光纤光栅应变传感器接口连接光纤光栅解调仪，同时连接安装测量应变配套软件的计算机，实时监控装配过程中各位置应变情况，分析连接件结合面的应力分布情况；
C、根据光纤光栅所获取的波长和材料杨氏模量E＝194GPa计算检测点的应变和应力，计算的结果用于监测应力分布是否均匀；
步骤二：利用安置好的加载监测装置计算相关系数；
A、计算螺栓间的弹性模量系数；将法兰螺栓间弹性相互作用系数矩阵
[A]n×n定义如下：
{Fvf}n×1-{Fv0}n×1＝[A]n×n({Ftf}n×1-{Ft0}n×1) (1)
弹性交互作用系数矩阵[A]n×n中，系数aij表示受影响螺栓的载荷变化量与加载螺栓加载量的比值，定义为同侧负影响系数；aij定义如下：
<mrow> <msub> <mi>a</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;Delta;f</mi> <mi>i</mi> </msub> </mrow> <mrow> <msub> <mi>&amp;Delta;f</mi> <mi>j</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>f</mi> <mi>i</mi> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>f</mi> <mi>i</mi> <mrow> <mi>v</mi> <mn>0</mn> </mrow> </msubsup> </mrow> <mrow> <msubsup> <mi>f</mi> <mi>j</mi> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>f</mi> <mi>j</mi> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
B、利用路径寻优及多次螺栓加载试验，在法兰螺栓加载一轮后算出最优路径，即给出下一轮的加载顺序，并将法兰螺栓加载轮数与紧密度建立关系；
步骤三：计算载荷，设定每一轮次的目标载荷值，将计算的第一轮的目标载荷设定为步骤四的目标载荷；螺栓的目标载荷值得取值范围，即在确定安装螺栓载荷后，就能确定工作垫片应力也就能预测其是否能到达相应的紧密度等级保留；并根据紧密度要求选择轮数；现以三轮次为例，按照目标载荷的百分比设定三轮次每一轮次的目标载荷；其中Gb，Gs，a为垫片长度，M为设计常数：
WP＝PD(A0m0+AH)＝A0max[Sm1,Sm2]+AHPD (3)
其中：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>G</mi> <mi>s</mi> </msub> <msup> <mrow> <mo>(</mo> <mn>18.02</mn> <msub> <mi>CP</mi> <mi>T</mi> </msub> <mo>)</mo> </mrow> <mi>a</mi> </msup> <mfrac> <mrow> <mi>lg</mi> <mi> </mi> <msub> <mi>T</mi> <mrow> <mi>p</mi> <mi>min</mi> </mrow> </msub> </mrow> <mrow> <mi>lg</mi> <mi> </mi> <msub> <mi>T</mi> <mi>p</mi> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>T</mi> <mrow> <mi>p</mi> <mi>min</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mi>D</mi> </msub> <mo>/</mo> <msup> <mi>P</mi> <mo>*</mo> </msup> <mo>)</mo> <mo>(</mo> <msubsup> <mi>L</mi> <mi>r</mi> <mo>*</mo> </msubsup> <mo>/</mo> <msub> <mi>L</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> </msqrt> <mo>=</mo> <msqrt> <mrow> <mo>(</mo> <mn>1</mn> <mo>/</mo> <mn>0.1013</mn> <mo>)</mo> <mo>(</mo> <mn>1</mn> <mo>/</mo> <mn>0.002</mn> <mo>&amp;times;</mo> <mn>150</mn> <mo>)</mo> </mrow> </msqrt> <msub> <mi>MP</mi> <mi>D</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mrow> <mi>m</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <msubsup> <mi>T</mi> <mrow> <mi>P</mi> <mi>T</mi> </mrow> <mi>a</mi> </msubsup> <mo>/</mo> <mn>1.5</mn> <mo>-</mo> <msub> <mi>P</mi> <mi>D</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>H</mi> </msub> <mo>/</mo> <msub> <mi>A</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤四：使法兰螺栓结构进入预装配状态，此时螺栓上不施加载荷；在法兰螺栓结构达到预装配状态后，执行步骤五到步骤七，分轮次对法兰螺栓进行拧紧，以使螺栓上预紧力达到目标值；
步骤五：执行第一轮次，首次均采用对称拧紧顺序施加载荷；第一轮轮次初始载荷为0，引入步骤二计算的同侧负影响系数，即：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mn>0</mn> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msubsup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中，为第一轮加载的目标载荷，弹性交互作用系数矩阵[A1]8×8为：
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>A</mi> <mn>1</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
代入计算，求得：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mrow> <mo>{</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>a</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>a</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>a</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mi>a</mi> <mo>)</mo> </mrow> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>,</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>1</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
步骤六：第二轮次螺栓的目标载荷为初始载荷为第一轮次加载的目标载荷，即：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mn>0</mn> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
根据装置显示的法兰螺栓加载最优路径，第二轮次采用正向顺序拧紧的加载顺序，其弹性交互作用系数矩阵[A2]8×8为：
<mrow> <msub> <mrow> <mo>&amp;lsqb;</mo> <msub> <mi>A</mi> <mn>2</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> </mtr> <mtr> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mi>a</mi> </mtd> </mtr> <mtr> <mtd> <mi>a</mi> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mi>a</mi> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
在第二轮次正向顺序加载过程中，螺栓一在加载前不受到其他螺栓影响，二到七位置的螺栓加载前分别受到一次同侧负影响，螺栓八受到两次同侧负影响；因此，有：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
代入计算，求得：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msubsup> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <mfrac> <mrow> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <mn>2</mn> <msubsup> <mi>aF</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>a</mi> <mo>)</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>aF</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>a</mi> <mo>)</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>aF</mi> <mn>1</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> <mtd> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>1</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> <mi>T</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤七：根据装置显示的法兰螺栓加载最优路径，第三轮次采用反向顺序拧紧的加载顺序，弹性交互作用系数矩阵[A3]8×8＝[A2]8×8，各螺栓受同侧负影响作用顺序与第二轮次相反，因此有：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mn>0</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>a</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>F</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
代入式计算，求得：
<mrow> <msub> <mrow> <mo>{</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>t</mi> <mi>f</mi> </mrow> </msubsup> <mo>}</mo> </mrow> <mrow> <mn>8</mn> <mo>&amp;times;</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <msubsup> <mfenced open = "{" close = "}"> <mtable> <mtr> <mtd> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mtd> <mtd> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>a</mi> <mo>)</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>aF</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mi>a</mi> <mo>)</mo> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>aF</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> <mtd> <mfrac> <mrow> <msubsup> <mi>F</mi> <mn>3</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> <mo>+</mo> <mn>2</mn> <msubsup> <mi>aF</mi> <mn>2</mn> <mrow> <mi>v</mi> <mi>f</mi> </mrow> </msubsup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>1</mn> <mo>&amp;times;</mo> <mn>8</mn> </mrow> <mi>T</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤八：通过步骤一构建的智能光纤光栅法兰螺栓加载监测装置监测法兰盘的应力分布，保证均匀性；
通过以上步骤，达到了用最小轮数实现法兰盘螺栓加载均匀的要求，解决螺栓装配过程中无规定流程，操作不确定性大，无法保证螺栓加载均匀性的问题。
2.根据权利要求1所述的一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法，其特征在于：在步骤二中所述的“路径寻优算法”，是利用动态A*算法进行路径寻优，其作法如下：将八螺栓法兰盘构建成一个8*8的二维矩阵，设一号螺栓为路径的起始点，要求需要遍历所有螺栓并回到一号螺栓，求得完成要求并保证最大载荷最小载荷之差较小的最优路径；其中，每一螺栓加载后载荷变化最大的螺栓设定为A*算法中的障碍物，根据弹性系数构建障碍物的动态方程；
所述的动态A*算法即：
从初始状态经由状态n到目标状态的代价估计＝在状态空间中从初始状态到状态n的实际代价+从状态n到目标状态的最佳路径的估计代价。
CN201710141266.1A 2017-03-10 2017-03-10 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法 Active CN107063648B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710141266.1A CN107063648B (zh) 2017-03-10 2017-03-10 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710141266.1A CN107063648B (zh) 2017-03-10 2017-03-10 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法

Publications (2)

Publication Number Publication Date
CN107063648A true CN107063648A (zh) 2017-08-18
CN107063648B CN107063648B (zh) 2019-02-19

Family

ID=59622370

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710141266.1A Active CN107063648B (zh) 2017-03-10 2017-03-10 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法

Country Status (1)

Country Link
CN (1) CN107063648B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109186822A (zh) * 2018-08-02 2019-01-11 大连理工大学 一种基于fbg传感器的螺栓紧固结合面面压检测方法
CN109186823A (zh) * 2018-08-02 2019-01-11 大连理工大学 基于fbg传感器检测螺栓连接结合面面压的标定方法
CN110909425A (zh) * 2019-11-04 2020-03-24 西安交通大学 一种可消除载荷超程的螺栓法兰连接工艺优化设计方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030128924A1 (en) * 2001-12-20 2003-07-10 Optinel Systems, Inc. Wavelength tunable filter device for fiber optic systems
KR20090046500A (ko) * 2007-11-06 2009-05-11 이승재 대공간 구조물의 상시 모니터링을 위한 광섬유 브래그 격자센서（ｆｂｇ）의 부착장치 및 ｆｂｇ를 부착하는 방법
CN102706544A (zh) * 2012-03-28 2012-10-03 上海市电力公司 一种增敏型光纤光栅法兰螺栓松动监测方法及监测装置
KR20130017820A (ko) * 2011-08-12 2013-02-20 (주)엠이엘 텔레콤 광섬유 브라그 격자 내장 파장 가변 모듈
CN104089729A (zh) * 2014-07-23 2014-10-08 上海市建筑科学研究院 一种光纤光栅式结构连接螺栓受力监测传感器及方法
US20140373636A1 (en) * 2013-06-19 2014-12-25 Chun-Chu Yang Synchronous pre-tensionable sensing screw with fiber bragg grating devices
CN104280169A (zh) * 2014-10-17 2015-01-14 武汉理工大学 一种圆环式光纤光栅测力装置及其应用
CN105043741A (zh) * 2015-08-11 2015-11-11 北京航空航天大学 一种考虑高温蠕变的螺栓预紧力退化过程实时检测装置
CN105447263A (zh) * 2015-12-14 2016-03-30 中国航空工业集团公司成都飞机设计研究所 一种用于复杂连接结构非线性分析的螺栓预紧力施加方法
CN106446459A (zh) * 2016-10-26 2017-02-22 中国特种设备检测研究院 一种螺栓连接组件的预紧优化方法
CN106425404A (zh) * 2016-11-24 2017-02-22 芜湖市元山机械制造有限公司 一种法兰螺栓压合工装及其工作方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030128924A1 (en) * 2001-12-20 2003-07-10 Optinel Systems, Inc. Wavelength tunable filter device for fiber optic systems
KR20090046500A (ko) * 2007-11-06 2009-05-11 이승재 대공간 구조물의 상시 모니터링을 위한 광섬유 브래그 격자센서（ｆｂｇ）의 부착장치 및 ｆｂｇ를 부착하는 방법
KR20130017820A (ko) * 2011-08-12 2013-02-20 (주)엠이엘 텔레콤 광섬유 브라그 격자 내장 파장 가변 모듈
CN102706544A (zh) * 2012-03-28 2012-10-03 上海市电力公司 一种增敏型光纤光栅法兰螺栓松动监测方法及监测装置
US20140373636A1 (en) * 2013-06-19 2014-12-25 Chun-Chu Yang Synchronous pre-tensionable sensing screw with fiber bragg grating devices
CN104089729A (zh) * 2014-07-23 2014-10-08 上海市建筑科学研究院 一种光纤光栅式结构连接螺栓受力监测传感器及方法
CN104280169A (zh) * 2014-10-17 2015-01-14 武汉理工大学 一种圆环式光纤光栅测力装置及其应用
CN105043741A (zh) * 2015-08-11 2015-11-11 北京航空航天大学 一种考虑高温蠕变的螺栓预紧力退化过程实时检测装置
CN105447263A (zh) * 2015-12-14 2016-03-30 中国航空工业集团公司成都飞机设计研究所 一种用于复杂连接结构非线性分析的螺栓预紧力施加方法
CN106446459A (zh) * 2016-10-26 2017-02-22 中国特种设备检测研究院 一种螺栓连接组件的预紧优化方法
CN106425404A (zh) * 2016-11-24 2017-02-22 芜湖市元山机械制造有限公司 一种法兰螺栓压合工装及其工作方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109186822A (zh) * 2018-08-02 2019-01-11 大连理工大学 一种基于fbg传感器的螺栓紧固结合面面压检测方法
CN109186823A (zh) * 2018-08-02 2019-01-11 大连理工大学 基于fbg传感器检测螺栓连接结合面面压的标定方法
CN110909425A (zh) * 2019-11-04 2020-03-24 西安交通大学 一种可消除载荷超程的螺栓法兰连接工艺优化设计方法
CN110909425B (zh) * 2019-11-04 2021-09-07 西安交通大学 一种可消除载荷超程的螺栓法兰连接工艺优化设计方法

Also Published As

Publication number Publication date
CN107063648B (zh) 2019-02-19

Similar Documents

Publication Publication Date Title
CN107063648A (zh) 一种智能光纤光栅法兰螺栓加载监测装置的加载监测方法
AU2009330323B2 (en) Structural shape for wind tower members
CN103712865A (zh) 高温部件的蠕变损伤评价方法以及损伤评价系统
JPS6389233A (ja) 接合ボルトおよびその締付力の調整方法
CN106446459A (zh) 一种螺栓连接组件的预紧优化方法
CN110991111B (zh) 一种基于摩擦接触的风电齿轮箱行星架的疲劳计算方法
CN108488254A (zh) 湿式离合器压力响应迟滞补偿方法
CN103399491A (zh) 光伏发电系统光伏组件机理模型参数辨识方法
CN110593969B (zh) 燃气轮机气缸的密封法兰及其设计方法
US20070078627A1 (en) Screw tightening designing method and screw tightening designing apparatus
CN110705140A (zh) 组合荷载作用下长条形基座螺栓群受力分布的确定方法
Husson Friction connections with slotted holes for wind towers
CN110083961B (zh) 汽轮机法兰强度和气密性获取系统
CN210686071U (zh) 燃气轮机气缸的密封法兰
CN206038211U (zh) 一种测量螺栓预紧力的装置
CN112632863A (zh) 多路阀仿真建模方法及多路阀仿真模型
CN109583126B (zh) 一种与时间相关的高温螺栓法兰接头泄漏率预测方法
CN203979745U (zh) 一种膨胀节自平衡式结构
CN103016866A (zh) 一种非匀质垫片的螺栓法兰连接结构设计方法
CN107061893A (zh) 一种特殊大型设备法兰
CN112989662A (zh) 一种风电齿轮箱柔性销结构的有限元计算方法
CN113408024A (zh) 装配式地下结构注浆式榫槽接头抗弯承载能力计算方法
CN207930148U (zh) 一种压力扩散焊接工装
CN110768569A (zh) 一种基于割率的压电陶瓷抗迟滞驱动方法
CN103454102A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211210

Address after: 255400 room 306, 3rd floor, bay hoisting building, No. 1, Fenghuangshan Road, Linzi District, Zibo City, Shandong Province

Patentee after: CETC North semiconductor technology (Shandong) Co.,Ltd.

Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road

Patentee before: BEIHANG University