Incremental deformation analysis of shell and corrugated diaphragm based on arbitrary configuration

With respect to an arbitrary configuration of a deformed structure, two sets of incremental equations are proposed for the deformation analysis of revolution shells and diaphragms loaded by both lateral pressures and the initial stresses produced in manufacturing. These general equations can be reduced to the simplified Koiter's Reissner-Meissner-Reissner (RMR) equations and the simplified Reissner's equations, when the initial stresses are set to zero. They can also be deduced to the total Lagrange form or the updated Lagrange form, respectively, as the structure is specified as the un-deformed or the former-deformed configurations. These incremental equations can be easily transformed into finite difference forms and solved by common numerical solvers of ordinary differential equations. Some numerical examples are presented to show the applications of the incremental equations to the deep shell of revolution and the corrugated diaphragms used in microelectronical mechanical system (MEMS). The results are in good agreement with those from finite element method (FEM). The project supported by the National Natural Science Foundation of China (10125211) and the 973 Program (G1999033108) The English text was polished by Keren Wang.     

Development of a Generalized Version of the Poisson– Nernst–Planck Equations Using the Hybrid Mixture Theory: Presentation of 2D Numerical Examples

A numerical scheme for the transient solution of a generalized version of the Poisson–Nernst–Planck (PNP) equations is presented. The finite element method is used to establish the coupled non-linear matrix system of equations capable of solving the present problem iteratively. The PNP equations represent a set of diffusion equations for charged species, i.e. dissolved ions, present in the pore solution of a rigid porous material in which the surface charge can be assumed neglectable. These equations are coupled to the ‘internally’ induced electrical field and to the velocity field of the fluid. The Nernst–Planck equations describing the diffusion of the ionic species and Gauss’ law in use are, however, coupled in both directions. The governing set of equations is derived from a simplified version of the so-called hybrid mixture theory (HMT). The simplifications used here mainly concerns ignoring the deformation and stresses in the porous material in which the ionic diffusion occurs. The HMT is a special version of the more ‘classical’ continuum mixture theories in the sense that it works with averaged equations at macroscale and that it includes the volume fractions of phases in its structure. The background to the PNP equations can by the HMT approach be described by using the postulates of mass conservation of constituents together with Gauss’ law used together with consistent constitutive laws. The HMT theory includes the constituent forms of the quasistatic version of Maxwell’s equations making it suitable for analyses of the kind addressed in this work. Within the framework of HTM, constitutive equations have been derived using the postulate of entropy inequality together with the technique of identifying properties by Lagrange multipliers. These results will be used in obtaining a closed set of equations for the present problem.     

An enhanced single-layer variational formulation for the effect of transverse shear on laminated orthotropic plates

A novel mixed formulation is derived by means of Reissner's variational approach-based on Castigliano's principle of least work in conjunction with a Lagrange multiplier method for the calculus of variations. The governing equations present an alternative theory for modeling the important three-dimensional structural aspects of plates in a two-dimensional form. By integrating the classical Cauchy's equilibrium equations with respect to the thickness co-ordinate, and enforcing continuity of shear and normal stresses at each ply interface, condenses the effect of the thickness. A reduced system of partial differential equations of sixth-order in one variable, is also proposed, which contains differential correction factors that formally modify the classical constitutive equations for composite laminates. The theory degenerates to classical composite plate analysis for thin configurations. Significant deviations from classical plate theory are observed when the thickness becomes comparable with the in-plane dimensions. A variety of case studies are presented and solutions are compared with other models available in the literature and with finite element analysis.     

浅曲梁大挠度分析的余能有限元方法

以高玉臣提出的弹性大变形余能原理为基础,利用Lagrange乘子,放松平衡方程和力边界条件对余能泛函的约束,推导出广义的余能原理.根据极分解定理,将变形分为刚性转动和纯变形两部分,则余能也包含相应的两部分,一部分与刚性转动有关,而另一部分与纯变形有关.使用线弹性本构关系,建立了可用于几何非线性计算的有限元模型.应用更新的Lagrange列式法,给出了增量形式的有限元公式.数值计算结果表明,该方法可用于浅曲粱的几何大变形计算.     

Buckling and large deflection behaviors of radially functionally graded ring-stiffened circular plates with various boundary conditions

