Exact integration of the von Mises elastoplasticity model with combined linear isotropic-kinematic hardening

The main purpose of this work is to present two semi-analytical solutions for the von Mises elastoplasticity model governed by combined linear isotropic-kinematic hardening. The first solution (SOLε) corresponds to strain-driven problems with constant strain rate assumption, whereas the second one (SOLσ) is proposed for stress-driven problems using constant stress rate assumption. The formulas are derived within the small strain theory Besides the new analytical solutions, a new discretized integration scheme (AMε) based on the time-continuous SOLε is also presented and the corresponding algorithmically consistent tangent tensor is provided. A main advantage of the discretized stress updating algorithm is its accuracy; it renders the exact solution if constant strain rate is assumed during the strain increment, which is a commonly adopted assumption in the standard finite element calculations. The improved accuracy of the new method (AMε) compared with the well-known radial return method (RRM) is demonstrated by evaluating two simple examples characterized by generic nonlinear strain paths.     

考虑热渗效应的有限长圆柱热固结解

在作者已建立的饱和多孔介质耦合非线性热弹性理论基础上,考虑热渗效应,建立了饱和多孔介质耦合热弹性固结方程,并推导了有限长圆柱热固结问题的解析解,进而以温控三轴试验的试样为例进行了算例分析,同时利用COMSOL软件进行了数值模拟,并将解析结果和数值结果进行对比．结果表明：在不排水条件下,影响试样最终孔压大小的参数是：土的泊松比>弹性模量>水的体膨胀系数,渗透系数对孔压变化影响不大,考虑与不考虑水土压缩性不影响孔压的计算;相对线弹性情况,考虑非线弹性膨胀系数随温度变化时的孔压有所下降,轴向应变变化很小．     

Fibre-reinforced composites with fibre-bending stiffness under azimuthal shear – Comparison of simulation results with analytical solutions

In Ref. [1], Spencer and Soldatos proposed an enhanced modelling approach for fibre-reinforced composites which accounts for the fibre-bending stiffness in addition to the directional dependency induced by the fibres. Although analytical solutions for simple geometries have been derived over the past years, often subject to specific assumptions such as small deformation kinematics, the application to more general and non-academic boundary value problems is desirable. Motivated by the latter, the numerical solution of the general system of partial differential equations by means of a multi-field finite element approach is proposed in Ref. [2] and the principal model properties are studied for a specific form of the elastic energy potential. In the present contribution a comparison of the numerical solution by means of the multi-field finite element approach against the analytical solution is presented for the azimuthal shear deformation of a tube-like structure. To this end, the general deformation pattern and especially the distribution of the stress and couple stress tensor are taken into account. We find that, although the analytical solution is derived subject to the assumption of small deformations, whereas the numerical solution is based on the finite strain counterpart of the theory, the simulation results are quasi identical, which verifies the numerical framework proposed.     

一维大应变固结理论的一类解析解

简述了一维固结理论的发展，研究了不同描述方法下固结系数的定义，运用变换群的方法求得了某些条件下一维大应变固结方程的完整解析解答，在此基础上与相同条件下的小应变固结理论作了比较，并通过对实例的计算，分析了不同假设条件下一维固结理论之间的差异．计算结果表明，在固结过程中，由大应变固结理论所得出的沉降量要大于由小应变理论所得出的沉降量，而两种固结理论所得出的最终沉降量相同．     

A CONSOLIDATION MODEL OF GROUP WELLS ON COMPOSITE FOUNDATION1)

针对复合地基固结问题,考虑井阻和土体径竖向渗流的影响,建立 以"桩-土-桩"为单元的复合地基群井计算模型,采用解析解法,推导出瞬时载荷作用下复合地基群井固结解析解,讨论了复合地基固结特点和孔压等的变化规律。研究表明,复合地基群井固结模型为复合地基固结理论的研究提供了一种新思路,在井径比大于5 时,基于中心桩、边桩流量相等假设,复合地基群井固结理论所求固结度和传统固结理论所求固结度相差很小,可忽略不计;桩体渗透系数、内外扰动区土体径向渗透系数、土体竖向渗透系数越大,复合地基固结速率越快;堆载越大,最终达到稳定时的沉降量和竖向应变越大;地基表层土体的孔压消散速率和竖向应变速率大于地基深部土体。     

