土与结构接触面弹塑性损伤模型用于单桩与地基相互作用分析

基于粗粒土与结构接触面弹塑性损伤静动力统一模型(称作EPDI模型)建立了可用于有限元分析的弹塑性损伤接触面单元。对接触面试验进行了模拟,采用不同的接触面本构模型及参数对单调和循环荷载作用下的单桩基础的侧摩阻力和桩顶位移进行了有限元分析。结果表明:包括剪应力应变关系和剪胀特性在内的接触面力学特性对桩土相互作用分析有重要影响,需要合理地加以描述。基于试验结果建立的弹塑性损伤接触面单元能够有效地用于土体与结构物相互作用分析,并能够合理地反映土与结构接触面的包括体应变及其与剪应变耦合特性在内的接触面主要静动力学特性。     

Transient dynamic crack analysis in piecewise homogeneous, anisotropic and linear elastic composites by a spatial symmetric Galerkin-BEM

Transient elastodynamic analysis of two-dimensional, piecewise homogeneous, anisotropic and linear elastic solids containing interior and interface cracks is presented in this paper. To solve the initial boundary value problem, a spatial symmetric time-domain boundary element method is developed. Stationary cracks subjected to impact loading conditions are considered. Elastodynamic fundamental solutions for homogenous, anisotropic and linear elastic solids are implemented. The piecewise homogeneous, anisotropic and linear elastic solids are modeled by the multi-domain technique. The spatial discretization is performed by a symmetric Galerkin-method, while a collocation method is utilized for the temporal discretization. An explicit time-stepping scheme is obtained for computing the unknown boundary data. Numerical examples are presented and discussed to show the effects of the interface cracks, the material anisotropy, the material combination and the dynamic loading on the dynamic stress intensity factors.     

SSI analysis of a massive concrete structure based on a novel convolution/deconvolution technique

Investigations conducted based on seismic soil-structure interaction analysis of a massive concrete structure supported on a raft foundation are presented in this paper. Linear transient dynamic analysis is carried out using finite element method and imposing transmitting boundary conditions at far field of layered elastic half-space. Analysis is conducted in two phases, namely: (i) free-field analysis of the layered half-space and (ii) seismic analysis of the structure by including soilstructure interaction effects. In the first phase, a simple and novel technique is used to establish free-field excitation at a depth in the half-space. In the second phase, seismic soil-structure interaction analysis of the structure is carried out for the free-field excitation determined in phase-I. Stress resultants experienced by the raft and the stresses at the interface between the rock and raft are evaluated. Critical examination of the results indicates tensile stresses of considerable magnitude at few locations in the rock-raft interface. Typical stress responses at the interface are presented and discussed in the paper.     

流固耦合问题的描述方法及分类简化准则

在流固耦合的工程实际问题中的大多数情况下,弹性薄壁构件的变形为几何非线性,再加上流体方程的非线性,将导致流体和弹性体相互作用界面上的强非线性。在界面上便可以结合拉格朗日法和欧拉法建立方程和接触条件。其方法主要有单一拉格朗日法、单一欧拉法、相容拉格朗日-欧拉法和任意拉格朗日-欧拉法四种方法。引入描述弹性体变形特征的数值m、n、k和描述流体运动特征的数值λ、ν,可将流固耦合问题进行分类。在流体弹性力学理论的基础上,介绍了流固耦合问题界面相互作用的描述方法,并根据诺沃日洛夫BB(НовожиловВВ)在非线性弹性力学中,从几何非线性方面对弹性力学问题的分类方法出发,对流体弹性力学中的流固耦合问题进行分类,由此,可为按类别对运动学条件、动力学条件及界面上的接触条件进行相应的简化提供可靠的依据。     

Quantifying moisture damage at mastic–aggregate interface

This study quantifies moisture damage at the asphalt mastic–aggregate interfaces using finite element method modelling technique in ABAQUS. A model aggregate surrounded by a layer of mastic is subjected to static loads of varying magnitudes and patterns. Using dynamic shear and elastic moduli of wet and dry mastic and aggregate as model inputs, moisture-induced damage is quantified through parameters such as contact stresses at interface, load to initiate damage and de-bonding at the interface. Results show that contact stresses are significantly higher in dry samples than wet samples. It is revealed that damage initiates at a smaller magnitude of deformation (0.0508 mm) in the wet sample than that (0.508 mm) in the dry sample. That is, a stiffer dry sample carries more loads and deforms less before damage initiation than a softer wet sample. In addition, approximately 6.8% interface de-bonding occurs in the dry sample, whereas 49.1% de-bonding occurs in the wet sample.     

