Showing papers on "Tangent stiffness matrix published in 2019"

Dual-support smoothed particle hydrodynamics in solid: variational principle and implicit formulation

Abstract: We derive the dual-support smoothed particle hydrodynamics (DS-SPH) in solid within the framework of variational principle. The tangent stiffness matrix of SPH can be obtained with ease, and can be served as the basis for the present implicit SPH. We propose an hourglass energy functional, which allows the direct derivation of hourglass force and hourglass tangent stiffness matrix. The dual-support is involved in all derivations based on variational principles and is automatically satisfied in the assembling of stiffness matrix. The implementation of stiffness matrix comprises with two steps, the nodal assembly based on deformation gradient and global assembly on all nodes. Several numerical examples are presented to validate the method.

On non-Fourier flux in nonlinear stretching flow of hyperbolic tangent material

Abstract: The present study explores the features of hyperbolic tangent material due to a nonlinear stretched sheet with variable sheet thickness. Non-Fourier flux theory is implemented for the development of energy expression. Such consideration accounts for the contribution by thermal relaxation. The resulting nonlinear differential system has been determined for the convergent series expressions of velocity and temperature. The solutions are demonstrated and analyzed through plots. Presented results indicate that velocity decays via larger material power law index and Weissenberg number. Temperature is the decreasing function of Prandtl number and thermal relaxation time.

A nonlinear quadrilateral thin flat layered shell element for the modeling of reinforced concrete wall structures

Abstract: In this article, a simple and accurate quadrilateral thin flat layered shell element formulation for the nonlinear analysis of reinforced concrete (RC) wall systems under static and cycling loads is presented. The 4 node shell element, with 6 degree of freedom (DOF) per node (3 displacements and 3 rotations) is created by superposing the quadrilateral layered membrane element with drilling degrees of freedom (12 DOF, 2 displacement and 1 rotation per node) developed by Rojas et al. (Eng Struct 124:521–538, 2016), and the Discrete Kirchhoff Quadrilateral Element (12 DOF, 1 displacement and 2 rotations) formulated by Batoz and Tahar (Int J Numer Methods Eng 18(11):1655–1677, 1982), to model the in-plane and the out of plane bending behavior of the shell element, respectively. In addition, to model the complex behavior and coupling of the axial, flexural and shear behavior, observed in complex RC wall structures, the transversal section of the shell element consists of a layered system of fully bonded, smeared steel reinforcement and smeared orthotropic concrete material with the rotating angle formulation. The formulation used a tangent stiffness matrix approach, which include the coupling of membrane and bending effects. For verification, the shell element formulation is used to model a set of experimental results for T-shaped RC walls that are available in the literature. The proposed element is robust, simple to implement, and it can predict the global results (load vs. displacement and maximum capacity) and also the local behavior (vertical strain at the base level along the web and the flange) observed in RC wall structures.

Prestress Design of Tensegrity Structures Using Semidefinite Programming

Abstract: Finding appropriate prestresses which can stabilize the system is a key step in the design of tensegrity structures. A semidefinite programming- (SDP-) based approach is developed in this paper to determine appropriate prestresses for tensegrity structures. Three different stability criteria of tensegrity structures are considered in the proposed approach. Besides, the unilateral property of members and the evenness of internal forces are taken into account. The stiffness of the whole system can also be optimized by maximizing the minimum eigenvalue of the tangent stiffness matrix. Deterministic algorithms are used to solve the semidefinite programming problem in polynomial time. The applicability of the proposed approach is verified by three typical examples. Compared to previous stochastic-based approaches, the global optimality of the solution of the proposed approach is theoretically guaranteed and the solution is exactly reproducible.

Analytical Formulation of Friction Contact Elements for Frequency-Domain Analysis of Nonlinear Vibrations of Structures With High-Energy Rubs

