Design and Optimization of FBG Implantable Flexible Morphological Sensor to Realize the Intellisense for Displacement

The measurement accuracy of the intelligent flexible morphological sensor based on fiber Bragg grating (FBG) structure was limited in the application of geotechnical engineering and other fields. In order to improve the precision of intellisense for displacement, an FBG implantable flexible morphological sensor was designed in this study, and the classification morphological correction method based on conjugate gradient method and extreme learning machine (ELM) algorithm was proposed. This study utilized finite element simulations and experiments, in order to analyze the feasibility of the proposed method. Then, following the corrections, the results indicated that the maximum relative error percentages of the displacements at measuring points in different bending shapes were determined to be 6.39% (Type 1), 7.04% (Type 2), and 7.02% (Type 3), respectively. Therefore, it was confirmed that the proposed correction method was feasible, and could effectively improve the abilities of sensors for displacement intellisense. In this paper, the designed intelligent sensor was characterized by temperature self-compensation, bending shape self-classification, and displacement error self-correction, which could be used for real-time monitoring of deformation field in rock, subgrade, bridge, and other geotechnical engineering, presenting the vital significance and application promotion value.


Introduction
Fiber Bragg grating (FBG) is a promising sensing element for the fabrication of flexible sensors owing to its advantageous characteristics, such as small size, light weight, water-and moisture-proofing abilities, and multiplex networking capacities. By embedding FBG into flexible substrates with properties of high flexibility, high ductility, and free bending, flexible sensors were developed and extensively utilized in civil engineering, medical engineering, aerospace, and robotics etc. for the measurement of morphology, displacement, strain, and other measurands [1][2][3][4]. Recently, morphological sensing with FBG embedded flexible sensors has become a hot topic [5][6][7][8]. Payo et al. [9] measured the deflections of a robotic arm based on FBG sensors and interpolation methods. Xu et al. and Bhamber et al. [10,11] adopted orthogonal FBG arrays to measure bidirectional curvatures, then reconstructed shapes using a curve fitting algorithm. Yi et al. [12] studied a shape reconstruction method based on spatial movable coordinates, and realized the shape reconstruction of an aircraft frames. Kang et al. [13] modeled the relationship of wavelength-strain-displacement and reconstructed the entire structural deformations and local displacements of aluminum and acrylic beam specimens. Also, Todd et al. [14] employed a material-adapted reference frame together with a local linearization approach, for the purpose of reconstructing the shapes of slender. In the field of geotechnical engineering, researchers have studied the morphological sensing methods for measuring displacement profiles in critical areas of side slopes, subgrades, bridges, and other objects based on the FBG embedded flexible sensors. Li et al. [15] established a deflection-strain relationship for FBG flexible rods, which was then used in physical model tests to measure the displacements of the underground caves. Guo et al. [16] calculated the displacement profiles of beams using the strains and rotation angles measured by the FBG sensing points. Kim N.S. et al. [17] measured the deflection curves using a regression analysis method, and then applied the method for the measurement of bridge displacements. Xu et al. [18] estimated the bending deformations of the beams based on the curvature functions of cross-sections.
In practical applications, several factors may affect the measurement accuracy of morphology, including the layout intervals of FBG, the strain transferring rate from the substrate to FBG, the cumulative errors induced by algorithm etc. Therefore, it is critical to correct the reconstructed shapes using effective methods, and to improve the shape measurement accuracy. In terms of morphological sensing error analyses and corrections, Zhang et al. [19] proposed a dynamic error analysis method based on a least mean square algorithm. This method has been used for the morphological reconstructions and error analyses of space plate structures. Sun et al. [20] experimentally calibrated the relationship between the wavelength and bending curvatures, and improved the reconstruction precision of polyimide film deformations through the linear interpolation of the curvatures. Wang et al. [21,22] proposed an in situ calibrated deformation reconstruction method for FBG embedded geogrid, and improved the accuracy effectively. In the aforementioned research, unified correction coefficients were mainly applied to calibrate and correct the entire morphology of the structure to be measured, as well as to improve the precision of the morphological sensing. However, it has been found that, in practice, due to the diversity of the measured structural morphologies, the unified coefficients have been difficult to be applicable to a variety of morphologies, leading to major differences of the correction results for different shapes. While the shape correction methods based on the different morphological classifications for the FBG embedded flexible sensors have not yet been reported.
In this paper, a morphological sensing method, which is capable of correcting the reconstructed shape based on the morphological classification, was proposed to further improve the measurement accuracy of FBG implantable flexible morphological sensors. First, an FBG implantable flexible morphological sensor with temperature self-compensation capabilities was developed, and an arc curve fitting based on morphological sensing algorithm was introduced. Then, temperature calibration experiment indicated that the temperature self-compensation could be realized by the differences in the central wavelength variations of the two sensing points in the detecting unit. In the morphological calibration experiment, a conjugate gradient method was adopted to define the correction coefficients k of each morphology. A morphological sensing method which is capable of correcting the measured shapes based on different morphological classifications was proposed. The classification of different bending shapes was performed by using the extreme learning machine (ELM) algorithm. Finally, numerical simulations and experimental analyses were performed to verify the feasibility and effectiveness of the proposed method. The results indicate that the proposed method can improve the morphological measurement accuracy of FBG implantable flexible morphological sensors effectively, and it is of vital significance for the measurement of displacement fields in subgrades, bridges, and other geotechnical engineering applications.