Thermal contact conductance of pressed contacts at low temperatures

The influence of variations of interface temperature in the range 50-300 K on the thermal contact conductance between aluminium and stainless steel joints was determined. Predictions were done by modeling the deformation at the interface for different values of surface finish and contact pressure over the range of interface temperatures. Both elastic and plastic deformation was considered. Experiments were carried out in a closed loop cryostat and the results were shown to compare well with the predictions. A reduction of the interface temperature resulted in a smaller value of thermal contact conductance. Interfacial pressure variation had much lower influence at the smaller value of temperatures. The role of surface roughness at the contact was also seen to be less significant at lower interface temperatures and the zone of hysteresis was smaller. A correlation was developed for estimating thermal contact conductance at joints over this temperature range. An explicit dependence of contact conductance on temperature was not seen to be necessary as long as the changes in the hardness and thermal conductivity of the material with temperature are incorporated in the correlation.     

Transformation of sliding motion to rolling during spheres collision

In this research, we have investigated the three-dimensional elastic collision of two balls, based on friction in the tangential plane. Our aim is to offer analytical closed form relations for post collision parameters such as linear and angular velocities, collision time and tangential and normal impulse in three dimensions. To simplify the problem, stick regime is not considered. In other words, balls have a low tangential coefficient of restitution. Sliding, sliding then rolling, and rolling at the beginning of contact are three cases that can occur during impact which have been considered in our research. The normal interaction force is described by the Hertz contact force and dimensionless analysis is used for investigating normal interaction force; furthermore, Coulomb friction is considered during sliding. Experimental data for collisions show when sliding exists through the impact, tangential impulses can be taken as frictional impulses using the Coulomb law if the dynamic regime is not stick regime. To identify transformation of sliding motion to rolling or sticking during the impact process, linear and trigonometric functions are considered as an approximation for the normal interaction force. Afterwards, we have obtained the condition for the possibility of this transformation; moreover, we can estimate the duration of sliding and rolling or sticking. We have obtained an analytical solution for maximum force and deformation, collision time, impulses and post-collision linear and angular velocities in three dimensions.     

粘弹性饱和土中深埋圆形隧道衬砌-土相互作用

基于Biot理论,采用渗流-力学耦合模型研究分析了粘弹性饱和土体中深埋圆形隧道衬砌-土相互作用问题,并考虑了衬砌材料的徐变特性。假定衬砌和土体完全接触,运用积分变换理论在Laplace变换域中得到了衬砌和土体接触面上的应力、变形及衬砌的薄膜力和弯矩。利用Laplace数值逆变换得到时域中的解,并由此分析讨论了衬砌和土体的相对刚度、阻尼比、以及衬砌的相对厚度对隧道衬砌与土体相互作用的影响。数值结果同时也反映了应力、变形、薄膜力和弯矩的随时间变化过程。     

Adaptive methods for thermomechanical coupled contact problems

In this paper an adaptive method for the analysis of thermomechanical coupled multi‐body contact problems is presented. The method is applied to non‐linear elastic solids undergoing finite (thermal) deformations. The contact model considers non‐linear pressure‐dependent heat flux as well as frictional heating in the interface. A time–space‐finite element discretization of the governing equations is formulated including unilateral constraints due to contact. A staggered solution algorithm has been constructed that allows an independent spatial discretization of the coupled subproblems. A posteriori projection‐based error estimators, which enforce implicitly the special boundary conditions due to thermal contact, are used to control the spatial discretization as well as the adaptive time stepping. Numerical examples are presented to corroborate the applicability of the adaptive algorithm to the considered problem type. Copyright © 2004 John Wiley & Sons, Ltd.     

土-结构三维动力分析的线性-非线性混合子结构法

为提高土-结构三维非线性动力分析问题的计算效率,根据结构存在局部塑性区域的特点,提出了适合于非线性子结构与多个线性子结构存在边界耦合情况的线性-非线性混合约束模态综合法。该方法将整体体系划分为线性和非线性两类子结构,对线性子结构依照势能判据截断准则进行自由度的缩减,并最终与非线性子结构进行综合,进而获得非线性体系的动态响应。进一步将该方法应用到土-高层建筑相互作用体系的三维动力分析算例中,针对两种不同的土体边界条件分别进行了多组地震动激励作用下的动力时程分析,验证了该方法的精确性与有效性。最后,针对非线性土体区域范围的多种不同取值进行了多组数值试验,结果表明:对于该文算例,当非线性土体区域的范围取为上部结构水平尺寸的3倍时即可获得满意的计算结果,进一步说明了该方法的有效性。     

