The strain hardening rotating hollow shaft subject to a positive temperature gradient

Summary The modeling of deformation in flexible multibody systems is still under intensive investigation. While the floating frame of reference formulation has become a standard for the modeling of deformable moving bodies, formulations based on absolute coordinates are comparatively new. The recently developed absolute nodal coordinate formulation uses solely nodal position and slopes as degrees of freedom for structural elements. The numerical treatment is similar to the absolute coordinate formulation, which is investigated in the following. An efficient formulation based on absolute coordinates with a reduced strain tensor (similar to corotational formulations) has been derived recently, and the analogy to the floating frame of reference formulation has been shown. The efficiency of this formulation is based on a co-rotated stiffness matrix which is factorized only once, however, only linear constraints have been treated up to now. The present paper treats the extension of the advantages of the co-rotated formulation with respect to nonlinear constraints, damping and contact. Constraints are discussed for the special cases of linear and nonlinear dependence on the deformation degrees of freedom. The (implicit) Newmark scheme is used to illustrate the time stepping procedure within the present method. In the case of a small number of nonlinear constraints, a single time step can be performed by solving only a small system of nonlinear equations by means of Newton's method. A numerical example shows a method to reduce the number of constraints and illustrates the computational advantages of the proposed method.     

A new iterative integral formulation for semilinear equations based on the generalized quasilinearization theory

This paper presents a new iterative integral approach for solving semilinear equations. The integral formulation is derived based on the generalized quasilinearization theory in which nonlinear equations are replaced by a set of iterative linear equations. An advantage of the new formulation is that its convergence is guaranteed under a given condition and the convergence rate can be quadratic. The effectiveness of the new approach has been demonstrated on several examples of the nonlinear Poisson type. Comparisons with some existing methods and a study of the convergence rate have also been conducted in this work.     

A nonlinear thick plate formulation based on the modified strain gradient theory

A geometric nonlinear first-order shear deformation theory-based formulation is presented to analyze microplates. The formulations derived herein are based on a modified strain gradient theory and the von Karman nonlinear strains. The modified strain gradient theory includes five material length scale parameters capable to capture the size effects in small scales. The governing equations of motion and the most general form of boundary conditions of an arbitrary-shaped plate are derived using the principle of virtual displacements. The analysis is general and can be reduced to the modified couple stress plate model or the classical plate model.     

Geometrically nonlinear formulation for the axisymmetric shell elements

A geometrically nonlinear formulation using total Lagrangian approach is presented for the axisymmetric shell elements. The basic element is formulated using the co-ordinates of the mid-surface nodes and the mid-surface nodal point normals. An important aspect of the formulation presented here is that the restriction on the magnitude of the nodal rotations is eliminated. This is accomplished by retaining true nonlinear nodal rotation terms in the definition of the displacement field and the consistent derivation of the element properties based on this displacement field. The element properties are derived and presented in detail. Numerical examples are also presented to demonstrate the element behaviour and the accuracy.     

Variational sensitivity analysis of a nonlinear solid shell element

The paper is concerned with variational sensitivity analysis of a nonlinear solid shell element, which is based on the Hu–Washizu variational principle. The sensitivity information is derived on the continuous level and discretized to yield the analytical expressions on the computational level. Especially, the pseudo load matrix and the sensitivity matrix, which dominate design sensitivity analysis of shape optimization problems, are derived. Because of the mixed formulation, condensation of the pseudo load matrix on the element level is performed to compute the sensitivity matrix. An illustrative example from the field of geometry‐based shape optimization demonstrates the possible application of the presented formulation. Copyright © 2013 John Wiley & Sons, Ltd.     

An integral-equation formulation of nonlinear deformation in a stack of buffered plates

An integral-equation formulation has been derived for nonlinear deformation in a stack of buffered Kirchhoff plates. The plates are assumed to follow a nonlinear bending moment-curvature law and the buffer material to follow the generalized Hooke’s law. By employing the recently derived special Green’s function for multilayers with interfacial membrane and flexural rigidities as the kernel, the integral-equation formulation only involves the surface loading area (for application to an indentation problem) and the portion of plates undergoing nonlinear deformation. Based on the integral equation, an efficient and accurate boundary element method has been derived to numerically solve the cylindrical indentation problem of the material with a bilinear flexural bending law for the plates. Numerical examples are presented to show a progressive damage process of yielding across a stack of plates as well as to demonstrate the validity and accuracy of the present integral-equation formulation.     

A peridynamic Reissner-Mindlin shell theory

In this work, we have developed a state-based peridynamics theory for nonlinear Reissener-Mindlin shells to model and predict large deformation of shell structures with thick wall. The nonlocal peridynamic theory of solids offers an integral formulation that is an alternative to traditional local continuum mechanics models based on partial differential equations. This formulation is applicable for solving the material failure problems involved in discontinuous displacement fields. The governing equations of the state-based peridynamic shell theory are derived based on the nonlocal balance laws by adopting the kinematic assumption of the Reissner and Mindlin plate and shell theories. In the numerical calculations, the stress points are employed to ensure the numerical stability. Several numerical examples are conducted to validate the nonlocal structure mechanics model and to verify the accuracy as well as the convergence of the proposed shell theory.     

