A new method for solving Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system

A new method is developed to solve Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system. Based on the governing equations of Biot＇s consolidation and the technique of Laplace transform, Fourier expansions and Hankel transform with respect to time t, coordinate θ and coordinate r, respectively, a relationship of displacements, stresses, excess pore water pressure and flux is established between the ground surface （z = 0） and an arbitrary depth z in the Laplace and Hankel transform domain. By referring to proper boundary conditions of the finite soil layer, the solutions for displacements, stresses, excess pore water pressure and flux of any point in the transform domain can be obtained. The actual solutions in the physical domain can be acquired by inverting the Laplace and the Hankel transforms.     

Modeling of the role of defects in sintered FeCrAlY foams

The metal sintering approach offers a cost- effective means for the mass-production of open-cell foams from a range of materials, including high-temperature steel alloys, which offer novel mechanical and acoustic properties. In a separate experimental study, the mechanical properties of open-celled steel alloy (FeCrAlY) foams have been characterized under uniaxial compression and shear loading. Compared to predictions from established models, a significant knockdown in material properties was observed. This knockdown was attributed to the presence of defects throughout the microstructure that result from the unique fabrication process. In the present paper, the microstructure of sintered FeCrAlY foams was modeled by using a finite element (FE) model. In particular, microstructural variations were introduced to a base lattice, and the effects on the strength and stiffness calculated. A range of defects identified under scanning electronic microscope (SEM) imaging were considered including broken ligaments, thickness variations, and pore blockages, which are the three primary imperfections observed in sintered foams. The corresponding levels of defect present in the material were subsequently input into the FE model, with the resulting predictions correlating well with experimental data. The project supported by the National Basic Research Program of China (2006CB601202), the National Natural Science Foundation of China (10328203, 10572111, 10632060), the National 111 Project of China (B06024), and the US Office of Naval Research (N000140210117).     

ANALYTICAL SOLUTION OF THERMAL CONSOLIDATION OF SATURATED SOILS WITH TWO PARALLEL CYLINDRICAL HEAT SOURCES1）

基于应力平衡条件、渗流连续方程、能量守恒方程,考虑土颗粒和孔隙水热膨胀系数的不同,建立考虑热水力耦合的饱和土体三维热固结控制方程。利用傅里叶变换和拉普拉斯变换导出变换域上的控制方程,解得点热源在变换域和实数域上的解析解,再利用区域积分给出两平行圆柱形热源热固结土体温度、孔压、位移的解析解,并对其进行分析,发现径距比增大会导致两热源温度相互影响程度减弱,热固结系数减小会导致孔压和位移的峰值增大。     

Prolongation of two phases in the model of fluid displacement through a capillaryProlongement de deux phases dans le modèle du déplacement d'un fluide dans un capillaire

The problem of a piston-like displacement of a fluid by another in a capillary is examined. It is suggested that each fluid is prolonged into the domain occupied by the other fluid. This enables the replacement of the two-phase flow problem by a transient single-phase flow problem, with discontinuity in velocity and pressure on a film interface. The problems related to the triple point are solved by introducing a limit fluid near the pore wall. The demonstration of the Washburn equation contributes to the physical justification of our model. To cite this article: Y. Lucas et al., C. R. Mecanique 334 (2006).     

A sublinear-scaling approach to density-functional-theory analysis of crystal defects

We develop a sublinear-scaling method, referred to as MacroDFT, for the study of crystal defects using ab-initio Density Functional Theory (DFT). The sublinear scaling is achieved using a combination of the Linear Scaling Spectral Gauss Quadrature method (LSSGQ) and a Coarse-Graining approach (CG) based on the quasi-continuum method. LSSGQ reformulates DFT and evaluates the electron density without computing individual orbitals. This direct evaluation is possible by recourse to Gaussian quadrature over the spectrum of the linearized Hamiltonian operator. Furthermore, the nodes and weights of the quadrature can be computed independently for each point in the domain. This property is exploited in CG, where fields of interest are computed at selected nodes and interpolated elsewhere. In this paper, we present the MacroDFT method, its parallel implementation and an assessment of convergence and performance by means of test cases concerned with point defects in magnesium.     

铁电陶瓷电致疲劳失效的研究进展

铁电陶瓷是具备力电转换功能的典型高技术材料.本文概述铁电陶瓷电致疲劳失效的研究进展.首先介绍电致疲劳的定义和特点,然后讨论电致疲劳失效在不同尺度下的表现行为,包括宏观尺度下裂纹的疲劳扩展;细观尺度下裂纹的萌生;微观尺度下点缺陷在循环电场下的积聚.随即阐述了铁电陶瓷在循环电场下缺陷汇聚的理论分析,运用微结构演化方法计算了单个孔洞随畴界的移动距离,推导了循环电场下铁电陶瓷内点缺陷浓度的演化方程,给出了点缺陷浓度与其汇聚程度之间的定量关系,从而提出了贯通不同尺度的铁电陶瓷电致疲劳失效机理.      

