The occurrence of earthquakes is closely related to the crustal geotectonic movement and the migration of mass,which consequently cause changes in gravity.The Gravity Recovery And Climate Experiment(GRACE)satellite da...The occurrence of earthquakes is closely related to the crustal geotectonic movement and the migration of mass,which consequently cause changes in gravity.The Gravity Recovery And Climate Experiment(GRACE)satellite data can be used to detect gravity changes associated with large earthquakes.However,previous GRACE satellite-based seismic gravity-change studies have focused more on coseismic gravity changes than on preseismic gravity changes.Moreover,the noise of the north–south stripe in GRACE data is difficult to eliminate,thereby resulting in the loss of some gravity information related to tectonic activities.To explore the preseismic gravity anomalies in a more refined way,we first propose a method of characterizing gravity variation based on the maximum shear strain of gravity,inspired by the concept of crustal strain.The offset index method is then adopted to describe the gravity anomalies,and the spatial and temporal characteristics of gravity anomalies before earthquakes are analyzed at the scales of the fault zone and plate,respectively.In this work,experiments are carried out on the Tibetan Plateau and its surrounding areas,and the following findings are obtained:First,from the observation scale of the fault zone,we detect the occurrence of large-area gravity anomalies near the epicenter,oftentimes about half a year before an earthquake,and these anomalies were distributed along the fault zone.Second,from the observation scale of the plate,we find that when an earthquake occurred on the Tibetan Plateau,a large number of gravity anomalies also occurred at the boundary of the Tibetan Plateau and the Indian Plate.Moreover,the aforementioned experiments confirm that the proposed method can successfully capture the preseismic gravity anomalies of large earthquakes with a magnitude of less than 8,which suggests a new idea for the application of gravity satellite data to earthquake research.展开更多
Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict te...Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.展开更多
The maximum principal stress, von Mises equivalent stress, equivalent creep strain, stress triaxiality in dissimilar metal welded joints between austenitic（HR3C） and martensitic heat-resistant steel（T91） are simul...The maximum principal stress, von Mises equivalent stress, equivalent creep strain, stress triaxiality in dissimilar metal welded joints between austenitic（HR3C） and martensitic heat-resistant steel（T91） are simulated by FEM at 873 K and under inner pressure of 42.26 MPa. The results show that the maximum principal stress and von Mises equivalent stress are quite high in the vicinity of weld/T91 interface, creep cavities are easy to form and expand in the weld/T91 interface. There are two peaks of equivalent creep strains in welded joint, and the maximum equivalent creep strain is in the place 27-32 mm away from the weld/T91 interface, and there exists creep constrain region in the vicinity of weld/T91 interface. The high stress triaxiality peak is located exactly at the weld/T91 interface. Accordingly, the weld/T91 interface is the weakest site of welded joint. Therefore, using stress triaxiality to describe creep cavity nucleation and expansion and crack development is reasonable for the dissimilar metal welded joint between austenitic and martensitic steel.展开更多
Critical plane method is one of the most promising approaches to predict the fatigue life when the structure is subjected to the multiaxial loading.The stress-strain status and the critical plane position for smooth s...Critical plane method is one of the most promising approaches to predict the fatigue life when the structure is subjected to the multiaxial loading.The stress-strain status and the critical plane position for smooth specimens are calculated using theoretical approaches when the loading mode is a continuous function.However,because of the existence of stress concentration and the characteristic of multiaxial non-proportion,it is difficult to calculate the stress-strain status and the critical plane position of geometric discontinuity structure by theory method.In this paper,a new numerical simulation method is proposed to determine the critical plane of geometric discontinuity structure under multiaxial loading.Firstly,the strain status of dangerous point is analyzed by finite element method.Secondly,the maximum shear strain amplitude of arbitrary plane is calculated using coordinate transformation principle.Finally,the plane which has the maximum shear strain amplitude is defined as the critical plane.The critical plane positions are analyzed when loading mode and notch parameters are different.Meanwhile,the relationship between notch depth and associated parameters on critical plane as well as that between loading amplitude and associated parameters on critical plane are given quantitatively.展开更多
Based on the multiple-term horizontal velocity solutions of 230 GPS monitoring sites in Tianshan and its adjacent region, the GPS site velocity fields and crustal horizontal strain fields in the area have been obtaine...Based on the multiple-term horizontal velocity solutions of 230 GPS monitoring sites in Tianshan and its adjacent region, the GPS site velocity fields and crustal horizontal strain fields in the area have been obtained. The results show that the crustal shortening rate of Tianshan, with the longitude （77°±1°）E as the boundary, gradually decreased towards two sides, from the south to the north, indicating that the pushing force of plate becomes weaker along with the fold deformation decreasing of the Tianshan. The direction of principal compressive strain of Tianshan and its adjacent area, nearly NNW, is basically perpendicular to the Tianshan cordillera trend, suggesting the distribution and variation of maximum principal compressive stress in Tianshan and its adjacent region resulted from collision and extrusion of Indian Plate. This paper indicates that the maximum shear strain field mainly con- centrates on two areas, one is Isyk lake of North Tianshan, Kyrgyzstan, and the other is the juncture of Jiashi （South Tianshan） and Pamir arc faults. In the above areas, it can be shown from the epicentral distribution that the strong earthquakes mostly occurs at the high shearing strain accumulation filed or its edge.展开更多
On January 21, 2016, a strong earthquake with a magnitude of Ms6.4 happened at Menyuan, Qinghai Province of China. In almost the same place, there was another strong earthquake happened in 1986, with similar magnitude...On January 21, 2016, a strong earthquake with a magnitude of Ms6.4 happened at Menyuan, Qinghai Province of China. In almost the same place, there was another strong earthquake happened in 1986, with similar magnitude and focal mechanism. In this paper, we analyze the characteristics of regional crustal deformation before the 2016 Menyuan Ms6.4 earth- quake by using the data from 10 continuous Global Positioning System （GPS） stations and 74 campaign-mode GPS stations within 200 km of this event： （a） Based on the velocity field from over ten years GPS observations, a regional strain rate field is calculated. The results indicate that the crustal strain rate and seismic moment accumulation rate of the Qilian- Haiyuan active fault, which is the seismogenic tectonics of the event, are significantly higher than the surrounding regions. In a 20 km~ 20 km area around the seismogenic region, the maximum and minimum principal strain rates are 21.5 nanostrain/a （NW-SE extension） and -46.6 nanostrain/a （NE-SW compression）, respectively, and the seismic moment accumulation rates is 17.4 Nm/a. The direction of principal compression is consistent with the focal mechanism of this event. （b） Based on the position time series of the continuous GPS stations for a time-span of about 6 years before the event, we calculate the strain time series. The results show that the dilatation of the seismogenic region is continuously reduced with a ＂non-linear＂ trend since 2010, which means the seismogenic region has been in a state of compression. However, about 2-3 months before the event, both the dilatation and maximum shear strain show significant inverse trends. These abnormal changes of crustal deformation may reflect the non-linear adjustment of the stress-strain accumulation of the seismogenic region, when the accumulation is approaching the critical value of rupture.展开更多
基金supported by the National Key Research and Development Program of China(Grant No.2019YFC1509202)the National Natural Science Foundation of China(Grant Nos.41772350,61371189,and 41701513).
文摘The occurrence of earthquakes is closely related to the crustal geotectonic movement and the migration of mass,which consequently cause changes in gravity.The Gravity Recovery And Climate Experiment(GRACE)satellite data can be used to detect gravity changes associated with large earthquakes.However,previous GRACE satellite-based seismic gravity-change studies have focused more on coseismic gravity changes than on preseismic gravity changes.Moreover,the noise of the north–south stripe in GRACE data is difficult to eliminate,thereby resulting in the loss of some gravity information related to tectonic activities.To explore the preseismic gravity anomalies in a more refined way,we first propose a method of characterizing gravity variation based on the maximum shear strain of gravity,inspired by the concept of crustal strain.The offset index method is then adopted to describe the gravity anomalies,and the spatial and temporal characteristics of gravity anomalies before earthquakes are analyzed at the scales of the fault zone and plate,respectively.In this work,experiments are carried out on the Tibetan Plateau and its surrounding areas,and the following findings are obtained:First,from the observation scale of the fault zone,we detect the occurrence of large-area gravity anomalies near the epicenter,oftentimes about half a year before an earthquake,and these anomalies were distributed along the fault zone.Second,from the observation scale of the plate,we find that when an earthquake occurred on the Tibetan Plateau,a large number of gravity anomalies also occurred at the boundary of the Tibetan Plateau and the Indian Plate.Moreover,the aforementioned experiments confirm that the proposed method can successfully capture the preseismic gravity anomalies of large earthquakes with a magnitude of less than 8,which suggests a new idea for the application of gravity satellite data to earthquake research.
基金he authors thank the financial support fromthe Key Program of National Natural Science Foundation of China(GrantNo.52039007)the Youth Science and Technology Innovation Research Team Fund of Sichuan Province(Grant No.2020JDTD0001).
文摘Evaluating the fracture resistance of rocks is essential for predicting and preventing catastrophic failure of cracked structures in rock engineering.This investigation developed a brittle fracture model to predict tensile mode(mode I)failure loads of cracked rocks.The basic principle of the model is to estimate the reference crack corresponding to the fracture process zone(FPZ)based on the maximum normal strain(MNSN)ahead of the crack tip,and then use the effective crack to calculate the fracture toughness.We emphasize that the non-singular stress/strain terms should be considered in the description of the MNSN.In this way,the FPZ,non-singular terms and the biaxial stress state at the crack tip are simul-taneously considered.The principle of the model is explicit and easy to apply.To verify the proposed model,laboratory experiments were performed on a rock material using six groups of specimens.The model predicted the specimen geometry dependence of the measured fracture toughness well.More-over,the potential of the model in analyzing the size effect of apparent fracture toughness was discussed and validated through experimental data reported in the literature.The model was demonstrated su-perior to some commonly used fracture models and is an excellent tool for the safety assessment of cracked rock structures.
基金Funded by the National Natural Science Foundation of China(No.51374154)
文摘The maximum principal stress, von Mises equivalent stress, equivalent creep strain, stress triaxiality in dissimilar metal welded joints between austenitic（HR3C） and martensitic heat-resistant steel（T91） are simulated by FEM at 873 K and under inner pressure of 42.26 MPa. The results show that the maximum principal stress and von Mises equivalent stress are quite high in the vicinity of weld/T91 interface, creep cavities are easy to form and expand in the weld/T91 interface. There are two peaks of equivalent creep strains in welded joint, and the maximum equivalent creep strain is in the place 27-32 mm away from the weld/T91 interface, and there exists creep constrain region in the vicinity of weld/T91 interface. The high stress triaxiality peak is located exactly at the weld/T91 interface. Accordingly, the weld/T91 interface is the weakest site of welded joint. Therefore, using stress triaxiality to describe creep cavity nucleation and expansion and crack development is reasonable for the dissimilar metal welded joint between austenitic and martensitic steel.
基金the National Natural Science Foundation of China(Nos.51778273 and 51605212)the Universities and Colleges Innovation Ability Improvement Project of Gansu(No.2019A-225)
文摘Critical plane method is one of the most promising approaches to predict the fatigue life when the structure is subjected to the multiaxial loading.The stress-strain status and the critical plane position for smooth specimens are calculated using theoretical approaches when the loading mode is a continuous function.However,because of the existence of stress concentration and the characteristic of multiaxial non-proportion,it is difficult to calculate the stress-strain status and the critical plane position of geometric discontinuity structure by theory method.In this paper,a new numerical simulation method is proposed to determine the critical plane of geometric discontinuity structure under multiaxial loading.Firstly,the strain status of dangerous point is analyzed by finite element method.Secondly,the maximum shear strain amplitude of arbitrary plane is calculated using coordinate transformation principle.Finally,the plane which has the maximum shear strain amplitude is defined as the critical plane.The critical plane positions are analyzed when loading mode and notch parameters are different.Meanwhile,the relationship between notch depth and associated parameters on critical plane as well as that between loading amplitude and associated parameters on critical plane are given quantitatively.
基金National Natural Science Foundation of China (40074024) and Natural Science Foundation of Xinjiang Uygur Autonomous Region (200321101).
文摘Based on the multiple-term horizontal velocity solutions of 230 GPS monitoring sites in Tianshan and its adjacent region, the GPS site velocity fields and crustal horizontal strain fields in the area have been obtained. The results show that the crustal shortening rate of Tianshan, with the longitude （77°±1°）E as the boundary, gradually decreased towards two sides, from the south to the north, indicating that the pushing force of plate becomes weaker along with the fold deformation decreasing of the Tianshan. The direction of principal compressive strain of Tianshan and its adjacent area, nearly NNW, is basically perpendicular to the Tianshan cordillera trend, suggesting the distribution and variation of maximum principal compressive stress in Tianshan and its adjacent region resulted from collision and extrusion of Indian Plate. This paper indicates that the maximum shear strain field mainly con- centrates on two areas, one is Isyk lake of North Tianshan, Kyrgyzstan, and the other is the juncture of Jiashi （South Tianshan） and Pamir arc faults. In the above areas, it can be shown from the epicentral distribution that the strong earthquakes mostly occurs at the high shearing strain accumulation filed or its edge.
基金supported by the National Science Foundation of China(41474090)Science for Earthquake Resilience(XH14063)the State Key Laboratory of Earthquake Dynamics(LED2013A02)
文摘On January 21, 2016, a strong earthquake with a magnitude of Ms6.4 happened at Menyuan, Qinghai Province of China. In almost the same place, there was another strong earthquake happened in 1986, with similar magnitude and focal mechanism. In this paper, we analyze the characteristics of regional crustal deformation before the 2016 Menyuan Ms6.4 earth- quake by using the data from 10 continuous Global Positioning System （GPS） stations and 74 campaign-mode GPS stations within 200 km of this event： （a） Based on the velocity field from over ten years GPS observations, a regional strain rate field is calculated. The results indicate that the crustal strain rate and seismic moment accumulation rate of the Qilian- Haiyuan active fault, which is the seismogenic tectonics of the event, are significantly higher than the surrounding regions. In a 20 km~ 20 km area around the seismogenic region, the maximum and minimum principal strain rates are 21.5 nanostrain/a （NW-SE extension） and -46.6 nanostrain/a （NE-SW compression）, respectively, and the seismic moment accumulation rates is 17.4 Nm/a. The direction of principal compression is consistent with the focal mechanism of this event. （b） Based on the position time series of the continuous GPS stations for a time-span of about 6 years before the event, we calculate the strain time series. The results show that the dilatation of the seismogenic region is continuously reduced with a ＂non-linear＂ trend since 2010, which means the seismogenic region has been in a state of compression. However, about 2-3 months before the event, both the dilatation and maximum shear strain show significant inverse trends. These abnormal changes of crustal deformation may reflect the non-linear adjustment of the stress-strain accumulation of the seismogenic region, when the accumulation is approaching the critical value of rupture.