基于互补模型的弹塑性桁架结构灵敏度分析过程

提出了一个用于弹塑性结构灵敏度分析的方法,称为NCP函数法。通过引入NCP函数,弹塑性结构分析的互补模型被等价为方程组形式;利用直接法对等价的方程组模型关于设计变量直接求导,建立了用于弹塑性灵敏度分析的线性方程组模型。新的灵敏度分析模型精确刻画了材料性质发生变化时灵敏度系数的不连续性,并避免了通过迭代过程来确定灵敏度系数的跳跃量。针对桁架结构的数值试验显示了算法的特点及其有效性。     

A physically and geometrically nonlinear scaled‐boundary‐based finite element formulation for fracture in elastomers

This paper is devoted to the formulation of a plane scaled boundary finite element with initially constant thickness for physically and geometrically nonlinear material behavior. Special two‐dimensional element shape functions are derived by using the analytical displacement solution of the standard scaled boundary finite element method, which is originally based on linear material behavior and small strains. These 2D shape functions can be constructed for an arbitrary number of element nodes and allow to capture singularities (e.g., at a plane crack tip) analytically, without extensive mesh refinement. Mapping these proposed 2D shape functions to the 3D case, a formulation that is compatible with standard finite elements is obtained. The resulting physically and geometrically nonlinear scaled boundary finite element formulation is implemented into the framework of the finite element method for bounded plane domains with and without geometrical singularities. The numerical realization is shown in detail. To represent the physically and geometrically nonlinear material and structural behavior of elastomer specimens, the extended tube model and the Yeoh model are used. Numerical studies on the convergence behavior and comparisons with standard Q1P0 finite elements demonstrate the correct implementation and the advantages of the developed scaled boundary finite element. Copyright © 2014 John Wiley & Sons, Ltd.     

Satellite attitude stabilization using fluid rings

The three-dimensional attitude stabilization of a satellite using fluid rings is proposed in the present paper. The general formulation of the system that comprises a satellite and three fluid rings, one on each axis, is obtained through Euler’s moment equations. The linearized system model is derived and the stability conditions are obtained. Linear control laws and nonlinear control laws based on sliding mode control techniques are developed for fluid control torques. The numerical simulation of the governing nonlinear equations of motion of the system along with stability analysis establishes the feasibility of achieving desired attitude stabilization of the satellite. The fluid controllers are successful in stabilizing the attitude of the satellite in presence of high attitude disturbance torques and intermittent actuators’ failures. In the case of sudden attitude disturbance torques, the nonlinear fluid controllers outperform the linear fluid controllers with virtually no effect on the satellite attitude response. The proposed attitude stabilization method may find applications in future satellite missions.     

A two-level mesh repartitioning scheme for the displacement-based lower-order finite element methods in volumetric locking-free analyses

We present an approach for repartitioning existing lower-order finite element mesh based on quadrilateral or triangular elements for the linear and nonlinear volumetric locking-free analysis. This approach contains two levels of mesh repartitioning. The first-level mesh re-partitioning is an h-adaptive mesh refinement for the generation of a refined mesh needed in the second-level mesh coarsening. The second-level mesh coarsening involves a gradient smoothing scheme performed on each pair of adjacent elements selected based on the first-level refined mesh. With the repartitioned mesh and smoothed gradient, the equivalence between the mixed finite element formulation and the displacement-based finite element formulation is established. The extension to nonlinear finite element formulation is also considered. Several linear and non-linear numerical benchmarks are solved and numerical inf-sup tests are conducted to demonstrate the accuracy and stability of the proposed formulation in the nearly incompressible applications.     

Geometrically nonlinear formulation for the curved shell elements

This paper presents a geometrically nonlinear formulation using total lagrangian approach for the three-dimensional curved shell elements. The basic element geometry is constructed using the coordinates of the middle surface nodes and the mid-surface nodal point normals. The element displacement field is described using three translations of the mid-surface nodes and the two rotations about the local axes. The existing shell element formulations are restricted to small nodal rotations between two successive load increments. The element formulation presented here removes such restrictions. This is accomplished by retaining nonlinear nodal rotation terms in the definition of the displacement field and the consistent derivation of the element properties. The formulation presented here is very general and yet can be made specific by selecting proper nonlinear functions representing the effects of nodal rotations. The element properties are derived and presented in detail. Numerical examples are also presented to demonstrate the behaviour and the accuracy of the elements.     

A nonlinear Timoshenko beam formulation based on the modified couple stress theory

This paper presents a nonlinear size-dependent Timoshenko beam model based on the modified couple stress theory, a non-classical continuum theory capable of capturing the size effects. The nonlinear behavior of the new model is due to the present of induced mid-plane stretching, a prevalent phenomenon in beams with two immovable supports. The Hamilton principle is employed to determine the governing partial differential equations as well as the boundary conditions. A hinged–hinged beam is chosen as an example to delineate the nonlinear size-dependent static and free-vibration behaviors of the derived formulation. The solution for the static bending is obtained numerically. The solution for the free-vibration is presented analytically utilizing the method of multiple scales, one of the perturbation techniques.     

