Micromechanical viscoelasto‐plastic models and finite element implementation for rate‐independent and rate‐dependent permanent deformation of stone‐based materials

This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.     

Rapid Identification and Emergency Investigation of Surface Ruptures and Geohazards Induced by the M_s 7.1 Yushu Earthquake

<正>The rapid identification based on InSAR technology was proved to be effective in our emergency investigation of surface ruptures and geohazards induced by the Yushu earthquake.The earthquake-generating fault of the Yushu earthquake is the Yushu section of the Garze-Yushu faults zone.It strikes NWW-NW,23 km long near the Yushu County seat,dominated by left-lateral strike slip,and appearing as a surface rupture zone.The macroscopic epicenter is positioned at Guo-yang-yan -song-duo of Gyegu Town(33°03'11N,96°51'26E),where the co-seismic horizontal offset measured is 1.75 m.Geohazards induced by the Yushu earthquake are mainly rockfalls,landslides,debris flows, and unstable slopes.They are controlled by the earthquake-generating fault and are mostly distributed along it.There are several geohazard chains having been established,such as earthquake,canal damage,soil liquefying,landslide-debris flow,earthquake,soil liquefying,roadbed deformation,etc.In order to prevent seismic hazards,generally,where there is a visible surface rupture induced by the Yushu earthquake,reconstruction should be at least beyond 20 m,on each side,from it.Sufficient attention should also be given to potential geohazards or geohazard chains induced by the earthquake.     

附加系统参数的多时相InSAR时空建模和形变估计

轨道误差和长波大气延迟组成的系统误差是影响InSAR形变监测精度的重要因素之一.传统方法在空间域对干涉图的系统误差建模,容易导致长波形变和系统误差相混淆.本文在时空域利用附加系统参数对系统误差建模,同时根据观测值质量对差分相位观测值定权,采用附加系统参数的加权最小二乘法估计形变参数和系统误差,实现了长波形变和系统误差的分离.模拟实验结果表明,在形变与系统误差的空间变化特性完全一致的极端情况下,本文方法能实现两者的有效分离,估计的形变速率均方根误差比传统方法降低了98.8%.ASAR数据实验显示当形变尺度较小且分散分布时,本文方法和传统方法得到的结果相似;当形变在研究区内表现为长波变化时,本文方法比传统方法估计的形变结果更为稳健.     

气象固定资产管理存在问题的探讨

通过对气象固定资产管理中存在的问题及成因分析,提出加强气象固定资产管理的一些建议。     

Inflation along Kilauea's Southwest Rift Zone in 2006

We report on InSAR and GPS results showing the first crustal inflation along the southwest rift zone at Kilauea volcano in over 20 years. Two independent interferograms (May 2–August 2, 2006 and June 22–Nov 7, 2006) from the ALOS PALSAR instrument reveal domal uplift located southwest of the main caldera. The uplift is bounded on the northeast by the caldera and follows the southwest rift zone for about 12 km. It is approximately 8 km wide. We use data derived from permanent GPS stations to calibrate the InSAR displacement data and estimate uplift of 7.7 cm during the first interferogram and 8.9 cm during the second with line-of-sight volumes of 2.8 × 10⁶ m³ and 3.0 × 10⁶ m³ respectively. The earthquake record for the periods before, during, and after inflation shows that a swarm of shallow earthquakes (z < 5 km) signaled the beginning of the uplift and that elevated levels of shallow seismicity along the rift zones occurred throughout the uplift period. GPS data indicate that the inflation occurred steadily over nine months between mid-January and mid-October, 2006 making injection of a sill unlikely. We attribute the inflation to recharge of a shallow ductile area under the SWRZ.     

根据2001年Mw 78可可西里强震InSAR同震测量结果反演东昆仑断裂两侧地壳弹性介质差异