Fabrication of the FBG Implantable Flexible Morphological Sensor
As detailed in Figure 1, for the proposed FBG implantable flexible morphological sensor, an acrylonitrile butadiene styrene (ABS) rod with high strength, corrosion resistance, and high temperature resistance as the flexible substrate. The flexible rod length was 900 mm, and the diameter was 5 mm. Two rectangular grooves (both depth and width: 1 mm) were slotted along the axis of the rod surface, and the interval of the two grooves was 180 • . Plastic welding technology was used to implant two FBG sensing arrays into the two grooves for fixation purposes. FBG sensing arrays were fabricated using the phase mask method by Suzhou NanZee Sensing Technology Co., Ltd., Suzhou, China, of which each array had 9 FBG sensing points, and the spacing of two adjacent sensing points was 100 mm. The detailed parameters were shown in Table 1. An FBG demodulation instrument SM130 (MOI) was used to collect the central wavelength of the FBG sensing arrays, whose wavelength measurement ranged from 1510 nm to 1590 nm, and the demodulation precision was 1 pm. The full spectrum of the FBG array on one side in the free state was shown in Figure 2.

Fabrication of the FBG Implantable Flexible Morphological Sensor
As detailed in Figure 1, for the proposed FBG implantable flexible morphological sensor, an acrylonitrile butadiene styrene (ABS) rod with high strength, corrosion resistance, and high temperature resistance as the flexible substrate. The flexible rod length was 900 mm, and the diameter was 5 mm. Two rectangular grooves (both depth and width: 1 mm) were slotted along the axis of the rod surface, and the interval of the two grooves was 180°. Plastic welding technology was used to implant two FBG sensing arrays into the two grooves for fixation purposes. FBG sensing arrays were fabricated using the phase mask method by Suzhou NanZee Sensing Technology Co., Ltd., Suzhou, China, of which each array had 9 FBG sensing points, and the spacing of two adjacent sensing points was 100 mm. The detailed parameters were shown in Table 1. An FBG demodulation instrument SM130 (MOI) was used to collect the central wavelength of the FBG sensing arrays, whose wavelength measurement ranged from 1510 nm to 1590 nm, and the demodulation precision was 1 pm. The full spectrum of the FBG array on one side in the free state was shown in Figure 2.

Morphological Sensing Model Based on Arc Curve Fitting Method
The FBG central wavelength is sensitive to axial strain and temperature, simultaneously [23,24]. The relationships between the central wavelength variations Δλ and strain Δε, as well as the temperature ΔT are shown as follows: where λ is the initial wavelength of FBG at free state, and αf and ξ is the thermal expansion coefficient and the thermal-optic coefficient of optical fiber, respectively. Pe is the effective photo-elastic coefficient and is equal to 0.22. The FBG sensing points divided the flexible rod into several microsegments, and each segment was taken as a detecting unit. The measuring points were located at the end of the detecting units, as detailed in Figure 1b. Then, in accordance with the principle of material mechanics [25,26], when the sensor was bent, the curvature Rx of the detecting unit was where Z denoted the distance between the sensing points and the neutral axis, and Δεx represented the strain at the sensing points. In regard to the FBG sensing points arranged on the A and B symmetric sides of the detecting units, when the sensor rotated clockwise, sensing points 1 and 2 underwent tensile and compressive strain, respectively, as shown in Figure 3. Due to the fact that the two sensing points were equal distances from the neutral axis, the bending-induced wavelength variations were also equal in the opposite direction. Meanwhile, the temperature-induced wavelength variations displayed the same direction. When the detecting unit was bent, the strain of the detecting unit could be expressed as follows: x e λ λ ε ε ε P λ λ 12 12 12 As can be seen in Equations (2) and (3), the relationship between the curvature radius ρ of the detecting unit, and the central wavelength variation of the FBG were as follows:

