Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

Porphyroblast inclusion-trail orientation data: eppure non son girate*!

Extensive examination of large numbers of spatially orientated thin sections of orientated samples from orogens of all ages around the world has demonstrated that porphyroblasts do not rotate relative to geographical coordinates during highly non-coaxial ductile deformation of the matrix subsequent to their growth. This has been demonstrated for all tectonic environments so far investigated. The work also has provided new insights and data on metamorphic, structural and tectonic processes including: (1) the intimate control of deformation partitioning on metamorphic reactions; (2) solutions to the lack of correlation between lineations that indicate the direction of movement within thrusts and shear zones, and relative plate motion; and (3) a possible technique for determining the direction of relative plate motion that caused orogenesis in ancient orogens.     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

H II regions and radio visibility of quasar nuclei

The core (injector) and the jet (relativistic plasma outflow) of AGN objects are surrounded by an ionized medium, an H II region observed in emission lines. The synchrotron radiation from the core and the jet is observed through a thin screen that cocoons the structure under consideration. The screen transparency depends on wavelength and distance from the injector. We consider the objects 3C 345 and 1803+784 whose core emission at decimeter wavelengths is absorption by more than 25 dB. The visible bright compact component is the bright nearby portion of the jet that extends outside the dense part of the screen. We explore the possibility of measuring the screen transparency from absorption in Hα recombination lines with different quantum numbers at centimeter wavelengths.     

网箱养殖罗非鱼综合症的病理组织学研究

于1989至1992年对水库网箱养殖罗非鱼综合症各发病期、各症状型的病鱼进行了病理组织学研究。研究结果表明:网箱养殖罗非鱼综合症的病理基础是肝肾功能障碍,病鱼的肝、肾、脾、肠、鳃、性腺、脑均有不同程度的坏死或变性。本文还对细菌继发性感染,各症状型的病理学变化及防治对策作了讨论。     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

青藏高原东部多年冻土退化对环境的影响

退化是本区多年冻土变化的基本趋势，不同区段退化幅度有差异，影响退化的因素除了全球气候转暖外，还有局地气候和地下水迳流。在退化敏感的地带，工程建设，牧区生态平衡都有不可忽视的影响。     

东营城区高精度三维地震采集方法研究

根据东营城区内建筑物较多，环境干扰较大，地下构造复杂，断裂非常发育，目的层埋藏深，资料信噪比较低的特点，以地震老资料和钻井资料为基础，建立地震地质模型，利用射线追踪技术选取合适的观测系统参数。充分利用卫星数字地图，没计了灵活多变的观测系统，采用可控震源与炸药震源配合使用方法，及时对城区内地下面元的覆盖次数、方位角和炮检距的分布进行质量监控，利用现场处理系统对资料进行及时分析，提高了地震资料品质。所获得的东营城区高精度三维地震采集资料显示，其浅中层分辨率和中深层信噪比都有了明显改善，取得了良好的地质效果。