Sliding frictional contact between a rigid punch and a laterally graded elastic medium

Analytical and computational methods are developed for contact mechanics analysis of functionally graded materials (FGMs) that possess elastic gradation in the lateral direction. In the analytical formulation, the problem of a laterally graded half-plane in sliding frictional contact with a rigid punch of an arbitrary profile is considered. The governing partial differential equations and the boundary conditions of the problem are satisfied through the use of Fourier transformation. The problem is then reduced to a singular integral equation of the second kind which is solved numerically by using an expansion–collocation technique. Computational studies of the sliding contact problems of laterally graded materials are conducted by means of the finite element method. In the finite element analyses, the laterally graded half-plane is discretized by quadratic finite elements for which the material parameters are specified at the centroids. Flat and triangular punch profiles are considered in the parametric analyses. The comparisons of the results generated by the analytical technique to those computed by the finite element method demonstrate the high level of accuracy attained by both methods. The presented numerical results illustrate the influences of the lateral nonhomogeneity and the coefficient of friction on the contact stresses.     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。     

A Newton method for three-dimensional fretting problems

The present paper concerns the numerical treatment of fretting problems using a finite element analysis. The governing equations resulting from a formal finite element discretization of an elastic body with a potential contact surface are considered in a quasi-static setting. The constitutive equations of the potential contact surface are Signorinis contact conditions, Coulombs law of friction and Archards law of wear. Using a backward Euler time discretization and an approach based on projections, the governing equations are written as an augmented Lagrangian formulation which is implemented and solved using a Newton algorithm for three-dimensional fretting problems of didactic nature. Details concerning the implementation are provided.     

线性区间有限元静力控制方程的组合解法

区间有限元的静力控制方程常被归结为区间方程组来求解。但实际上两者并不等价。本文根据不确定结构有限元分析的力学背景，直接从问题的基本参量的不确定性出发，将基本区间参量的边界组合与求解区间方程组的有关解法相结合，提出了线性区间有限元静力控制方程的两种组合解法－参量边界全组合法和组合迭代法。可以以较小的计算量获得或逼近位移和应力区间的准确界限。且不受基本参量变化范围的限制。算例分析表明文中方法是实用和可行的。     

Finite elasto-plastic deformation—I: Theory and numerical examples

It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate.Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admiting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.     

FUZZY ARITHMETIC AND SOLVING OF THE STATIC GOVERNING EQUATIONS OF FUZZY FINITE ELEMENT METHOD

The key component of finite element analysis of structures with fuzzy parameters, which is associated with handling of some fuzzy information and arithmetic relation of fuzzy variables , was the solving of the governing equations of fuzzy finite element method. Based on a given interval representation of fuzzy numbers, some arithmetic rules of fuzzy numbers and fuzzy variables were developed in terms of the properties of interval arithmetic. According to the rules and by the theory of interval finite element method, procedures for solving the static governing equations of fuzzy finite element method of structures were presented. By the proposed procedure, the possibility distributions of responses of fuzzy structures can be generated in terms of the membership functions of the input fuzzy numbers. It is shown by a numerical example that the computational burden of the presented procedures is low and easy to implement. The effectiveness and usefulness of the presented procedures are also illustrated.     

From Single-Phase to Compositional Flow: Applicability of Mixed Finite Elements

In this paper we discuss the formulation of the governing equations that describe flow of fluids in porous media. Various types of fluid flow, ranging from single-phase flow to compositional flow, are considered. It is shown that all the differential equations governing these types of flow can be effectively rewritten in a fractional flow formulation; i.e., in terms of a global pressure and saturation (or saturations), and that mixed finite element methods can be accurately exploited to solve the pressure equation. Numerical results are presented to see the performance of the mixed methods for the flow equations in three space dimensions.     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

摩擦接触弹塑性分析的数学规划法

本文提出摩擦接触弹塑性分析的数学规划法,是弹性接触问题解法的推广,使之能处理比例加载时的塑性效应。利用所导出接触问题总刚度阵的逆,结合参变量线性互补公式,仍可消去由接触单元而引入的惩罚因子。于是,可以获得一个弹塑性问题、弹性接触和弹塑性接触问题的有限元参数二次规划法模型。所选的几个实例说明了方法的实用性。     

Effect of imperfect bonding on the dynamic response of a pre-stressed sandwich plate-strip with elastic layers and a piezoelectric core

In this paper,we consider the dynamic response of a pre-stressed sandwich plate-strip with a piezoelectric core and elastic layers under the action of a time-harmonic force resting on a rigid foundation.The investigation is carried out within the framework of the piecewise homogeneous body model by utilizing the exact equations of motion and relations of the linear theory of electro-elasticity.It is assumed that there is a shear-spring-type imperfect contact between the layers,but a complete contact between the plate-strip and the rigid foundation.A mathematical model of the problem is constructed,and the governing equations of motion are solved by employing the finite element method(FEM).Numerical results illustrating the influence of a change in the value of the shear-spring parameter on the dynamic response of the plate-strip are then presented.     

自由单元法及其在结构分析中的应用

