Special Issue

## New Trends in Numerical Simulation and Data Analysis

View this Special Issue

Research Article | Open Access

Volume 2020 |Article ID 6701696 | https://doi.org/10.1155/2020/6701696

Jinwei Qiao, Na Liu, "Numerical Simulation of Bending Stiffness Analysis for Spring Linkage Applied to In-Pipe Robot", Mathematical Problems in Engineering, vol. 2020, Article ID 6701696, 11 pages, 2020. https://doi.org/10.1155/2020/6701696

# Numerical Simulation of Bending Stiffness Analysis for Spring Linkage Applied to In-Pipe Robot

Accepted07 Aug 2020
Published27 Aug 2020

#### Abstract

Spring linkage can be applied to in-pipe robots for connecting different modules together and can make it pass through elbows more easily. However, its stiffness cannot be set to be too hard or too soft. This paper tries to make a balance between the compressive stiffness and the bending stiffness of the spring. After a brief introduction to the construction mechanism and some assumptions, the mathematical representation of the spring bending stiffness was deduced based on the Kirchhoff theory which describes the spatial curve with displacement rather than time. Then, some simulations aiming at verifying the correctness of the deduced bending stiffness expression were carried out. Finally, the relationship between the two rigidities was found out, which helps to find a way to decrease the bending stiffness of spring while keeping its compressive stiffness strong enough.

#### 1. Introduction

Many in-pipe robots have been widely applied to pipeline surveillance, detection, and maintenance . In most cases, the in-pipe robot with only one module is always not sufficient enough to carry heavy payload or to pass through bent parts of a pipeline [5,6]. The linkage for connecting modules is incredibly necessary . However, few concerns have been arisen. Reported linkages can be divided into two classes, the active and the passive [8, 9]. As the name suggests the active one can change direction actively . It has steering mechanism that usually consists of motor, transmission system, and the executor. For example, in , the Explorer Robot adopts motor actuated gear transmission box as its linkage to transmit not only force and motion but also torque, developed by the Robotics Institute at Carnegie Mellon University. Similar steering mechanisms appeared in . The steering mechanism adds the weight of the robot body and is always quite complicated, too large to be used in closed pipeline. In view of this point, the passive one is better . As the typical one, the universal joint has obtained the widest applications  such as the robots described in [19, 20]. It is found that the use of traditional joints often leads to some inconvenience, such as the enlarged volume, raised weight, and blocked transportation vessel among modules. In order to get through the blocked transformation, a specially designed universal joint was processed in . It is quite inconvenient.

One alternative is the spring that was firstly introduced in . With the addition of a ball in the updated version, it obtained the ability to pass though elbows. However, there would be pressure acting directly on the pipe wall, caused by bending connecting spring, as the in-pipe robot tries to pass through elbows. The pressure then turns to be resistant friction for moving robot. Therefore, the connecting spring cannot be designed to be too soft, which means large bending stiffness may push the two ends of spring against the pipeline hardly and create horrible friction [23, 24]. At the same time, the compressive stiffness of spring cannot be too small . Effective transmitting of force and motion relies on strong compressive stiffness; otherwise, the transmitted force will be decreased and the transmitted motion will be delayed. It needs a balance.

The numerical simulation is the main researching method for bending stiffness of spring linkage [26, 27]. There are many numerical methods in recent years, which mainly include finite element method (FEM) [28, 29], boundary element method (BEM) [30, 31], and meshless method . The FEM is the most widely used numerical simulation method, but in the simulation of nonlinear problems, the meshless method and the BEM have more advantages than the finite element method [38, 39]. In this paper, the FEM is used to simulate bending stiffness of spring linkage.

The rest of this paper is mainly organized as follows. Firstly, it introduces the basic concept of the linkage mechanism and explains its advantages. Then, the expression for spring stiffness is deduced, followed by some simulations and discussion. The conclusion obtained from this work is given at last.

#### 2. Mechanism Description

As shown in Figure 1, two separate modules of in-pipe robots with end screw threads are connected using a compressible round coil spring. The screw pitches l1 of the two modules share the same size, larger than that of the spring, termed as l2, while the bottom diameter of the screw thread and the inner diameter of the spring share the same size. Therefore, the spring pitch could be elongated when engaged with screw thread, so that the tensile force from lengthened spring will hold the screwed joint stably.