利用2001年 Mw 78 可可西里强震InSAR同震测量结果，反演了青藏高原北部东昆仑断裂两侧地壳弹性介质差异.InSAR测量结果显示断层南侧的同震位移比北侧的大20％～30%.根据人工地震反射剖面建立岩石圈模型，以断层两侧杨氏模量差异和震源破裂深度为反演变量，通过有限元方法模拟实测得到的同震位移剖面.反演得到最佳断层破裂深度为20～22km，断层南侧杨氏模量相对北侧比值为81%～92%.结果表明，断层两侧弹性介质性质存在明显差异，由于构造运动作用，断层南部地壳不及北部地壳坚硬.前人利用地震层析成像和大地电磁测深等手段推断青藏高原内昆仑山断裂以南可可西里－羌塘地块地壳内广泛发育低速高导层，我们通过形变场力学分析得到与此相一致的结果.     

遥感技术在地震研究中的应用进展

本文介绍了遥感技术在地震监测和研究中的应用和进展,分析了卫星热红外遥感、InSAR技术和气体遥感方法在地震前后地表热红外增温、形变和气体含量改变的研究现状,初步探讨了该方法的成因机制.本文同时分析了该技术在地震预报中存在的问题,提出了未来的研究重点.     

Centrifuge modeling of buried continuous pipelines subjected to normal faulting

Seismic ground faulting is the greatest hazard for continuous buried pipelines.Over the years,researchers have attempted to understand pipeline behavior mostly via numerical modeling such as the finite element method.The lack of well-documented field case histories of pipeline failure from seismic ground faulting and the cost and complicated facilities needed for full-scale experimental simulation mean that a centrifuge-based method to determine the behavior of pipelines subjected to faulting is best to verify numerical approaches.This paper presents results from three centrifuge tests designed to investigate continuous buried steel pipeline behavior subjected to normal faulting.The experimental setup and procedure are described and the recorded axial and bending strains induced in a pipeline are presented and compared to those obtained via analytical methods.The influence of factors such as faulting offset,burial depth and pipe diameter on the axial and bending strains of pipes and on ground soil failure and pipeline deformation patterns are also investigated.Finally,the tensile rupture of a pipeline due to normal faulting is investigated.     

干涉雷达人工角反射器技术

人工角反射器是众多雷达目标中的典型散射结构.本文介绍了人工角反射器技术应用在雷达干涉测量(CD-InSAR)的基本概念,描述了其中常用的二面角、三角锥形三面角、正方体三面角和长方体锥形三面角等四种典型的人工角反射器特点和最大雷达散射截面的计算机模拟,最后对人工角反射器制作和布设中的一些注意事项及主要应用领域进行了探讨.     

InSAR analysis of a blind thrust rupture and related active folding: the 1999 Ain Temouchent earthquake (M w 5.7, Algeria) case study

We study the surface deformation associated with the 22 December 1999 earthquake, a moderate sized but damaging event at Ain Temouchent (northwestern Algeria) using Interferometric Satellite Aperture Radar images (InSAR). The mainshock focal mechanism solution indicates reverse faulting with a NE–SW trending rupture comparable to other major seismic events of this section of the Africa–Eurasia plate boundary. Previously, the earthquake fault parameters were, however, poorly known because no aftershocks were precisely determined and no coseismic surface ruptures were observed in the field. Using a pair of ERS data with small baseline and short temporal separation in the ascending orbit we obtained an interferogram that shows the coseismic surface displacement field despite poor coherence. The interferogram measures four fringes and displays an ellipse-shaped lobe with ∼11 cm peak line-of-sight displacement. The elastic modeling using a boundary element method (Poly3Dinv) indicate coseismic slip reaching up to 1 m at 5 km depth on the N 57° E trending, dipping 32° NW Tafna thrust fault. The geodetic estimate of seismic moment is 4.7 × 10¹⁷ N m. (Mw 5.7) in is good agreement with seismological results. The elliptical shape of the surface displacement field coincides with the NE–SW trending Berdani fault-related fold. The consistency between the geological observations and InSAR solution shed light on the precise earthquake location and related Tafna fault parameters.