横观各向同性饱和层状地基的三维稳态动力响应

首先引入状态向量,将直角坐标系下横观各向同性饱和土的Biot波动方程转化为一组状态方程,然后基于双重Fourier变换,求解了状态方程并得到传递矩阵。进而利用传递矩阵,并结合饱和地基的边界条件、排水条件及层间接触和连续条件,求解了横观各向同性饱和层状地基的稳态动力响应问题。数值算例表明采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的动力特性。     

An elliptical disc anchor in a damage‐susceptible poroelastic medium

A flat rigid elliptical anchorage located in a damage‐susceptible fluid‐saturated poroelastic medium is subjected to an in‐plane load, which induces a pure translation in the plane of the anchor. This paper develops computational estimates for the time‐dependent displacement of the disc anchor for the classical problem that involves Biot consolidation and compares the results with situations where the porous skeleton can experience micro‐mechanical damage that leads to an alteration in both its elasticity and fluid transport characteristics. Copyright © 2005 John Wiley & Sons, Ltd.     

A fully coupled finite element analysis of water‐table fluctuation and land deformation in partially saturated soils due to surface loading

A fully coupled numerical model is presented for the water‐table fluctuation and land deformation in partially saturated soils due to surface loading. This numerical model is developed based on the poroelastic governing equations for groundwater flow in deforming variably saturated porous media and the Galerkin finite element method. The numerical model is verified and validated against a one‐dimensional consolidation problem concerning surface loading on a soil column which has six different initial water‐table elevations. The numerical model is then applied to a two‐dimensional consolidation problem of surface loading on a partially saturated soil at a construction site. Results from the numerical simulations of both problems show that the water table fluctuates in the partially saturated soils, and the unsaturated zone above the water table has significant effects on the consolidation behaviour of the partially saturated soils under surface loading. Such effects are caused by the permanent absorption of a portion of the mechanical loading stress and the weak hydromechanical coupling between the solid skeleton deformation field and the groundwater flow field in the unsaturated zone due to its partial saturation. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamic behavior of saturated poroviscoelastic media

Summary The Vardoulakis-Beskos model for the dynamics of nearly and fully saturated poroelastic soils and the Aifantis-Beskos model for the dynamics of fully saturated, fissured, poroelastic rocks are here extended to include viscoelastic material behavior. Linear hereditary isotropic viscoelasticity of the relaxation type for the solid phase is considered. Correspondence principles are established for both models in the Laplace transformed domain. The one-dimensional dynamic column problem associated with both of the above proroviscoelastic soil and rock models is solved by an analytical-numerical procedure to illustrate the theory and assess the effect of viscoelasticity on the response.     

A coupled discrete element‐finite difference approach for modeling mechanical response of fluid‐saturated porous materials

A model of fluid‐saturated poroelastic medium was developed based on a combination of the discrete element method and grid method. The developed model adequately accounts for the deformation, fracture, and multiscale internal structure of a porous solid skeleton. The multiscale porous structure is taken into account implicitly by assigning the porosity and permeability values for the enclosing skeleton, which determine the rate of filtration of a fluid. Macroscopic pores and voids are taken into account explicitly by specifying the computational domain geometry. The relationship between the stress–strain state of the solid skeleton and pore fluid pressure is described in the approximations of simply deformable discrete element and Biot's model of poroelasticity. The developed model was applied to study the mechanical response of fluid‐saturated samples of brittle material. Based on simulation results, we constructed a generalized logistic dependence of uniaxial compressive strength on loading rate, mechanical properties of fluid and enclosing skeleton, and on sample dimensions. The logistic form of the generalized dependence of strength of fluid‐saturated elastic–brittle porous materials is due to the competition of two interrelated processes, such as pore fluid pressure increase under solid skeleton compression and fluid outflow from the enclosing skeleton to the environment. Copyright © 2015 John Wiley & Sons, Ltd.     

Detection and Modelling of Nonlinear Elastic Response in Damaged Composite Structures

