A CRITICAL EVALUATION OF VARIOUS NONLINEAR BEAM FINITE ELEMENTS

2012 ◽  
Vol 09 (04) ◽  
pp. 1250045
Author(s):  
A. LAULUSA ◽  
J. N. REDDY

The characteristics of interdependent interpolation and mixed interpolation nonlinear beam finite elements are investigated in comparison with the equal-order interpolation element with uniform reduced integration. The stiffness matrix of the 3-noded and 4-noded equal order interpolation elements is identical to that of the 2-noded interdependent interpolation element if the internal nodal degrees-of-freedom are eliminated. The extension of the latter to include nonlinear kinematics by approximating the extensional displacement and the twist rotation with quadratic and cubic Lagrange polynomials yields unsatisfactory results. The 2-noded, 3-noded, and 4-noded mixed interpolation elements using one-, two-, and three-point quadrature rules, respectively, are shown to be equivalent to the corresponding uniform interpolation elements employing the same quadrature rules. The equivalence is established in the framework of nonlinear kinematics and anisotropic couplings.

Author(s):  
W. V. Nack

Abstract An approximate analysis method for brake squeal is presented. Using MSC/NASTRAN a geometric nonlinear solution is run using a friction stiffness matrix to model the contact between the pad and rotor. The friction coefficient can be pressure dependent. Next, linearized complex modes are found where the interface is set in a slip condition. Since the entire interface is set sliding, it produces the maximum friction work possible during the vibration. It is a conservative measure for stability evaluation. An averaged friction coefficient is measured and used during squeal. Dynamically unstable modes are found during squeal. They are due to friction coupling of neighboring modes. When these modes are decoupled, they are stabilized and squeal is eliminated. Good correlation with experimental results is shown. It will be shown that the complex modes baseline solution is insensitive to the type of variations in pressure and velocity that occur in a test schedule. This is due to the conservative nature of the approximation. Convective mass effects have not been included.


2018 ◽  
Vol 141 (3) ◽  
Author(s):  
Shengye Lin ◽  
Shuyun Jiang

This paper studies the stiffness characteristics of preloaded duplex angular contact ball bearings. First, a five degrees-of-freedom (5DOF) quasi-static model of the preloaded duplex angular contact ball bearing is established based on the Jones bearing model. Three bearing configurations (face-to-face, back-to-back, and tandem arrangements) and two preload mechanisms (constant pressure preload and fixed position preload) are included in the proposed model. Subsequently, the five-dimensional stiffness matrix of the preloaded duplex angular contact ball bearing is derived analytically. Then, an experimental setup is developed to measure the radial stiffness and the angular stiffness of duplex angular contact ball bearings. The simulated results match well with those from experiments, which prove the validity of the proposed model. Finally, the effects of bearing configuration, preload mechanism, and unloaded contact angle on the angular stiffness and the cross-coupling are studied systematically.


2010 ◽  
Vol 77 (6) ◽  
Author(s):  
M. Jafari ◽  
M. J. Mahjoob

In this paper, the exact stiffness matrix of curved beams with nonuniform cross section is derived using direct method. The considered element has two nodes and 12 degrees of freedom, with three forces and three moments applied at each node. The noncoincidence effect of shear center and center of area is also considered in this element. The deformations of the beam are due to bending, torsion, tensile, and shear loads. The line passing through center of area is a general three-dimensional curve and the cross section properties may change arbitrarily along it. The method is extended to deal with distributed loads on the curved beams. The stiffness matrix of some selected types of beams is determined by this method. The results are compared (where possible) with previously published results, simple beam finite element analysis and analytic solution. It is shown that the determined stiffness matrix is exact and that any type of beam can be analyzed by this method.


2000 ◽  
Vol 42 (1) ◽  
pp. 141-150 ◽  
Author(s):  
Mark Ainsworth ◽  
Bill McLean ◽  
Thanh Tran

AbstractA boundary integral equation of the first kind is discretised using Galerkin's method with piecewise-constant trial functions. We show how the condition number of the stiffness matrix depends on the number of degrees of freedom and on the global mesh ratio. We also show that diagonal scaling eliminates the latter dependence. Numerical experiments confirm the theory, and demonstrate that in practical computations involving strong local mesh refinement, diagonal scaling dramatically improves the conditioning of the Galerkin equations.


2019 ◽  
Vol 64 (3) ◽  
pp. 1-10
Author(s):  
Matteo Filippi ◽  
Alfonso Pagani ◽  
Erasmo Carrera

This paper proposes a geometrically nonlinear three-dimensional formalism for the static and dynamic study of rotor blades. The structures are modeled using high-order beam finite elements whose kinematics are input parameters of the analysis. The displacement fields are written using two-dimensional Taylor- and Lagrange-like expansions of the cross-sectional coordinates. As far as the Taylor-like polynomials are concerned, the linear case is similar to the first-order shear deformation theory, whereas the higher-order expansions include additional contributions that describe the warping of the cross section. The Lagrange-type kinematics instead utilizes the displacements of certain physical points as degrees of freedom. The inherent three-dimensional nature of the Carrera unified formulation enables one to include all Green–Lagrange strain components as well as all coupling effects due to the geometrical features and the three-dimensional constitutive law. A number of test cases are considered to compare the current solutions with experimental and theoretical results reported in terms of large deflections/rotations and frequencies related to small amplitude vibrations.


Mathematics ◽  
2020 ◽  
Vol 8 (8) ◽  
pp. 1358
Author(s):  
Shumin Man ◽  
Qiang Gao ◽  
Wanxie Zhong

Variational integrators for dynamic systems with holonomic constraints are proposed based on Hamilton’s principle. The variational principle is discretized by approximating the generalized coordinates and Lagrange multipliers by Lagrange polynomials, by approximating the integrals by quadrature rules. Meanwhile, constraint points are defined in order to discrete the holonomic constraints. The functional of the variational principle is divided into two parts, i.e., the action of the unconstrained term and the constrained term and the actions of the unconstrained term and the constrained term are integrated separately using different numerical quadrature rules. The influence of interpolation points, quadrature rule and constraint points on the accuracy of the algorithms is analyzed exhaustively. Properties of the proposed algorithms are investigated using examples. Numerical results show that the proposed algorithms have arbitrary high order, satisfy the holonomic constraints with high precision and provide good performance for long-time integration.


1977 ◽  
Vol 28 (2) ◽  
pp. 123-141 ◽  
Author(s):  
P J Holt ◽  
J P H Webber

SummaryThe formulation of curved finite elements to represent a two-dimensional circular sandwich ring with honeycomb core and laminated faces is investigated. Assumed stress hybrid and equilibrium methods are found to be easier to employ in this case than the displacement approach. Using these methods, an element stiffness matrix is developed. The approximations of membrane faces and an infinite core normal stiffness are then used to develop simpler elements. Test cases show that these assumptions may become invalid, but that they are adequate for most practical cases where the core thickness to radius ratio and the face thickness to core thickness ratio are both low.


Sign in / Sign up

Export Citation Format

Share Document