CN110887448B - 一种基于光纤应变测量的梁结构形态重构方法 - Google Patents
一种基于光纤应变测量的梁结构形态重构方法 Download PDFInfo
- Publication number
- CN110887448B CN110887448B CN201911192971.XA CN201911192971A CN110887448B CN 110887448 B CN110887448 B CN 110887448B CN 201911192971 A CN201911192971 A CN 201911192971A CN 110887448 B CN110887448 B CN 110887448B
- Authority
- CN
- China
- Prior art keywords
- strain
- optical fiber
- section
- axis
- bending
- 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
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 60
- 238000005259 measurement Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005452 bending Methods 0.000 claims abstract description 39
- 238000006073 displacement reaction Methods 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 239000000835 fiber Substances 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 19
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Optical Transform (AREA)
Abstract
本发明公开了一种基于光纤应变测量的梁结构形态重构方法,该方法包括以下步骤:步骤1:获取圆形截面梁末端截面扭转角;并沿圆形截面梁轴线方向设置多个测量点以获取沿梁轴线方向的表面应变;步骤2:计算测量点处对应的梁截面扭转角,结合测量点处所获取的表面应变,进行应变分解,得到梁结构绕垂直于轴向的两个正交方向上弯曲变形产生的表面应变;步骤3:根据得到的正交应变信息,采用Ko位移理论分别计算梁在两个正交方向上的弯曲变形挠度,完成三维形态重构。该方法相比于其他的空间三维曲线重构方法,考虑了扭转变形的影响,在小变形情况下具有计算速度快、精度高的优点。
Description
技术领域
本发明涉及柔性结构的形变测量与形态重构技术研究领域,具体涉及一种基于光纤应变测量的梁结构形态重构方法。
背景技术
柔性结构的形变测量与形态重构是结构受力分析与振动控制的重要基础。在航空航天、土木工程等领域,柔性梁结构是一种比较典型的结构形式,由于其复杂的工作环境,极易发生空间弯曲、扭转等组合变形或振动,从而影响整个***的性能,因此对其进行高精度形变测量和形态重构具有重要意义。在测量中,光纤传感由于其质量轻、体积小、抗电磁干扰、耐腐蚀、对结构影响小等优点,得到了广泛的关注与应用。
目前,基于光纤测量的重构技术大部分基于的是结构的曲率信息,利用坐标变换及曲率递推实现平面或空间曲线的重构。一般而言,此类方法涉及大量的坐标转换,计算量较大,因此在实时监测中比较适合于静态变形或慢变柔性结构,如地下管线、医用内窥镜的形态感知等。而对于存在振动的柔性梁结构,变形监测的实时性要求更高,需要有更快的计算速度,此类方法则往往难以满足要求。另外,对于基材截面是圆形的柔性结构,多采用周向间隔90°或120°排布的多条光纤来获取结构的正交曲率信息,但是在计算中大多不考虑结构的扭转变形对光纤测量的影响,这给重构结果带来了一定的误差;或者虽考虑扭转影响但计算过程比较复杂,结合重构方法本身特点,计算量较大。因此,对于存在空间弯扭组合变形的梁结构,有必要研究一种计算速度快、重构精度高的形态重构方法。
发明内容
本发明的目的是提出一种基于光纤应变测量的梁结构形态重构方法,用于圆形截面梁结构在空间弯曲及扭转组合变形下的静态或动态三维形态重构,采用光纤应变传感器测量梁表面两个正交方向的应变信息,同时通过其他测量手段获取梁的截面扭转角,在此基础上根据应变分解及Ko位移理论实现梁的三维形态重构。
具体包括以下步骤:
步骤1:获取圆形截面梁末端截面扭转角;并沿圆形截面梁轴线方向设置多个测量点以获取沿梁轴线方向的表面应变;
步骤2:计算测量点处对应的梁截面扭转角,结合测量点处所获取的表面应变,进行应变分解,得到梁结构绕垂直于轴向的两个正交方向上弯曲变形产生的表面应变;
步骤3:根据得到的正交应变信息,采用Ko位移理论分别计算梁在两个正交方向上的弯曲变形挠度,完成三维形态重构。
上述的基于光纤应变测量的梁结构形态重构方法,其中,所述多个测量点通过光纤应变传感器a和光纤应变传感器b进行测量;所述光纤应变传感器 a和光纤应变传感器b沿梁轴线方向粘贴于梁结构表面,并在梁截面圆周方向上间隔90°排布。
上述的基于光纤应变测量的梁结构形态重构方法,其中,所述步骤2中,首先通过梁端部的截面扭转角插值计算内部各测量点处的截面扭转角。对于只在端部受扭转外载荷作用的梁结构,一般可通过线性插值计算得到
其中,φ0、φf分别为梁起始端和末端的截面扭转角,φi为光纤应变传感器上第i个测量点处的梁截面扭转角,L为梁的总长度,xi为第i个测量点在梁轴线方向上的坐标;
其次,利用光纤应变传感器上第i个测量点处所获取的表面应变εa.i、εb.i和截面扭转角φi,进行应变分解,获取梁绕垂直于轴线方向的y轴和z轴弯曲产生的表面应变εy.i、εz.i:
其中,
εay.i、εaz.i分别为光纤应变传感器a处绕y轴和z轴的弯曲应变,εby.i、εbz.i分别为光纤应变传感器b处绕y轴和z轴的弯曲应变,εt.i为扭转应变,tanγ为梁的表面扭转角正切,r为圆形截面梁的截面半径。
上述的基于光纤应变测量的梁结构形态重构方法,其中,所述步骤3中,根据Ko位移理论,梁在两个正交方向的弯曲挠度分别计算如下:
y方向:
z方向:
其中,Δl为测量点间距,tanθz.i=dy/dx为梁在oxy平面内的斜率, tanθy.i=dz/dx为梁在oxz平面内的斜率,yi、zi分别为梁在y、z方向上的挠度;
利用各个测量点处的三轴坐标(xi,yi,zi)即可重构出梁结构的三维形态。
所述方法中采用的光纤应变传感器可以用其他应变传感器如压电应变片等来代替,只需保证获取相同测量点相同方向的应变信息即可。
相对于现有技术,本发明具有以下有益效果:
(1)本发明在圆形截面梁的形态重构中考虑了扭转变形的影响,计算精度较高。现有的基于光纤测量的梁形态重构方法,往往直接基于光纤测量应变或曲率进行计算,这必须以结构不发生扭转变形为前提。当结构存在扭转时,光纤位置发生偏转,所测信号中包含了扭转变形的影响。因此,传统的重构计算方法会产生一定误差。本发明对其进行了修正,可以获得更高的重构精度。
(2)本发明采用Ko理论进行重构计算,计算速度较快,可用于结构的在线实时重构。相比于曲率递推等曲线重构方法,Ko理论不涉及复杂的坐标转换,计算量较小,可以直接应用于结构的小弯曲变形情况。而对于梁的大弯曲挠度情况,通常在结构局部范围内也可以进行小变形假设,可通过分段的方式应用Ko理论,以较低的计算量进行重构,保证计算速度。
附图说明
图1为圆形截面梁上的光纤传感器布置方式(变形前);
图2为圆形截面梁上的光纤传感器布置方式(变形后);
图3为圆形截面梁的横截面视图;
图4为圆形截面梁表面扭转角的计算示意图;
图5为数值仿真算例给出的梁结构重构挠度曲线;
图6为梁结构重构挠度曲线的末端局部放大图;
图7为两种方法重构结果的挠度误差曲线对比图。
具体实施方式
以下结合附图通过具体实施例对本发明作进一步的描述,这些实施例仅用于说明本发明,并不是对本发明保护范围的限制。
本发明提出了一种基于光纤应变测量的梁结构形态重构方法,用于圆形截面梁结构在空间弯曲及扭转组合变形下的静态或动态三维形态重构,采用光纤应变传感器测量梁表面两个正交方向的应变信息,同时通过其他测量手段获取梁的截面扭转角,在此基础上根据应变分解及Ko位移理论实现梁的三维形态重构;具体包括以下实施步骤:
步骤1:利用圆形截面梁表面粘贴的光纤应变传感器,获取沿梁轴线方向的表面应变;同时利用布置在梁末端的其他测量装置获取末端截面扭转角。
所采用的光纤应变传感器沿梁轴线方向粘贴于结构表面,共计a、b两根,如图1-图3所示。其中图1、图2给出了两根光纤应变传感器在梁结构上的整体布置方式,沿轴线方向粘贴;图3为梁的横截面图,可以看出光纤应变传感器a、光纤应变传感器b在圆周方向上间隔90°排布,在结构发生扭转变形前分别位于A′、B′点,变形后分别位于A、B点。
梁结构末端截面扭转角φf则可以通过其他多种测角手段获取。
步骤2:根据测量得到的末端截面扭转角插值得到应变测量点处对应的梁截面扭转角,结合光纤应变信息,进行应变分解,得到梁结构绕垂直于轴向的两个正交方向上弯曲变形产生的表面应变。
首先,对于只在端部受到扭转外载荷作用的梁结构,一般可通过线性插值计算各应变测量点处的截面扭转角
其中,φ0、φf分别为测量得到的梁起始端和末端的截面扭转角,一般φ0=0,φi为光纤上第i个应变测量点处的梁截面扭转角,L为梁的总长度,xi为第i个应变测量点在梁轴线方向上的坐标。
其次,利用第i个测量点处的光纤测量应变εa.i、εb.i和截面扭转角φi,进行应变分解。
如图3所示,梁结构发生变形后,光纤应变传感器a、光纤应变传感器b 分别位于截面的A、B点处,所测应变中同时包含了弯曲和扭转产生的应变,为了计算变形挠度,需要从中提取出弯曲应变。在图3中,N-N′为梁的中性面与横截面的交线,表示梁绕N-N′方向发生弯曲,其上弯曲应变为0。曲率矢量沿N-N′方向,可以将其在y轴和z轴两个正交方向上进行分解,得到和分别表示绕y轴和z轴的弯曲变形。因此,光纤所测应变可以进行如下分解
其中εa.i、εb.i为光纤直接测量应变,εay.i、εaz.i分别为光纤应变传感器a处(A 点)绕y轴和z轴的弯曲应变,εby.i、εbz.i分别为光纤应变传感器b处(B点) 绕y轴和z轴的弯曲应变,εt.i为扭转应变。若要求解εay.i、εaz.i、εby.i、εbz.i,需要进一步补充其他方程。
注意,对于圆形截面梁,弯曲应变满足如下关系
其中r为圆形截面梁的截面半径,ρy.i、ρz.i分别为梁绕y轴和绕z轴弯曲的曲率半径。由上式可以推出
另外,对于扭转应变εt.i,可以进行如下估算。如图4所示,取某测量点附近的一小段梁结构,该单元同时存在弯曲与扭转变形,其中由扭转造成的表面扭转角为γ。假设只发生弯曲时,光纤为长度为发生扭转后,光纤为长度为则在小弯曲变形下,扭转应变可以表示为
由于梁结构只在端部受到扭转外载荷作用,表面扭转角γ为常值,故有
联立式(2)和式(4),可以解得
即为光纤应变传感器a、光纤应变传感器b绕y轴和z轴的弯曲应变。
由于计算弯曲挠度需要采用截面上A′和C点的弯曲应变εy.i和εz.i(见图 3),所以还需要进行进一步处理。由于εy.i=r/ρy.i,εz.i=r/ρz.i,故有
步骤3:根据得到的正交应变信息,采用Ko位移理论分别计算梁在两个正交方向上的弯曲变形挠度,完成三维形态重构。
梁在y方向上的变形挠度可以由C点绕z轴的弯曲应变εz.i计算得到,梁在z方向上的变形挠度可以由A′点绕y轴的弯曲应变εy.i计算得到。根据Ko 位移理论,梁在y、z两个方向上的弯曲挠度分别计算如下:
y方向:
z方向:
其中,Δl为测量点间距,tanθz.i=dy/dx为梁在oxy平面内的斜率, tanθy.i=dz/dx为梁在oxz平面内的斜率,yi、zi分别为梁在y、z方向上的挠度。
利用各个测量点处的三轴坐标(xi,yi,zi)即可重构出梁结构的三维形态。
所述方法中采用的光纤应变传感器可以用其他应变传感器如压电应变片等来代替,只要能保证获取相同测量点处相同方向的应变信息即可。
下面以仿真算例来说明本发明的重构效果。
选取长3m,截面直径为0.1m的圆形截面悬臂梁结构,为使其出现弯扭组合变形,分别在x=1m,x=2m,x=3m处施加y向和z向的不同力载荷,并在x=3m处施加绕x方向的扭转力矩载荷。利用Ansys有限元分析软件获取结构变形后的三维形态,作为参考基准。同时,依据节点坐标及环面螺旋线计算公式来模拟光纤测量应变。
分别采用本发明方法与直接采用光纤测量应变(不考虑扭转影响)两种方法进行重构,得到梁结构的挠度曲线,如图5-图6所示。图中实线为本发明方法获得的挠度曲线,虚点线为直接采用光纤测量应变重构得到的挠度曲线,虚划线为Ansys给出的挠度曲线。图6对挠度曲线的末端进行了局部放大,可以明显看出,采用本发明的计算结果与Ansys结果更为接近,而用光纤应变直接计算的结果则存在较大误差。图7给出了两种方法的挠度误差曲线,采用本发明方法的最大绝对误差为1.2×10-3m,最大相对误差为0.53%;而用测量应变直接重构的最大绝对误差为3.7×10-3m,最大相对误差为1.62%。这里的最大相对误差定义为最大绝对误差与最大变形挠度的比值。
以上结果充分表明,当梁结构同时存在弯曲和扭转变形时,在形态重构计算中考虑扭转影响非常必要,验证了本发明方法的有效性。
综上所述,本发明针对同时存在空间弯曲-扭转组合变形形式的圆形截面梁结构提出了一种基于应变测量的三维形态重构方法。利用在梁结构表面布设周向相隔90°的两条轴向应变测量光纤,并结合结构末端转角,可以获取梁结构在两个正交方向上的应变信息,再基于Ko位移理论,可以对梁的三维形态进行精确重构。该方法相比于其他的空间三维曲线重构方法,考虑了扭转变形的影响,在小变形情况下具有计算速度快、精度高的优点。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (3)
1.一种基于光纤应变测量的梁结构形态重构方法,用于圆形截面梁结构在空间弯曲及扭转组合变形下的静态或动态三维形态重构,其特征在于,包括以下步骤:
步骤1:获取圆形截面梁末端截面扭转角;并沿圆形截面梁轴线方向设置多个测量点以获取沿梁轴线方向的表面应变;
所述多个测量点通过光纤应变传感器a和光纤应变传感器b进行测量;所述光纤应变传感器a和光纤应变传感器b沿梁轴线方向粘贴于梁结构表面,并在梁截面圆周方向上间隔90°排布;
步骤2:计算测量点处对应的梁截面扭转角,结合测量点处所获取的表面应变,进行应变分解,得到梁结构绕垂直于轴向的两个正交方向上弯曲变形产生的表面应变;
步骤3:根据得到的正交应变信息,采用Ko位移理论分别计算梁在两个正交方向上的弯曲变形挠度,完成三维形态重构;所述步骤3中,首先通过线性插值计算各测量点处的截面扭转角:
其中,φ0、φf分别为梁起始端和末端的截面扭转角,φi为光纤应变传感器上第i个测量点处的梁截面扭转角,L为梁的总长度,xi为第i个测量点在梁轴线方向上的坐标;
其次,利用光纤应变传感器上第i个测量点处所获取的表面应变εa.i、εb.i和截面扭转角φi,进行应变分解,获取梁绕垂直于轴线方向的y轴和z轴弯曲产生的表面应变εy.i、εz.i:
其中,
εay.i、εaz.i分别为光纤应变传感器a处绕y轴和z轴的弯曲应变,εby.i、εbz.i分别为光纤应变传感器b处绕y轴和z轴的弯曲应变,εt.i为扭转应变,tanγ为梁的表面扭转角正切,r为圆形截面梁的截面半径。
3.如权利要求1所述的基于光纤应变测量的梁结构形态重构方法,其特征在于,所述多个测量点通过应变传感器测量;所述应变传感器包括但不限于光纤应变传感器、压电应变片。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911192971.XA CN110887448B (zh) | 2019-11-28 | 2019-11-28 | 一种基于光纤应变测量的梁结构形态重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911192971.XA CN110887448B (zh) | 2019-11-28 | 2019-11-28 | 一种基于光纤应变测量的梁结构形态重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110887448A CN110887448A (zh) | 2020-03-17 |
CN110887448B true CN110887448B (zh) | 2021-08-20 |
Family
ID=69749283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911192971.XA Active CN110887448B (zh) | 2019-11-28 | 2019-11-28 | 一种基于光纤应变测量的梁结构形态重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110887448B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112287454A (zh) * | 2020-10-28 | 2021-01-29 | 北京天骥空间科技有限公司 | 一种基于应力应变的箭体外形重构技术 |
CN113283086B (zh) * | 2021-05-27 | 2022-04-12 | 华中科技大学 | 一种基于拓展Frenet框架的三维曲线重构方法和*** |
CN113916123B (zh) * | 2021-11-22 | 2022-05-27 | 北京航空航天大学 | 基于应变的大柔性梁式结构三维空间构型确定方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN202101649U (zh) * | 2011-06-23 | 2012-01-04 | 重庆交通大学 | 二维微小扭转角测量*** |
CN103926085B (zh) * | 2014-03-31 | 2016-08-24 | 宝山钢铁股份有限公司 | 一种h型大梁框架弯曲扭转力学性能测试方法及装置 |
CN105136359B (zh) * | 2015-09-11 | 2018-06-01 | 北京航空航天大学 | 基于光纤传感器测得的梁表面应变值计算工作载荷的方法 |
CN105547235B (zh) * | 2016-01-20 | 2019-01-04 | 武汉理工大学 | 变截面梁结构拉压与弯曲复合变形场的测量方法 |
CN108413860B (zh) * | 2018-01-30 | 2019-10-29 | 大连理工大学 | 一种闭口截面薄壁梁约束扭转变形的实时监测方法 |
CN108413861B (zh) * | 2018-01-30 | 2019-10-29 | 大连理工大学 | 一种开口截面薄壁梁约束扭转变形的实时监测方法 |
CN110069832B (zh) * | 2019-04-01 | 2021-07-27 | 南京航空航天大学 | 基于几何坐标转换的非连续板结构变形反演与拼接方法 |
-
2019
- 2019-11-28 CN CN201911192971.XA patent/CN110887448B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110887448A (zh) | 2020-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110887448B (zh) | 一种基于光纤应变测量的梁结构形态重构方法 | |
Floris et al. | Fiber optic shape sensors: A comprehensive review | |
CN108519175B (zh) | 基于布拉格光纤光栅的可变量程的土体压力测量方法 | |
Tang et al. | An improved pendulum method for the determination of the center of gravity and inertia tensor for irregular-shaped bodies | |
CN105403344B (zh) | 管道实时应力的获取方法 | |
CN105043285A (zh) | 一种基于光纤传感器的测量变截面梁变形的方法 | |
CN109100073B (zh) | 一种基于应变反演的六维力传感器及其测量方法 | |
CN107389284B (zh) | 一种基于应变的框架结构弹性变形的测量方法 | |
Floris et al. | Experimental study of the influence of FBG length on optical shape sensor performance | |
CN114754691B (zh) | 直升机桨叶弯扭形态分布式光纤监测与反演方法 | |
CN108692668B (zh) | 基于光纤光栅传感的三维形状检测***及方法 | |
CN110069832A (zh) | 基于几何坐标转换算法的非连续板结构变形反演与拼接方法 | |
JP5602655B2 (ja) | 亀裂サイズ推定方法 | |
CN111879279A (zh) | 高堆石坝的心墙变形监测方法 | |
Todd et al. | A local material basis solution approach to reconstructing the three-dimensional displacement of rod-like structures from strain measurements | |
CN111232239B (zh) | 曲面挠变位移场重构方法、装置及设备 | |
CN104748716B (zh) | 观测三维应变张量变化的方法 | |
CN109724778B (zh) | 基于三芯光纤扭曲补偿的三维位姿恢复方法 | |
CN105547235A (zh) | 变截面梁结构拉压与弯曲复合变形场的测量方法 | |
Yi et al. | An orthogonal curvature fiber bragg grating sensor array for shape reconstruction | |
Jin et al. | Analysis and correction method of axial strain error in multi-core fiber shape sensing | |
CN112632831A (zh) | 一种基于光纤光栅传感器的箭体结构性能识别方法 | |
CN114112358B (zh) | 机械臂的刚度测量方法和刚度测量*** | |
JP4013986B1 (ja) | 固定構造体の曲げ応力測定方法、記録媒体及びコンピュータ | |
CN116164659A (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 |