The buckling and large deflection behaviors of axis-symmetric radially functionally graded (RFG) ring-stiffened circular plates are investigated by the dynamic relaxation (DR) method combined with the finite difference discretization technique. The material properties of the constituent components of the RFG plate are assumed to vary continuously according to the Mori-Tanaka distribution along the radial direction. The nonlinear governing equations are obtained in the incremental form based on the first-order shear deformation plate theory (FSDT) and the von Karman relations for large deflection. In the buckling analysis, an external in-plane load is applied to the plate incrementally so that, in each load-step, the incremental form of the governing equations can be solved by a numerical code prepared based on the DR method. After converging the DR code in the first increment, the latter load-step is added to the previous one, and the program is repeated again. The critical buckling load is determined from the compressive load-displacement curve obtained by solving the incremental form of the governing equa- tions. Based on the present incremental form of formulation, a bending analysis can also be conducted if the whole load is applied simultaneously. Finally, a detailed parametric study is carried out to investigate the influences of various boundary conditions, grading indices, thickness-to-radius ratios, stiffener’s positions and depths on the critical buckling load, and displacements and stresses resulted from the bending analysis. It is observed that the effect of the stiffener on the results is much greater in the functionally graded plate with higher material grading indices. The results also reveal that, by increasing the depth of the stiffer, the values of ascending the critical buckling load are approximately identical for both simply supported and clamped boundary conditions.     

Linear and higher order displacement theories for adhesively bonded lap joints

This work presents an adhesive model for stress analysis of bonded lap joints, which can be applied to model thin and thick adhesive layers. In this theory, linear variations of displacement components along the adhesive thickness are firstly assumed, and the longitudinal strain and the Poisson's effect of the adhesive are modeled. A differential form of the equilibrium equations for the adherends is analytically solved by means of compatible relations of the adhesive deformation. The derived shear and peel stresses are compared with the classical adhesive model of continuous springs with constant shear and peel stresses, and validated with two-dimensional finite element results of the geometrically nonlinear analysis using a commercial package. The numerical results show that the present linear displacement theory can be applied to both thin and moderately thick adhesive layers. The present formulation of the linear displacement theory is then extended to the higher order displacement theory for stress analysis of a thick adhesive, whose numerical results are also compared with those of the finite element computation.     

A class of controllable wave propagations in finite elasticity

Governing equations of axisymmetric finite dynamic deformations of an incompressible, isotropic and elastic cylindrical shell made of Neo-Hookean materials are derived. The non-linear partial differential equations are simplified for the cases where all deformation variations along the thickness of the tube may be neglected. The simplified non-linear equations are then solved exactly to arrive at traveling wave solutions along the axis. These wave solutions are called controllable because they can be maintained by prescribable surface stresses, bounded amplitude and frequency of excitations alone.     

空间刚性杆——弹簧组合结构轨道、姿态耦合动力学分析

以空间太阳帆塔在轨运行中遇到的强耦合动力学问题为研究背景,建立了空间刚性杆-- 弹簧组合结构轨道与姿态耦合 问题的动力学模型,采用辛 (几何) 算法研究了其轨道与姿态耦合的动力学行为,研究结果可以从系统的能量保持情况间接得到验 证. 首先,基于变分原理,通过引入对偶变量将描述空间刚性杆-- 弹簧组合结构动力学行为的拉格朗日方程导入哈 密尔顿体系,建立简化模型的正则控制方程;随后,采用辛龙格库塔方法模拟分析了地球非球摄动对轨道、姿态的影响及系统能 量的数值偏差问题. 数值模拟结果显示:随着初始姿态角速度增大,轨道半径的扰动 增大,轨道与姿态之间的耦合效应加剧; 带谐摄动对空间刚性杆-- 弹簧组合结构模型的轨道、姿态产生的影响比田谐摄动要高出至少两个数量级;同时辛龙格库塔方法能更好 地快速模拟地球非球摄动影响下空间刚性杆-- 弹簧组合结构的动力学行为,并能够长时间保持系统的总能量,有望为 超大空间结构实时反馈控制提供实时动力学响应结果.      

Elasto-plastic postbuckling of damaged orthotropic plates

