The Stability of a Movable High-Strength Inverted-Triangular Steel Bridge
The overall stability of a movable high-strength inverted-triangular steel bridge is worth studying because of its new truss structure. In this study, an approach was proposed based on the stiffness equivalence principle in which the inverted-triangle truss structure was modeled as a thin-walled triangular beam. On this basis, the calculation of the critical load of elastic stability of a movable high-strength inverted-triangular steel bridge with variable rigidity at both ends and locally uniformly distributed load was carried out based on the energy theory, which was in good agreement with existing theories. A material performance test at BS700 was carried out to establish the material properties, and then a finite element model of the bridge was established, the results of which were compared with those of the experimental load test, in order to verify the accuracy of the finite element model. Considering material nonlinearity and geometric nonlinearity, nonlinear buckling analysis of the bridge was conducted and the factors influencing the bridge’s ultimate bearing capacity were analyzed.
High-strength steel with a yield strength of 460–690 MPa has been commonly used in construction projects in Japan, the United States, and some other countries [1–3]. In addition, the 1100 MPa-level high-strength steel has been used in some military bridges, such as FB48. Over the past few years, an increasing number of high-strength steel buildings have been built in China [4–6], with Q420 and Q460 as the main two types of high-strength steel being used. According to the research, the application of high-strength steel can effectively contribute to self-weight reduction and cost saving , which has made it more suitable for the design of movable emergency bridges. This is the reason why high-strength steel with a yield strength greater than 700 MPa (≥ 700 MPa) is used in some movable bridges . Movable high-strength inverted-triangular steel bridges developed with BS700 high-strength steel (= 700 MPa) are a new type of long-span bridge whose length and deck load can reach 51 m and 60 ton, respectively. This kind of bridge consists of two side bridge segments with variable cross-sections and four identical central bridge segments.
Apart from BS700 high-strength steel, a new inverted-triangular truss structure was used in this bridge. Two pieces of inverted-triangular truss structures made up two lanes, which were connected and integrated in between on the top using cross beams. Central bridge segments are shown in Figure 1. Normally, a rectangular truss system is used in a truss bridge, yet research findings suggested that a triangular treadway bridge performs better in saving materials and reducing weight than a rectangular truss bridge . Meanwhile, the triangular truss was usually employed on a limited basis for a number of structures such as roof truss girders, transmission towers, and crane booms . In some composite bridges the triangular trusses were also widely used . The mentioned triangular trusses were usually vertically triangular, but an inverted-triangular truss system emerged with the development of long-span spatial structures . Today the inverted-triangular truss is more commonly applied in composite bridges [12, 13]. Inverted-triangular trusses were used in the movable bridge designed in this study. To make it convenient for mechanized bridging, there was no connection to the bottom chord, and the top chord was also weakly connected, which offers a promising area for further studies on the bridge’s overall stability.
Considerable research on the stability of triangular or inverted-triangular truss structures has been done. The eigenvalue stability of truss structures was studied through critical point theory [14, 15]. Some laboratory-scale tests and full-scale testing of old steel truss tests were performed to get the load-carrying capacity of truss structures [16, 17]. A full-scale triangular composite truss structure specimen was tested and compared with finite element analysis , and the finite element model (FEM) was validated with experimental results. Recently the truss bridge’s linear buckling analysis and nonlinear stability analysis, which can consider material nonlinearity and initial geometric defects, were both done by the finite element method [18–20]. The stability ultimate load influence factors of truss structures were also studied. It was found that the joint rigid connection has more influence on the ultimate bearing capacity of triangular space-truss structure . Height-width ratio has more influence on the integral stability and the proposed height-width ratio is 0.75~0.875 . In comparison with the paper that has already been published, the main purpose of this work is to investigate the stability of a novel high-strength inverted-triangular steel truss bridge that has varied sections in its side bridge segments. In this study, the elastic stability of the movable inverted-triangular steel bridge was analyzed theoretically, and a calculation equation was obtained to predict its elastic load-carrying capacity. The material performance test of BS700 and the field test of load-carrying capacity were conducted. The bridge’s elastic stability load-carrying capacity was calculated using the finite element method and compared with the theoretical result. Both material nonlinearity and geometric nonlinearity were considered in analyzing the nonlinear stability of the bridge by the finite element method. Moreover, factors influencing the stability load-carrying capacity of the movable high-strength inverted-triangular steel bridge were analyzed.
2. Analysis of Elastic Stability Load-Carrying Capacity
Analysis of the elastic stability load-carrying capacity lays the foundation for analyzing the nonlinear stability load-carrying capacity. To perform the theoretical analysis, the truss structure is usually modeled as a beam, following the principle of equivalent stiffness . In terms of the truss structure of a movable high-strength inverted-triangular steel bridge, the stiffness parameters are mainly axial rigidity , in-plane flexural rigidity , out-plane flexural rigidity , and torsional rigidity .
2.1. Simplification to Thin-Walled Triangular Beam Based on Equivalent Stiffness
The simplified diagram of the cross section of the movable high-strength inverted-triangular steel bridge is shown in Figure 2.
There are two different sections on the top chord whose areas are and . The area of the bottom chord is . The position of a section centroid is determined according to =0. Supposing that point is the centroid and and are the distances between the centroid and the top and bottom chords, thenThe equivalent area isNeglecting the contribution of the component’s cross section to the overall moment of inertia moment of inertia, the equivalent moment of inertia isTherefore, the axial rigidity , in-plane equivalent flexural rigidity , and out-plane flexural rigidity could be calculated, respectively. In terms of torsional rigidity , since the web member connecting the chord primarily bears the shear force, its contribution to torsion could not be ignored. There are two types of web members used in the movable high-strength inverted-triangular steel bridge: one is connected to the bottom and top chord I and the other is connected to the bottom chord and top chord II (Figure 1). For simplification, the middle chord (top chord II) and web member are evenly distributed to both sides, forming two sections of trusses on the left and right sides (Figure 3).
One section of the movable high-strength inverted-triangular steel bridge is taken as a calculation unit and converted into a thin plate whose thickness is te at each side. The whole section to be calculated is converted into a thin-walled triangular beam, with its torsion strength calculated. If the intersegment length is and height is , according to the equation of strain energy density, the strain energy of each piece of thin plate isThe thickness can be determined based on the condition that strain energy and of a truss with a length of d are equal. Suppose the total shear force that results from shear flow acting upon the truss plane is . A segment of the bridge is taken as the research object, as shown in Figure 4.
The internal force generated by a slanted web member with the action of shear force isIf the length and area of the web member are and , and it bears a constant force , its strain energy isThe chord’s axial force varies linearly within the length of . Assuming the maximum force , the chord’s strain energy isThe strain energy of each piece of truss unit with a length of d can be obtained according to (6) and (7):
Top chord:where .
Slanted web member:where and is the area of central slanted web member.
Bottom chord:where .
Vertical web member:where is the area of the slanted web member.
The equation for the piece of truss unit could be calculated by summing up the strain energy of the aforementioned four chords and web members:The following equation is obtained according to the principle of equal strain energy:Then, the thickness of equivalent thin-walled beam isAt this point, the converted thicknesses of the triangle’s three sides, , , and , are calculated. The equivalent torsion moment of inertia is obtained according to the constant equation of free torsion of a thin-walled closed section of arbitrary shape .where is the area of the closed section formed by the medians of the thin-walled triangle; , , and are the widths of three thin-walled plates, i.e., lengths of the three sides of the triangle’s cross section; and , , and are the converted thicknesses of three sides. , the torsional rigidity of the thin-walled triangular beam simplified from the movable high-strength inverted-triangular steel bridge, can be obtained based on (14) and (15).
The equivalent rigidity can be calculated using MATLAB according to this theoretical derivation . Due to the fact that the rigidities of both end units of the movable high-strength inverted-triangular steel bridge are variable while the central rigidity remains unchanged, two programs are compiled to calculate their respective equivalent rigidities, which can be calculated by entering the bridge’s height and width and the sizes of various chords and bars into the programs. With this calculation, the moments of area of the sections where central rigidity remained unchanged were , , and .
2.2. Critical Load of Linear-Elastic Stability
It is difficult to find the solution to the torsional rigidity if the movable high-strength inverted-triangle steel bridge is simplified into a thin-walled triangular beam, since the bridge’s rigidity varies at the two ends while the central rigidity stays the same. To this end, it is solved using the energy method . The axial compressive strain energy and shear strain energy are too low to be taken into account. In addition, a study made by Liu  suggests that the warping strain energy of a thin-walled triangular beam is negligible. This strain energy contains lateral bending strain energy , in-plane bending strain energy , and pure torsion strain energy .
Lateral bending strain energy can be obtained according to the strain energy equation .In-plane bending strain energy isFree torsion strain energy iswhere is the lateral displacement, is the vertical displacement, and is the torsion angle.
Therefore, total strain energy isChen  reported that the flexural bars’ elastic flexural-tensional buckling critical flexural moment was increased to some extent by the effect of large counter-arch due to flexural deformation. Therefore, ignoring the impact of in-plane flexure would result in a slightly smaller critical load and a more conservative result. Equation (19) is thus reduced toWithout considering residual stress, the potential energy of external force can be calculated with the following equation when any section is subjected to the combined action of axial force and flexural moment .The potential energy of the external force should be (21) combined with the strain energy generated by the changes of the load’s action spot. Supposing the distance between the load’s action spot and shearing center is (Figure 5), the potential generated isHence, the total potential energy of the external force iswhere is the section asymmetry coefficient and because the section is symmetrical. is the distance between the shearing center and the centroid. According to our calculation, .
Total potential energy isThe integral of the total potential energy should be taken in segments as the height of the movable high-strength inverted-triangular steel bridge varies linearly at both ends and the uniformly distributed load acts locally upon the mid-span (Figure 6). Due to the linear variation of height, supposing the flexural rigidity and torsional rigidity are also subject to linear variations, the maximum rigidities are and , respectively, and the minimum rigidities are and , respectively. Assume the width of the local uniformly distributed load, the length of the segment with variable rigidity, the intensity of uniformly distributed load, and the span of the simply supported bridge are , , , and , respectively.
The expressions of segmented load and rigidity are shown as follows.Therefore, the integral of total potential energy is performed on five intervals along the beam’s lengthwise direction; namely,The Rayleigh–Ritz method is used to solve the critical flexural moment . Suppose the shape function that matches the boundary conditions isTotal potential energy isAfter the segment integral, total potential energy iswhere are the parameters related to , , and .
According to the minimal potential energy principle, we haveThree parameters are introduced:The buckling condition of the equivalent thin-walled triangular beam isTo separate the load from the equation above, several coefficients, i.e., , and , are introduced:The critical load is obtained:To verify the correctness of this expression, a special circumstance, where load is uniformly distributed along the full span and rigidity remained unchanged along the beam’s axial direction, was modeled. In this case, the critical load is as follows according to (42):According to Chen , the elastic critical flexural moment with full-span uniform load distribution isWithout considering warping the moment of inertia, the result of (43) is consistent with that of (44). Not considering the self-weight, the elastic critical load of the movable high-strength inverted-triangular steel bridge was 2300.6 kN, which was calculated by substituting the data obtained based on the strain energy principle into (43).
3. Finite Element Analysis of Elastic Stability Load-Carrying Capacity
To further verify the accuracy of the equation for the elastic stability ultimate load, a model was built using the finite element method to perform elastic buckling analysis.
3.1. Establishment and Verification of the Finite Element Model
The beam element 188 of the ANSYS software, which can import a user’s section data, was adopted in the finite element model. In order to analyze many modes used for studying the influence factors on stability ultimate load capacity, the model uses many parameters which can be easily changed in an APDL document (Analysis Parameter Design Language) . The whole finite element model contains a basic model, eigenvalue buckling analysis, and the results of picking-up. The material model made use of the multilinear hardening model.
In order to get the material model data, uniaxial tension tests were carried out as shown in Figure 7(a). No obvious plastic flow or necking was observed in the tests. The failure planes were nearly placed along the 45° cross section of the specimens as shown in Figure 7(b). Shear failure is dominant in the material.
As shown in Figure 7, the complete stress-strain curve of the BS700 steel is gently strain-hardening. According to the test data, the elastic modulus of the material is 210.9 GPa. The nominal yield stress is 734 MPa and the nominal yield strain is 0.00426. For the material model we chose the multilinear isotropic hardening (MISO) model. The MISO model needs points, as shown in Table 1.
3.2. Finite Element Model
In order to verify the basic finite element model validity, the load capacity test of the bridge was performed. A 51 m-long bridge, where the bearing length at both ends was 2 m and clearance height under the span was 0.6 m, was constructed. Two sections of planks were placed at the middle of bridge span to simulate two decks. Steel plates weighing 25.42 t→48.97 t→58.42 t→63.15 t→67.87 t→71.17 t→74.56 t→80.84 t (1 t≈10 kN) were loaded to and piled up on the planks. There were three deformation points at L/6, 5L/6, and L/2 (L is 51 m). Measuring points and the loading diagram of the test and field test images are shown in Figures 8 and 9, respectively.
The test results compared with finite element analysis (FEA) results are shown in Table 2.
According to Table 2 it can be declared that the finite element model can simulate the bridge accurately; the errors were below 5%, excepting a few points.
3.3. Eigenvalue Buckling Analysis
The typical equation for eigenvalue buckling analysis was where is the matrix of elastic rigidity; is the matrix of stress rigidity; is the characteristic displacement vector; and is the eigenvalue.
Without considering self-weight, the eigenvalue obtained from eigenvalue buckling analysis multiplied by the load applied is the eigenvalue buckling load. In fact, the load kept changing, so thatWhen the self-weight was not considered, the first-order eigenvalue buckling mode was obtained, as shown in Figure 10. The first-order eigenvalue buckling load was 2405.4 kN, which was 4.6% larger than the 2300.6 kN load that was derived by the thin-walled beam equivalent theory in this study. This suggests a highly accurate value of the eigenvalue buckling load of the movable high-strength inverted-triangular steel bridge was derived in this study.
4. Nonlinear Ultimate Bearing Capacity of the Bridge
4.1. Theories on Nonlinear Analysis
Values of bridge’s elastic ultimate bearing capacity obtained by theoretical derivation and eigenvalue buckling analysis were usually so large that they could only be used for qualitative analysis of the structure’s load-carrying capacity instead of genuinely reflecting it. Therefore, nonlinear buckling analysis was required. Considering factors such as the structure’s geometric defect and the nonlinearity of materials, a nonlinear buckling analysis was able to more truthfully present the structure’s load-carrying capacity. In consideration of the material-geometry dual nonlinearity, the basic buckling equation could be converted intowhere is the small-displacement elastic stiffness matrix, is the initial stress stiffness matrix, is the large displacement elastic stiffness matrix, also known as the initial strain matrix, and and are nodal displacement and load vector, respectively. Solutions were found by means of the iteration method by gradually applying load increments. Structure stiffness changed as load rose. The structure reached its maximum load-carrying capacity when . Ultimate bearing capacity was tracked using the arc-length method.
4.2. Ultimate Bearing Capacity of the Bridge
The geometric deflection was introduced to the analysis model by the consistent-deflection-mode method, which simulates deflection distribution by the lowest buckling mode shape. The maximal displacement of buckling analysis was multiplied by a coefficient to amplify to the 1/1000 span of bridge.
Based on the above finite element model and considering material and geometric dual nonlinearity factors, the stability ultimate load of the bridge can be obtained by using the arc-length method . Firstly, an elastic buckling analysis was performed when the load was 600 kN, and then the elastic buckling load, which was the top limit load of bridge, was found. Secondly, setting the elastic buckling load as the initial load, the accumulated loads were then introduced gradually until the solution was emanative. The load-displacement curve can be traced through time post-process (post 26). The highest point of the load-displacement curve reveals the stability ultimate load capacity of a movable high-strength inverted-triangular steel truss bridge. The load-displacement curves are shown in Figure 11.
From Figure 11 we can obtain the stability ultimate load capacity () of the bridge, which is 1637.6 kN. The load-displacement curve is approximately linear, and the vertical displacement is bigger than the lateral displacement before the load arrives at the ultimate load capacity. Therefore, the bridge’s deformation is mainly in the moment plane; it is a kind of bending buckling when the load is small. When the load arrives at the stability ultimate load capacity of the bridge, the lateral deformation becomes bigger, and the bridge exhibits true bending-torsion buckling.
5. Stability Ultimate Load Influence Factors Analysis
There are many factors which can affect the stability ultimate load capacity of a movable high-strength inverted-triangular steel truss bridge, such as bridge span, bridge height, track width, and member’s section. In order to study the influence of these factors on the stability ultimate capacity of the bridge, many FEA (finite element analysis) models with different parameters were analyzed.
5.1. Influence of Height and Track Width
The height and track width are two important parameters for design. When the analysis was performed, one parameter was changed each time. Through FEA analysis the stability ultimate load capacity was obtained. Based on the actual bridge’s ultimate load capacity (which is 1637.6 kN) the change percentage can be determined . The results are shown in Figure 12.
From Figure 12 a conclusion can be made that the height has more influence on the stability ultimate load capacity than the track width. The reason for this is that the bridge’s height has more influence on section bending stiffness than track width does. Meanwhile, the track width changes but the whole bridge’s width does not change; the lateral torsion stiffness changes a little. Thus, track width has less influence on the stability ultimate load capacity of the bridge.
5.2. Influence of Chord
The top and bottom chord are the important members in the bridge. How a chord influences the stability ultimate load was also studied in this paper. The top chord has two different box sections, which were labeled chord I and chord II. The top chord area changed at the same time of the width and height changes of the section. The multiplied area of the top chord can be found through adjusting the width and height of section gradually. The bottom chord is a kind of complex section. The multiplied area of the bottom chord can be found through adjusting the main dimensions of the section gradually. The results are shown in Figure 13.
From Figure 13 it can be concluded that the top chord has more influence on the stability load capacity than the bottom chord. The top chord I has more influence than chord II. The bottom chord has a little influence on the stability load capacity that changes below 5%. Even top chord I has more influence than height when both two parameters change by the same percentage, as shown by comparing Figure 12 with Figure 13.
5.3. Influence of Web Member
The web member mainly bears shear force in the bridge structure but also provides some bending stiffness. There are two different web members shown in Figure 1, the sections of which are box shapes. Thus, the influence of the web member on the stability ultimate load was studied using the same method as the chord. The results are shown in Figure 14.
According to Figure 14, the web member has a little influence on the stability ultimate load capacity of the bridge. The influence of the upright web member is more than the slanted web member.
This paper investigated the stability of a movable high-strength inverted-triangular steel bridge by the theoretical, experimental, and finite element methods. From this study the following conclusions can be drawn:
The movable high-strength inverted-triangular steel bridge can be modeled as equivalent to a thin-walled triangular beam according to the stiffness equivalence principle. The equivalent rigidities of the end units which are variable and the central units rigidity which remains unchanged were calculated.
The equation which can calculate the elastic stability critical load of the inverted-triangular steel bridge with variable rigidities at both ends was put forward. This equation was well agreed with existing theoretical calculation. The calculated elastic stability load of the bridge was 2300.6 kN by substituting equivalent rigidity into the equation.
The material performance test of BS700 (= 700 MPa) high-strength steel was carried out, and the material properties used for finite element analysis were obtained. The load capacity test of the movable high-strength inverted-triangular steel bridge whose length is 51m was performed.
A finite element model for the elastic buckling analysis of the movable high-strength inverted-triangular steel bridge was established. The finite element simulated results well agreed with field test results. The elastic buckling critical load of the bridge is 2405.4 kN, which further verifies the accuracy of the equation of the bridge’s elastic stability critical load proposed in this study.
The nonlinear buckling analysis model was put forward with consideration of material and geometric nonlinearity. The material model made use of the multilinear hardening model. The geometric deflection was introduced to the analysis model by the consistent-deflection-mode method. According to calculations, the bridge’s ultimate bearing capacity is 1637.6 kN which is lower than the elastic stability critical load.
Analysis of factors influencing the bridge’s ultimate bearing capacity suggests that top chord and bridge height are two major factors, of which top chord 1 exerts the dominant influence. Track width, bottom chord, and web member have less influence on the bridge’s ultimate bearing capacity.
|:||In-plane flexural rigidity|
|:||Out-plane flexural rigidity|
|:||The top chord I area|
|:||The top chord II area|
|:||The bottom chord area|
|:||The distances between the centroid and the top chord|
|:||The distances between the centroid and the bottom chord|
|:||The thickness of equivalent thin-walled beam|
|:||The intersegment length|
|:||The intersegment height|
|:||The strain energy of each piece of thin plate|
|:||The total shear force that results from shear flow acting upon the truss plane|
|:||The top chord strain energy of each piece of truss unit|
|:||The bottom chord strain energy of each piece of truss unit|
|:||The vertical web member strain energy of each piece of truss unit|
|:||The total strain energy of the truss unit|
|:||Lateral bending strain energy|
|:||In-plane bending strain energy|
|:||Pure torsion strain energy|
|U:||The total strain energy|
|:||Total potential energy|
|:||The lateral displacement|
|:||The vertical displacement|
|:||The torsion angle|
|:||The potential energy of external force|
|:||The section asymmetry coefficient|
|:||The distance between the shearing center and the top side of equivalent thin-walled beam|
|:||The width of the local uniformly distributed load|
|:||The length of the segment with variable rigidity|
|:||The intensity of uniformly distributed load|
|:||The span of the simply supported bridge|
|:||The parameters related to , , and|
|:||The shape function parameters of the lateral displacement and the torsion angle|
|:||The parameters related to , , , ,,,and|
|:||The parameters related to , , , and|
|:||The critical load|
|:||The small-displacement elastic stiffness matrix|
|:||The initial stress stiffness matrix|
|:||The initial strain matrix|
|:||Nodal displacement vector|
|:||Nodal load vector|
|:||The stability ultimate load capacity of the bridge.|
The bridge’s height, width, sizes of various chords, and bars used to support the moments of area of the sections of this study are available from the corresponding author upon request. The parameters related to L, B, and m used to support the elastic critical flexural moment of this study are available from the corresponding author upon request. The material model used to support the stress-strain curve of BS700 steel of this study is currently under embargo while the research findings are commercialized. Requests for data, [ months] after publication of this article, will be considered by the corresponding author. 1/1000 span of bridge is the data used to support maximal displacement of buckling analysis of this study and is included within the article. It can be found in [13, 17].
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The authors wish to acknowledge the financial support from the project of the Major State Basic Research Development of China (973 Program, No. 2014CB046801) and China Postdoctoral Science Foundation (Grant No. 2017M623403).
G. A. Rosier and J. E. Groll, “High strength quenched and tempered steels in structures,” Steel Construction, vol. 21, no. 3, p. 13, 1987.View at: Google Scholar
G. Pocock, “High strength steel use in Australia, Japan and the US,” Structural Engineering International, vol. 84, no. 21, pp. 27–30, 2006.View at: Google Scholar
Z. M. Chen, Y. L. Zhang, and M. X. Peng, “Application of high-strength steel and thick steel plates to CCTV new site building,” Steel Construction, vol. 24, no. 2, pp. 34–38, 2009.View at: Google Scholar
S. H. Zhou, Z. Y. Zhu, and W. H. Qi, “Structural design on the project of Phoenix International Media Center,” Building Structure, vol. 41, pp. 56–62, 2001.View at: Google Scholar
G. Shi, H. Y. Ban, and Y. J. Shi, “Overview of research progress for high strength steel structures,” Engineering Mechanics, vol. 30, no. 1, pp. 1–13, 2013.View at: Google Scholar
M. K. Gou and L. Tao, “The application of high strength steel (=700 MPa) to the movable bridge,” Steel Structure, vol. 17, no. 5, p. 6, 2002.View at: Google Scholar
G. E. Blandford, “Progressive failure analysis of inelastic space truss structures,” Computer Structure, vol. 458, no. 4, pp. 981–990, 1996.View at: Google Scholar
G. H. Li, Bridge Stability and Vibration, Chinese Railway Press, Beijing, China, 1st edition, 2002.
Y. P. Tan and S. Q. Wang, “Buckling analysis of a long-span fast launching simple-supported bridge,” Journal of Shenzhen University science & engineering, vol. 18, no. 1, pp. 64–70, 2001.View at: Google Scholar
M. Tong, F. Mao, and H. Qiu, “Structural stability analysis for truss bridge,” in Proceedings of the International Workshop on Automobile, Power and Energy Engineering, APEE 2011, pp. 546–553, China, April 2011.View at: Google Scholar
F. N. Sun, Research on the integral stability of large span truss with inverted triangle section, Gui Zhou University, Guiyang, 2008.
Y. S. Liu, Analysis for out-plane global stability of spatially latticed rigid frame with triangular section, Harbin Institute of Technology, Harbin, 2000.
J. Chen, Stability of Steel Structures Theory and Design, Science Press, Beijing, China, 5th edition, 2011.
X. M. Zhang, Matlab Language and Application Case, South-east University Press, Nanjing, China, 1st edition, 2010.
C. H. Zhang, ANSYS 16.1 for the Analysis of Structure Engineering, China Machinery Press, Beijing, China, 4th edition, 2016.