In this paper the nonlinear material response of damaged composite structures under periodic excitation is experimentally and numerically investigated. In particular, the nonlinear wave propagation problem was numerically analysed through a finite element model able to predict the nonlinear interaction of acoustic/ultrasonic waves with damage precursors and micro-cracks. Such a constitutive model is based on the Landau’s semi-analytical approach to account for anharmonic effects of the medium, and is able to provide an understanding of nonlinear elastic phenomena such as the second harmonic generation. Moreover, Kelvin tensorial formulation was used to extend the wave propagation problem in orthotropic materials to the 3D Cartesian space. In this manner, the interaction of the stress waves with the 3D crack could be analysed. This numerical model was then experimentally validated on a composite plate undergone to impact loading. Good agreement between the experimental and numerical second harmonic response was found, showing that this material model can be used as a simple and useful tool for future structural diagnostic applications.     

Damage effects of adhesives in modern glass façades: a micro-mechanically motivated volumetric damage model for poro-hyperelastic materials

This paper presents a micro-mechanically motivated volumetric damage model accounting for cavitation effects in modern glass connections, e.g. laminated glass connections. The volumetric part of an arbitrary Helmholtz free energy function is equipped with an isotropic damage formulation. To develop a micro-mechanical damage model, the porous micro-structure of a transparent structural silicone adhesive is analyzed numerically applying hydrostatic loading conditions. Based on the structural responses of different types of cubic representative volume elements incorporating an initial void fraction, three damage parameters are fitted utilizing the Levenberg–Marquard algorithm. The present volumetric damage model is implemented into ANSYS FE Code using a UserMat subroutine, where the algorithmic setting is described in detail in the present paper. To compare the structural responses of cubic equivalent homogeneous materials with representative volume elements, benchmark tests under hydrostatic loading are performed. The results indicate that the novel damage model accounts adequately for volumetric damage due to the cavitation effect. A special form of the pancake test is described briefly. The test allows for visualizing the cavitation effect during experimental testing. The experimental results of the pancake test are compared with numerical results, where the pancake test is simulated incorporating the micro-mechanical damage model. The micro-mechanically motivated scalar, internal damage variable is equipped with the obtained damage parameters from the structural response of the representative volume elements. The results show an adequate approximation of the experiment through the simulation. However, to optimize the results of the simulation, an optimization study on the damage parameters is conducted utilizing the Downhill-Simplex algorithm. Using the optimized damage parameters, the simulation of the pancake tests is further improved. Hence, it is shown that the novel micro-mechanically motivated volumetric damage model is excellently suited to represent the cavitation effect in poro-hyperelastic materials.     

Leakage through a permeable capillary tube into a poroelastic tumor interstitium

Fluid flow through a permeable circular tube embedded in an infinite poroelastic ambient medium is studied as a model of blood flow through the vasculature of a solid tumor. The flow through the interstitium is described by Darcy's law for an isotropic porous medium with a pressure-dependent permeability, the flow along the tube is described by Poiseuille's law, and the extravasation flux across the tube surface is described by Starling's law involving the transmural pressure. Kirchhoff's transformation is applied to derive Laplace's equation for a modified interstitial pressure. Given the arterial, venous, and ambient pressures, the problem is formulated in terms of a coupled system of integral, differential, and algebraic equations for the vascular and interstitial pressures. The overall hydrodynamics is described in terms of hydraulic conductivity coefficients for the arterial, venous, and extravasation flow rates. Solutions obtained by a boundary-element method confirm that interstitium dilatation promotes the rate of extravasation.     

BEM modeling of saturated porous media susceptible to damage

This paper deals with the numerical analysis of saturated porous media, taking into account the damage phenomena on the solid skeleton. The porous media is taken into poro-elastic framework, in full-saturated condition, based on Biot’s Theory. A scalar damage model is assumed for this analysis. An implicit boundary element method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and two-dimensional elastostatic problems. The integration over boundary elements is evaluated using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the relevant domain integrals. The non-linear problem is solved by a Newton-Raphson procedure. Numerical examples are presented, in order to validate the implemented formulation and to illustrate its efficacy.     

Porous asymmetric SiO2-g-PMMA nanoparticles produced by phase inversion

A new kind of asymmetric organic–inorganic porous structure has been proposed. Asymmetric lattices of polymer grafted silica nanoparticles were manufactured by casting and phase inversion in water. Silica nanoparticles were first functionalized with 3-(dimethylethoxysilyl)propyl-2-bromoisobutyrate, followed by grafting of poly(methylmethacrylate) (PMMA) segments, performed by atom-transfer radical polymerization. Mechanically stable self-standing films were prepared by casting a dispersion of functionalized nanoparticles in different solvents and immersion in water. The resulting asymmetrically porous morphology and nanoparticle assembly was characterized by scanning electron and atomic force microscopy. The PMMA functionalized SiO₂ hybrid material in acetone or acetone/dioxane led to the best-assembled structures. Porous asymmetric membranes were prepared by adding free PMMA and PMMA terminated with hydrophilic hydroxyl group. Nitrogen flow of 2800 L m^?2 h^?1 was measured at 1.3 bar demonstrating the porosity and potential application for membrane technology.     