Transient analysis of three-dimensional dynamic interaction between multilayered soil and rigid foundation

The study of dynamic soil-structure interaction is significant to civil engineering applications, such as machine foundation vibration, traffic-induced vibration, and seismic dynamic response. The scaled boundary finite element method (SBFEM) is a semi-analytical algorithm, which is used to solve the dynamic response of a three-dimensional infinite soil. It can automatically satisfy the radiation boundary condition at infinity. Based on the dynamic stiffness matrix equation obtained by the modified SBFEM, a continued fraction algorithm is proposed to solve the dynamic stiffness matrix of layered soil in the frequency-domain. Then, the SBFEM was coupled with the finite element method (FEM) at the interface to solve the dynamic stiffness matrices of the rigid surface/buried foundation. Finally, the mixed-variable algorithm was used to solve the three-dimensional transient dynamic response of the foundation in the time domain. Numerical examples were performed to verify the accuracy of the proposed algorithm in solving the dynamic stiffness matrix of the infinite domain in the frequency domain and the dynamic transient displacement response of the foundation in the time domain. Compared with the previous numerical integration technique, the dynamic stiffness matrix in the frequency domain calculated by using the proposed algorithm has higher accuracy and higher efficiency.     

Variational approach to the analysis of steady‐state thermo‐elastic wear regimes

A transient wear process on frictional interface of two thermo‐elastic bodies in a relative steady sliding motion induces shape evolution of contact interface and tends to a steady state for which the wear process occurs at fixed contact stress and strain distribution. The temperature field generated by frictional and wear dissipation on the contact surface is assumed to reach a steady state. This state is assumed to correspond to minimum of the wear dissipation power and the temperature field corresponds to maximum of the heat entropy production. The stationarity conditions of the response functionals provide the contact pressure distribution and the corresponding temperature field. The present approach extends the authors previous analyses of optimal or steady‐state contact shapes by accounting for coupled wear and thermal distortion effects. The wear rule is assumed as a non‐linear relation of wear rate to shear stress and relative sliding velocity. The analysis of disk and drum brakes is presented with account for thermal distortion effect. It is shown that the contact shape in a steady thermo‐elastic state essentially differs from that specified for mechanical loading with neglect of thermal effects. The thermal instability regimes are not considered in the paper. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of 3D time domain boundary element formulation to wave propagation in poroelastic solids

The dynamic responses of fluid-saturated semi-infinite porous continua to transient excitations such as seismic waves or ground vibrations are important in the design of soil-structure systems. Biot's theory of porous media governs the wave propagation in a porous elastic solid infiltrated with fluid. The significant difference to an elastic solid is the appearance of the so-called slow compressional wave. The most powerful methodology to tackle wave propagation in a semi-infinite homogeneous poroelastic domain is the boundary element method (BEM). To model the dynamic behavior of a poroelastic material in the time domain, the time domain fundamental solution is needed. Such solution however does not exist in closed form. The recently developed ‘convolution quadrature method’, proposed by Lubich, utilizes the existing Laplace transformed fundamental solution and makes it possible to work in the time domain. Hence, applying this quadrature formula to the time dependent boundary integral equation, a time-stepping procedure is obtained based only on the Laplace domain fundamental solution and a linear multistep method. Finally, two examples show both the accuracy of the proposed time-stepping procedure and the appearance of the slow compressional wave, additionally to the other waves known from elastodynamics.     

3D non-linear time domain FEM–BEM approach to soil–structure interaction problems

Dynamic soil–structure interaction is concerned with the study of structures supported on flexible soils and subjected to dynamic actions. Methods combining the finite element method (FEM) and the boundary element method (BEM) are well suited to address dynamic soil–structure interaction problems. Hence, FEM–BEM models have been widely used. However, non-linear contact conditions and non-linear behavior of the structures have not usually been considered in the analyses. This paper presents a 3D non-linear time domain FEM–BEM numerical model designed to address soil–structure interaction problems. The BEM formulation, based on element subdivision and the constant velocity approach, was improved by using interpolation matrices. The FEM approach was based on implicit Green's functions and non-linear contact was considered at the FEM–BEM interface. Two engineering problems were studied with the proposed methodology: the propagation of waves in an elastic foundation and the dynamic response of a structure to an incident wave field.     

土与结构接触面土体软/硬化本构模型及数值实现