Exact three-dimensional analysis for static torsion of piezoelectric rods

Based on the theory of elasticity, exact analytical and numerical solutions of piezoelectric rods under static torsion are studied. In this paper, direct solution method is used. The main scope is to check the extension of validity of assumptions in previous papers that had been made based on linear distribution of electric potential through the cross section and their influences on deflection and the angle of rotation. Stress and electric induction functions are employed to obtain the exact solution of the static and electrostatic equilibrium equations under torsional loading. It is shown that previous assumptions are valid only in some types of piezoelectric materials, while in other types these assumptions lead to considerable deviations from accurate modeling. The present analytical solutions are compared with three-dimensional finite element analysis results and absolute agreements are found. At the end of this article, torsional rigidity, shape-effects on induced piezoelectric deformation and the range of valid region for linear distribution of electric potential assumption have been studied.     

Finite amplitude one-dimensional wave propagation in a thermoelastic half-space

The problem of finite wave propagation in a nonlinearly thermoelastic half-space is considered. The surface of the half-space is subjected to a time-dependent thermal and normal mechanical loading. The solution is obtained by a numerical procedure, which is shown to furnish accurate results, and linear dynamic thermoelastic problems are obtained as special cases. The accuracy of the results is checked by comparison with some known analytical solutions which can be obtained in some special cases of both the linear and the nonlinear problems. In those cases where the solution contains shocks, it is shown that the numerical results satisfy the necessary jumps conditions which need to hold across such discontinuities.     

功能梯度矩形板的近似理论与解析解

提出了压电功能梯度矩形板在竖向载荷作用下的近似理论与解析解. 引入了板理论的Kirchhoff假设、Reissner-Mindlin假设和提出的补充假设, 并假设材料常数在板厚方向按指数规律变化. 推导了板在周边简支同时又接地情况下中性层法线转角的解和用Fourier级数表示的电势解. 该解在形式上比精确解简单得多, 进行数值计算时也相当方便与快捷. 计算结果与ANSYS软件用三维实体单元的有限元计算结果进行了比较, 证实了该方法即使在厚板情况下仍然具有很高的精度.      

Analytical solution to one-dimensional consolidation in unsaturated soils

This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. The analytical solution is for Fredlund＇s one-dimensional consolidation equation in unsaturated soils. The transfer relationship between the state vectors at top surface and any depth is obtained by using the Laplace transform and Cayley-Hamilton mathematical methods to the governing equations of water and air, Darcy＇s law and Fick＇s law. Excess pore-air pressure, excess pore-water pressure and settlement in the Laplace-transformed domain are obtained by using the Laplace transform with the initial conditions and boundary conditions. By performing inverse Laplace transforms, the analytical solutions are obtained in the time domain. A typical example illustrates the consolidation characteristics of unsaturated soil from an- alytical results. Finally, comparisons between the analytical solutions and results of the finite difference method indicate that the analytical solution is correct.     

Mode I stress intensity factor solutions for spot welds in lap-shear specimens

The analytical solutions of the mode I stress intensity factor for spot welds in lap-shear specimens are investigated based on the classical Kirchhoff plate theory for linear elastic materials. First, closed-form solutions for an infinite plate containing a rigid inclusion under counter bending conditions are derived. The development of the closed-form solutions is then used as a guide to develop approximate closed-form solutions for a finite square plate containing a rigid inclusion under counter bending conditions. Based on the J integral, the closed-form solutions are used to develop the analytical solutions of the mode I stress intensity factor for spot welds in lap-shear specimens of large and finite sizes. The analytical solutions of the mode I stress intensity factor based on the solutions for infinite and finite square plates with an inclusion are compared with the results of the three-dimensional finite element computations of lap-shear specimens with various ratios of the specimen half width to the nugget radius. The results indicate that the mode I stress intensity factor solution based on the finite square plate model with an inclusion agrees well with the computational results for lap-shear specimens for the ratio of the half specimen width to the nugget radius between 4 and 15. Finally, a set of the closed-form stress intensity factor solutions for lap-shear specimens at the critical locations are proposed for future applications.     

关节软骨两相多孔介质非线性模型的有限元方法