Based on the elasto-plastic mechanics and continuum damage theory, a yield criterion related to spherical tensor of stress is proposed to describe the mixed hardening of damaged orthotropic materials. Its dimensionless form is isomorphic with the Mises criterion for isotropic materials. Furthermore, the incremental elasto-plastic damage constitutive equations and damage evolution equations are established. Based on the classical nonlinear plate theory, the incremental nonlinear equilibrium equations of orthotropic thin plates considering damage effect are obtained, and solved with the finite difference and iteration methods. In the numerical examples, the effects of damage evolution and initial deflection on the elasto-plastic postbuckling of orthotropic plates are discussed in detail.     

Finite element implementation of a generalised non-local rate-dependent crystallographic formulation for finite strains

This work describes the finite element implementation of a generalised strain gradient and rate-dependent crystallographic formulation for finite strains and general anisothermal conditions based on a multiplicative decomposition of the deformation gradient. The implementation involved the development of both a novel finite element formulation to determine the spatial slip rate gradients at each material point, and an implicit numerical integration scheme at the constitutive level to update the stresses and solution dependent variables. The time-integration procedure uses a Newton–Raphson scheme with a single level of iteration to solve the incremental non-linear equations associated with the non-local constitutive formulation. Closed-form solutions for the relevant fourth-order Jacobian tensors are given. The proposed numerical scheme is formulated in a general form and hence should be applicable to most existing crystallographic models. The crystallographic formulation is then used to investigate the effect of the morphology and volume fraction of the reinforcing phase of a two-phase single crystal on its macroscopic behaviour.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。     

A generalized variational formulation for convective heat transfer

The Scope of this paper is to develop the basic equations for a variational formulation which can be used to solve problems related to convection and/or diffusion dominated flows. The formulation is based on the introduction of a generalized quantity defined as the hear displacement. The governing equation is expressed in terms of this quantity and a variational formulation is developed which leads to a system of equations similar in form to Lagrange's equations of mechanics. These equations can be used for obtaining approximate solutions, though they are of particular interest for application of the finite element method. As an example of the formulation two finite element models are derived for solving convectiondiffusion boundary value problems. The performance of the two models is investigated and numerical results are given for different cases of convection and diffusion with two types of boundary conditions. The applications of the developed formulations are not limited to convection-diffusion problems but can also be applied to other types of problems such as mass transfer, hydrodynamics and wave propagation.     

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.     

Analytical modelling of thermal stresses in anisotropic composites

This paper represents a continuation of the author's previous work which deals with an analytical model of thermal stresses which originate during a cooling process of an anisotropic solid elastic continuum. This continuum consists of anisotropic spherical particles which are periodically distributed in an anisotropic infinite matrix. The infinite matrix is imaginarily divided into identical cubic cells with central particles. This multi-particle–matrix system represents a model system which is applicable to two-component materials of the precipitate–matrix type. The thermal stresses, which originate due to different thermal expansion coefficients of components of the model system, are determined within the cubic cell. The analytical modelling results from fundamental equations of continuum mechanics for solid elastic continuum (Cauchy's, compatibility and equilibrium equations, Hooke's law). This paper presents suitable mathematical procedures which are applied to the fundamental equations. These mathematical procedures lead to such final formulae for the thermal stresses which are relatively simple in comparison with the final formulae presented in the author's previous work which are extremely extensive. Using these new final formulae, the numerical determination of the thermal stresses in real two-component materials with anisotropic components is not time-consuming.     

基于子结构解耦的连接特性识别新方法

连接部件动态特性的准确辨识对预测装配式机械结构的动力学行为有重要意义. 针对传统基于子结构解耦的连接结构动力学特性识别方法难以直接利用可测量数据进行辨识及辨识结果受噪声影响显著等问题, 本文提出了一种新方法. 首先, 提取子结构解耦基本方程在测试自由度上的分量, 并经矩阵变换得到显含连接动刚度矩阵的形式, 而后由真实连接动刚度矩阵分解为已知的初始矩阵与待求的增量矩阵, 推导了具有收敛性质的增量型方程以增强界面自由度较多时辨识的数值稳定性, 并采用多项式拟合动刚度将其转化为了拟合系数的频域估计方程, 按给定准则选取合适的频率点联立方程后, 得到了只需装配体测试自由度上的频响函数来辨识连接特性的迭代公式. 最后, 以若干算例说明了算法的具体流程. 对10自由度弹簧?质量块系统进行了数值仿真, 验证了所提方法的正确性及抗噪性; 对包含一处胶接连接的T形梁结构和包含两处螺栓连接的L形梁结构进行了试验, 所辨识连接结构与残余结构重组的装配体有限元模型计算的频响函数与测量值在较宽频带内吻合较好, 表明了该方法能有效识别实际装配体结构中的连接特性.      

