scholarly journals Probing afterslip and viscoelastic relaxation following the 2015 Mw 7.8 Gorkha earthquake based on the 3-D Finite Element Model

2020 ◽  
Author(s):  
Lina Su ◽  
Fuqiang Shi ◽  
Weijun Gan ◽  
Xiaoning Su ◽  
Junyi Yan

Abstract We analyzed daily displacement time series from 34 continuous GPS stations in Nepal and 5 continuous GPS stations in South Tibet, China, and extracted the first 4.8 years of postseismic motion after the 2015 Mw7.8 Gorkha earthquake. With the longer duration GPS observations, we find that postseismic displacements mainly exhibit southward and uplift motion. To study the postseismic afterslip and viscoelastic relaxation, we built a 3-D spherical finite element model (FEM) with heterogeneous material properties and surface topography across the Himalayan range, accounting for the strong variations in material properties and surface elevation along central Himalayan arc. On the basis of the FEM, we reveal that the predicted viscoelastic relaxation of cm level moves southward to the north of the Gorkha earthquake rupture, but in an opposite direction to the observed postseismic deformation in the south; the postseismic deformation with viscoelastic relaxation is well explained by afterslip downdip of the coseismic rupture, which is still dominant 4.8 years after the 2015 Mw7.8 Gorkha earthquake. The lack of slip on a shallow portion and western segment of the MHT during and after the 2015 Gorkha earthquake implies continued seismic hazard in the future.

2020 ◽  
Author(s):  
Lina Su ◽  
Fuqiang Shi ◽  
Weijun Gan ◽  
Xiaoning Su ◽  
Junyi Yan

Abstract We analyzed daily displacement time series from 34 continuous GPS stations in Nepal and 5 continuous GPS stations in South Tibet, China, and extracted the first 4.8 years of postseismic motion after the 2015 Mw7.8 Gorkha earthquake. With the longer duration GPS observations, we find that postseismic displacements mainly exhibit southward and uplift motion. To study the postseismic afterslip and viscoelastic relaxation, we built a 3-D spherical finite element model (FEM) with heterogeneous material properties and surface topography across the Himalayan range, accounting for the strong variations in material properties and surface elevation along central Himalayan arc. On the basis of the FEM, we reveal that the predicted viscoelastic relaxation of cm level moves southward to the north of the Gorkha earthquake rupture, but in an opposite direction to the observed postseismic deformation in the south; the postseismic deformation excluding viscoelastic relaxation is well explained by afterslip downdip of the coseismic rupture. The afterslip is dominant during 4.8 years after the 2015 Mw7.8 Gorkha earthquake; the contribution by the viscoelastic relaxation gradually increases slightly. The lack of slip on a shallow portion and western segment of the MHT during and after the 2015 Gorkha earthquake implies continued seismic hazard in the future.


2020 ◽  
Vol 72 (1) ◽  
Author(s):  
Lina Su ◽  
Fuqiang Shi ◽  
Weijun Gan ◽  
Xiaoning Su ◽  
Junyi Yan

Abstract We analyzed daily displacement time series from 34 continuous GPS stations in Nepal and 5 continuous GPS stations in South Tibet, China, and extracted the first 4.8 years of postseismic motion after the 2015 Mw7.8 Gorkha earthquake. With the longer duration GPS observations, we find that postseismic displacements mainly exhibit southward and uplift motion. To study the postseismic afterslip and viscoelastic relaxation, we built a 3-D spherical finite-element model (FEM) with heterogeneous material properties and surface topography across the Himalayan range, accounting for the strong variations in material properties and surface elevation along the central Himalayan arc. On the basis of the FEM, we reveal that the predicted viscoelastic relaxation of cm level moves southward to the north of the Gorkha earthquake rupture, but in an opposite direction to the observed postseismic deformation in the south; the postseismic deformation excluding viscoelastic relaxation is well explained by afterslip downdip of the coseismic rupture. The afterslip is dominant during 4.8 years after the 2015 Mw7.8 Gorkha earthquake; the contribution by the viscoelastic relaxation gradually increases slightly. The lack of slip on a shallow portion and western segment of the MHT during and after the 2015 Gorkha earthquake implies continued seismic hazard in the future.


2020 ◽  
Author(s):  
Lina Su ◽  
Fuqiang Shi ◽  
Weijun Gan ◽  
Xiaoning Su ◽  
Junyi Yan