Morphological Sensing Model Based on Arc Curve Fitting Method
The FBG central wavelength is sensitive to axial strain and temperature, simultaneously [23,24]. The relationships between the central wavelength variations ∆λ and strain ∆ε, as well as the temperature ∆T are shown as follows: where λ is the initial wavelength of FBG at free state, and α f and ξ is the thermal expansion coefficient and the thermal-optic coefficient of optical fiber, respectively. P e is the effective photo-elastic coefficient and is equal to 0.22. The FBG sensing points divided the flexible rod into several microsegments, and each segment was taken as a detecting unit. The measuring points were located at the end of the detecting units, as detailed in Figure 1b. Then, in accordance with the principle of material mechanics [25,26], when the sensor was bent, the curvature R x of the detecting unit was where Z denoted the distance between the sensing points and the neutral axis, and ∆ε x represented the strain at the sensing points. In regard to the FBG sensing points arranged on the A and B symmetric sides of the detecting units, when the sensor rotated clockwise, sensing points 1 and 2 underwent tensile and compressive strain, respectively, as shown in Figure 3. Due to the fact that the two sensing points were equal distances from the neutral axis, the bending-induced wavelength variations were also equal in the opposite direction. Meanwhile, the temperature-induced wavelength variations displayed the same direction. When the detecting unit was bent, the strain of the detecting unit could be expressed as follows: Sensors 2018, 18, 2342 5 of 20 As can be seen in Equations (2) and (3), the relationship between the curvature radius ρ of the detecting unit, and the central wavelength variation of the FBG were as follows: Therefore, the curvature radius of the detecting unit could be obtained by acquiring the central wavelength variation of the FBG.
Therefore, the curvature radius of the detecting unit could be obtained by acquiring the central wavelength variation of the FBG. When the morphological sensor was bent, each detecting unit was continuous, in which the curvature radius ρ was greater than the arc length L. According to the differential principle, when the detecting unit length is small enough, the FBG implantable flexible morphological sensor can be regarded as being composed of many detecting units. Therefore, the curvature of the detecting unit at the FBG sensing point can be used to replace the curvature of the entire arc curve. For the convenience of calculation, it was stipulated in this study that when the curve was bent counterclockwise, then the curvature of the detecting unit was positive. Furthermore, when bending clockwise, the curvature was negative. Then, by assuming the starting point of the FBG implantable flexible morphological sensor was fixed with the coordinate of O (0, 0), each detecting unit could be considered as an arc microsegment. Following this, Oi, Pi, ρi, and θi represented the center, end point, bending radius, and center angle of the i-th arc segment, respectively. The angle βi of Pi on the i-th arc curve was the angle between the PiOi and the horizontal line. When this was positive, the sensor was bent counterclockwise, while it was bent clockwise when it was negative, as shown in Figure 4.  When the morphological sensor was bent, each detecting unit was continuous, in which the curvature radius ρ was greater than the arc length L. According to the differential principle, when the detecting unit length is small enough, the FBG implantable flexible morphological sensor can be regarded as being composed of many detecting units. Therefore, the curvature of the detecting unit at the FBG sensing point can be used to replace the curvature of the entire arc curve. For the convenience of calculation, it was stipulated in this study that when the curve was bent counterclockwise, then the curvature of the detecting unit was positive. Furthermore, when bending clockwise, the curvature was negative. Then, by assuming the starting point of the FBG implantable flexible morphological sensor was fixed with the coordinate of O (0, 0), each detecting unit could be considered as an arc microsegment. Following this, O i , P i , ρ i , and θ i represented the center, end point, bending radius, and center angle of the i-th arc segment, respectively. The angle β i of P i on the i-th arc curve was the angle between the P i O i and the horizontal line. When this was positive, the sensor was bent counterclockwise, while it was bent clockwise when it was negative, as shown in Figure 4.
Therefore, the curvature radius of the detecting unit could be obtained by acquiring the central wavelength variation of the FBG. When the morphological sensor was bent, each detecting unit was continuous, in which the curvature radius ρ was greater than the arc length L. According to the differential principle, when the detecting unit length is small enough, the FBG implantable flexible morphological sensor can be regarded as being composed of many detecting units. Therefore, the curvature of the detecting unit at the FBG sensing point can be used to replace the curvature of the entire arc curve. For the convenience of calculation, it was stipulated in this study that when the curve was bent counterclockwise, then the curvature of the detecting unit was positive. Furthermore, when bending clockwise, the curvature was negative. Then, by assuming the starting point of the FBG implantable flexible morphological sensor was fixed with the coordinate of O (0, 0), each detecting unit could be considered as an arc microsegment. Following this, Oi, Pi, ρi, and θi represented the center, end point, bending radius, and center angle of the i-th arc segment, respectively. The angle βi of Pi on the i-th arc curve was the angle between the PiOi and the horizontal line. When this was positive, the sensor was bent counterclockwise, while it was bent clockwise when it was negative, as shown in Figure 4.  In regard to the first arc curve OP 1 , the coordinate of the center O 1 (O 1x , O 1y ) was (0, ρ 1 ), and the arc OP 1 could be represented in rectangular coordinates as follows: When t = θ 1 − π/2, the coordinate of point P 1 (P 1x , P 1y ) could be expressed as follows: When i ≥ 2, the arc curve P i−1 P i could be represented as follows: Due to the tangency of the arc curves, the center O i of arc curve P i−1 P i , center O i+1 of P i P i+1 , and tangent point P i existed along the same straight line. When transiting from the curve P i−1 P i to P i P i+1 , the curvature changed from positive to negative. At this moment, the bending radius ρ i+1 and center angle θ i+1 of the curve P i P i+1 were negative. Then, the equation of the curve P i P i+1 could be expressed as follows: The angle β i of point P i on the arc curve P i−1 P i was i ∑ n = 1 θ n − π/2, and the angle β i+1 on curve Then, in accordance with Equations (6)-(8), the measuring point coordinate P i (P ix , P iy ) at the end of each arc curve was as follows:

Temperature Calibration Experiment
In order to avoid the influences of temperature on the morphological sensing precision, a temperature calibration experiment was first carried out to verify the temperature characteristics and self-compensation effects of the FBG implantable flexible morphology sensor. The FBG implantable flexible morphology sensor was placed in the calorstat, with the temperature gradually adjusted from −20 • C to 40 • C. The temperature was raised by 10 • C each time, and the calorstat temperature measurement precision was 0.01 • C. SM130 was used to collect the central wavelengths of the two FBG arrays.
The temperature response coefficients of the two FBG arrays were shown in Table 2. It can be seen in Table 2 that the temperature response coefficients of the sensing points in the same detecting unit were basically the same. The temperature response curves of two sensing points in one of the detecting units are shown in Figure 5. Sensing point 1 and sensing point 2 were implanted on the sides of A and B, respectively, in which the relationships between the central wavelength variation and the temperature variation were λ 1 = 0.0110∆T + 0.2181 and λ 2 = 0.0109∆T + 0.2171, respectively. Therefore, the temperature self-compensation could be realized by the differences in the central wavelength variations of the two sensing points in the detecting unit. detecting units are shown in Figure 5. Sensing point 1 and sensing point 2 were implanted on the sides of A and B, respectively, in which the relationships between the central wavelength variation and the temperature variation were λ1 = 0.0110ΔT + 0.2181 and λ2 = 0.0109ΔT + 0.2171, respectively. Therefore, the temperature self-compensation could be realized by the differences in the central wavelength variations of the two sensing points in the detecting unit.

