A numerical method of viscoelastic finite element coupled with spring-block model to simulate the interaction of faults and crustal blocks

A numerical method of viscoelastic finite element coupled with spring-block model is developed to study temporal processes from the slow tectonic motion of large-scale crust to the rapid failure of small-scale faults. Our modeling demonstrates that the motion of crustal blocks is driven by forces from tectonic plate boundaries, and the deformation is distributed on faults for the stress accumulating. The coupling model generates earthquake sequences that display a magnitude-frequency scaling consistent with Gutenberg-Richter law. The frictional heterogeneities affeci earthquakes occurrence and stresses distribution of crustal blocks. Rupture of earthquakes starts at the nucleation node, and propagates bilaterally along faults with the stress triggering, release and redistribution. The failure of faults releases part of crustal stresses, the stress state of crustal blocks near fault is affected by the rupture of local segments on the fault, and the stress state of crustal blocks far away from the fault is controlled by the seismic activity of the whole fault.     

Numerical simulations of deformation and movement of blocks within North China in response to 1976 Tangshan earthquake

Using methods of discontinuous deformation analysis and finite element (DDA+FEM), this paper simulates dynamic processes of the Tangshan earthquake of 1976, which occurred in the northern North China where its internal blocks apparently interacted. Studies focus upon both the movement and deformation of the blocks, in particular, the Ordos block, and variations of stress states on the boundary faults. The Tangshan earthquake was composed of three events: slipping motions of NNE-striking major fault, NE-striking fault near the northeastern end of the NNE-striking fault, and NW-striking fault on the southeastern side of the NNE-striking fault. Compared with previous studies, our model yields a result that is more agreeable with the configuration of aftershock distributions. A number of data are presented, such as the principle stress field during the earthquake, contours of the maximum shear stress, the strike-slip deformation between blocks near the earthquake focus, time-dependent variations of slips of earthquake-triggered faulting, the maximum slip distance, and stress drops. These results are in accord with the earthquake source mechanism, basic parameters from earthquake wave study, macro-isoseismic line, observed horizontal displacement vectors, etc. The Tangshan earthquake exerted different influences on the adjacent blocks and boundary faults between them, thus resulting in differential movement and deformation. The Ordos block seems to have experienced the small-scale counterclockwise rotation and deformation, but its northeast part, bounded on the east by the Taihangshan and on the north by the Yanshan and Yinshan belts, underwent relatively stronger deformation. The Tangshan earthquake also changed the stress state of boundary faults of the North China, leading to an increase in shear stress and a decrease in normal stress in the NW-trending Zhangjiakou-Penglai fault through Tangshan City and the northern border faults of the Ordos block, and therefore raises the potential risk of earthquake occurrence. This result is supported by the facts that a series of Ms ≥ 6 earthquakes took place at the northern margin of the Ordos block after the Tangshan earthquake.     

Active blocks and strong seismic activity in North China region

The active North China block consists of three second-order blocks: Ordos, North China Plain, and East Shandong-Huanghai Sea blocks. Two active tectonic zones, the Anyang-Heze-Linyi and Tangshan-Cixian zones, exist in the active North China Plain block and have separated the active block into 3 third-order active blocks, Taihangshan, Hebei-Shandong, and Henan-Huai blocks. The 3 third-order active blocks are characterized by their entire motion and are clearly different in their Cenozoic structures and deep structures. The active boundary tectonic zones between the third-order active blocks are less than those between the first- and second-order active blocks in their movement strength, extent, and seismic activity. The density of M· ·6 earthquakes in the boundary zones between active blocks is higher than that within the blocks by 9–22 times in the North China region, up to one order of magnitude on average. M· · 7 earthquakes occurred basically in the boundary zones between active blocks. The difference is not occasional, but reflects the nature of intraplate movement and the characteristics of strong seismic activity and is the powerful evidence for hypothesis of active blocks.     

Fine structure of Pn velocity beneath Sichuan-Yunnan region

We use 23298 Pn arrival-time data from Chinese national and provincial earthquake bulletins to invert fine structure of Pn velocity and anisotropy at the top of the mantle beneath the Sichuan-Yunnan and its adjacent region. The results suggest that the Pn velocity in this region shows significant lateral variation; the Pn velocity varies from 7.7 to 8.3 km/s. The Pn-velocity variation correlates well with the tectonic activity and heat flow of the region. Low Pn velocity is observed in southwest Yunnan, Tengchong volcano area, and the Panxi tectonic area. These areas have very active seismicity and tectonic activity with high surface heat flow. On the other hand, high Pn velocity is observed in some stable regions, such as the central region of the Yangtze Platform; the most pronounced high velocity area is located in the Sichuan Basin, south of Chengdu. Pn anisotropy shows a complex pattern of regional deformation. The Pn fast direction shows a prominent clockwise rotation pattern from east of the Tibetan block to the Sichuan-Yunnan diamond block to southwest Yunnan, which may be related to southeastward escape of the Tibetan Plateau material due to the collision of the Indian Plate to the Eurasia Plate. Thus there appears to be strong correlation between the crustal deformation and the upper mantle structure in the region. The delay times of events and stations show that the crust thickness decreases from the Tibetan Plateau to eastern China, which is consistent with the results from deep seismic sounding.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 2: Complex modes