On the equivalent static load method for flexible multibody systems described with a nonlinear finite element formalism

The equivalent static load (ESL) method is a powerful approach to solve dynamic response structural optimization problems. The method transforms the dynamic response optimization into a static response optimization under multiple load cases. The ESL cases are defined based on the transient analysis response whereupon all the standard techniques of static response optimization can be used. In the last decade, the ESL method has been applied to perform the structural optimization of flexible components of mechanical systems modeled as multibody systems (MBS). The ESL evaluation strongly depends on the adopted formulation to describe the MBS and has been initially derived based on a floating frame of reference formulation. In this paper, we propose a method to derive the ESL adapted to a nonlinear finite element approach based on a Lie group formalism for two main reasons. Firstly, the finite element approach is completely general to analyze complex MBS and is suitable to perform more advanced optimization problems like topology optimization. Secondly, the selected Lie group formalism leads to a formulation of the equations of motion in the local frame, which turns out to be a strong practical advantage for the ESL evaluation. Examples are provided to validate the proposed method. Copyright © 2016 John Wiley & Sons, Ltd.     

Electricity market equilibrium analysis based on nonlinear interior point algorithm with complementarity constraints

This paper presents an interior point algorithm based on a.c. network model for determining the Nash supply function equilibrium (SFE) of bid-based electricity markets. The SFE problem is considered as a bi-level game. At the first level, the problem begins with the formulation of an optimal power flow (OPF) interior point-based algorithm to handle the independent system operator (ISO) problem for maximising social welfare. This algorithm is based on the OPF with a.c. network transmission model taking into account all the operating aspects such as the generation capacity limits, bus voltage limits, transmission line constraints, network losses and especially the effect of the reactive power. The resulting Karush-Kuhn-Tucker (KKT) conditions of the problem at the first level are then reformulated using nonlinear complementarity constraints and incorporated as equality constraints in the second-level formulation for maximising the individual profit for each strategic generating firm. By employing a special nonlinear complementarity function, the complementarity constraints are then transformed into nonlinear algebraic expressions, thus the KKT conditions of the resulting combined problem can be derived. The final problem is then solved iteratively based on the solution techniques of the interior point algorithm. Numerical examples of a three-bus system, the IEEE 14-bus system and the IEEE 30-bus system, show that the algorithm can successfully determine the electricity market equilibrium with the a.c. network model.     

有限元线法求解非线性模型问题——Ⅰ.薄膜大挠度

本研究将新近发展起来的有限元线法应用于非线性问题,分析求解了若干具有代表性的模型问题,探讨了统一的求解模式及相应的处理手段。作为这一系列工作的首例,本文将该法应用于薄膜大挠度这一几何非线性模型问题,对任意形状的薄膜作了理论公式推导,通过对几种典型形状薄膜的具体数值计算,揭示了该类问题存在极限变形状态这一重要特性。数值算例的精确性与可靠性以及求解的高效性表明,本法是求解这类几何非线性问题的高效能的方法。     

A nonlinear viscoelastic bushing element in multibody dynamics

This paper presents a formulation for incorporating nonlinear viscoelastic bushing elements into multibody systems. The formulation is based on the assumption that the relaxation function can be expressed as a sum of functions which are nonlinear in deformation and exponentially decreasing in time. These forces can represent elastomeric mounts or bushings in automotive suspension systems. The numerical implementation of the nonlinear viscoelastic bushing model into a general purpose rigid multibody dynamics code is described, and an extension of the formulation is also presented wherein component flexibility is included. Model validation was performed by comparing experimental data to simulation results obtained using the nonlinear viscoelastic model and a nonlinear elastic model. The experimental data were obtained at the Center's facilities by testing an automotive lower control arm/bushing system, subjected to a simulated road load event. The comparison demonstrates the better load prediction capability of the viscoelastic bushing model compared to the conventional model.     

热冲击下的复合材料壳刚柔耦合动力学研究

研究了在热冲击下任意形状(仅一个方向有曲率)复合材料壳的非线性刚柔耦合动力学响应。根据Mindlin理论,建立了任意形状的复合材料壳的非线性应变-位移关系。借助于数学理论以及几何关系,描述了壳上任意点的变曲率。用虚功原理建立了动力学变分方程,并采用等参单元对壳的连续动力学方程进行离散,建立了中心刚体-复合材料壳的刚-柔耦合动力学方程。用高斯积分计算常值阵,为了提高计算效率,采用广义-α法结合Newton-Raph-son迭代法对动力学方程进行积分。将采用该方法计算得到的频率与ANSYS软件计算得到的作对比,验证了模型的正确性。通过算例分析了在热冲击作用下复合材料壳的线性、非线性的动力学特性,以及曲率、材料特性对动力学响应的影响。