Morphological Calibration Experiment
A morphological bending experiment was carried out in this study in order to verify the feasibility of the morphological classification correction method, as shown in Figure 6. One end of the FBG implantable flexible morphological sensor was fixed on the calibration platform with coordinate paper (cell: mm 2 ). The coordinate of the fixed point on side B was O (0, 0), and the coordinates of each measuring point on side B were Pj (100*j, 0), with j = 1, 2, … 9 (see Figure 7). On the other end of the sensor, the tail fiber was connected to the SM130, and the wavelength data of the two FBG arrays obtained by the SM130 demodulation were transmitted to the computer via a network cable. Then, displacements were applied to the different locations of the sensor for the formation of the different bending shapes, and the displacements of each measuring point were read by a high-precision displacement sensor (precision: 0.01 mm).

Morphological Calibration Experiment
A morphological bending experiment was carried out in this study in order to verify the feasibility of the morphological classification correction method, as shown in Figure 6. One end of the FBG implantable flexible morphological sensor was fixed on the calibration platform with coordinate paper (cell: mm 2 ). The coordinate of the fixed point on side B was O (0, 0), and the coordinates of each measuring point on side B were P j (100*j, 0), with j = 1, 2, . . . 9 (see Figure 7). On the other end of the sensor, the tail fiber was connected to the SM130, and the wavelength data of the two FBG arrays obtained by the SM130 demodulation were transmitted to the computer via a network cable. Then, displacements were applied to the different locations of the sensor for the formation of the different bending shapes, and the displacements of each measuring point were read by a high-precision displacement sensor (precision: 0.01 mm).
The experiment was conducted on three typical bending shapes of the FBG implantable flexible morphological sensors, as shown in Figure 7a-c, and six different sizes of displacement were applied to each type. For Type 1 (Figure 7a), six different sizes of displacement were applied to the sensor at (900, 0) along the y direction, at 80 mm, 142 mm, 215 mm, 282 mm, 347 mm, and 400 mm, respectively. For Type 2 (Figure 7b), the displacements were applied along the y-axis at (500, 0) at 70 mm, 123 mm, 142 mm, 157 mm, 181 mm, and 211 mm, respectively. For Type 3 (Figure 7c), the displacements were applied to the sensor at (500, 0) and (900, 0) along the y direction, respectively, and the displacement values in the six cases were 50 mm/17 mm, 97 mm/49 mm, 96 mm/42 mm, 121 mm/63 mm, 151 mm/71 mm, and 190 mm/165 mm, respectively. In each case, ten repetitions were performed. Then, SM130 was used to collect the central wavelength of the FBG sensing points. When the displacement was determined to be stable, the real displacement values of the measuring points were read by a high-precision displacement sensor. By using the actual displacement values of the measuring points as the benchmark, the optimization goal was to minimize the mean relative errors (MREs) between the reconstructed displacement and real displacement at each measuring point.
The calculation formula of the MRE was 1 n n ∑ i = 1 |y Ci (k) − y Ai |, where n = 9, y Ci (k) was the calculation displacement of each measuring point based on corrected strain value k *∆ε x of each detecting unit and arc curve fitting method, and y Ai was the real displacement of each measuring point. Then, a conjugate gradient method was adopted to define the correction coefficients k of each morphology. Conjugate gradient method has the characteristics of simple iterative format, small amounts of storage required, high stability and fast convergence speed, which can be used to solve optimization problems. A detailed formula description of conjugate gradient method was shown as follows. The experiment was conducted on three typical bending shapes of the FBG implantable flexible morphological sensors, as shown in Figure 7a-c, and six different sizes of displacement were applied to each type. For Type 1 (Figure 7a), six different sizes of displacement were applied to the sensor at (900, 0) along the y direction, at 80 mm, 142 mm, 215 mm, 282 mm, 347 mm, and 400 mm, respectively. For Type 2 (Figure 7b), the displacements were applied along the y-axis at (500, 0) at 70 mm, 123 mm, 142 mm, 157 mm, 181 mm, and 211 mm, respectively. For Type 3 (Figure 7c), the displacements were applied to the sensor at (500, 0) and (900, 0) along the y direction, respectively, and the displacement values in the six cases were 50 mm/17 mm, 97 mm/49 mm, 96 mm/42 mm, 121 mm/63 mm, 151 mm/71 mm, and 190 mm/165 mm, respectively. In each case, ten repetitions were performed. Then, SM130 was used to collect the central wavelength of the FBG sensing points. When the displacement was determined to be stable, the real displacement values of the measuring points were read by a highprecision displacement sensor. By using the actual displacement values of the measuring points as The experiment was conducted on three typical bending shapes of the FBG implantable flexible morphological sensors, as shown in Figure 7a-c, and six different sizes of displacement were applied to each type. For Type 1 (Figure 7a), six different sizes of displacement were applied to the sensor at (900, 0) along the y direction, at 80 mm, 142 mm, 215 mm, 282 mm, 347 mm, and 400 mm, respectively. For Type 2 (Figure 7b), the displacements were applied along the y-axis at (500, 0) at 70 mm, 123 mm,  The optimization problem was The sequence {k m } could be obtained by the following iterative format, where d m could be expressed as follows: where α m was search step size, which obtained by satisfying the following weak Wolfe conditions, where 0 < δ ≤ σ < 1. β m was the scalar of the representation algorithm and PRP method was selected in this paper, which could be expressed as follows: The correction coefficient k of each morphology was obtained by the conjugate gradient method (see Table 3). Correction coefficients for Types 1 and 2 were determined to be k 1 = 1.45 and k 2 = 1.