Abstract: In gas-turbine engines and other rotating machinery structures rubbing contact interactions can occur when the contacting components have large relative motion between components: such as in rotating bladed disc-casing rubbing contacts, rubbing in rotor bearing and labyrinth seals, etc. The analysis of vibrations of structures with rubbing contacts requires the development of a mathematical model and special friction contact elements that would allow for the prescribed relative motion of rubbing surfaces in addition to the motion due to vibrations of the contacting components. In the proposed paper, the formulation of the friction contact elements is developed which includes the effects of the prescribed relative motion on the friction stick-slip transitions and, therefore, on the contact interaction forces. For a first time, the formulation is made for the frequency domain analysis of coupled rubbing and vibrational motion, using the multiharmonic representation of the vibration displacements. The formulation is made fully analytically to express the multiharmonic contact interaction forces and multiharmonic tangent stiffness matrix in an explicit analytical form allowing their calculation accurately and fast. The dependency of the friction and contact stiffness coefficients on the energy dissipated during high-energy rubbing contacts and, hence, on the corresponding increase of the contact interface temperature is included in the formulation. The efficiency of the developed friction elements is demonstrated on a set of test cases including simple models and a large-scale realistic blade.

Effect of Shear and Distortion Deformations on Lateral Buckling Resistance of Box Elements in the Framework of Eurocode 3

Abstract: This paper treats the distortional and shear deformation effects on the elastic lateral torsional buckling of thin-walled box beam elements, under combined bending and axial forces. For the purpose, a nonlinear kinematic model based on higher order theory is used applicable to both short and long thin-walled box beams. Because in the kinematic model of the higher order theory integrates additional flexibility terms related to shear, distortion and warping effects, it accurately predicts the lateral torsional buckling of the straight box beams. Ritz’s method is adopted as solution strategy in order to obtain the nonlinear governing equilibrium equations, then the buckling loads are computed by solving the eigenvalue problem basing on the singularity of the tangential stiffness matrix. Owing to flexural–torsional and distortional couplings, new matrices are obtained in both geometric and initial stress parts of the tangent stiffness matrix. The proposed method with the new stiffness terms, is efficient and accurate in lateral torsional buckling predictions, when compared with the commercial FEM code ABAQUS results. Based on the existing European guidelines EC3, an extensive numerical investigation is performed to demonstrate the effects of both shear and distortional deformations on the moment carrying capacity. The convenience of the model is outlined and the limit of models developed without shear and distortion deformation effects on lateral buckling loads evaluation is discussed.

Nonlinear Analysis of Spatial Cable of Long-Span Cable-Stayed Bridge considering Rigid Connection

Abstract: To determine the cable sag effect and solve the rigid connection problem at cable ends of a long-span cable-stayed bridge, a two-node spatial catenary cable element with arbitrary rigid arms is developed. Using the finite rotation formula of the space vector and a differential method, the incremental relation between displacement and force at both ends of the rigid arm is given. Then, explicit expression of the tangent stiffness matrix of the element with arbitrary rigid arms is derived based on the catenary equations. Two numerical examples are provided to verify the validity of the new element. A long-span cable-stayed bridge application model is established, and the cables are simulated using three methods. The results show that the rigid end effect has influence on displacements, bending moments and rigidity and should not be ignored. The catenary cable element with arbitrary rigid arms can be used to simulate the geometric nonlinear mechanical behavior of the cables and can well solve the rigid connection problem.

A fiber beam element model for elastic-plastic analysis of girders with shear lag effects

Abstract: This paper proposes a one-dimensional fiber beam element model taking account of materially non-linear behavior, benefiting the highly efficient elastic-plastic analysis of girders with shear-lag effects. Based on the displacement-based fiber beam-column element, two additional degrees of freedom (DOFs) are added into the proposed model to consider the shear-lag warping deformations of the slabs. The new finite element (FE) formulations of the tangent stiffness matrix and resisting force vector are deduced with the variational principle of the minimum potential energy. Then the proposed element is implemented in the OpenSees computational framework as a newly developed element, and the full Newton iteration method is adopted for an iterative solution. The typical materially non-linear behaviors, including the cracking and crushing of concrete, as well as the plasticity of the reinforcement and steel girder, are all considered in the model. The proposed model is applied to several test cases under elastic or plastic loading states and compared with the solutions of theoretical models, tests, and shell/solid refined FE models. The results of these comparisons indicate the accuracy and applicability of the proposed model for the analysis of both concrete box girders and steel-concrete composite girders, under either elastic or plastic states.

The tangent stiffness matrix for an absolute interface coordinates floating frame of reference formulation