A new analytical solution to axisymmetric Biot＇s consolidation of a finite soil layer

A new analytical method is presented to study the axisymmetric Biot＇s consolidation of a finite soil layer. Starting from the governing equations of axisymmetric Blot＇s consolidation, and based on the property of Laplace transform, the relation of basic variables for a point of a finite soil layer is established between the ground surface （z= 0） and the depth z in the Laplace and Hankel transform domains. Combined with the boundary conditions of the finite soil layer, the analytical solution of any point in the transform domain can be obtained. The actual solution in the physical domain can be obtained by inverse Laplace and Hankel transforms. A numerical analysis for the axisymmetric consolidation of a finite soil layer is carried out.     

Response of saturated porous media subjected to local thermal loading on the surface of semi-infinite space

Heat source function method is adopted in the present paper to derive elementary solutions of coupled thermo-hydro-mechanical consolidation for saturated porous media under conjunct actions of instantaneous point heat source, instantaneous point fluid source and constant volume force. By using the so-called fictitious heat source method and images method, the solutions of a semi-infinite saturated porous medium subjected to a local heat source with time-varied intensity on its free surface are developed from elementary solutions. The numerical integral methods for calculating the unsteady temperature, pore pressure and displacement fields are given. The thermomechanical response are analyzed for the case of a circular planar heat source. Besides, the thermal consolidation characteristics of a saturated porous medium subjected to a harmonic thermal loading are also given, and the fluctuation processes of the field variables located below the center of heat source are analyzed. The project supported by the National Natural Science Foundation of China (50578008) The English text was polished byYunming Chen.     

Ultrasonic characterization of point defects induced by cross slip under pure shear

Longitudinal wave velocity is used to characterize the point defects in crystalline solids. High purity Al single crystal was selected for both the finite element analysis and experimental work. Since the jog motions of dislocations caused by intersected slides such as cross slips induce point defects, the total amount of cross slips was calculated instead of calculating directly from the point defects. The effect of crystal orientations on total amount of cross slips under pure shear was also investigated via the finite element method. The result suggest that if the initial shear stress direction is located at the inner side of stereographic triangle, only single slip activities occurred at the beginning of plastic deformation and no effects due to point defects were present. However, as the shear stress direction rotates along the slip direction, point defects are induced by cross slips between primary and secondary slip systems due to work-hardening. This phenomenon was then examined by measuring longitudinal wave velocity changes propagating in Al single crystal subjected to the combination loads of equi-biaxial tension and compression (a pure shear state). Good qualitative agreement between the finite element result and measured data suggest that the longitudinal wave velocity can be used as an index to characterize point defects in crystalline materials.     

SPONTANEOUS POLARIZATION AND THE DISTRIBUTION OF INTERNAL STRESS ANDELECTRIC FIELDS AROUND A TRIPLE POINT

In this article, distributions of internal stress and internal electric fields around a triple point of ferroelectric polycrystals generated by the spontaneous deformation and spontaneous polarization were investigated. It was found that when all three grains consist of a single domain, the internal stresses and the internal electric fields do not vanish. Though it may be determined according to the principle of energy, the spontaneous configuration will not be unique without involving other conditions due to the symmetry of the crystal structure. Sponsored by National Natural Science Foundation of China(No. 19672053) and Special Funds for Doctoral Programs (No. 96061305).     

Meshless method based on collocation with consistent compactly supported radial basis functions

Based on our previously study, the accuracy of derivatives of interpolating functions are usually very poor near the boundary of domain when Compactly Supported Radial Basis Functions (CSRBFs) are used, so that it could result in significant error in solving partial differential equations with Neumann boundary conditions. To overcome this drawback, the Consistent Compactly Supported Radial Basis Functions (CCSRBFs) are developed, which satisfy the predetermined consistency conditions. Meshless method based on point collocation with CCSRBFs is developed for solving partial differential equations. Numerical studies show that the proposed method improves the accuracy of approximation significantly. The project supported by the National Natural Science Foundation of China (10172052)     

Numerical analysis of Biot's consolidation process by radial point interpolation method