通过吸收有限元与无网格法的优点,提出了一种新的数值方法------自由单元法.此方法在离散方面,采用有限元法中的等参单元,表征几何形状和进行物理量的插值;在算法方面,采用单元配点技术,逐点产生系统方程.主要特点是,在每个配置点只需要一个和周围自由选择的节点而形成的一个独立的等参单元,因而不需要考虑物理量在单元之间的相互连接关系与导数连续性问题. 本文介绍强形式与弱形式两种自由单元法,前者直接由控制方程和边界条件直接产生系统方程,后者通过在自由单元上建立控制方程的加权余量式产生弱形式积分式,并通过像传统有限元法中的积分过程建立系统方程组.本文提出的方法是一种单元配点法,对于域内点为了获得较高的导数精度,需要采用至少具有一个内部点的等参单元,为此除了可使用各阶次的拉格朗日四边形单元外, 还 给出了七节点三角形等参单元,用于模拟较为复杂的几何形状问题.      

Damage analysis and dynamic response of elasto-plastic laminated composite shallow spherical shell under low velocity impact

Based on the elasto-plastic mechanics, the damage analysis and dynamic response of an elasto-plastic laminated composite shallow spherical shell under low velocity impact are carried out in this paper. Firstly, a yielding criterion related to spherical tensor of stress is proposed to model the mixed hardening orthotropic material, and accordingly an incremental elasto-plastic damage constitutive relation for the laminated shallow spherical shell is founded when a strain-based Hashin failure criterion is applied to assess the damage initiation and propagation. Secondly, using the presented constitutive relations and the classical nonlinear shell theory, a series of incremental nonlinear motion equations of orthotropic moderately thick laminated shallow spherical shell are obtained. The questions are solved by using the orthogonal collocation point method, Newmark method and iterative method synthetically. Finally, a modified elasto-plastic contact law is developed to determine the normal contact force and the effect of damage, geometrical parameters, elasto-plastic contact and boundary conditions on the contact force and the dynamic response of the structure under low velocity impact are investigated.     

Formulations of a hydromechanical interface element

A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.     

Applying the Virtual Fields Method to the identification of elasto-plastic constitutive parameters

This paper deals with the identification of constitutive parameters governing elasto-plastic constitutive equations. The procedure used is the virtual fields method. It consists in writing the principle of virtual work with particular virtual fields. A cost function is then deduced. Minimizing this cost function provides the unknown parameters. An important feature is the fact that no finite element analysis is required with this approach. Moreover, it is not sensitive to the distribution of the loading if suitable virtual fields are used. Some numerical examples illustrate the relevancy of the method and its sensitivity to noisy data.     

BUCKLING AND POSTBUCKLING ANALYSIS OF ELASTO-PLASTIC FIBER METAL LAMINATES

The elasto-plastic buckling and postbuckling of fiber metal laminates （FML） are studied in this research. Considering the geometric nonlinearity of the structure and the elasto- plastic deformation of the metal layers, the incremental Von Karman geometric relation of the FML with initial deflection is established. Moreover, an incremental elasto-plastic constitutive relation adopting the mixed hardening rule is introduced to depict the stress-strain relationship of the metal layers. Subsequently, the incremental nonlinear governing equations of the FML subjected to in-plane compressive loads are derived, and the whole problem is solved by the iterative method according to the finite difference method. In numerical examples, the effects of the initial deflection, the loading state, and the geometric parameters on the elasto-plastic buckling and postbuckling of FML are investigated, respectively.     

圆柱壳撞水时的弹塑性动力屈曲研究

本文研究圆柱壳在流固中冲击载荷下的弹塑性动力曲问题。建立了液体－气体－固体三相的数学模型。其中结构部分控制方程由弹塑性力学中关于加速度的最小原理获得，本构关系采用增量理论，液体假设为不可压缩，空气层认为最等熵压缩。分别讨论了不同冲击高度时压力变化规律，屈曲对载荷的影响。屈曲沿壳长的分布及发展规律等。     

A dynamic finite element method for inhomogeneous deformation and electromechanical instability of dielectric elastomer transducers

We present a three-dimensional nonlinear finite element formulation for dielectric elastomers. The mechanical and electrical governing equations are solved monolithically using an implicit time integrator, where the governing finite element equations are given for both static and dynamic cases. By accounting for inertial terms in conjunction with the Arruda–Boyce rubber hyperelastic constitutive model, we demonstrate the ability to capture the various modes of inhomogeneous deformation, including pull-in instability and wrinkling, that may result in dielectric elastomers that are subject to various forms of electrostatic loading. The formulation presented here forms the basis for needed computational tools that can elucidate the electromechanical behavior and properties of dielectric elastomers that are used for engineering applications.     

Shear buckling of infinite plates resting on tensionless elastic foundations

The buckling problem of an infinite thin plate resting on a tensionless Winkler foundation and subjected to shearing loads is investigated. The infinite plate is simplified to a one-dimensional mechanical model by assuming a lateral buckling mode function and a borderline function between contact and non-contact regions. After the governing differential equations for the plate sections in the contact and non-contact regions have been solved, the problem reduces to two nonlinear algebraic equations. Buckling coefficients for plates with simply supported edges and clamped edges are determined for a range of relative foundation stiffness factors. Comparison of the results with existing theory and finite element analyses shows good agreement.     

Solving material nonlinear structural problems by a finite analytic method in local coordinates

The purpose of this paper is to develop a finite analytic (FA) numerical solution for the elasto-plastic problem of the total theory. Schemes for the FA method in local coordinates for solving non-linear governing equations in the form of Navier equations are derived, which can be utilized to solve the problem in a domain of arbitrary geometry. Numerical illustration shows that the schemes are effective and practical.