Correction Methods
The morphological sensing algorithm which was based on arc curve fitting has the ability to evaluate the shape of the FBG implantable flexible morphological sensor. Thus, it is able to evaluate the overall displacement profiles of the structures to be measured. However, the following factors may cause large morphological sensing errors and affect the precision of the displacement measurements. First, during the sensor fabrication processes, curvature measurement errors at the FBG sensing points may easily be caused due to the deviations when implanting the FBG into the substrate. Second, the curvatures measured by the flexible sensor are discontinuous, which is caused by a certain interval between the adjacent FBG sensing points. The curvatures of the FBG sensing points are used to replace the curvatures of the detecting units in the calculation process of the arc curve fitting algorithm, which may induce the measurement errors. Third, the morphological sensing method is the integration of the piece-wise detecting units, and the displacement errors of each segment will gradually accumulate to generate accumulated errors at the end of the sensor.
To further improve the measurement accuracy of the FBG implantable flexible morphological sensor, the different bending shapes of the sensor were automatically categorized using the intelligent classification algorithm. By calibrating the typical bending shapes of the sensor, the correction coefficients of different shapes were determined, and the corresponding correction coefficients of different bending shapes was then selected automatically to correct the reconstructed shapes for the purpose of minimizing the errors.
ELM algorithm was adopted to intelligently classify the bending shapes of the FBG implantable flexible morphological sensor in this research study. ELM is a novel intelligence optimization algorithm for single-hidden layer feedforward neural networks (SLFNs), which was proposed in 2004 [27]. This algorithm randomly selects hidden nodes and analytically determines the output weights of SLFNs. Compared with conventional popular learning algorithms, this algorithm has the characteristics of extremely fast learning speed, good generalization performance, and tuning of free parameters. A detailed formula description of ELM was presented in Ref. [28].
Define that there are N arbitrary distinct training samples (x i , y i ) ∈ R n × R m , where x are the training inputs and y are the training targets. The function of ELM classifier with L hidden layer neurons and a sigmoid function g(x) is mathematically modeled as follows: where β j is output weight matrix, a j is input weight matrix, b j is bias connecting the input neurons. According to any continuous probability distribution, a j ,b j can be randomly generated. Therefore, the Equation (18) can be written in the form of the following matrix. where where H is the hidden layer output matrix of ELM. These training algorithms need to adjust the input weights and the hidden layer biases, and the output weights can be given as follows: where H † is the Moore-Penrose generalized inverse of matrix H.
When H T H is nonsingular, or when HH T is nonsingular, Therefore, the output function of ELM classifier can be shown as follows: The steps of classification of the morphological correction method were presented as follows: Step 1: The central wavelengths of the two FBG arrays were obtained using an FBG demodulator, and Equation (4) was used to convert the wavelength variations at the sensing point to the curvature of detecting unit. Then, based on an arc curve fitting method, the bending shapes of the FBG implantable flexible morphological sensor were first reconstructed.
Step 2: The sensor was calibrated in order to obtain the actual displacement of each measuring point in the typical bending shape. Then, a conjugate gradient method was applied to determine the correction coefficients k of the different bending shapes, with the actual displacement as the standard.
Step 3: The experimental process was completed on different bending shapes, and then Step 1 was repeated. The displacement of the measuring point P i along the y-axis was taken as the input. Meanwhile, different bending shape types were regarded as the output in order to train the ELM classifier model.
Step 4: An ELM classifier was used to classify the bending shapes first reconstructed in Step 1. Then, the corrected strain value k *∆ε x of each detecting unit determined in Step 2 was applied to reconstruct the bending shapes of the FBG implantable flexible morphological sensor in order to obtain an accurate bending shape. The specific process was shown in Figure 8.  Figure 8. Flow chart of classification morphological correction method based on ELM intelligent algorithm.

Simulation Verification
In this study, a finite element method was applied to simulate the bending shapes of the FBG implantable flexible morphological sensor, and to compare the morphological sensing effects which had been obtained by the unified coefficient and classification coefficient corrections. First, a cylinder model with a length of 900 mm and a diameter of 5 mm was established, as shown in Figure  9a1,b1,c1,d1,e1. The initial end center coordinates of the cylinder model was O (0, 0, 0), and the cylinder model was divided into nine segments along the x-axis. Each segment was taken as a detecting unit of the proposed FBG implantable flexible morphological sensor. The column material was ABS, and its mechanical property parameters were shown in Table 4. In this study's simulation, the cylinder model had a total of 4379 meshes with minimum unit masses of 0.351. Then, with consideration given to the geometrical nonlinearity of the material used in the simulation, an elastic model was selected.  The above method was adopted for bending shape classification correction of the proposed FBG implantable flexible morphological sensor, which had the ability to effectively reduce the displacement errors of the measuring points caused by the sensor preparation errors, FBG network capacity limitations, and other factors. Therefore, the classification morphological correction method proposed in this study displayed the ability to improve the intellisense displacement abilities of the FBG flexible morphological sensor.