A method, based on the Hilbert–Huang spectral analysis, has been proposed by the authors to identify linear structures in which normal modes exist (i.e., real eigenvalues and eigenvectors). Frequently, all the eigenvalues and eigenvectors of linear structures are complex. In this paper, the method is extended further to identify general linear structures with complex modes using the free vibration response data polluted by noise. Measured response signals are first decomposed into modal responses using the method of Empirical Mode Decomposition with intermittency criteria. Each modal response contains the contribution of a complex conjugate pair of modes with a unique frequency and a damping ratio. Then, each modal response is decomposed in the frequency–time domain to yield instantaneous phase angle and amplitude using the Hilbert transform. Based on a single measurement of the impulse response time history at one appropriate location, the complex eigenvalues of the linear structure can be identified using a simple analysis procedure. When the response time histories are measured at all locations, the proposed methodology is capable of identifying the complex mode shapes as well as the mass, damping and stiffness matrices of the structure. The effectiveness and accuracy of the method presented are illustrated through numerical simulations. It is demonstrated that dynamic characteristics of linear structures with complex modes can be identified effectively using the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

Predictive instantaneous optimal control of elastic structures during earthquakes

A predictive instantaneous optimal control (PIOC) algorithm is proposed for controlling the seismic responses of elastic structures. This algorithm compensates for the time delay that happens in practical control applications by predicting the structural response over a period that equals the time delay, and by substituting the predicted response in the instantaneous optimal control (IOC) algorithm. The unique feature of this proposed PIOC algorithm is that it is simple and at the same time compensates for the time delay very effectively. Numerical examples of single degree of freedom structures are presented to compare the performance of PIOC and IOC systems for various time delay magnitudes. Results show that a time delay always causes degradation of control efficiency, but PIOC can greatly reduce this degradation compared to IOC. The effects of the structure's natural periods and the choice of control gains on the degradation induced by the time delay are also analyzed. Results show that shorter natural periods and larger control gains are both more sensitive and more serious to the degradation of control efficiency. Finally, a practical application of PIOC is performed on a six‐story moment‐resisting steel frame. It is demonstrated that PIOC contributes significantly to maintain stability in multiple degree of freedom structures, and at the same time PIOC has a satisfactory control performance. Copyright © 2003 John Wiley & Sons, Ltd.     

Predictive instantaneous optimal control of inelastic structures during earthquakes

A new inelastic structural control algorithm is proposed by incorporating the force analogy method (FAM) with the predictive instantaneous optimal control (PIOC) algorithm. While PIOC is very effective in compensating for the time delay for elastic structures, the FAM is highly efficient in performing the inelastic analysis. Unlike conventional inelastic analysis methods of changing stiffness, the FAM analyzes structures by varying the structural displacement field, and therefore the state transition matrix needs to be computed only once. This greatly simplifies the computation and makes inelastic analysis readily applicable to the PIOC algorithm. The proposed algorithm compensates for the time delay that happens in practical control systems by predicting the inelastic structural response over a period that equals the magnitude of the time delay. A one‐story frame with both strain‐hardening and strain‐softening inelastic characteristics is analyzed using this algorithm. Results show that the proposed control algorithm is feasibile for any inelastic structures. While the control efficiency deteriorates with the increase in magnitude of the time delay, the PIOC maintains acceptable performance within a wide range of time delay magnitudes. Finally, a computer model of a six‐story moment‐resisting steel frame is analyzed to show that PIOC has good control results for real inelastic structures. Copyright © 2003 John Wiley & Sons, Ltd.     

冲绳海槽南部沉积层序的粒度特征

对取自冲绳海槽南部的A2 3孔经前处理后的沉积物样品进行了粒度测试,探讨了陆源沉积层序的粒度特征。沉积物平均粒径、分选系数、偏态、峰态等粒度参数的变化显示出A2 3孔的沉积层序以 4 0 0cm为界可分为上、下两段,下段各参数稳定,代表了比较稳定的水动力条件的沉积环境;上段则波动较大,小规模浊流沉积频繁发育,暗示了不稳定的沉积环境。浊流层内粗粒与细粒组分含量及粒度分布的变化具有明显的规律性。通过分析A2 3孔陆源沉积层序的各沉积参数特征,认为冲绳海槽南部的沉积环境和动力在时间尺度上有较大的变化,东海陆架物质向南部海槽的输运可能是以峡谷为通道的重力作用引起的床底沉积作用为主。另外,对陆源碎屑和全样粒度结果的对比,显示了前处理方法对于海洋沉积物粒度分析的重要性.     

鄱阳湖组(第四系)的修订及特征

根据对鄱阳湖地区 Zk0 1钻孔柱状剖面的研究 ,对鄱阳湖组的定义进行修订、补充 ,认为鄱阳湖组为一套以青灰色淤泥质粘土、黄褐色粘土、粉质粉土为主夹少量粉砂或细砂的湖泊—河口三角洲相沉积。年代地层为第四系全新统 ,底界年龄约 6 2 10± 12 5 a B.P.,开阔湖泊沉积始于 32 4 0± 180 a B.P.,且由北向南时代渐新 ,具有明显的穿时性