工程结构三维疲劳裂纹最大应力强度因子计算

介绍了裂纹的类型、裂纹尖端应力场的奇异性。以一维问题为例,推导论证了奇异单元能够很好的反映裂纹尖端应力场的奇异性。应力强度因子一般表达式表明应力强度因子与载荷呈线性关系,并依赖于物体和裂纹的几何形状和尺寸。本文借助大型通用有限元软件ANSYS,采用位移外插法计算了三维表面裂纹前沿不同位置处的应力强度因子,并与《应力强度因子手册》基于实验的理论公式计算结果相比较。结果表明:有限元结果与理论解误差较小,裂纹最深处应力强度因子最大。     

Ring foundation on elastic subgrade: an analytical solution for computer modelling using the Lagrangian multiplier method

In the practice of geotechnical engineering, the case of a ring footing carrying a set of concentrated point loads is a common problem. At times, the induced vertical and angular displacements for the ring footing need to be evaluated at a relatively precise level. By making use of the governing set of equations derived for the case of a general curved beam, expressions that can be easily implemented in modern computing software are derived for the vertical and angular displacements of a ring footing of rectangular cross section as functions of the radial position. The loading case considered is a vertical point load, and the soil is modelled as elastic. Estimates of the displacements have been shown for a common range of practical applications. The behaviour for a set of concentrated loads may be evaluated using the derived equations through direct superposition. Nonlinear finite element analysis is used to evaluate the vertical deflection and angular twist of the ring foundation. Numerical analysis performed for three ring foundations with different radii and cross sections is reported to validate the accuracy of the derived analytical solution. Copyright © 2016 John Wiley & Sons, Ltd.     

A constitutive model for mechanical and chemo‐mechanical compaction in sedimentary basins and finite element analysis

Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd.     

An approach to the modeling of viscoelastic damage. Application to the long‐term creep of gypsum rock materials

A three‐dimensional phenomenological model is developed to describe the long‐term creep of gypsum rock materials. The approach is based on the framework of continuum damage mechanics where coupling with viscoelasticity is adopted. Specifically, a local damage model based on the concept of yield surface is proposed and deeply investigated. Among the many possibilities, we choose in this work its coupling with a generalized Kelvin–Voigt rheological model to formulate the whole behavior. Long‐term as well as short‐term relaxation processes can be integrated in the model by means of as many as necessary viscoelastic processes. The numerical discretization is described for an easy integration within a finite element procedure. Finally, a set of numerical simulations is given to show the possibilities of the presented model. It shows good agreement with some experimental results found in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Palaeoflood activity and climate change over the last 1400 years recorded by lake sediments in the north‐west European Alps

A high‐resolution sedimentological and geochemical study of a high‐altitude proglacial lake (Lake Blanc, Aiguilles Rouges, 2352 m a.s.l.) revealed 195 turbidites, 190 of which are related to flood events over the last 1400 years. We used the coarsest sediment fraction of each turbidite as a proxy for the intensity of each flood event. Because most flood events at this locality are triggered by localized summer convective precipitation events, the reconstructed sedimentary record reveals changes in the frequency and intensity of such events over the last millennium. Comparisons with other temperature, palaeohydrological and glacier reconstructions in the region suggest that the most intense events occurred during the warmest periods, i.e. during the Medieval Climate Anomaly (AD 800–1300) and the current period of global warming. On a multi‐decadal time scale, almost all the flood frequency peaks seem to correspond to warmer periods, whereas multi‐centennial variations in flood frequency appear to follow the regional precipitation pattern. Consequently, this new Alpine flood record provides further evidence of a link between climate warming and an increase in the frequency and intensity of flooding on a multi‐decadal time scale, whereas the centennial variability in flood frequencies is related to regional precipitation patterns. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Coupled finite‐element simulation of injection well testing in unconsolidated oil sands reservoir

This paper presents a finite‐element (FE) model for simulating injection well testing in unconsolidated oil sands reservoir. In injection well testing, the bottom‐hole pressure (BHP) is monitored during the injection and shut‐in period. The flow characteristics of a reservoir can be determined from transient BHP data using conventional reservoir or well‐testing analysis. However, conventional reservoir or well‐testing analysis does not consider geomechanics coupling effects. This simplified assumption has limitations when applied to unconsolidated (uncemented) oil sands reservoirs because oil sands deform and dilate subjected to pressure variation. In addition, hydraulic fracturing may occur in unconsolidated oil sands when high water injection rate is used. This research is motivated in numerical modeling of injection well testing in unconsolidated oil sands reservoir considering the geomechanics coupling effects including hydraulic fracturing. To simulate the strong anisotropy in mechanical and hydraulic behaviour of unconsolidated oil sands induced by fluid injection in injection well testing, a nonlinear stress‐dependent poro‐elasto‐plastic constitutive model together with a strain‐induced anisotropic permeability model are formulated and implemented into a 3D FE simulator. The 3D FE model is used to history match the BHP response measured from an injection well in an oil sands reservoir. Copyright © 2013 John Wiley & Sons, Ltd.     

Analyse fréquentielle régionale des précipitations journalières maximales annuelles au Québec,Canada / Regional frequency analysis of annual maximum daily precipitation in Quebec,Canada

Abstract

Abstract A complete regional analysis of daily precipitations is carried out in the southern half of the province of Quebec, Canada. The first step of the regional estimation procedure consists of delineating the homogeneous regions within the area of study and testing for homogeneity within each region. The delineation of homogeneous regions is based on using L-moment ratios. A simulation-based testing of statistical homogeneity allows one to verify the inter-site variability. The second step of the procedure deals with the identification of the regional distribution and the estimation of its parameters. The General Extreme Value (GEV) distribution was identified as an appropriate parent distribution. This distribution has already been recommended by several previous research studies for regional frequency analysis of precipitation extremes. The parameters of the GEV distribution are estimated based on the computation of the regional L-CV, L-C_S and the mean of annual maximal daily precipitations. The third step consists of the estimation of precipitation quantiles corresponding to various return periods. The final procedure allows for the estimation of these quantiles at sites where no precipitation information is available. The use of a jack-knife resampling procedure with data from the province of Quebec allows one to demonstrate the robustness and efficiency of the regional estimation procedure. Values of the root mean square error were below 10% for a return period of 20 years, and 20% for a return period of 100 years.     

Dynamics of unsaturated poroelastic solids at finite strain

We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called u^s ? p^w ? p^a formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element analysis of the effect of corrosion on the behavior of reinforced earth walls

A simple model for the corrosion‐induced loss of stiffness and strength of the steel strips of earth‐reinforced walls was introduced in a finite element simulation of the long‐term behavior of the wall, in which the backfill‐strips interactions are taken into account by means of a generalized homogenization procedure (called a multiphase model). The results show an initial phase of slow displacements induced by the loss of stiffness, followed after a few decades by a steep acceleration of the displacements, leading to wall failure. The influences of the parameter controlling corrosion, the backfill cohesion and the heterogeneity of the corrosion process are discussed. Results are used to discuss a strategy for reinforced earth wall surveillance. Copyright © 2011 John Wiley & Sons, Ltd.