Contact dynamics of elasto-plastic thin beams simulated via absolute nodal coordinate formulation

Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation (ANCF). The inter-nal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation. A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient con-tact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previ-ously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the pro-posed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.     

Steady-state dynamics of a 3D tensegrity structure: Simulations and experiments

This paper considers a modeling and analysis approach for the investigation of the linear and nonlinear steady-state dynamics of a base excited 3D tensegrity module carrying a top mass. The tensegrity module contains three compressive members, which may buckle and six cables (tendons). First, a dynamic model of the system is derived using Lagrange’s equation with constraints. The buckling modeling of the compressive members is based on the assumed-mode method with a single mode discretization. The tendons are modeled as piecewise linear springs, which can only take tensile forces. This research focusses on the dynamic stability of the tensegrity structure by defining the geometrical and material properties in such a way that the system is just below the static stability boundary. Static and linear dynamic analysis is performed. In the nonlinear steady-state analysis, frequency-amplitude plots, power spectral density plots, bifurcation point continuation diagrams, and Poincaré maps are presented. A tensegrity structure is designed and manufactured and an experimental set-up is realized in order to validate the model by comparing experimentally and numerically obtained responses. In the validation stage, the numerical results are based on an amplifier-shaker-tensegrity structure model. It can be concluded that the numerical results match partly quantitatively and partly qualitatively with the experimentally obtained responses.     

可变体系的几何稳定平衡状态及受力分析

索和撑杆组成的张拉集成体系，膜结构及悬索体系，一般都要施加预应力使其成为能承受荷载的几何不变体系，整体受力表现为几何非线性。本文将其受力状态分解为两部分，一是几何稳定平衡状态，二是弹性稳定平衡状态。文中提出几何稳定平衡状态的分析计算方法，证明增量迭代法的收敛性。结果表明，所提出的方法收敛快，精度高。     

直接基于单元平衡的梁非线性分析方法

以单元力平衡条件为基础建立的梁单元在近年来受到广泛关注,与传统的以单元位移相容条件为基础建立的梁单元相比,基于平衡的梁单元可获得较高的精度。针对基于平衡的梁单元,给出了一种几何非线性分析方法,使用该方法进行几何非线性分析时不需进行单元内部位移场的构造,避免了使用现有位移场构造方法时可能引起的振荡问题,同时降低了有限元程...     

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.     

A novel method of determining the sole configuration of tensegrity structures

A novel analysis method is presented for form-finding of tensegrity structures. The spectral decomposition of the force density matrix and the singular value decomposition of the equilibrium matrix are performed iteratively to find the feasible sets of nodal coordinates and force densities. An algorithm of determining the sole configuration of free-form tensegrities is provided by specifying an independent set of nodal coordinates, which indicates the geometrical and mechanical properties of the structures can be at least partly controlled by the proposed method. Several illustrative examples are presented to demonstrate the efficiency and robustness in finding self-equilibrium configurations of tensegrity structures.     

Using co-rotational method for cracked frame analysis

A cracked element is formulated using the two-filed Hellinger–Reissner functional. Due to utilization of the linear elastic fracture mechanics, only geometrical nonlinearities can be considered for the cracked element. A clear step-by-step algorithm for the element state determination is also presented. The element flexibility matrix is derived in a basic coordinate system. Co-rotational approach is used to transform the element stiffness matrix and the resisting force vector from the basic system to the global one. The suggested element is applicable for static and dynamic analysis, as well as, the stress intensity factor calculation, and also inverse crack detection. Various numerical problems verify accuracy of the proposed element for linear and nonlinear structural analysis.     

预应力混凝土平面杆系结构的有限元方法

建立了基于有限元方法的考虑材料和几何非线性的任意截面预应力混凝土平面杆系结构的数值分析模型,可用于模拟预应力混凝土大跨度梁、单向偏压细长柱等的非线性全过程结构响应。引入修正的Rodriguez截面模型确定截面切线刚度,其中混凝土的贡献通过截面边界顶点定义的梯形单元来实现;在此基础上利用传统的平面非线性杆单元导出了标准有限元公式。通过两个算例验证了该模型的可靠性和适用性。     

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

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

Stability conditions of prestressed pin-jointed structures

Pin-jointed structures are first classified to trusses, tensile structures, and tensegrity structures in view of their respective stability properties. A sufficient condition for stability of an equilibrium state is derived for tensegrity structures. The condition is based on the bilinear forms of the linear and geometrical stiffness matrices considering the flexibility of members. The stability is defined by the positive definiteness of the tangent stiffness matrix, whereas the definition of prestress-stability is based on the geometrical stiffness matrix and the infinitesimal mechanisms. Numerical examples verify that the so-called super-stability condition might not be satisfied by a stable tensegrity structure, and that a prestress-stable structure can be unstable if the prestresses are moderately large.     

Dynamic analysis of an offshore pipe laying operation using the reel method

A system designed for a rigid and flexible pipe laying purposes is presented in the paper.Mathematical and numerical models are developed by using the rigid finite element method(RFEM).The RFEM is an efficient solution in the time domain.Static and dynamic problems related to pipe installation are solved by taking the advantage of simple interpretation and implementation of the method.Large deformations of the pipe during spooling and when it is reeled out at sea are considered.A material model implemented is used to take into consideration nonlinear material properties.In particular,the full elasto-plastic material characteristics with hardening and Bauschinger effect are included.Dynamic analyses are performed and the results attached in this work demonstrates how the sea conditions influence the machinery and pipeline,assuming a passive reel drive system. The influence of several other operational parameters on dynamic loads is verified.An active system,implemented as a part of the mathematical model,improves the system performance.Some results are presented as well.     