Abstract We analyzed GPS coordinate time series from 34 continuous GPS stations in Nepal and 5 continuous GPS stations in South Tibet of China, and extracted the first 4.8 years postseismic displacements after the 2015 Mw7.8 Gorkha earthquake. With the longer duration GPS observations, we found that postseismic displacements mainly exhibit the southward and uplift movement at the epcientral area. To study the postseismic afterslip and viscoelastic relaxation, we then built 3-D spherical finite element model (FEM) with heterogeneous material properties and surface topography across the Himalayan range, accounting for the strong variations of material properties and surface elevation along central Himalayan arc. The sophisticated FEM is more realistic and perform better than the flat layered earth model. On the basis of it, we reveal that the predicted viscoelastic relaxation of cm level is opposite to the observed postseismic deformation; the postseismic deformation with viscoelastic relaxation deducted is well explained by afterslip downdip of coseismic rupture, which indicates the afterslip is still dominant during 4.8 years postseismic deformation after the 2015 Mw7.8 Gorkha earthquake; The lack of slip on a shallow portion and western segment of the MHT during and after the 2015 Gorkha earthquake implies continued seismic hazard in the future.


2011 ◽  
Vol 9 (67) ◽  
pp. 351-361 ◽  
Author(s):  
O. Panagiotopoulou ◽  
S. D. Wilshin ◽  
E. J. Rayfield ◽  
S. J. Shefelbine ◽  
J. R. Hutchinson

Finite element modelling is well entrenched in comparative vertebrate biomechanics as a tool to assess the mechanical design of skeletal structures and to better comprehend the complex interaction of their form–function relationships. But what makes a reliable subject-specific finite element model? To approach this question, we here present a set of convergence and sensitivity analyses and a validation study as an example, for finite element analysis (FEA) in general, of ways to ensure a reliable model. We detail how choices of element size, type and material properties in FEA influence the results of simulations. We also present an empirical model for estimating heterogeneous material properties throughout an elephant femur (but of broad applicability to FEA). We then use an ex vivo experimental validation test of a cadaveric femur to check our FEA results and find that the heterogeneous model matches the experimental results extremely well, and far better than the homogeneous model. We emphasize how considering heterogeneous material properties in FEA may be critical, so this should become standard practice in comparative FEA studies along with convergence analyses, consideration of element size, type and experimental validation. These steps may be required to obtain accurate models and derive reliable conclusions from them.


Author(s):  
Sean M. Finley ◽  
J. Harley Astin ◽  
Evan Joyce ◽  
Andrew T. Dailey ◽  
Douglas L. Brockmeyer ◽  
...  

OBJECTIVE The underlying biomechanical differences between the pediatric and adult cervical spine are incompletely understood. Computational spine modeling can address that knowledge gap. Using a computational method known as finite element modeling, the authors describe the creation and evaluation of a complete pediatric cervical spine model. METHODS Using a thin-slice CT scan of the cervical spine from a 5-year-old boy, a 3D model was created for finite element analysis. The material properties and boundary and loading conditions were created and model analysis performed using open-source software. Because the precise material properties of the pediatric cervical spine are not known, a published parametric approach of scaling adult properties by 50%, 25%, and 10% was used. Each scaled finite element model (FEM) underwent two types of simulations for pediatric cadaver testing (axial tension and cardinal ranges of motion [ROMs]) to assess axial stiffness, ROM, and facet joint force (FJF). The authors evaluated the axial stiffness and flexion-extension ROM predicted by the model using previously published experimental measurements obtained from pediatric cadaveric tissues. RESULTS In the axial tension simulation, the model with 50% adult ligamentous and annulus material properties predicted an axial stiffness of 49 N/mm, which corresponded with previously published data from similarly aged cadavers (46.1 ± 9.6 N/mm). In the flexion-extension simulation, the same 50% model predicted an ROM that was within the range of the similarly aged cohort of cadavers. The subaxial FJFs predicted by the model in extension, lateral bending, and axial rotation were in the range of 1–4 N and, as expected, tended to increase as the ligament and disc material properties decreased. CONCLUSIONS A pediatric cervical spine FEM was created that accurately predicts axial tension and flexion-extension ROM when ligamentous and annulus material properties are reduced to 50% of published adult properties. This model shows promise for use in surgical simulation procedures and as a normal comparison for disease-specific FEMs.


Author(s):  
X. G. Tan ◽  
R. Kannan ◽  
Andrzej J. Przekwas