Simulation Verification
In this study, a finite element method was applied to simulate the bending shapes of the FBG implantable flexible morphological sensor, and to compare the morphological sensing effects which had been obtained by the unified coefficient and classification coefficient corrections. First, a cylinder model with a length of 900 mm and a diameter of 5 mm was established, as shown in Figure 9a 1 ,b 1 ,c 1 ,d 1 ,e 1 . The initial end center coordinates of the cylinder model was O (0, 0, 0), and the cylinder model was divided into nine segments along the x-axis. Each segment was taken as a detecting unit of the proposed FBG implantable flexible morphological sensor. The column material was ABS, and its mechanical property parameters were shown in Table 4. In this study's simulation, the cylinder model had a total of 4379 meshes with minimum unit masses of 0.351. Then, with consideration given to the geometrical nonlinearity of the material used in the simulation, an elastic model was selected. A displacement constraint was fixed on one side of the cylinder model, and five different types of displacement were applied in order to simulate five different typical bending shapes. For Type 1 (Figure 9a 1 ), a 140 mm displacement was applied along the y-axis at (900, −2.5, 0) in the cylinder model. For Type 2 (Figure 9b 1 ), a 50 mm displacement was applied along the y-axis at (500, −2.5, 0). For Type 3 (Figure 9c 1 ), 92 mm and 30 mm displacements were applied along the y-axis at (500, −2.5, 0) and (900, −2.5, 0), respectively. For Type 4 (Figure 9d 1 ), 30 mm, 9.5 mm, and 40 mm displacements were applied along the y-axis at (300, −2.5, 0), (600, −2.5, 0), and (900, −2.5, 0), respectively. For Type 5 (Figure 9e 1 ), 30.7 mm, 11 mm, 38 mm, and 18.2 mm displacements were applied along the y-axis at (300, −2.5, 0), (500, −2.5, 0), (700, −2.5, 0), and (900, −2.5, 0), respectively. For the five typical bending shapes, the strain value at i = 0, 1, ... 8 in the coordinate points of S i (50 + 100*i, −2.5, 0) were extracted from the cylinder model. Meanwhile, the displacement at j = 1, 2, ... 9 in the coordinate points of P j (100*j, −2.5, 0) were extracted as the standard displacement values. The strain values which had been extracted were used for the first reconstruction of the different bending shapes of the cylinder model by arc curve fitting method, in order to obtain the displacement of the measuring points under different bending shapes, as shown in Figure 9a 2 ,b 2 ,c 2 ,d 2 ,e 2 . It was observed that, having been influenced by multiple factors (such as the accumulated errors of detecting units), there were certain errors between the first reconstructed displacement and the simulated displacement. For Types 1 and 2, the relative errors of the measuring points away from the fixed end had gradually increased (with maximum errors of 12.36 mm and 16.66 mm, respectively), and both were located at the ninth measuring point. This was determined to be due to the fact that the cylinder model had undergone uniaxial stress conditions for Types 1 and 2, and each measuring point displayed the same displacement error direction. The displacement errors of the measuring points tended to accumulate point by point. However, for Types 3, 4, and 5, the measuring points with the maximum relative errors were located in the fifth, third, and seventh points, with relative error values of 13.81 mm, 4.93 mm, and 3.79 mm, respectively. These results may be due to the fact that the cylinder model had undergone multidirectional stress conditions for Types 3, 4, and 5, and the displacement errors of the measuring points displayed the phenomena of positive and negative error offsets. in the coordinate points of Pj (100 * j, −2.5, 0) were extracted as the standard displacement values. The strain values which had been extracted were used for the first reconstruction of the different bending shapes of the cylinder model by arc curve fitting method, in order to obtain the displacement of the measuring points under different bending shapes, as shown in Figure 9a2,b2,c2,d2,e2. It was observed that, having been influenced by multiple factors (such as the accumulated errors of detecting units), there were certain errors between the first reconstructed displacement and the simulated displacement. For Types 1 and 2, the relative errors of the measuring points away from the fixed end had gradually increased (with maximum errors of 12.36 mm and 16.66 mm, respectively), and both were located at the ninth measuring point. This was determined to be due to the fact that the cylinder model had undergone uniaxial stress conditions for Types 1 and 2, and each measuring point displayed the same displacement error direction. The displacement errors of the measuring points tended to accumulate point by point. However, for Types 3, 4, and 5, the measuring points with the maximum relative errors were located in the fifth, third, and seventh points, with relative error values of 13.81 mm, 4.93 mm, and 3.79 mm, respectively. These results may be due to the fact that the cylinder model had undergone multidirectional stress conditions for Types 3, 4, and 5, and the displacement errors of the measuring points displayed the phenomena of positive and negative error offsets.
During this study's experimental process, based on the actual applied displacement as the benchmark, weighted corrections were conducted on the first reconstructed morphology. Then, a During this study's experimental process, based on the actual applied displacement as the benchmark, weighted corrections were conducted on the first reconstructed morphology. Then, a conjugate gradient method was adopted to define the correction coefficient of each morphology. As detailed in Table 5, correction coefficients for Types 1 and 2 were determined to be k 1 = 1. respectively. The unified correction coefficient was the average of k 1 , k 2 , k 3 + , k 3 − , k 4 + , k 4 − , k 5 + , and k 5 − , and denoted as k= 1.18. Figure 10 shows the MREs of the sensing displacement of the cylindrical model under the five bending shapes based on different correction methods. It was observed that under the different bending shapes, the MREs of the measuring points were significantly different when the strain values in the center of the detecting units were corrected using the unified coefficients. For Types 1, 2, 3, and 4, the MREs of the corrected measuring points were determined to decline, of which that of Type 3 displayed the minimum decline of 0.84 mm, while that of Type 1 displayed the maximum decline of 4.27 mm. For Type 5, the MREs had increased following the corrections of the unified coefficients. However, for the different bending shapes, the errors in the measuring points which had been corrected by the different coefficients had obviously decreased. The MREs had been reduced by 4.33 mm (Type 1), 7.78 mm (Type 2), 6.47 mm (Type 3), 1.18 mm (Type 4), and 1.03 mm (Type 5), respectively, after the classification corrections. These findings indicated that, when compared with the unified coefficient corrections, the classification corrections of the different bending shapes had improved the measurement precision of the displacements, which confirmed that it was necessary to use different coefficient corrections for the various bending shapes. corrected by the different coefficients had obviously decreased. The MREs had been reduced by 4.33 mm (Type 1), 7.78 mm (Type 2), 6.47 mm (Type 3), 1.18 mm (Type 4), and 1.03 mm (Type 5), respectively, after the classification corrections. These findings indicated that, when compared with the unified coefficient corrections, the classification corrections of the different bending shapes had improved the measurement precision of the displacements, which confirmed that it was necessary to use different coefficient corrections for the various bending shapes.