An algorithm is proposed to solve Biot's consolidation problem using meshless method called a radial point interpolation method (radial PIM). The radial PIM is advantageous over the meshless methods based on moving least-square (MLS) method in implementation of essential boundary condition and over the original PIM with polynomial basis in avoiding singularity when shape functions are constructed. Two variables in Biot's consolidation theory, displacement and excess pore water pressure, are spatially approximated by the same shape functions through the radial PIM technique. Fully implicit integration scheme is proposed in time domain to avoid spurious ripple effect. Some examples with structured and unstructured nodes are studied and compared with closed-form solution or finite element method solutions.     

The singularity expansion method as applied to the elastodynamic scattering problem

A basic problem in non-destructive testing of materials is the classification of defects with respect to size and geometry. Utilizing elastrodynamic scattering, ultrasonic test methods reduce this problem either to the application of reconstruction algorithms or to the parametrization of experimental data in terms of quantities being related to the physical scattering mechanisms. The so-called Singularity Expansion Method (SEM), originally developed for broadband electromagnetic scattering by arbitrarily shaped targets proves essentially useful to predict and understand impulsive scattering of ultrasound; in addition, even though not yet fully solved in practice, it seems possible to parametrize experimentally obtained time records in a sense which is physically intuitive.

SEM starts either with the eigenfunction expansion in the complex Laplace domain for canonical objects such as a sphere or spheroids, or with a corresponding integral equation formulation for more arbitrarily shaped defects. The essential point is then an expansion in terms of the singularities of the scattered field in the Laplace domain similar to the expansion of transfer functions in linear system theory. It turns out that the location of these singularities is characteristic for the geometry of the scattering body; therefore, it might be a useful tool to parametrize size and shape of the defects.

Several theoretically derived singularity patterns are presented for various body shapes and material compositions, which yield a thorough and physically intuitive interpretation in terms of distinct creeping wave modes. They are compared with first experimental results for a spherical void in a steel test specimen.     

条形荷载下梯度非均匀非饱和土的动力响应分析

基于多孔介质混合物理论, 建立了梯度非均匀非饱和土地基模型, 研究了条形荷载作用下梯度非均匀非饱和土地基的动力响应问题. 通过傅里叶积分变换和Helmholtz矢量分解原理, 获得频域内非饱和土地基动力响应问题的通解, 结合回传射线矩阵法和边界条件, 求解获得了非均匀非饱和土层中位移、应力以及孔隙压力的计算列式. 假设沿深度方向梯度非均匀非饱和土的物理力学性质按幂函数连续变化, 通过数值傅里叶逆变换得到了非均匀非饱和土地基中的应力、位移以及孔隙压力等物理量的数值解, 分析讨论了土体非均匀性对非饱和土介质动力响应的影响规律. 结果表明: 土体非均匀性显著改变了非饱和土中竖向位移、正应力和孔隙压力在其深度方向上的振动模态, 其中孔隙气压在其深度方向的振动频率随着梯度因子的增加而不断增大, 波峰值不断靠近地表处附近; 竖向位移随着梯度因子的增大不断减小; 正应力和孔隙水压随着梯度因子的增大先增大后减小, 并且土体非均匀程度越高, 正应力与孔隙水压的幅值越大.      

A generalized self-consistent Mori-Tanaka model for predicting the effective moduli of the coated inclusion based composite materials

The weak point of the generalized self-consistent method (GSCM) is that its solution for the effective shear moduli involves determining the complicated displacement and strain fields in constitutents. Furthermore, the effective moduli estimated by GSCM cannot be expressed in an explicit form. Instead of following the procedure of GSCM, in this paper a generalized self-consistent Mori-Tanaka method (GSCMTM) is developed by means of Hill's interface condition and the assumption that the strain in the inclusion is uniform. A comparison with the existing theoretical and experimental results shows that the present GSCMTM is sufficiently accurate to predict the effective moduli of the coated inclusion-based composite materials. Moreover, it is interesting to find that the application of Hill's interface condition in volumetric domain is equivalent to the Mori-Tanaka average field approximation. This project was supported by the National Natural Science Foundation of China and China Postdoctoral Science Foundation.     

Multiscale diffusion method for simulations of long-time defect evolution with application to dislocation climb

In many problems of interest to materials scientists and engineers, the evolution of crystalline extended defects (dislocations, cracks, grain boundaries, interfaces, voids, precipitates) is controlled by the flow of point defects (interstitial/substitutional atoms and/or vacancies) through the crystal into the extended defect. Precise modeling of this behavior requires fully atomistic methods in and around the extended defect, but the flow of point defects entering the defect region can be treated by coarse-grained methods. Here, a multiscale algorithm is presented to provide this coupling. Specifically, direct accelerated molecular dynamics (AMD) of extended defect evolution is coupled to a diffusing point defect concentration field that captures the long spatial and temporal scales of point defect motion in the presence of the internal stress fields generated by the evolving defect. The algorithm is applied to study vacancy absorption into an edge dislocation in aluminum where vacancy accumulation in the core leads to nucleation of a double-jog that then operates as a sink for additional vacancies; this corresponds to the initial stages of dislocation climb modeled with explicit atomistic resolution. The method is general and so can be applied to many other problems associated with nucleation, growth, and reaction due to accumulation of point defects in crystalline materials.     