Reservoir water level effects on nonlinear dynamic response of arch dams

In this study, reservoir water level effects on nonlinear dynamic response of arch dams are investigated. For this purpose, the nonlinear behaviour of the dam concrete is idealized as elasto-plastic using the Drucker–Prager model based on the associated flow rule assumption. Water in the reservoir is represented by the Lagrangian (displacement-based) fluid finite elements. The program NONSAP is modified for elasto-plastic analysis of fluid–structure systems and employed in the response calculations. Nonlinear dynamic analysis of an arch dam subjected to earthquake ground motion is performed for five different water levels. The El-Centro N–S component of the Imperial Valley earthquake, on May 18, 1940, has been used as the ground motion. The crest displacements, the maximum tensile stresses on the upstream and downstream faces of the dam and the time history of the yield function of an element in the dam body are presented. The results obtained from nonlinear analyses for different water levels are compared with each other. It is apparent that the reservoir water level effects must be considered in the elasto-plastic analyses of arch dams to earthquake ground motion.     

基于非线性梁理论的有限质点法

有限质点法是以向量式力学为基础,用有限数量的质点来模拟结构的变形行为,质点的运动由牛顿运动定律来计算。在有限质点法中,质点通过构件相连,构件约束着质点的运动,并且其内力由质点的运动变量来描述。基于向量式力学的基本思想和非线性梁理论,提出了一种新的有限质点法,该方法在共旋单元坐标系中描述梁的非线性变形。以空间梁系结构为例,推导了计算构件内力的非线性公式,并考虑了弯扭耦合变形。通过两个连续欧拉角的变换公式得到共旋坐标系的旋转矩阵。与传统的有限质点法相比,本文提出的方法避免了刚体虚转动分析。通过四个结构的数值求解,验证了本文方法在计算结构大变形响应时具有较高的精度。     

Geometric and material nonlinear analysis of tensegrity structures

A numerical method is presented for the large deflection in elastic analysis of tensegrity structures including both geometric and material nonlinearities.The geometric nonlinearity is considered based on both total Lagrangian and updated Lagrangian formulations,while the material nonlinearity is treated through elastoplastic stress-strain relationship.The nonlinear equilibrium equations are solved using an incremental-iterative scheme in conjunction with the modified Newton-Raphson method.A computer program is developed to predict the mechanical responses of tensegrity systems under tensile,compressive and flexural loadings.Numerical results obtained are compared with those reported in the literature to demonstrate the accuracy and efficiency of the proposed program.The flexural behavior of the double layer quadruplex tensegrity grid is sufficiently good for lightweight large-span structural applications.On the other hand,its bending strength capacity is not sensitive to the self-stress level.     

Adaptive force density method for form-finding problem of tensegrity structures

A numerical method is presented for form-finding of tensegrity structures. Eigenvalue analysis and spectral decomposition are carried out iteratively to find the feasible set of force densities that satisfies the requirement on rank deficiency of the equilibrium matrix with respect to the nodal coordinates. The equilibrium matrix is shown to correspond to the geometrical stiffness matrix in the conventional finite element formulation. A unique and non-degenerate configuration of the structure can then be obtained by specifying an independent set of nodal coordinates. A simple explanation is given for the required rank deficiency of the equilibrium matrix that leads to a non-degenerate structure. Several numerical examples are presented to illustrate the robustness as well as the strong ability of searching new configurations of the proposed method.     

ANALYSIS OF ELASTO-PLASTIC POSTBUCKLING AND ENERGY RELEASE RATE FOR DELAMINATED FIBER METAL LAMINATED BEAMS IN THERMAL ENVIRONMENT

The elasto-plastic postbuckling of fiber metal laminated beams with delamination and the energy release rate along the delamination front are discussed in this paper. Considering geometrical nonlinearity, thermal environment and geometrical initial imperfection, the incremental nonlinear equilibrium equations of delaminated fiber metal laminated beams are established,which are solved using the differential quadrature method and iterative method. Based on these,according to the J-integral theory, the elasto-plastic energy release rate is studied. The effects of some important parameters on the elasto-plastic postbuckling behavior and energy release rate of the aramid reinforced aluminum laminated beams are discussed in details.     

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

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

Finite element based form-finding algorithm for tensegrity structures

This paper presents a novel form-finding algorithm for tensegrity structures that is based on the finite element method. The required data for the form-finding is the topology of the structure, undeformed bar lengths, total cable length, prestress of cables and stiffness of bars. The form-finding is done by modifying the single cable lengths such that the total cable length is preserved and the potential energy of the system is minimized. Two- and three-dimensional examples are presented that demonstrate the excellent performance of the proposed algorithm.     

Non-linear dynamic analysis of tensegrity structures using a co-rotational method

A new formulation is presented for the non-linear dynamic analysis of space truss structures. The formulation is based on the dynamics of 3D co-rotational rods. In the co-rotation method, the rigid body modes are assumed to be separated from the total deformations at the local element level. In this paper a new co-rotational formulation is proposed based on the direct derivation of the inertia force vector and the tangent dynamic matrix. A closed-form equation is derived for the calculation of the inertia force, the tangent dynamic matrix, the mass matrix and the gyroscopic matrix. The new formulation is used to perform dynamic analysis of example tensegrity structures. The developed formulation is applicable to tensegrity structures with non-linear effects due to internal mechanisms or geometric non-linearities, and is applied to two numerical examples. The efficiency of the proposed approach is compared to the conventional Lagrangian method, and savings in computation of about 55%, 54% and 37% were achieved.