基于土与结构接触面变形特性分析,将接触面土体的剪切滑动面与单元体三维应力状态下的八面体面相对应,通过土的三维弹塑性本构模型在八面体面上的剪切应力-应变关系,建立了接触面土体剪切应力-应变关系;将接触面土体法向压缩变形与侧限压缩条件相对应,通过侧限压缩条件下的荷载变形关系,建立了接触面土体法向应力-应变关系;进一步将接触面土体切向与法向耦合,建立了接触面土体本构模型,模型只有4个材料参数,参数物理意义明确,可由等向压缩试验和常规三轴压缩试验确定。与接触面土体试验结果的对比分析表明,所建立的本构模型可较好地描述接触面土体切向软/硬化特性与法向变形规律。结合有限元软件ABAQUS,编制了FRIC模型子程序,通过模拟土与结构界面剪切滑移过程表明,编制的FRIC子程序可较好地模拟土与结构界面接触的非线性力学行为。     

Fractional-derivative viscoelastic model of the shock interaction of a rigid body with a plate

The impact of a rigid body upon an elastic isotropic plate is investigated for the case when the equations of motion take rotary inertia and shear deformation into account. The impactor is considered as a mass point, and the contact between it and the plate is established through a buffer involving a linear-spring–fractional-derivative dashpot combination, i.e., the viscoelastic features of the buffer are described by the fractional-derivative Maxwell model. It is assumed that a transient wave of transverse shear is generated in the plate, and that the reflected wave has insufficient time to return to the location of the spring’s contact with the plate before the impact process is completed. To determine the desired values behind the transverse-shear wave front, one-term ray expansions are used, as well as the equations of motion of the impactor and the contact region. As a result, we are led to a set of two linear differential equations for the displacements of the spring’s upper and lower points. The solution of these equations is found analytically by the Laplace-transform method, and the time-dependence of the contact force is obtained. Numerical analysis shows that the maximum of the contact force increases, tending to the maximal contact force when the fractional parameter is equal to unity.     

土-结构相互作用对结构主动控制的影响研究

利用弹性半空间假定,推导了结构考虑土与结构相互作用(SSI)体系的主动控制方程.分析了土与结构相互作用对结构主动控制的影响.并分析研究了土及结构各种主要参数对结构控制效果的影响.研究结果表明,土与结构相互作用对结构主动控制效果有一定影响;地基土的性质、结构刚度、基础埋深等在一定范围内对结构的主动控制有着显著的影响.     

Exterior interface cracks

Continuing a recent investigation of interface cracks, attention is paid to the exterior crack. Two elastic solids bonded over a finite segment of their boundary and capable of transmitting shear and tensile tractions are considered. It is found that one of the edge cracks remains completely closed under shear alone, and opens gradually as the level of tension is increased. Both crack tips, however, must remain closed at least over a small interval. Stress intensity factors and bond and contact stresses are given for a specific example.     

Dynamic analysis of interfacial crack problems in anisotropic bi-materials by a time-domain BEM

A time-domain boundary element method (BEM) together with the sub-domain technique is applied to study dynamic interfacial crack problems in two-dimensional (2D), piecewise homogeneous, anisotropic and linear elastic bi-materials. The bi-material system is divided into two homogeneous sub-domains along the interface and the traditional displacement boundary integral equations (BIEs) are applied on the boundary of each sub-domain. The present time-domain BEM uses a quadrature formula for the temporal discretization to approximate the convolution integrals and a collocation method for the spatial discretization. Quadratic quarter-point elements are implemented at the tips of the interface cracks. A displacement extrapolation technique is used to determine the complex dynamic stress intensity factors (SIFs). Numerical examples for computing the complex dynamic SIFs are presented and discussed to demonstrate the accuracy and the efficiency of the present time-domain BEM.     

Level set topology optimization of structural problems with interface cohesion

This paper presents a finite element topology optimization framework for the design of two‐phase structural systems considering contact and cohesion phenomena along the interface. The geometry of the material interface is described by an explicit level set method, and the structural response is predicted by the extended finite element method. In this work, the interface condition is described by a bilinear cohesive zone model on the basis of the traction‐separation constitutive relation. The non‐penetration condition in the presence of compressive interface forces is enforced by a stabilized Lagrange multiplier method. The mechanical model assumes a linear elastic isotropic material, infinitesimal strain theory, and a quasi‐static response. The optimization problem is solved by a nonlinear programming method, and the design sensitivities are computed by the adjoint method. The performance of the presented method is evaluated by 2D and 3D numerical examples. The results obtained from topology optimization reveal distinct design characteristics for the various interface phenomena considered. In addition, 3D examples demonstrate optimal geometries that cannot be fully captured by reduced dimensionality. The optimization framework presented is limited to two‐phase structural systems where the material interface is coincident in the undeformed configuration, and to structural responses that remain valid considering small strain kinematics. Copyright © 2017 John Wiley & Sons, Ltd.