Morphological Classification Based on ELM Algorithm
In Section 3.2, morphological calibration experiment was conducted on three typical bending shapes of the FBG implantable flexible morphological sensors. Six different sizes of displacement were applied to each type, and ten repetitions were performed in each case. An arc curve fitting algorithm was used to calculate the initial morphology of the FBG implantable flexible morphological

Morphological Classification Based on ELM Algorithm
In Section 3.2, morphological calibration experiment was conducted on three typical bending shapes of the FBG implantable flexible morphological sensors. Six different sizes of displacement were applied to each type, and ten repetitions were performed in each case. An arc curve fitting algorithm was used to calculate the initial morphology of the FBG implantable flexible morphological sensor. On this basis, an ELM algorithm was adopted to train the displacement data in the initial morphology along the y direction. Also, the different bending shapes were intelligently classified. To be specific, for the six circumstances of each bending shape, the repetitive experimental data of the previous seven times of each circumstance were taken as the training samples (training samples totaled 126 groups). The displacement data of the training sample along the y-axis were used as the input, while the bending shape Types (Types 1, 2, and 3) were adopted as the output, in order to establish an ELM model. The repetitive experimental data of the later three times were taken as the test samples (test samples totaled 54 groups) in order to verify the classification effects, and the results are shown in Figure 11. The classification intervals of Types 1, 2, and 3 were [0.75, 1.25], [1.75, 2.25], and [2.75, 3.25], respectively. It was observed that for the 54 groups of test sample data, the ELM algorithm was able to effectively realize the intelligent classification of the three typical bending shapes of the FBG implantable flexible morphological sensor. The ELM algorithm was implemented based on MATLA programming software, with an entire program running time of 477 ms (computer configuration: Inter(R) Core(TM) i5-4210M CPU and 4G RAM). The results indicated that the algorithm had a fast calculation speed with good generalization ability, and was able to effectively ensure the real-time performance of the morphological sensing.
algorithm was able to effectively realize the intelligent classification of the three typical bending shapes of the FBG implantable flexible morphological sensor. The ELM algorithm was implemented based on MATLA programming software, with an entire program running time of 477 ms (computer configuration: Inter(R) Core(TM) i5-4210M CPU and 4G RAM). The results indicated that the algorithm had a fast calculation speed with good generalization ability, and was able to effectively ensure the real-time performance of the morphological sensing.

Comparison of Morphological Sensing Based on Different Correction Methods
In this study, for each type of bending shape, the repetitive experimental data of the previous three times in one experimental circumstance were selected to compare the effects from the unified coefficient correction and the classification coefficient correction on the displacement sensing of the FBG implantable flexible morphological sensor. For each group of experimental data, the classification morphological correction method and an arc curve fitting algorithm were used to calculate the displacement of each measuring point, respectively. The standard deviations of the displacements at each measuring point for each bending shape are shown in Table 6. For each bending shape, the average of repetitive experimental data of the previous three times was taken. Figure 12 shown the contrast between the displacements of the measuring points before Figure 11. Classification results of bending shapes based on testing samples.