A Contribution to Time-Dependent Damage Modeling of Composite Structures

The paper presents a new damage model for predicting stiffness loss due to creep loading and cyclic fatigue. The model, developed within a continuum damage mechanics framework, is based on the idea of a time-dependent damage spectrum, some elements of which occur rapidly and others slowly. The use of this spectrum allows a single damage kinematic to model creep and fatigue damage and to take into account the effect of stress amplitude, R ratio, and frequency. The evolution equations are based on similar equation than the one describing the viscoelasticity model and are relatively easy to implement. The new model is compared to the experimental results on carbon fiber/epoxy tubes. Quasi-static, creep and fatigue tests are performed on filament-wound tubular specimens to characterize the elastic, viscoelastic and plastic behavior of the composite material. Varying amounts of damage are observed and discussed depending on stress level and R ratio. The experimental work aims to develop and validate the damage model for predicting stiffness loss due to creep loading and cyclic fatigue.     

Porous copper template from partially spark plasma-sintered Cu-Zn aggregate via dezincification

Present work deals with the preparation of spark plasma-sintered Cu-Zn aggregate (5, 10 and 20 wt% Zn) with interfacial bonding only starting from elemental powders of Cu and Zn (99.9% purity) and subsequently making of porous template of Cu by dezincification. Sintering is done so as to achieve only interfacial bonding with the aim to maintain maximum potential difference between the Cu and Zn particles during dezincification process in various solutions, viz. 1 N HCl and 3.5 wt% NaCl solutions. X-ray diffraction, optical microscopy and SEM-EDS are carried out to examine microstructural evolution and subsequent changes in hardness with sintering temperatures and different Zn percentages. Dezincification and pore formation are conducted on sintered 0.5 mm thick 12 mm diameter disc samples. The size, distribution and nature of pores in porous templates of Cu are then investigated using optical microscopy and SEM-EDS analysis.     

各向异性双重孔隙介质的应力与油水两相渗流耦合理论模型

针对天然裂缝性油藏的特性,建立了描述双重孔隙介质中油水两相流体流动特性的流固耦合理论模型。该模型不仅考虑了渗透率的各向异性,而且考虑了岩石固体骨架变形的各向异性。渗流方程是依据双重孔隙的概念建立起来的,而固体骨架变形控制方程则是根据Biot 的等温、线性孔隙弹性理论建立起来的。同时,给出了横向各向同性及结构各向异性、固体材料各向同性时的双重孔隙介质的应力与油水两相渗流耦合理论模型。对该模型进行了简化,并将其简化后模型与单相流的各项同性和各向异性双重孔隙介质流固耦合理论模型进行了比较。     

A Spherically Anisotropic Microstructure Model for Fluid-saturated Poroelastic Isotropic Materials

The microstructure model proposed in this work is an assemblage of hollow spheres saturated by a fluid. The solid phase of each sphere is linearly elastic and spherically anisotropic. On the basis of this microstructure model, the effective bulk modulus, Biot’s coefficient and porosity variation are determined. It is shown that local anisotropy has important effects on the macroscopic isotropic poroelastic properties via a dimensionless material parameter which characterizes the degree of anisotropy and exponentially affects the porosity.     

Heat transfer in photonic mirrors

The use of secondary mirrors in solar energy concentration is common. However, high concentrated solar radiation heats these mirrors thereby degrading their physical properties. In particular, aluminum mirrors melt because of high temperature due to storage by high radiative heat transfer. In contradistinction photonic crystals could present “perfect reflection” and they can be fabricated using porous silicon which has a higher melting point than aluminum (porous silicon has a melting point higher than 900 K). Porous silicon is a nanostructured semiconductor material which can be fabricated with different porosities and refractive indices. Multilayers of alternating periodic refractive index conform the structure of these photonic crystals. The light that propagates in these structures interacts with its periodic refractive index that generates wavelength gaps of forbidden transmission and so these multilayers conform a mirror. Even these photonic structures are heated when they are exposed to high concentrated solar radiation. In this work we experimentally analyze this heating process and model it using an effective medium approach to explain the increasing temperature behavior.     