Until today the modeling of human body biomechanics poses many great challenges because of the complex geometry and the substantial heterogeneity of human body. We developed a detailed human body finite element model in which the human body is represented realistically in both the geometry and the material properties. The model includes the detailed head (face, skull, brain, and spinal cord), the skeleton, and air cavities (including the lung). Hence it can be used to accurately acquire the stress wave propagation in the human body under various loading conditions. The blast loading on the human surface was generated from the simulated C4 blast explosions, via a novel combination of 1-D and 3-D numerical formulations. We used the explicit finite element solver in the multi-physics code CoBi for the human body biomechanics. This is capable of solving the resulting large system containing millions of unknowns in an extremely scalable fashion. The meshes generated for these simulations are of good quality. This enables us to employ relatively large time step sizes, without resorting to the artificial time scaling treatment. In order to study the human body dynamic response under the blast loading, we also developed an interface to apply the blast pressure loading on the external human body surface. These newly developed models were used to conduct parametric simulations to find out the brain biomechanical response when the blasts impact the human body. Under the same blast loading we also show the differences of brain response when having different material properties for the skeleton, the existence of other body parts such as torso.


2017 ◽  
Vol 89 (2) ◽  
pp. 274-279
Author(s):  
Thomas Wright ◽  
Imran Hyder ◽  
Mitchell Daniels ◽  
David Kim ◽  
John P. Parmigiani

Purpose The purpose of this paper is to determine which of the ten material properties of the Hashin progressive damage model significantly affect the maximum load-carrying ability of center-notched carbon fiber panels under in-plane tension and out-of-plane bending. Design/methodology/approach The approach used is to calculate the maximum load using a finite element model for a range of material property values as specified by a fraction factorial design. The finite element model used has been experimentally validated in prior work. Findings Results showed that for the laminates considered, at most three and as few as one of the ten Hashin material properties significantly affected the magnitude of the maximum load. Practical implications While the results of this paper only specifically apply to the laminates included in the study, the results suggest that, in general, only a small number of the Hashin material properties affect laminate load-carrying ability. Originality/value Knowing which properties are significant is of value in selecting materials to optimize performance and also in determining which properties need to be known to a high accuracy.


2016 ◽  
Vol 17 (4) ◽  
pp. 497-503 ◽  
Author(s):  
Rinchen Phuntsok ◽  
Marcus D. Mazur ◽  
Benjamin J. Ellis ◽  
Vijay M. Ravindra ◽  
Douglas L. Brockmeyer

OBJECT There is a significant deficiency in understanding the biomechanics of the pediatric craniocervical junction (CCJ) (occiput–C2), primarily because of a lack of human pediatric cadaveric tissue and the relatively small number of treated patients. To overcome this deficiency, a finite element model (FEM) of the pediatric CCJ was created using pediatric geometry and parameterized adult material properties. The model was evaluated under the physiological range of motion (ROM) for flexion-extension, axial rotation, and lateral bending and under tensile loading. METHODS This research utilizes the FEM method, which is a numerical solution technique for discretizing and analyzing systems. The FEM method has been widely used in the field of biomechanics. A CT scan of a 13-month-old female patient was used to create the 3D geometry and surfaces of the FEM model, and an open-source FEM software suite was used to apply the material properties and boundary and loading conditions and analyze the model. The published adult ligament properties were reduced to 50%, 25%, and 10% of the original stiffness in various iterations of the model, and the resulting ROMs for flexion-extension, axial rotation, and lateral bending were compared. The flexion-extension ROMs and tensile stiffness that were predicted by the model were evaluated using previously published experimental measurements from pediatric cadaveric tissues. RESULTS The model predicted a ROM within 1 standard deviation of the published pediatric ROM data for flexion-extension at 10% of adult ligament stiffness. The model's response in terms of axial tension also coincided well with published experimental tension characterization data. The model behaved relatively stiffer in extension than in flexion. The axial rotation and lateral bending results showed symmetric ROM, but there are currently no published pediatric experimental data available for comparison. The model predicts a relatively stiffer ROM in both axial rotation and lateral bending in comparison with flexion-extension. As expected, the flexion-extension, axial rotation, and lateral bending ROMs increased with the decrease in ligament stiffness. CONCLUSIONS An FEM of the pediatric CCJ was created that accurately predicts flexion-extension ROM and axial force displacement of occiput–C2 when the ligament material properties are reduced to 10% of the published adult ligament properties. This model gives a reasonable prediction of pediatric cervical spine ligament stiffness, the relationship between flexion-extension ROM, and ligament stiffness at the CCJ. The creation of this model using open-source software means that other researchers will be able to use the model as a starting point for research.


Sign in / Sign up

Export Citation Format

Share Document