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.     

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.     

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.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 John Wiley & Sons, Ltd.     

区域地质灾害危害程度等级划分探讨

区域性地质灾害危害程度等级尚无比较科学的划分方法。通过福建省区域地质灾害危害程度等级划分的探讨,采用单位面积地质灾害造成的人员死亡数、财产损失值作为“地质灾害模数”,对区域地质灾害危害程度等级进行划分,判别不同区域地质灾害严重程度,为区域地质灾害防治决策提供科学依据。     

福建省滨海火电厂地质灾害问题及风险控制探讨

滨海火力发电厂工程主要包括厂区建筑、码头、管道、取排水、填海和贮灰场等工程。其主要面临着福建省海岸带构造运动、断裂及地震活动、港湾淤积、海底滑坡、软土地基、海底活动地貌、基岩不均匀风化以及人类工程活动等主要的灾害性地质因素。通过对这些因素潜在的致灾特点分析，提出了滨海火电厂地质灾害风险控制应包括选址阶段地质灾害风险回避、设计施工阶段地质灾害风险处理及运行阶段地质灾害风险监控等3方面。地质灾害风险评估是滨海火电厂地质灾害风险控制的首要任务。针对滨海电厂工程的特点，评估内容应着重于地质灾害危险性评估及易损性评估。选址阶段地质灾害风险回避主要是对构造不稳定的回避。地质灾害风险处理主要是电厂工程的基础处理及管道抗冲刷处理。电厂运行阶段地质灾害风险监控主要是对建筑物基础稳定性及海域冲淤变化的监控。     

Nonlinear earthquake response analysis and energy calculation for seismic slit shear wall structures

Based on the concept of structural passive control, a new type of slit shear wall, with improved seismic performance when compared to an ordinary solid shear wall, was proposed by the authors in 1996. The idea has been verified by a series of pseudo-static and dynamic tests. In this paper a macro numerical model is developed for the wall element and the energy dissipation device. Then, nonlinear time history analysis is carried out for a 10-story slit shear wall model tested on a shaking table. Furthermore, the seismic input energy and the individual energy dissipated by the components are calculated by a method based on Newmark-β assumptions for this shear wall model, and the advantages of this shear wall are further demonstrated by the calculation results from the viewpoint of energy. Finally, according to the seismic damage criterion on the basis of plastic accumulative energy and maximum response, the optimal analysis is carried out to select design parameters for the energy dissipation device.     

基于应变模态法识别刚架桥梁的损伤

基于应变模态方法，对刚架桥的损伤识别进行了研究。通过对某刚架桥在不同损伤工况下的数值仿真计算，探讨了应变模态方法用于刚架桥损伤识别的能力。计算结果表明：利用应变模态差曲线能比较准确地识别出刚架桥的损伤位置；应变模态差曲线在刚架桥损伤单元处的跳跃幅值随单元损伤程度的增加而增大，依此可定性地识别出刚架桥的损伤程度。     

高密度电法勘探在岩溶查找中的应用

介绍了武汉市汤逊湖污水处理厂厂址的地质概况和地球物理特征，以及高密度电法勘探的施工方案，分析了视电阻率异常特征，并说明了其异常验证情况。高密度电法勘探结果表明该区存在溶洞，为施工提供了决策依据，建议在有溶洞的地段采取相应措施，以防止工程建设后发生地质灾害。     

高耸塔架结构节点损伤基于神经网络的两步诊断法

本文针对高耸钢塔架结构的损伤特点，建立了具有节点损伤的有限元模型，提出了一种分层神经网络两步诊断法：第一步，由基于区域残余力理论的第一层神经网络进行结构损伤区域的初诊；第二步，由基于应变模态理论的第二层神经网络进行损伤区域内的具体损伤节点位置和程度的诊断。对一平面塔架结构的数值仿真分析表明：本文提出的损伤诊断方法的结果是令人满意的。