Compared with the traditionally rigid gimbals, the spring linker has quite small bending stiffness and of course still has the ability to transfer force and motion from one module of the robot to another one. It has couple of advantages illustrated as follows:(a)No Singular Points. The traditional rigid universal joint can be applied in transmitting force and torque from one direction to another direction in 3D space. Nevertheless, there are usually some dead points on the boundary of its reachable space that constrains the range of the transmission angle. This singularity can be avoided in flexible spring linkages transmitting.(b)Hollow Conformation. This nature helps the robots find their ways to pass through pneumatic lines and electric and signal cables.(c)Length Shortened. The length can be eliminated to a minimum size due to their high flexibility and compact assembly.(d)Crash Cushion. Pipes always have cracks, corruption pots, impurity, and something else in it. Therefore, there must be a lot of stochastic crashes that occurred to moving robot perform a great harm to the machine mechanism. The linker spring relaxes these strikes and makes the motion of robot more smooth.(e)Energy Storage. The linker spring works as an energy storage unit if the action of telescopic mechanism could be properly controlled. The energy storage function improved the energy utilization ratio and synthesis characteristic of the robot.

With nominal diameter of 2R and helix angle of θ, the spring has the same circular cross sections from head to foot. r is the radius of the spring wire whose original length is l. Before analysis, we assume the following. Firstly, the material is well distributed so that its elastic modulus and shear modulus remain constant everywhere . Secondly, the deformation along the axis of spring wire is ignored for its tiny influence. Thirdly, the gravity is neglected and no contact happens. Lastly, in the deforming process, the rigid cross section keeps perpendicular to the axis line of spring wire all the time. All the assumptions made herein are reasonable, further meeting the hypothesis of Kirchhoff theory [40, 41]. Additionally, note that the compressive stiffness of a straight spring is equal to its tensile stiffness, so only the compressive stiffness will be referred to in the following description.

#### 3. Bending Stiffness

##### 3.1. Coordinate System Built

In order to describe the deformation, some reference frames have been built on the vertical posed spring. As shown in Figure 2, the general Descartes reference system is located at the start point at the bottom. Its ζ-axis is set to be parallel with the central axis of the spring and the ξ-axis set horizontally left forward. An arc coordinate is built up from the start point so that every point along the axis of spring wire can be expressed as , where s is the accumulated arc length. Translating from start point to point yields local coordinate system . The rotated around ζ-axis by α is , which is then rotated around ξ′-axis by the helix angle θ forming the principle coordinate system of the circular cross section. As the cross section rotates around the tangential direction to produce a new system , the rotated angle is denoted as β. As a result, the three Euler angles α, θ, and β, which are all functions of the arc length s, together define the orientation of every cross section. Figure 3 shows the conversion process of different frames from beginning to end.

##### 3.2. Initial State Description

Based on the three Euler angles, the curvature vector for a spatial curve has the following expression :

For a spring shown in Figure 1 with circular cross section, the β is zero . Thus, the initial state of every piece of cross sections can be simplified toin which not only , but also . That is because the relative torsion angle and its derivative are both equal to zero as the spring stays free. Moreover, according to the parameters defined above, we have the relation shown as

##### 3.3. Applied Moment Description

Figure 4 shows the loading condition of the linkage spring as the in-pipe robot is moving inside elbows. Two ends of the connecting spring would be under the action of two equally opposite moments M, whose direction is along the ξ-axis of the general reference system .

In order to transfer the moment from the general coordinate system to the principle one, the relationship between these two is written down according to their definitions illustrated above:in which, and are the coordinate basis vectors of general and principle reference systems, respectively . Thus, the acting moment M in ξ-axis direction can be transmitted to

##### 3.4. Bent Spring Description

Under the action of moment M, the connecting spring gets bent around ξ-axis and two end faces turn to form an included angle γ. It is worth noticing that the helix angle θ cannot remain as a constant any more due to the flexural deformation of the spring. As a result, the curvature vector of the bent spring can be expressed under the principle coordinate system as

Substituting (4) into (6) yields

##### 3.5. Solution

According to the Kirchhoff equation, for each piece of cross sections, we havein which, E and G are the elastic modulus and shear modulus of the spring material, respectively. Ix and Iy are two inertia moments of the related cross section around x- and y-axis, while Iz is its polar moment of inertia around z-axis. For the spring with circular cross sections, 2Ix = 2Iy = Iz = πr4/2.

By substituting (2), (5), and (7) into (8) simultaneously, the following holds:

By combining the last two equations in (9), the following relation holds:

Then, we solve (10) by recurring to (3) and taking into consideration. The result isin which, the second and higher order small quality are both neglected.

On the other hand, the substitution of into the first equation of (9) results in

That is,

Taking (13) back into (11) leads to

Multiplied by ds on the two sides, (14) takes integration along the whole spring wire:in which, l is the whole length of the spring coil. It can be expressed as

Remember that the two ends of the spring are both fixed to the robot body, which means no variation or variation rate of the Euler angles has happened at those sections. Therefore, the integration of the second item laid on the left side of the equal sign in (15) equals to zero. Besides, the integration on the right hand is the sum of rotation of all the pieces of the cross sections, which is γ. Thus,

In this work, we define the bending stiffness Kb of connecting spring as the ratio of the applied moment M at its two ends to the resulting rotation angle γ. By taking into consideration, the expression of result is

By making use of (20) and (16), the bending stiffness of the connecting spring expressed in (18) can be rewritten down asin which the following relation was applied:

It can be seen from (18) that the bending stiffness Kb is an absolutely nonlinear function, in terms of material and construction parameters. Meanwhile, it indicates that the bending stiffness Kb of spring decreases with the raised initial helix angle or the elongated spring wire but increases as the radius of the spring wire grows thicker.

In particular, for , spring changes to straight rod with bending stiffness of

The result shown in (21) is similar to but more accurate than the traditional one expressed as .

It is worth remembering that the compressive stiffness of a spring has the definition 

#### 4. Simulation and Verification

It can be seen from (19) that the bending stiffness is a function of n, R, r, and θ. In order to verify its rightness and find out its scope of application, this section makes some simulations that were performed using the following steps:Step 1: select a set of parameters as an example to make model in SolidWorksStep 2: import the model from SolidWorks into Workbench to get meshed and loadedStep 3: make simulation on compressive stiffness and then compare the result with that calculated in (22) to ensure the validity of the simulation itselfStep 4: make simulation on bending stiffness and compare its result with that calculated in (19) to verify the correctness of (19)Step 5: repeat operation from step 1 to step 4 with varying construction parameters n, R, r, and θ, in order to find out its applicable scope

##### 4.1. Data Selection

The parameter values used for the example simulation are listed in Table 1.

 Item Value Elastic modulus (EPa) 2 × 1011 Poisson ratio ν 0.3 Shear modulus (GPa) 0.77 × 1011 Radius of spring wire r (mm) 0.5 Normal radius of spring R (mm) 5 Active coils number n 8 Pitch l2 (mm) 10 Spring length L (mm) 80 Length of spring wire l (mm) 263.753 Helix angle θ0 (deg) 17.66
##### 4.2. Simulation Performance

Based on the data listed in Table 1, a modified spring model was made as shown in Figure 5 which has two extra solid connectors for the sake of easy loading. Then, the model was imported from SolidWorks into Workbench for next simulation.

As verification of the simulation itself, the imported spring model was firstly simulated on compressive stiffness and then on bending stiffness. Figure 5 shows the result of spring under compressive force of 10 N that was acting along the spring axial. The result is 8.8581 mm. Thus, the compressive stiffness of the spring is 1.13 N/mm, while the result of (22) is 1.20 N/mm. The similarity of the two results successfully ensured the validity of simulation.

Once the validity of simulation got verified by the comparison of calculation and simulation results, we began to make bending simulation. One of the results is shown in Figure 6, with two moments of 10 Nmm acting in the way shown in Figure 4. In order to find out the rotated angle of the end face, we made use of the maximum and minimum displacement in the axial direction. They are 1.8787 mm and −1.879 mm, respectively, in Figure 6.

By considering the geometry relationship shown in Figure 7, the following equation holds:in which, h1 + h2 is the spring outside diameter, equal to 2(R + r). As a result, the rotated angle γ is 0.329, and then the bending rigidity is 30.379 Nmm/rad. On the other hand, the corresponding result can be obtained by substituting all needed values shown in Table 1 into (19). It is 31.2 Nmm/rad, very close to that of simulation result.

#### 5. Results and Discussion

After a demonstrative simulation, the parameters that appeared in (19) were all changed to reconstruct the spring and redo the simulation, one at a time. Firstly, the active coils number n was chosen to be changed, from 4 to 18. Figure 8 shows the results of both simulation result represented by star and the calculation result represented by line.