Abstract: In this work, a full and complete development of the tangent stiffness matrix is presented, suitable for the use in an absolute interface coordinates floating frame of reference formulation. For simulation of flexible multibody systems, the floating frame formulation is used for its advantage to describe local elastic deformation by means of a body’s linear finite element model. Consequently, the powerful Craig–Bampton method can be applied for model order reduction. By establishing a coordinate transformation from the absolute floating frame coordinates and local interface coordinates corresponding to the Craig–Bampton modes to absolute interface coordinates, it is possible to satisfy kinematic constraints without the use of Lagrange multipliers. In this way, the floating frame does not need to be located at an interface point and can be positioned close to the body’s center of mass, without requiring an interface point at the center of mass. This improves simulation accuracy. In this work, the expression for the new method’s tangent stiffness matrix is obtained by taking the variation of the equation of equilibrium. The global tangent stiffness matrix is expressed as a local tangent stiffness matrix, consisting of both material stiffness and geometric stiffness terms, transformed to the global frame by the rotation matrix of the floating frame. Simulations of static and dynamic validation problems are performed. These simulations show the importance of including the tangent stiffness matrix for both convergence and simulation efficiency.

Accelerated complex-step finite difference for expedient deformable simulation

Abstract: In deformable simulation, an important computing task is to calculate the gradient and derivative of the strain energy function in order to infer the corresponding internal force and tangent stiffness matrix. The standard numerical routine is the finite difference method, which evaluates the target function multiple times under a small real-valued perturbation. Unfortunately, the subtractive cancellation prevents us from setting this perturbation sufficiently small, and the regular finite difference is doomed for computing problems requiring a high-accuracy derivative evaluation. In this paper, we graft a new finite difference scheme, namely the complex-step finite difference (CSFD), with physics-based animation. CSFD is based on the complex Taylor series expansion, which avoids subtractions in first-order derivative approximation. As a result, one can use a very small perturbation to calculate the numerical derivative that is as accurate as its analytic counterpart. We accelerate the original CSFD method so that it is also as efficient as the analytic derivative. This is achieved by discarding high-order error terms, decoupling real and imaginary calculations, replacing costly functions based on the theory of equivalent infinitesimal, and isolating the propagation of the perturbation in composite/nesting functions. CSFD can be further augmented with multicomplex Taylor expansion and Cauchy-Riemann formula to handle higher-order derivatives and tensor-valued functions. We demonstrate the accuracy, convenience, and efficiency of this new numerical routine in the context of deformable simulation - one can easily deploy a robust simulator for general hyperelastic materials, including user-crafted ones to cater specific needs in different applications. Higher-order derivatives of the energy can be readily computed to construct modal derivative bases for reduced real-time simulation. Inverse simulation problems can also be conveniently solved using gradient/Hessian-based optimization procedures.

Higher order nonlocal operator method

Abstract: We extend the nonlocal operator method to higher order scheme by using a higher order Taylor series expansion of the unknown field. Such a higher order scheme improves the original nonlocal operator method proposed by the authors in [A nonlocal operator method for solving partial differential equations], which can only achieve one-order convergence. The higher order nonlocal operator method obtains all partial derivatives with specified maximal order simultaneously without resorting to shape functions. The functional based on the nonlocal operators converts the construction of residual and stiffness matrix into a series of matrix multiplication on the nonlocal operator matrix. Several numerical examples solved by strong form or weak form are presented to show the capabilities of this method.

Linearized Dynamic Analysis of Weightless Sagging Planar Elastic Cables

Abstract: Nonlinear dynamic analysis of elastic structures is known to be much more complex than their linear analysis. There are many sources of nonlinearity of the structural response of elastic cables, viz., physical nonlinearity due to nonlinear tension–extension relations, geometric nonlinearity associated with finite elastic displacements and nonlinearity of nodal load–displacement relations due to the presence of self-weight. Incremental second-order differential equations of motion are used to predict the vibration amplitudes relative to the equilibrium state caused by additional dynamic forces. Generally, the tangent stiffness matrices are determined by adding the tangent elastic and geometric stiffness matrices. Many a time, an approximate linearized dynamic analysis is attempted. In this paper, the initial tangent stiffness matrix corresponding to the equilibrium state is used in the second-order linear differential equation of motion. The dynamic response relative to the equilibrium state of the structure subjected to additional dynamic loads is predicted. The predictions of linearized dynamic analysis are generally considered acceptable for small elastic displacements from the equilibrium state. The validity of such linearized dynamic analysis for elasto-flexible cables obeying third-order differential equation of motion is explored.