Comparison of Morphological Sensing Based on Different Correction Methods
In this study, for each type of bending shape, the repetitive experimental data of the previous three times in one experimental circumstance were selected to compare the effects from the unified coefficient correction and the classification coefficient correction on the displacement sensing of the FBG implantable flexible morphological sensor. For each group of experimental data, the classification morphological correction method and an arc curve fitting algorithm were used to calculate the displacement of each measuring point, respectively. The standard deviations of the displacements at each measuring point for each bending shape are shown in Table 6. For each bending shape, the average of repetitive experimental data of the previous three times was taken. Figure 12 shown the contrast between the displacements of the measuring points before and after the corrections and the real displacements Figure 12a 1 ,b 1 ,c 1 under the different bending shapes, as well as the relative errors of the measuring points before and after the corrections Figure 12a 2 ,b 2 ,c 2 . It can be seen in Figure 12a 2 ,b 2 ,c 2 that there were certain errors between the first reconstructed displacements and the real displacements. It was observed that among these, with the increasing distances between the measuring points and fixed point, the relative errors of Types 1, 2, and 3 had gradually increased. The maximum relative errors had all occurred to the ninth measuring point (for example: 38.57 mm, 67.15 mm, and 53.06 mm). It can be seen that the first reconstructed displacement of the measuring point had a larger error. Following the corrections, the errors of the three bending shapes had been significantly reduced. The MREs obtained by the first reconstruction, along with the different correction methods, were shown in Figure 13. After the central strain values of detecting units were corrected by the unified coefficients, the MREs of the measuring points displayed larger differences. For example, for Type 1, the MRE of the measuring points after correction was the minimum at 3.84 mm; for Type 2, the MRE was 15.79 mm; for Type 3, the MRE was the maximum at 17.95 mm. However, for the different bending shapes with different coefficient corrections, it could be seen that the MREs of measuring points had obviously decreased. The MREs of Types 1, 2, and 3 after the classification corrections were 4.8 mm, 1.67 mm, and 2.4 mm, respectively. The MREs had been reduced by 13.16 mm (Type 1), 29.90 mm (Type 2), and 30.56 mm (Type 3), respectively, following the classification corrections. When compared with the unified coefficient corrections, the classification coefficient corrections displayed a good correction ability for the displacement errors of the measuring points with different bending shapes. point (for example: 38.57 mm, 67.15 mm, and 53.06 mm). It can be seen that the first reconstructed displacement of the measuring point had a larger error. Following the corrections, the errors of the three bending shapes had been significantly reduced. The MREs obtained by the first reconstruction, along with the different correction methods, were shown in Figure 13. After the central strain values of detecting units were corrected by the unified coefficients, the MREs of the measuring points displayed larger differences. For example, for Type 1, the MRE of the measuring points after correction was the minimum at 3.84 mm; for Type 2, the MRE was 15.79 mm; for Type 3, the MRE was the maximum at 17.95 mm. However, for the different bending shapes with different coefficient corrections, it could be seen that the MREs of measuring points had obviously decreased. The MREs of Types 1, 2, and 3 after the classification corrections were 4.8 mm, 1.67 mm, and 2.4 mm, respectively. The MREs had been reduced by 13.16 mm (Type 1), 29.90 mm (Type 2), and 30.56 mm (Type 3), respectively, following the classification corrections. When compared with the unified coefficient corrections, the classification coefficient corrections displayed a good correction ability for the displacement errors of the measuring points with different bending shapes.     Table 7 shows the relative error percentages at each measuring point under the three bending shapes following the classification corrections. It could be seen that under the different bending shapes, the maximum relative error percentages of each measuring point of the FBG implantable flexible morphological sensor were 6.39% (Type 1), 7.04% (Type 2), and 7.02% (Type 3), respectively. Following the corrections of the different bending shapes, the relative error percentages of each measuring point were observed to be relatively low. These findings indicated that the proposed classification morphological correction method was feasible, and the method was found to effective in improving the precision of the FBG implantable flexible morphological sensor when sensing deformation fields in geotechnical engineering related applications.  Table 7 shows the relative error percentages at each measuring point under the three bending shapes following the classification corrections. It could be seen that under the different bending shapes, the maximum relative error percentages of each measuring point of the FBG implantable flexible morphological sensor were 6.39% (Type 1), 7.04% (Type 2), and 7.02% (Type 3), respectively. Following the corrections of the different bending shapes, the relative error percentages of each measuring point were observed to be relatively low. These findings indicated that the proposed classification morphological correction method was feasible, and the method was found to effective in improving the precision of the FBG implantable flexible morphological sensor when sensing deformation fields in geotechnical engineering related applications.

Conclusions
In this research, an FBG implantable flexible morphological sensor with the capacity of temperature self-compensation, bending shape self-classification, and displacement error self-correction was designed by internally implanting a 180 • spacing FBG array in an ABS rod. On this basis, the causes of the morphological sensing errors were analyzed for an FBG implantable flexible morphological sensor. Meanwhile, the key focus was to propose a morphological sensing method for an FBG implantable flexible sensor based on classification morphological corrections. A temperature calibration experiment was performed to verify that the influences of temperature could be eliminated by using the central wavelength variations of two FBG sensing points in each detecting unit. Although limited by such factors as the yielding strength of the flexible rod substrate and various experimental conditions, this study's experiments were carried out on three typical bending shapes of the FBG implantable flexible morphological sensor, and the correction coefficients k of each morphology were obtained by the conjugate gradient method. A finite element simulation was adopted for five typical bending shapes of a cylinder model, and the necessity and feasibility of the proposed method were analyzed. The intelligent classification of the different bending shapes were realized using an ELM classifier. During this study's experimental process, the MREs of the measuring points were observed to be significantly reduced following the coefficient corrections of the different bending shapes. The MREs were found to be reduced by 13.16 mm (Type 1), 29.90 mm (Type 2), and 30.56 mm (Type 3), respectively, following the classification corrections. When compared with the unified coefficient corrections, the classification coefficient corrections displayed a good correction ability for the displacement errors at the measuring points in the different bending shapes. Following the corrections by the classification morphological correction method, the maximum relative error percentages of the displacements at measuring points in different bending shapes were determined to be 6.39% (Type 1), 7.04% (Type 2), and 7.02% (Type 3), respectively. Following the corrections of the different bending shapes, the relative error percentages of each of the measuring points were relatively low. These findings confirmed that the proposed classification correction method was feasible and effective in improving the intelligent displacement sensing ability of the FBG implantable flexible morphological sensor. However, due to the limitations of the experimental conditions and the properties of the substrate material, the experimental validation was only carried out on three morphologies in this study. In the future, the correction performances of various bending shapes will be considered in detail, and the effects of the proposed classification morphological correction method under three-dimensional coordinates will be further examined. This study's designed FBG implantable flexible morphological sensor can be potentially used for the real-time monitoring of deformation fields in rock masses, subgrades, bridges, as well as other geotechnical engineering processes, and presents vital significance and application promotion value.