基于混合物理论的两相多孔介质模型可以准确描述关节软骨的力学行为，关节软骨的渗透率与固体相体积应变相关。本文研究这一模型的非线性有限元法。具体采用伽辽金加权残值法得到有限元平衡方程，编制了有限元程序，进而对关节软骨围限压缩蠕变和应力松弛行为进行了数值模拟。与视渗透率为常数的线性模型的计算结果比较表明，在变形较大时，渗透率随固体相体积应变变化这一非线性效应不容忽视。     

Iterative method for large deflection nonlinear problem of laminated composite shallow shells and plates

ＩｎｔｒｏｄｕｃｔｉｏｎＡｓｃｏｍｐｏｓｉｔｅｓａｒｅｅｘｔｅｎｓｉｖｅｌｙａｐｐｌｉｅｄｉｎｍｏｄｅｒｎｅｎｇｉｎｅｒｉｎｇｓｔｒｕｃｔｕｒｅｓａｎｄｍｅｃｈａｎｉｃｓｏｆｃｏｍｐｌｅｘｍａｔｅｒｉａｌｓｈａｓｂｅｎｈｉｇｈｌｙａｄｖａｎｃｅｄ,ｉ...     

A new exact integration method for the Drucker–Prager elastoplastic model with linear isotropic hardening

This paper presents the exact stress solution of the non-associative Drucker–Prager elastoplastic model governed by linear isotropic hardening rule. The stress integration is performed under constant strain-rate assumption and the derived formulas are valid in the setting of small strain elastoplasticity theory. Based on the time-continuous stress solution, a complete discretized stress updating algorithm is also presented providing the solutions for the special cases when the initial stress state is located in the apex and when the increment produces a stress path through the apex. Explicit expression for the algorithmically consistent tangent tensor is also derived. In addition, a fully analytical strain solution is also derived for the stress-driven case using constant stress-rate assumption. In order to get a deeper understanding of the features of these solutions, two numerical test examples are also presented.     

Plane strain consolidation of soil layer with anisotropic permeability

This paper presents an alternative analytical technique to study a plane strain consolidation of a poroelastic soil by taking into account the anisotropy of permeability. From the governing equations of a saturated poroelastic soil, the relationship of basic variables for a point of a soil layer is established between the ground surface （z=0） and the depth z in the Laplace-Fourier transform domain. Combined with the boundary conditions, an exact solution is derived for plane strain Biot＇s consolidation of a finite soil layer with anisotropic permeability in the transform domain. Numerical inversions of the Laplace transform and the Fourier transform are adopted to obtain the actual solution in the physical domain. Numerical results of plane strain Biot＇s consolidation for a single soil layer show that the anisotropic of permeability has a great influence on the consolidation behavior of the soils.     

Analytical and numerical solutions for vacuum preloading considering a radius related strain distribution

A system of PVDs combined with other preloading methods such as vacuum preloading and surcharge preloading is an effective and economical method which is widely used in the ground treatment. The consolidation theories for drain wells under equal strain condition are often used in the design of ground treatment by PVDs. A radius related strain distribution is proposed to get analytical solutions for the excess pore water pressure and settlement. A linear distributed vacuum pressure along the drain depth and the smear zone as well as well resistance are considered. The numerical results for vacuum loading process are obtained by developed FEM model to compare with the analytical solutions. The results indicate that the influence of the equal strain hypothesis cannot be neglected when n is larger than 10, where n denotes the ratio of diameter of the model to diameter of the drain. The analytical solutions proposed in this paper are more consistent with the numerical results than the analytical results obtained employing the equal strain condition.     

欧拉描述的大变形固结理论

以往大变形固结理论主要基于一般的固体力学模型，其控制方程忽视了固结过程中排水引起 的质量变化. 提出饱和土的连续介质模型，并基于连续介质力学的公理体系推导了反映质量 变化的欧拉描述的大变形固结控制方程. 发现传统固结理论中：(1)忽视了渗流速度对土体平衡条件的影响；(2)决定土体平衡的总应力张量只有在土体变形速度和渗流速度方向相同时才具有对称性等. 在忽略变质量效应等条件下，传统理论成为本文理论的特例. 通过算例 的有限元分析，比较了欧拉描述与两种物质描述方法的差别，得到初步结论：(1)欧拉描述 方法计算的地基沉降量要小于物质描述方法的结果；(2)欧拉描述方法计算的侧向位移偏大 于两种物质描述结果.     