Continuum effective-stress approach for high-rate plastic deformation of fluid-saturated geomaterials with application to shaped-charge jet penetration

A practical engineering approach for modeling the constitutive response of fluid-saturated porous geomaterials is developed and applied to shaped-charge jet penetration in wellbore completion. An analytical model of a saturated thick spherical shell provides valuable insight into the qualitative character of the elastic–plastic response with an evolving pore fluid pressure. However, intrinsic limitations of such a simplistic theory are discussed to motivate the more realistic semi-empirical model used in this work. The constitutive model is implemented into a material point method code that can accommodate extremely large deformations. Consistent with experimental observations, the simulations of wellbore perforation exhibit appropriate dependencies of depth of penetration on pore pressure and confining stress.     

饱和软粘土地基的损伤模型与震陷计算

基于各向同性弹塑性损伤和Prevost模型的基本理论,把弹塑性等向硬化、运动硬化和各向同性损伤结合起来,推导了循环荷载作用下不排水饱和软粘土弹塑性动力损伤本构模型.考虑到地震作用下土体应力的不规则性,对循环三轴试验获得的粘土残余应变经验计算公式进行了修正,最后将该残余应变引入到损伤演化方程中.通过对地基的震陷计算并与不考虑损伤的模型进行对比,结果表明,该模型能合理地考虑屈服面内应力循环对地基残余塑性应变的贡献.     

Effect of phenolic resin infiltration content on the structural and electrochemical properties of hierarchical porous carbons

Hierarchical porous (micro-, meso- and macro-porous) carbons (HPCs) are synthesized by a facile replica template method with phenolic resin (PR) as a carbon source and hollow mesoporous silica as a hard template. The morphology of the HPCs can be easily controlled by altering the mass ratio of PR to SiO₂ spheres. Structural characterizations reveal that a well-defined HPC with a large surface area of 1141 m² g^?1 is formed when PR/SiO₂ is 1:1. With further increasing PR infiltration content, macropores of carbons disappear while solid structures appear. A possible formation mechanism for the morphological transformation of HPCs is proposed. The effect of phenolic resin infiltration content on the electrochemical properties of HPCs-based materials as supercapacitor electrodes is further evaluated. The HPCs-1(PR/SiO₂ = 1) electrode exhibits the highest specific capacitance of 256 F g^?1 due to its faster diffusion of electrolyte ions and lower charge transfer resistance. The relationship between the morphology and the electrochemical behavior of HPCs is also elucidated.     

Anisotropic Gurson‐Tvergaard‐Needleman plasticity and damage model for finite element analysis of elastic‐plastic problems

Implementation and analysis of the anisotropic version of the Gurson‐Tvergaard‐Needleman (GTN) isotropic damage criterion are performed on the basis of Hill's quadratic anisotropic yield theory with the definition of an effective anisotropic coefficient to represent the elastic‐plastic behavior of ductile metals. This study aims to analyze the extension of the GTN model suitable for anisotropic porous metals and to investigate the GTN model extension. An anisotropic damage model is implemented using the user material subroutine in ABAQUS/standard finite element code. The implementation is verified and applied to simulate a uniaxial tensile test on a commercially produced aluminum sheet material for three‐dimensional and plane stress test cases. Spherical and ellipsoidal micro voids are considered in the matrix material, and their effects on the uniaxial stress‐strain response of the material are analyzed. Hill's quadratic anisotropic yield theory predicts substantially large damage evolution and a low stress‐strain curve compared with those predicted by the isotropic model. An approximate model for anisotropic materials is proposed to avoid increased damage evolution. In this approximate model, Hill's anisotropic constants are replaced with an effective anisotropy coefficient. All model‐generated stress‐strain predictions are compared with the experimental stress‐strain curve of AA6016‐T4 alloy.     

横观各向同性土体三维比奥固结有限层解法

建立了三维横观各向同性土体固结的有限层求解方法,并编制了相应的计算程序;通过2个算例的分析对比,验证了计算方法和程序的正确性;探讨了横观各向同性地基参数对三维地基固结的影响,给出了对工程有价值的相关计算图表和结论.