Influential factors of loess liquefaction and pore pressure development

The liquefaction of loess under dynamic loading is studied experimentally with a dynamic triaxial test system. The effects of over-consolidation ratio (OCR), saturation degree and the frequency of dynamic loading upon loess liquefaction are investigated. The development of pore pressure within loess samples is also discussed. Based on the experimental results, the empirical relationship between pore pressure ratio and loading cycle number ratio is established for normal consolidated saturated loess.The project supported by the National Natural Science Foundation of China (50178005)     

Mechanical characterization of the role of defects in sintered FeCrAlY foams

Open celled metal foams fabricated through metal sintering are a new class of material that offers novel mechanical and acoustic properties. Previously, polymer foams have been widely used as a means of absorbing acoustic energy. However, the structural applications of these foams are limited. The metal sintering approach offers a cost- effective means for the mass-production of open-cell foams from a range of materials, including high-temperature steel alloys. In this first part of two-paper series, the mechanical properties of open-celled steel alloy (FeCrAlY) foams were characterized under uniaxial compression and shear loading. Compared to predictions from established models, a significant knockdown in material properties was observed. This knockdown was attributed to the presence of defects throughout the microstructure that result from the unique fabrication process. Further in situ tests were carried out in a SEM (scanning electronic microscope) in order to investigate the effects of defects on the properties of the foams. Typically, the onset of plastic yielding was observed to occur at defect locations within the microstructure. At lower relative densities, ligament bending dominates, with the deformation initializing at defects. At higher relative densities, an additional deformation mechanism associated with membrane elements was observed. In the follow-up of this paper, a finite element model will be constructed to quantify the effects of defects on the mechanical performance of the open-cell foam. The project supported by the US Office of Naval Research (N000140210117), the National Basic Research Program of China (2006CB601202), the National Natural Science Foundation of China (10328203, 10572111, 10632060), and the National 111 Project of China (B06024).     

Special properties of Eshelby tensor for a regular polygonal inclusion

When studying the regular polygonal inclusion in 1997, Nozaki and Taya discovered numerically some remarkable properties of Eshelby tensor: Eshelby tensor at the center and the averaged Eshelby tensor over the inclusion domain are equal to that of a circular inclusion and independent of the orientation of the inclusion. Then Kawashita and Nozaki justified the properties mathematically. In the present paper, some other properties of a regular polygonal inclusion are discovered. We find that for an N-fold regular polygonal inclusion except for a square, the arithmetic mean of Eshelby tensors at N rotational symmetrical points in the inclusion is also equal to the Eshelby tensor for a circular inclusion and independent of the orientation of the inclusion. Furthermore, in two corollaries, we point out that Eshelby tensor at the center, the averaged Eshelby tensor over the inclusion domain, and the line integral average of Eshelby tensors along any concentric circle of the inclusion are all identical with the arithmetic mean.The project supported by the National Natural Science Foundation of China (10172003 and 10372003) The English text was polished by Keren Wang.     

The motion of a spherical particle in the stokes flow outside a circular orifice

For any study of a suspension entering a pore, the knowledge of the force and moment exerted on a solute particle in an arbitrary position outside the pore is essential. This paper for the first time presents approximate analytical expressions (in closed form) of all the twelve force and moment coefficients for a sphere outsied a circular orifice, on the basis of a number of discrete data computed by Yan et al. (1987). These coefficients are then applied to calculate the trajectory and angular velocity of a spherical particle approaching the pore at zero Reynolds number. The trajectory is in excellent agreement with the available experimental results. An analysis of the relative importance of the coefficients shows that the rotation effect cannot be neglected near the pore opening or near the wall, and that the lateral force effect must be taken into account in the neighborhood of the edge of the pore opening. It is due to neglecting these factors that previous theoretical results deviate from the experimental ones near the pore opening. The effects of the ratio of the particle to pore radii as well as the influences of the gravity/buoyance on the particle trajectory, velocity distribution and rotation are discussed in detail. It is pointed out that in the experiments of neutrally-buoyant suspensions, the restriction on the density of the particle is most demanding for a large particle size. The expressions of forces and moments presented herein are complete, relatively accurate and convenient, thus providing a good prerequisite for further studies of any problems involving the entrance of particles to a pore.Project supported by the National Natural Science Foundation of China.