It can be seen from Figure 8 that both the compressive and bending rigidities have their calculated and simulated results staying close to each other after n grows larger than 4. The active coils number n provides an inverse proportion effect on two rigidities.

The nominal radius R was the second one to be changed from 6 mm to 12 mm. The results shown in Figure 9 also indicate an inverse effect of normal radius R on rigidities.

The effect of coil radius is positive, according to the results of Figure 10. In particular, the square symbol “□” describes the results of Figures 5 and 6.

Notice that there is no helix angle θ0 involved in the expression of compressive rigidity shown in (22). This is a big problem. However, the simulation result of spring bending rigidity indeed changes with varying helix angles as shown in Figure 11(a), totally different from the horizontal line representing the calculation result. Figure 11(b) shows that the calculation and simulation results meet at the helix angle of 7°, although diverged quickly before and after. Therefore, we can obtain a conclusion here that the expression used for calculation spring bending rigidity is of the best utility when helix angle θ0 is set around 7.

#### 6. Compression and Bending Stiffness Balancing

All simulations made above have verified the correctness of the deduced bending expression shown in (19). In order to find out the relationship between K and Kb, we substitute (22) into (18), and the result isIt seems that to modify the helix angle θ0 may be an effective way to change Kb while leaving K unchanged, because there is no helix angle θ0 involved in K. However, the results shown in Figure 11 denied the possibility of this alternation. Besides the material parameter ν, there is only the nominal radius R left. Therefore, to modify R is an effective way to adjust Kb while keeping the K unchanged.

All the notations mentioned above are listed in Table 2, as well as their meanings.

 Notation Meaning E Material elastic modulus G Material shear modulus H1, h2 The only two components of the outer diameter, corresponding to the mix and max deformation in each cross section Coordinate basis vectors of principle reference systems Coordinate basis vectors of general reference systems Ix, Iy, Iz Inertia moment around x-, y-, and z-axis K Spring compression stiffness Kb Spring bending stiffness L Spring wire length l1 Screw pitch of robot module l2 Screw pitch of spring M Acting moment General Descartes reference system Start point at bottom of spring Point along the axis of spring wire Principle coordinate system Local coordinate system R Spring wire radius R Spring nominal radius S Accumulated arc length Curvature vector Α Rotation angle about ζ-axis Β Rotation angle of the cross section around the tangential direction Γ Included angle υ Poisson ratio Θ Helix angle θ 0 Initial value of θ

#### 7. Experiments

For in-pipe robots, those turning a corner, an elastic connector is always required . In order to achieve effective power transmission from one module to another, the spring compression stiffness should be large enough. Meanwhile, the bending elastic connector must produce normal pressure against the inner pipe wall when passing through a corner. The normal pressure will generate friction that hinders the forward movement of in-pipe robot. The above calculation and simulation results tell us that the balance between compression stiffness and bending stiffness could be achieved, according to (24).

In this part, we applied some springs with different shape parameters to assemble the same robot. Then we make velocity test when passing through a bending pipe. The maximum velocity will express the best combination of spring shape parameters and also the best balance between compression and bending stiffness. Figure 12 shows the basic components of an in-pipe robot whose elastic connector can move smoothly through the bent pipe thanks to the spring. Figure 13 shows the experiment platform for the velocity test. It has in-pipe robot drive circuit and control algorithm for driving, as well as bent pipe for running environment.

During the experiment process, we kept spring material and normal radius constant, R = 10 mm, and changed the wire radius, helix angle, and active coil number one by one.

It can be seen from the experiment results that the spring compression and bending stiffness both increase dramatically with bigger coil radius. The maximum speed comes at the 6th parameter set in Table 3. The bending passing velocity increased initially with bigger K; however, it got slower with lasting increasing K. It indicates the same trends with calculation and simulation result.

 No. Coil radius r (mm) n Helix angle θ0 (rad) Compression stiffness K (N/m) Bending stiffness Kb (N/m) Velocity (mm/s) 1 0.5 6 0.1 201 259 3.3323 2 0.5 6 0.2 201 254 4.1148 3 0.5 6 0.3 201 246 4.7221 4 1 10 0.1 1925 2484 4.7324 5 1 10 0.2 1925 2436 5.2970 6 1 10 0.3 1925 2357 5.6852 7 1.5 14 0.1 6961 8982 2.6500 8 1.5 14 0.2 6961 8808 3.1622 9 1.5 14 0.3 6961 8525 3.6225