New way to construct high order Hamiltonian variational integrators

This paper develops a new approach to construct variational integrators. A simplified unconventional Hamilton’s variational principle corresponding to initial value problems is proposed, which is convenient for applications. The displacement and momentum are approximated with the same Lagrange interpolation. After the numerical integration and variational operation, the original problems are expressed as algebraic equations with the displacement and momentum at the interpolation points as unknown variables. Some particular variational integrators are derived. An optimal scheme of choosing initial values for the Newton-Raphson method is presented for the nonlinear dynamic system. In addition, specific examples show that the proposed integrators are symplectic when the interpolation point coincides with the numerical integration point, and both are Gaussian quadrature points. Meanwhile, compared with the same order symplectic Runge-Kutta methods, although the accuracy of the two methods is almost the same, the proposed integrators are much simpler and less computationally expensive.     

Free vibration analysis of elastically supported functionally graded annular plates subjected to thermal environment

Free vibration analysis of functionally graded (FG) thin-to-moderately thick annular plates subjected to thermal environment and supported on two-parameter elastic foundation is investigated. The material properties are assumed to be temperature-dependent and graded in the thickness direction. The equations of motion and the related boundary conditions, which include the effects of initial thermal stresses, are derived using the Hamilton’s principle based on the first order shear deformation theory (FSDT). The initial thermal stresses are obtained by solving the thermoelastic equilibrium equations. Differential quadrature method (DQM) as an efficient and accurate numerical tool is adopted to solve the thermoelastic equilibrium equations and the equations of motion. The formulations are validated by comparing the results in the limit cases with the available solutions in the literature for isotropic and FG circular and annular plates. The effects of the temperature rise, elastic foundation coefficients, the material graded index and different geometrical parameters on the frequency parameters of the FG annular plates are investigated. The new results can be used as benchmark solutions for future researches.     

Analyse contrastive de modèles incrémentaux et hypoplastiques

In this paper we first establish two necessary and sufficient conditions in order that incremental constitutive equations expressing the strain rate tensor as a function of the Jaumann's derivative of the Cauchy's stress tensor can be inverted under the general form of hypoplastic models when the stress state is located inside the domain bounded by the limit state surface. We are then interested in the physical meaning of these conditions with regard to the incremental response of the material.     

形状记忆合金管接头空间轴对称有限元分析

本文采用形状记忆合金（ＳＭＡ）的三维本构方程和有限变形理论，考虑拉、压不同应力状态对相变点移动的规律，编制了ＳＭＡ轴对称大变形的有限元程序，与单向拉伸下解析所得的应力、应变曲线相比，证实程序的正确性．文末计算一ＳＭＡ管接头，并指出按空间轴对称计算的必要性．     

Fully-coupled nonlinear analysis of single lap adhesive joints

This paper presents novel closed-form and accurate solutions for the edge moment factor and adhesive stresses for single lap adhesive bonded joints. In the present analysis of single lap joints, both large deflections of adherends and adhesive shear and peel strains are taken into account in the formulation of two sets of nonlinear governing equations for both longitudinal and transverse deflections of adherends. Closed-form solutions for the edge moment factor and the adhesive stresses are obtained by solving the two sets of fully-coupled nonlinear governing equations. Simplified and accurate formula for the edge moment factor is also derived via an approximation process. A comprehensive numerical validation was conducted by comparing the present solutions and those developed by Goland and Reissner, Hart-Smith and Oplinger with the results of nonlinear finite element analyses. Numerical results demonstrate that the present solutions for the edge moment factor (including the simplified formula) and the adhesive stresses appear to be the best as they agree extremely well with the finite element analysis results for all ranges of material and geometrical parameters.