In order to obtain the optimal section for super-large cross-section loess tunnels with different overburdens, the ANSYS software is applied in this paper.
Based on the calculation nonconvergence criteria and Mohr-Coulomb criteria, the static stability of loess tunnel section is analyzed by the finite element static strength reduction method. According to the safety factors of rock mass surrounding tunnel in the case of critical failure, the safety factors of super-large cross-section loess tunnel is discussed with different section forms (rectangular cross-section, circular cross-section, horseshoe cross-section and curve wall cross-section) and different overburdens under gravity.
The results show that the safety factors of circular section and curve wall cross-section are bigger than the safety factors of the horseshoe cross-section and rectangular cross-section. The curve wall cross-section is considered to be the optimal section because of the fact that the force around the vault and the arch bottom is uniform and symmetrical, and the maximum vertical displacement of the lining is small. The horseshoe cross-section should be avoided in the loess tunnel, because the safety factor of horseshoe cross-section is so small, and stress concentration phenomenon is obvious in the case of deep buried tunnel.
Open Peer Review Details | |||
---|---|---|---|
Manuscript submitted on 02-06-2017 |
Original Manuscript | Optimal Sections of Tunnels’ Cross Sections with Different Overburdens |
According to the existing literature, many experts and scholars have put forward many kinds of optimization design methods for the forms of tunnel section, and verified a variety of optimization design methods for the theoretical and practical aspects. Bahmanikashkooli et al. [1A. Bahmanikashkooli, M. Zare, and B. Safarpour, "Application of particle swarm optimization algorithm for computing critical depth of horseshoe cross section tunnel", APCBEE Procedia, vol. 9, pp. 207-211, 2014.
[http://dx.doi.org/10.1016/j.apcbee.2014.01.037] ] calculated the critical depth of horseshoe cross-section tunnel by using particle swarm optimization. Xu Zemin et al. [2Z.M. Xu, R.Q. Hunag, and S.T. Wang, "Tunnel classifying in light of depth", Chin. J. Geol. Hazard Control, vol. 11, no. 4, pp. 8-13, 2000.] further divided the tunnel into shallow buried tunnel, deep buried tunnel and super deep buried tunnel according to the bearing capacity of soil mass surrounding tunnel, the deformation and failure position and the deformation mode, and considered the ratio of the vertical principal stress and the minimum horizontal principal stress. Wang et al. [3M.N. Wang, J. Guo, and L.S. Luo, "Study of critical buried depth of large cross-section loess tunnel for high speed railway", Yantu Lixue, vol. 31, no. 4, pp. 1157-1162, 2010.] preliminarily determined the boundary value for the shallow and deep overburden of the large section loess tunnel, and obtained that the shallow buried was less than 11m, and that the boundary value of shallow and deep buried was 40m - 60m, and that the deep overburden was more than 60m. Mobaraki and Vaghefi [4B. Mobaraki, and M. Vaghefi, "Numerical study of the depth and cross-sectional shape of tunnel under surface explosion", Tunn. Undergr. Space Technol., vol. 47, pp. 114-122, 2015.
[http://dx.doi.org/10.1016/j.tust.2015.01.003] ] studied the dynamic response of the tunnel when the overburden of the tunnel was 3m, 5m, 7m, 10.5m and 14m, and the tunnel with different cross-section forms was analyzed by using LS-DYNA. Lu et al. [5A.Z. Lu, H.Y. Chen, and Y. Qin, "Shape optimization of the support section of a tunnel at great depths", Comput. Geotech., vol. 61, pp. 190-197, 2014.
[http://dx.doi.org/10.1016/j.compgeo.2014.05.011] ] studied the optimal shape of the support section for deep buried tunnel, and obtained the optimal shape of the support section suggesting that the minimum value of the objective function based on the mixed compensation function method could satisfy a given constraint condition by using the conformal mapping method of plane elastic complex variable function to optimize the support section. Lv [6A.Z. Lv, "Study of critical buried depth of large cross-section loess tunnel for high speed railway", J. China coal soc., vol. 22, no. 5, pp. 1157-1162, 1997.] obtained optimal cross-section by comparing the damage degree of soil mass surrounding tunnel with the elastic assumption, and sought optimal support cross-section shape with the assumption of the lining and rock mass contact smoothly,in other words, the stress vector and normal displacement on
the contact surface satisfy the continuous condition, but the sliding between surrounding rock and lining may occur [7A.Z. Lv, H.Y. Chen, and Y. Tan, "Shape optimization of tunnel support section under contact condition of pure slip", Chin J. Jock Mech. Eng., vol. 33, no. 8, pp. 1563-1571, 2014.]. Theoretical study about stability of tunnel cross-section is carried out by Peila [8D. Peila, "A theoretical study of reinforcement influence on the stability of a tunnel face", Geotech. Geol. Eng., vol. 12, no. 3, pp. 145-168, 1994.
[http://dx.doi.org/10.1007/BF00426984] ], and the corresponding force, plastic zone and displacement were respectively analyzed. Cheng et al. [9X.S. Cheng, W. Shi, and R.R. Tian, "Safety assessment of lined loess tunnel under earthquake", Electron. J. Geotech. Eng., vol. 19, no. I, pp. 1903-1915, 2014. [J].-11X.S. Cheng, H.Y. Zhu, and R.R. Tian, "Soil mass stability analysis of circular unlined loess tunnel under earthquake action", International Symposium on Geomechanics and Geotechnics. From Micro to Macro,, pp. 1127-1134, 2010.
[http://dx.doi.org/10.1201/b10528-169] ] proposed two new methods of seismic safety factor and fracture surface to evaluate the seismic stability of tunnel. Zheng et al. [12Y.R. Zheng, C.Y. Qiu, and H. Zhang, "Exploration of stability analysis methods for surrounding rocks of soil tunnel", Chin. J. Jock Mech. Eng., vol. 27, no. 10, pp. 1968-1980, 2008., 13Y.R. Zheng, H.L. Ye, and Q. Xiao, "The analysis about the stability of seismic slope and tunnel based on full dynamic analysis method", J. Disaster Prev Mitigation Eng., vol. 30, no. S1, pp. 279-285, 2010.] obtained the critical tensile failure state of tunnel by reduction parameters of soils shear strength. Cheng et al. [14X.S. Cheng, R.R. Tian, and J.L. Wang, "Static and seismic stability analysis of soil mass surrounding tunnel about large-span unlined egg-shaped loess tunnel", Tumu Gongcheng Xuebao, vol. 43, no. S1, pp. 582-587, 2010.] analyzed static and dynamic stability of long-span egg shaped loess tunnel without lining by strength reduction method, and obtained the influence on loess tunnel span, the thickness of covering soil and seismic intensity on safety factor and stability. Li et al. [15H. Li, X.Y. Zhu, and Y.F. Xu, "Influence of tunnel cross-section shape on stress and displacement of surrounding rock", Tunnel Constr, vol. 29, no. 1, pp. 38-44, 2009.] investigated distribution regular of displacement and stress concentration for different cross-section forms excavated by using FLAC3D. Dong et al. [16S.M. Dong, Q.C. Xin, and S.S. Lu, "The influence of the tunnel cross-section shape on stability of surrounding rocks", China Rural Water Hydropower, vol. 01, no. 102, pp. 104-107, 2011.] analyzed and calculated different cross-sections of tunnel excavated by using ANSYS. All things considered, they believe egg shaped section is the most reasonable form.
In order to get a safe, economical and reasonable section form, based on the above research, this paper will analyze the stability of different section forms and different overburdens super-large cross-section loess tunnel by using the static strength reduction method and finite element analysis software ANSYS. This provides a theoretical reference for the practical engineering application of the optimal design for the loess tunnel section.
Zienkiewicz et al. [17O.C. Zienkiewicz, C. Humpheson, and R.W. Lewis, "Associated and non-associated visco-plasticity and plasticity in soil mechanics", Geotechnique, vol. 25, no. 4, pp. 671-689, 1975.
[http://dx.doi.org/10.1680/geot.1975.25.4.671] ] proposed the concept of shear strength reduction coefficient for the first time in the finite element analysis of elasto-plastic soil mass. The finite element strength reduction method constantly reduced the shear strength index of the loess during the finite element model analysis until it reaches the critical failure state, and so corresponding safety factor is obtained. The traditional stability limit equilibrium method is based on the Coulomb criterion, and the safety factor is calculated by the force balance. The safety coefficient is defined as the ratio of the shear strength of slip surface τ_{c} and actual shearτ. It can be expressed as
(1) |
where η is the safety coefficient, c and φ is respectively cohesion force and internal friction angle. Eq. (1) with both sides divided by the safety coefficient η, that is
(2) |
where,
From Eq. (2), the concept of finite element static strength reduction method can be developed. The method can reflect the distribution situation of plastic zone vividly.
The ideal elastic-plastic model is used in the constitutive model of rock and soil mass in the finite element strength reduction method. It is important to select a reasonable yield criterion. At present, Coulomb criterion is widely adopted to reduce the shear strength parameters (Cohesion force c and internal friction angle φ) of rock/soil mass surrounding tunnel until the calculation does not converge by using the static analysis model. So, the safety factor of the tunnel without lining under earthquake is obtained [18Y.R. Zheng, Z.Y. Chen, and G.X. Wang, Eng. Treat. Slope Landslide., China Communications Press, 2011.-21X.S. Cheng, and Y.R. Zheng, "Calculation discussion about safety factor of unlined loess tunnel wall rock structure under earthquake", Yantu Lixue, vol. 32, no. 3, pp. 761-766, 2011.]. That is:
(3) |
where I_{1,}J_{2,} and θ_{σ} are the first invariant of the stress tensor, the second invariant of the stress partial tensor and the Lode angle, respectively.
Due to the yield surface of Mohr Coulomb, yield criterion is an irregular hexagonal in the π plane, the spires, edges and corners bring inconvenience to numerical calculation. So the simplified method is adopted by finite element method by using generalized Mises yield criterion. This is D-P yield criterion, and the yield surface is a circle. It can be expressed as
(4) |
where α and k areconstant which is related to the rock and soil material cohesion force c and internal friction angle φ. It can be expressed as
(5) |
It was assumed that the soil mass surrounding tunnel of the whole super-large section loess tunnel structure is isotropic and ideal elastic-plastic, and the whole stress and deformation of the tunnel structure is plane strain problem and only affected by gravity. The soil mass yield criterion is based on the Mohr-Coulomb criterion. The lining concrete surrounding the tunnel is considered an elastic material. The index of physical and mechanical property for the super-large cross-section loess tunnel is shown in Table 1.
Optimal cross-section form is obtained by analyzing different overburdens because the pressure of the soil mass is affected by tunnel overburden. As an example, the depth of super shallow buried tunnel is 10m, shallow buried is 30m, deep buried is 80m. Optimal cross-section form of super-large cross-section with different tunnel overburdens is obtained by comparing the safety factor, the critical strain nephogram, the first principal stress nephogram and maximum vertical displacement nephogram of lining for super-large cross-section loess tunnel with super shallow buried, shallow buried and deep buried by using ANSYS.
The stability analysis of the super-large cross-section loess tunnel is studied with different tunnel overburdens. The cross-section form for the super-large cross-section loess tunnel is divided into the rectangular cross-section, the circular cross-section, the horseshoe cross-section and the curve wall cross-section. The model size of each cross-section of the tunnel is 15m×15m (width×height). The calculation range is 5 times of the tunnel diameter on both sides (left side and right side) downward side on super-large cross-section loess tunnel. The tunnel overburden is 30m. The boundary conditions are: the left and right boundaries are constrained by the horizontal direction, and the bottom boundary is constrained by the vertical direction, and the top is free surface, which is free from any restriction. The calculation model of different section forms of super large section loess tunnel is shown in Table 2.
The safety factor for loess tunnel under the action of gravity can be obtained by using static finite element strength reduction method. The soil mass thickness with 800mm outside initial lining is taken as the shear strength reduction zone and the same material parameters is used by the reduction zone and soil mass. The stability of different cross-section forms with the same conditions is analyzed by comparing the safety factor, and the first principal stress and the maximum vertical displacement of the tunnel lining in the case of super shallow buried tunnel. According to the boundary value of deep and shallow buried large cross-section loess tunnel, the overburden thickness of super shallow buried tunnel is 10m. The safety factor and critical strain nephogram in the case of super shallow buried loess tunnel with different cross-section is shown in Table 3.
From Table 3, the safety factor of super-large cross-section loess tunnel with different cross-sections forms can be obtained by strength reduction method. The safety factor of circular super-large cross-section is 1.96. The safety factor of horseshoe super-large cross-section is 1.74. The safety factor of rectangular super-large cross-section is 1.47. The safety factor of curve wall super-large cross-section is 1.75. This shows that the safety factors are ordered from large to small, namely, circular, curved wall, horseshoe and rectangle cross-section. The safety factors of the horseshoe and curve wall cross-section is almost the same, but they differ greatly from the safety factor of rectangular cross-section.
From critical strain nephogram, in the case of super shallow buried tunnel, strain in the arch feet and the vault are more concentrated, and the maximum critical strain occurs at the bottom of the arch, and the range of cross-section strain increases with the increase in the safety factor. The critical first principal stress nephogram in the case of super shallow buried loess tunnel with different cross-sections is shown in Table 4.
From Table 4, tensile stress appears in the vault and arch bottom, the maximum tensile stress occurs at the arch feet, and the range of stress nephogram increases with the increase in the safety factor. This shows that the force cross-section bore increases with the increase in the safety factor. The stress distribution of circular cross-section is more uniform, and the stress of rectangle and horseshoe is more concentrated. The maximum vertical displacement nephogram of lining in the case of super shallow buried loess tunnel with different cross-sections is shown in Table 4.
Table 5 shows that the maximum vertical displacement occurs in the vault. The maximum vertical displacement of lining for circular cross-section is 32.15mm. The maximum vertical displacement of lining for horseshoe cross-section is 36.52mm. The maximum vertical displacement of lining for rectangular cross-section is 32.6mm. The maximum vertical displacement of lining for curve wall cross-section is 22.27mm. The maximum vertical displacement of lining for super-large cross-sections loess tunnel with different cross-sections ranges from small to large, namely, curve wall, circular, rectangular, horseshoe cross-section. This shows that the optimal cross-section of the super-large cross-section loess tunnel is the circular and the curve wall cross-section in the case of super shallow buried tunnel.
According to the boundary value of deep and shallow buried large cross-section loess tunnel, the overburden thickness of shallow buried tunnel is 30m. The process of analysis is similar to the static stability analysis for super shallow buried tunnel. The safety factor and critical strain nephogram in the case of shallow buried tunnel are shown in Table 6.
Table 6 shows that the safety factor of super-large cross-section loess tunnel with different cross-section forms can be obtained by strength reduction method. The safety factor of circular super-large cross-section is 1.94. The safety factor of horseshoe super-large cross-section is 1.34. The safety factor of rectangular super-large cross-section is 1.31. The safety factor of curve wall super-large cross-section is 1.58. This shows that the safety factors are ordered from large to small, namely, circular, curved wall, horseshoe and rectangle cross-section. From critical strain nephogram, the strain in the arch feet and the vault is more concentrated, and the range of cross-section strain increases with the increase of the safety factor. The critical first principal stress nephogram in the case of shallow buried loess tunnel with different cross-sections is shown in Table 7.
Table 7 shows that the maximum stress appears in the vault and arch feet, and the range of stress nephogram increases with the increase of the safety factor. This shows that the bearing capacity of cross-section increases with the increase in the safety factor. The stress distribution of circular and curve wall cross-section is more uniform, and the stresses of rectangle and horseshoe cross-section are more concentrated. Compared with the critical first principal stress of the super shallow buried loess tunnel with the overburden of 10m, the critical first principal stress of the shallow buried loess tunnel is increased. The maximum vertical displacement nephogram of lining in the case of shallow buried loess tunnel with different cross-sections is shown in Table 8.
Table 8 shows that the maximum vertical displacement occurs in the vault. The maximum vertical displacement of lining for circular cross-section is 44.39mm. The maximum vertical displacement of lining for horseshoe cross-section is 42.09mm. The maximum vertical displacement of lining for rectangular cross-section is 46.85mm. The maximum vertical displacement of lining for curve wall cross-section is 32.78mm. The maximum vertical displacement of lining for super-large cross-sections loess tunnel with different cross-sections starting from large to small, namely, curve wall, horseshoe, circular, and rectangular cross-section. Compared with the super shallow buried tunnel, the maximum vertical displacement of the lining for the shallow buried tunnel is increased. This shows that the optimal cross-section of the super-large cross-section loess tunnel is the circular and the curve wall cross-section in the case of shallow buried tunnel. Meanwhile, it shows that the stress of the super-large cross-section loess tunnel in the case of shallow buried is greater than the stress in the case of super shallow buried.
According to the boundary value of deep and shallow buried large cross-section loess tunnel, the overburden thickness of deep-buried tunnel is 80m. The process of analysis is similar to the static stability analysis for super shallow buried tunnel. The safety factor and critical strain nephogram in the case of deep-buried tunnel are shown in Table 9.
From Table 9, the safety factor of super-large cross-section loess tunnel can be obtained by strength reduction method in the case of deep buried tunnel with different cross-sections forms. The safety factor of circular super-large cross-section is 1.63. The safety factor of horseshoe super-large cross-section is 1.04. The rectangular super-large cross-section does not converge. The safety factor of curve wall super-large cross-section is 1.80. This shows that the safety factors are ordered from large to small, namely curved wall, circular, horseshoe cross-section. From critical strain nephogram, strain occurs in the arch feet and the vault, and the range of cross-section strain increases with the increase of the safety factor. This shows that the safety factor of the super-large cross-section loess tunnel with different cross-sections in the case of super shallow buried, shallow-buried and deep-buried are ordered from large to small, namely circular, curve wall, horseshoe, and rectangular cross-section. The critical first principal stress nephogram in the case of deep-buried loess tunnel with different cross-sections is shown in Table 10.
Table 11. Shows that the maximum vertical displacement occurs in vault in the case of deep-buried tunnel. The maximum vertical displacement of lining for circular cross-section is 78.10mm. The maximum vertical displacement of lining for horseshoe cross-section is 73.99mm. The maximum vertical displacement of lining for rectangular cross-section is 88.28mm. The maximum vertical displacement of lining for curve wall cross-section is 60.03mm. The maximum vertical displacement of lining for super-large cross-section loess tunnel with different cross-sections occurs from small to large, namely, curve wall, horseshoe, circular, rectangular cross-section. The maximum first principal stress of tunnel appears at the bottom of arch, and the stress distribution of the circular and the curve wall cross-section is more uniform and symmetrical, but the stress distribution of the rectangular and the horseshoe cross-section is not uniform. Compared with the super shallow buried and shallow-buried loess tunnel, the maximum vertical displacement of the lining for the deep buried tunnel occurs in vault and increases with the increase of overburden. This shows that the pressure and the maximum vertical displacement of the lining for loess tunnel increase with the increase of the depth of loess tunnel.
Table 10 shows that the maximum first principal stress of the super-large cross-section loess tunnel occurs at the arch feet. The stress distribution of circular and curve wall cross-section is more uniform in vault and arch feet, and the stress of horseshoe cross-section is more concentrated in arch bottom. The stress distribution of the rectangular cross-section is around the tunnel, and the maximum stress occurs at the arch feet. Compared with the critical first principal stress of the super shallow buried and shallow buried loess tunnel, the critical first principal stress of the deep buried loess tunnel is gradually increased. This shows that the force cross-section bore increases with the increase of the thickness of cover layer. The maximum vertical displacement nephogram of lining in the case of deep-buried loess tunnel with different cross-sections is shown in Table 11.
Not applicable.
The authors declare no conflict of interest, financial or otherwise.
This paper is a part of the national natural science foundation of China (Grant number: 51478212), and a part of the education ministry doctoral tutor foundation of China (Grant number: 20136201110003), and a part of science & technology development project of China railway 21 group Co., LTD (Grant number: 14B-3).
[1] | A. Bahmanikashkooli, M. Zare, and B. Safarpour, "Application of particle swarm optimization algorithm for computing critical depth of horseshoe cross section tunnel", APCBEE Procedia, vol. 9, pp. 207-211, 2014. [http://dx.doi.org/10.1016/j.apcbee.2014.01.037] |
[2] | Z.M. Xu, R.Q. Hunag, and S.T. Wang, "Tunnel classifying in light of depth", Chin. J. Geol. Hazard Control, vol. 11, no. 4, pp. 8-13, 2000. |
[3] | M.N. Wang, J. Guo, and L.S. Luo, "Study of critical buried depth of large cross-section loess tunnel for high speed railway", Yantu Lixue, vol. 31, no. 4, pp. 1157-1162, 2010. |
[4] | B. Mobaraki, and M. Vaghefi, "Numerical study of the depth and cross-sectional shape of tunnel under surface explosion", Tunn. Undergr. Space Technol., vol. 47, pp. 114-122, 2015. [http://dx.doi.org/10.1016/j.tust.2015.01.003] |
[5] | A.Z. Lu, H.Y. Chen, and Y. Qin, "Shape optimization of the support section of a tunnel at great depths", Comput. Geotech., vol. 61, pp. 190-197, 2014. [http://dx.doi.org/10.1016/j.compgeo.2014.05.011] |
[6] | A.Z. Lv, "Study of critical buried depth of large cross-section loess tunnel for high speed railway", J. China coal soc., vol. 22, no. 5, pp. 1157-1162, 1997. |
[7] | A.Z. Lv, H.Y. Chen, and Y. Tan, "Shape optimization of tunnel support section under contact condition of pure slip", Chin J. Jock Mech. Eng., vol. 33, no. 8, pp. 1563-1571, 2014. |
[8] | D. Peila, "A theoretical study of reinforcement influence on the stability of a tunnel face", Geotech. Geol. Eng., vol. 12, no. 3, pp. 145-168, 1994. [http://dx.doi.org/10.1007/BF00426984] |
[9] | X.S. Cheng, W. Shi, and R.R. Tian, "Safety assessment of lined loess tunnel under earthquake", Electron. J. Geotech. Eng., vol. 19, no. I, pp. 1903-1915, 2014. [J]. |
[10] | X.S. Cheng, R.R. Tian, and J.L. Wang, "Parameter determination about loess tunnel analysis model under earthquake action", The 11^{th} International Symposium on Structural Engineering., pp. 34-38, 2010. |
[11] | X.S. Cheng, H.Y. Zhu, and R.R. Tian, "Soil mass stability analysis of circular unlined loess tunnel under earthquake action", International Symposium on Geomechanics and Geotechnics. From Micro to Macro,, pp. 1127-1134, 2010. [http://dx.doi.org/10.1201/b10528-169] |
[12] | Y.R. Zheng, C.Y. Qiu, and H. Zhang, "Exploration of stability analysis methods for surrounding rocks of soil tunnel", Chin. J. Jock Mech. Eng., vol. 27, no. 10, pp. 1968-1980, 2008. |
[13] | Y.R. Zheng, H.L. Ye, and Q. Xiao, "The analysis about the stability of seismic slope and tunnel based on full dynamic analysis method", J. Disaster Prev Mitigation Eng., vol. 30, no. S1, pp. 279-285, 2010. |
[14] | X.S. Cheng, R.R. Tian, and J.L. Wang, "Static and seismic stability analysis of soil mass surrounding tunnel about large-span unlined egg-shaped loess tunnel", Tumu Gongcheng Xuebao, vol. 43, no. S1, pp. 582-587, 2010. |
[15] | H. Li, X.Y. Zhu, and Y.F. Xu, "Influence of tunnel cross-section shape on stress and displacement of surrounding rock", Tunnel Constr, vol. 29, no. 1, pp. 38-44, 2009. |
[16] | S.M. Dong, Q.C. Xin, and S.S. Lu, "The influence of the tunnel cross-section shape on stability of surrounding rocks", China Rural Water Hydropower, vol. 01, no. 102, pp. 104-107, 2011. |
[17] | O.C. Zienkiewicz, C. Humpheson, and R.W. Lewis, "Associated and non-associated visco-plasticity and plasticity in soil mechanics", Geotechnique, vol. 25, no. 4, pp. 671-689, 1975. [http://dx.doi.org/10.1680/geot.1975.25.4.671] |
[18] | Y.R. Zheng, Z.Y. Chen, and G.X. Wang, Eng. Treat. Slope Landslide., China Communications Press, 2011. |
[19] | Y.R. Zheng, S.Y. Zhao, and C.J. Deng, "Development of finite element limit analysis method and its applications in geotechnical engineering", Eng. Sci., vol. 8, no. 12, pp. 39-61, 2006. |
[20] | Y.R. Zheng, S.Y. Zhao, and L.Y. Zhang, "Slope stability analysis by strength reduction FEM", Eng. Sci., vol. 4, no. 10, pp. 57-62, 2002. |
[21] | X.S. Cheng, and Y.R. Zheng, "Calculation discussion about safety factor of unlined loess tunnel wall rock structure under earthquake", Yantu Lixue, vol. 32, no. 3, pp. 761-766, 2011. |