#### 8. Conclusions

In order to improve the characteristic of connecting spring applied to in-pipe robot, a balance is obtained between the compressive stiffness and the bending stiffness of the spring:(1)The detail of the spring linkage mechanism was introduced at first. It was simple, reliable, and really suitable for the application in in-pipe robot.(2)In order to get a higher bending passing speed, the spring linkage should be neither too soft nor too hard. The compression and bending rigidity need an optimization design.(3)The expression for spring bending rigidity calculation was deduced. Combining the simulation and experiment results, we found the suitable method to optimize the spring. The most sensitive parameter was the helix angle θ0.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported by the Key Research and Development Plan of Shandong Province (grant no. 2019GGX104044) and the Major Scientific and Technological Innovation Projects of Shandong Province (grant no. 2019JZZY020323).

1. Y. S. Zhang, S. Y. Jiang, X. W. Zhang, X. Y. Ruan, and D. M. Guo, “A variable-diameter capsule robot based on multiple wedge effects,” IEEE/ASME Transactions on Mechatronics, vol. 16, no. 2, pp. 241–254, 2011. View at: Publisher Site | Google Scholar
2. M. Ono and S. Kato, “A study of an earthworm type inspection robot movable in long pipes,” International Journal of Advanced Robotic Systems, vol. 7, no. 1, pp. 85–90, 2010. View at: Publisher Site | Google Scholar
3. M. M. Moghaddam, M. Arbabtafti, and A. Hadi, “In-pipe inspection crawler adaptable to the pipe interior diameter,” International Journal of Robotics & Automation, vol. 26, no. 2, pp. 135–145, 2011. View at: Publisher Site | Google Scholar
4. J. Chen, T. Chen, and Z. Q. Deng, “Dynamic properties on double driving articulated in-pipe robot inspecting system,” Applied Mechanics and Materials, vol. 128-129, pp. 1430–1433, 2012. View at: Publisher Site | Google Scholar
5. Z. Hu, Y. Wang, S. Yu, and G. Cui, “A method to enhance fidelity of force feedback control in virtual and human-robot micro interaction cardiovascular intervention surgery,” The International Journal of Advanced Manufacturing Technology, vol. 105, no. 12, pp. 4883–4897, 2019. View at: Publisher Site | Google Scholar
6. H. Tourajizadeh, M. Rezaei, and A. H. Sedigh, “Optimal control of screw in-pipe inspection robot with controllable pitch rate,” Journal of Intelligent & Robotic Systems, vol. 90, no. 3-4, pp. 269–286, 2018. View at: Publisher Site | Google Scholar
7. T. Sudhawiyangkul, K. Yoshida, S. I. Eom, and J.-w. Kim, “A novel bending microactuator with integrated flexible electro-rheological microvalves using an alternating pressure source for multi-actuator systems,” Microsystem Technologies, vol. 26, no. 5, pp. 1507–1519, 2020. View at: Publisher Site | Google Scholar
8. Q. Liu, J. Zhao, H. Zhu, W. Zheng, and Y. Yang, “A novel double bevel support structure for downhole robot,” Arabian Journal for Science and Engineering, vol. 44, no. 2, pp. 1069–1079, 2019. View at: Publisher Site | Google Scholar
9. S. Yaqub, A. Ali, M. Usman et al., “A spiral curve gait design for a modular snake robot moving on a pipe,” International Journal of Control, Automation and Systems, vol. 17, no. 10, pp. 2565–2573, 2019. View at: Publisher Site | Google Scholar
10. Y. Shi, Z. Mu, M. Cai, H. Song, and Y. Wang, “Advances in motion control of gas pipeline detection robot,” Science China Technological Sciences, vol. 63, no. 5, pp. 877-878, 2020. View at: Publisher Site | Google Scholar
11. F. Moreira, A. Abundis, M. Aguirre, J. Castillo, and P. A. Bhounsule, “An inchworm-inspired robot based on modular body, electronics and passive friction pads performing the two-anchor crawl gait,” Journal of Bionic Engineering, vol. 15, no. 5, pp. 820–826, 2018. View at: Publisher Site | Google Scholar
12. Z. Wang, S. Guo, Q. Fu, and J. Guo, “Characteristic evaluation of a magnetic-actuated microrobot in pipe with screw jet motion,” Microsystem Technologies, vol. 25, no. 2, pp. 719–727, 2019. View at: Publisher Site | Google Scholar
13. H. Schempf, E. Mutschler, A. Gavaert, G. Skoptsov, and W. Crowley, “Visual and nondestructive evaluation inspection of live gas mains using the explorer (TM) family of pipe robots,” Journal of Field Robotics, vol. 27, no. 3, pp. 217–249, 2010. View at: Publisher Site | Google Scholar
14. O. Adria, H. Streich, and J. Hertzberg, “Dynamic replanning in uncertain environments for a sewer inspection robot,” International Journal of Advanced Robotic Systems, vol. 1, no. 1, pp. 33–38, 2008. View at: Publisher Site | Google Scholar
15. S. Wakimoto, J. Nakajima, M. Takata, T. Kanda, and K. Suzumori, A Micro Snake-like Robot for Small Pipe Inspection, Institute of Electrical and Electronics Engineers Inc, Nagoya, Japan, 2003.
16. S. Hirose, H. Ohno, T. Mitsui, and K. Suyama, “Design of in-pipe inspection vehicles for 25,50,150 pipes,” in Proceedings of the IEEE International Conference on Robotics and Automation, pp. 2309–2314, Detroit, MI, USA, May 1999. View at: Google Scholar
17. C. Lin, Y. Hu, Y. Wei, and Z. Cai, “Characteristics analysis of supporting and locking mechanism based on the non-circular gear compound transmission,” Journal Of Mechanical Science And Technology, vol. 34, no. 6, pp. 2561–2571, 2020. View at: Publisher Site | Google Scholar
18. H. Takeshima and T. Takayama, “Development of a steerable in-pipe locomotive device with six braided tubes,” ROBOMECH Journal, vol. 5, no. 1, 2018. View at: Publisher Site | Google Scholar
19. M. Muramatsu, N. Namiki, R. Koyama, and Y. Suga, Autonomous Mobile Robot in Pipe for Piping Operations, Institute of Electrical and Electronics Engineers Inc, Takamatsu, Japan, 2000.
20. S. Iwashina, I. Hayashi, N. Iwatsuki, and K. Nakamura, Development of In-Pipe Operation Micro Robots, IEEE, Nagoya, Japan, 1994.
21. C. Q. Xu, Research on Design and Motion Stability of a Novel Creeping Micro In-Pipe Robot, National University of Defense Technology, Changsha, China, 2010.
22. I. Hayashi, N. Iwatsuki, and S. Iwashina, Running Characteristics of a Screw-Principle Microrobot in a Small Bent Pipe, IEEE, Nagoya, Japan, 1995.
23. L. Brown, J. Carrasco, S. Watson, and B. Lennox, “Elbow detection in pipes for autonomous navigation of inspection robots,” Journal of Intelligent & Robotic Systems, vol. 95, no. 2, pp. 527–541, 2019. View at: Publisher Site | Google Scholar
24. Q. Tu, Q. Liu, T. Ren, and Y. Li, “Obstacle crossing and traction performance of active and passive screw pipeline robots,” Journal of Mechanical Science and Technology, vol. 33, no. 5, pp. 2417–2427, 2019. View at: Publisher Site | Google Scholar
25. Q. Liu, J. Zhao, H. Zhu, and W. Zhang, “Mechanical model of drilling robot driven by the differential pressure of drilling fluid,” Arabian Journal for Science and Engineering, vol. 44, no. 2, pp. 1447–1458, 2019. View at: Publisher Site | Google Scholar
26. S. Perreault, P. Cardou, and C. Gosselin, “Approximate static balancing of a planar parallel cable-driven mechanism based on four-bar linkages and springs,” Mechanism and Machine Theory, vol. 79, pp. 64–79, 2014. View at: Publisher Site | Google Scholar
27. F. B. Liu and Y. M. Cheng, “The improved element-free galerkin method based on the nonsingular weight functions for inhomogeneous swelling of polymer gels,” International Journal of Applied Mechanics, vol. 10, Article ID 18500474, 2018. View at: Publisher Site | Google Scholar
28. J. Cheng, “Data analysis of the factors influencing the industrial land leasing in Shanghai based on mathematical models,” Mathematical Problems in Engineering, vol. 2020, Article ID 9346863, 11 pages, 2020. View at: Publisher Site | Google Scholar
29. L. Chen, H. P. Ma, and Y. M. Cheng, “Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems,” Chinese Physics, vol. 22, no. 5, Article ID 050202, 2013. View at: Publisher Site | Google Scholar
30. M. Peng and Y. Cheng, “A boundary element-free method (BEFM) for two-dimensional potential problems,” Engineering Analysis with Boundary Elements, vol. 33, no. 1, pp. 77–82, 2009. View at: Publisher Site | Google Scholar
31. K. M. Liew and Y. Cheng, “Complex variable boundary element-free method for two-dimensional elastodynamic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 49-52, pp. 3925–3933, 2009. View at: Publisher Site | Google Scholar
32. Z. Liu, G. Wei, and Z. Wang, “The radial basis reproducing kernel particle method for geometrically nonlinear problem of functionally graded materials,” Applied Mathematical Modelling, vol. 85, pp. 244–272, 2020. View at: Publisher Site | Google Scholar
33. Z. Liu, H. F. Gao, G. F. Wei, and Z. M. Wang, “The meshfree analysis of elasticity problem utilizing radial basisreproducing kernel particle method,” Results in Physics, vol. 17, Article ID 103037, 2020. View at: Publisher Site | Google Scholar
34. Z. Liu, G. Wei, and Z. Wang, “Numerical solution of functionally graded materials based on radial basis reproducing kernel particle method,” Engineering Analysis with Boundary Elements, vol. 111, pp. 32–43, 2020. View at: Publisher Site | Google Scholar
35. Z. Liu, G. F. Wei, and Z. M. Wang, “Numerical analysis of functionally graded materials using reproducing kernel particle method,” International Journal of Applied Mechanics, vol. 11, no. 6, Article ID 1950060, 2019. View at: Publisher Site | Google Scholar
36. B. H. Wang, Y. Q. Ma, and Y. M. Cheng, “The improved complex variable element-free Galerkin method for bending problem of thin plate on elastic foundations,” International Journal of Applied Mechanics, vol. 11, no. 10, Article ID 1950105, 2019. View at: Publisher Site | Google Scholar
37. Q. Wu, F. B. Liu, and Y. M. Cheng, “The interpolating element-free Galerkin method for three-dimensional elastoplasticity problems,” Engineering Analysis with Boundary Elements, vol. 115, pp. 156–167, 2020. View at: Publisher Site | Google Scholar
38. F. B. Liu, Q. Wu, and Y. M. Cheng, “A meshless method based on the nonsingular weight functions for elastoplastic large deformation problems,” International Journal of Applied Mechanics, vol. 11, no. 1, Article ID 1950006, 2019. View at: Publisher Site | Google Scholar
39. P. P. Peng, Q. Wu, and Y. M. Cheng, “The dimension splitting reproducing kernel particle method for three-dimensional potential problems,” International Journal for Numerical Methods in Engineering, vol. 121, no. 1, pp. 146–164, 2020. View at: Publisher Site | Google Scholar
40. A. E. H. Love, A Treatise on the Mathematical Theory of Elasticity, Dover Publications, New York, NY, USA, 1927.
41. N. Eqra, S. Taghvaei, and R. Vatankhah, “Optimal kinematic design of a single-DOF planar grasper based on metaheuristic optimization,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 41, no. 10, 2019. View at: Publisher Site | Google Scholar
42. Y. Z. Liu, Nonlinear Mechanics of Thin Elastic Rod-Theoretical Basis of Mechanical Model of DNA, Qtinghua University Press, Beijing, China, 2006.
43. Y. Miyazaki and K. Kondo, “Analytical solution of spatial elastica and its application to kinking problem,” International Journal of Solids and Structures, vol. 34, no. 27, pp. 3619–3636, 1997. View at: Publisher Site | Google Scholar
44. R. Huang, S. J. Zheng, Z. S. Liu, and T. Y. Ng, “Recent advances of the constitutive models of smart materials-hydrogels and shape memory polymers,” International Journal of Applied Mechanics, vol. 12, no. 2, Article ID 2050014, 2020. View at: Publisher Site | Google Scholar
45. B. Arun Sundaram, K. Kesavan, and S. Parivallal, “Recent advances in health monitoring and assessment of in-service oil and gas buried pipelines,” Journal of The Institution of Engineers (India): Series A, vol. 99, no. 4, pp. 729–740, 2018. View at: Publisher Site | Google Scholar

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.