基于能量等效原理的应变局部化分析:Ⅱ.有限元解法

以非局部塑性理论为基础,应用状态空间理论,通过局部和非局部两个状态空间的塑性能量耗散率等效原理,提出了一种求解应变局部化问题的新方法,以得到与网格无关的数值解.针对二维问题的屈服函数和流动法则导出了求解非局部内变量的一般方程,并提出了在有限元环境中求解应变局部化问题的应力更新算法.为了验证所提出的方法,对1个一维拉杆和3个二维平面应变加载试件进行了有限元分析.数值结果表明,塑性应变的分布和载荷-位移曲线都随着网格的变小而稳定地收敛,应变局部化区域的尺寸只与材料内尺度有关,而对有限元网格的大小不敏感.对于一维问题,当有限元网格尺寸减小时,数值解收敛于解析解.对于二维剪切带局部化问题,数值解随着网格尺寸的减小而稳定地向唯一解收敛.当网格尺寸减小时,剪切带的宽度和方向基本上没有变化.而且得到的塑性应变分布和网格变形是平滑的.这说明,所提方法可以克服经典连续介质力学模型导致的网格相关性问题,从而获得具有物理意义的客观解.此模型只需要单元之间的位移插值函数具有C~0连续性,因而容易在现有的有限元程序中实现而无需对程序作大的修改.     

Analytical solutions to general anti-plane shear problems in finite elasticity

This paper presents a pure complementary energy variational method for solving a general anti-plane shear problem in finite elasticity. Based on the canonical duality–triality theory developed by the author, the nonlinear/nonconvex partial differential equations for the large deformation problem are converted into an algebraic equation in dual space, which can, in principle, be solved to obtain a complete set of stress solutions. Therefore, a general analytical solution form of the deformation is obtained subjected to a compatibility condition. Applications are illustrated by examples with both convex and nonconvex stored strain energies governed by quadratic-exponential and power-law material models, respectively. Results show that the nonconvex variational problem could have multiple solutions at each material point, the complementary gap function and the triality theory can be used to identify both global and local extremal solutions, while the popular convexity conditions (including rank-one condition) provide mainly local minimal criteria and the Legendre–Hadamard condition (i.e., the so-called strong ellipticity condition) does not guarantee uniqueness of solutions. This paper demonstrates again that the pure complementary energy principle and the triality theory play important roles in finite deformation theory and nonconvex analysis.     

用DQM求解成层地基一维非线性固结问题

由于多层地基的一维非线性固结问题求解的复杂性,其解析解很难求得。本文基于Davis和Raymond一维非线性固结理论,利用DQM(Differential Quadrature Method)导了初始有效应力沿深度变化、任意边界条件、任意荷载作用下成层地基一维非线性固结的统一表达式,求得了孔压、有效应力和平均固结度的解答。通过解的收敛性分析讨论了DQM解的有效性。由于DQM解对于固结间题各种复杂条件具有统一的矩阵表达式,更便于编程计算和工程应用。最后,用本文解答对三层地基一维非线性固结问题进行了讨论。     

CONSOLIDATION THEORY OF UNSATURATED SOIL BASED ON THE THEORY OF MIXTURE(Ⅱ)

The present paper uses the mathematics model for consolidation of unsaturatedsoil developed in ref.[1]to solve boundary value problems.The analytical solutionsfor one-dimensional consolidation problem are gained by making use of Laplacetransform and finite Fourier transform.The displacement and the pore water pressureas well as the pore gas pressure are found from governing equations simultaneously.The theoretical formulae of coefficient and degree of consolidation are also given inthe paper.With the help of the method of Galerkin Weighted Residuals,the finiteelement equations for two-dimensional consolidation problem are derived.A FORTRANprogram named CSU8 using8-node isoparameter element is designed.A plane strainconsolidation problem is solved using the program,and some distinguishing features onconsolidation of unsaturated soil and certain peculiarities on numerical analysis arerevealed.These achievements make it convenient to apply the theory proposed by theauthor in engineering practice.