Penalty Formulations for the Dynamic Analysis of Elastic Mechanisms

1989 ◽  
Vol 111 (3) ◽  
pp. 321-327 ◽  
Author(s):  
E. Bayo ◽  
M. A. Serna

A series of penalty methods are presented for the dynamic analysis of flexible mechanisms. The proposed methods formulate the equations of motion with respect to a floating frame that follows the rigid body motion of the links. The constraint conditions are not appended to the Lagrange’s equations in the form of algebraic or differential constraints, but inserted in them by means of a penalty formulation, and therefore the number of equations of the system does not increase. Furthermore, the discretization of the equations using the finite element method leads to a system of ordinary differential equations that can be solved using standard numerical algorithms. The proposed methods are valid for three dimensional analysis and can be very easily implemented in existing codes. Furthermore, they can be used to model any type of constraint conditions, either holonomic or nonholonomic, and with any degree of redundancy. A series of mechanisms composed of elastic members are analyzed. The results demonstrate the capabilities of the proposed methods for simulation analysis.

Author(s):  
Olivier A. Bauchau ◽  
Minghe Shan

Abstract The application of the finite element method to the modeling of Cosserat solids is investigated in detail. In two- and three-dimensional elasticity problems, the nodal unknowns are the components of the displacement vector, which form a linear field. In contrast, when dealing with Cosserat solids, the nodal unknowns form the special Euclidean group SE(3), a nonlinear manifold. This observation has numerous implications on the implementation of the finite element method and raises numerous questions: (1) What is the most suitable representation of this nonlinear manifold? (2) How is it interpolated over one element? (3) How is the associated strain field interpolated? (4) What is the most efficient way to obtain the discrete equations of motion? All these questions are, of course intertwined. This paper shows that reliable schemes are available for the interpolation of the motion and curvature fields. The interpolated fields depend on relative nodal motions only, and hence, are both objective and tensorial. Because these schemes depend on relative nodal motions only, only local parameterization is required, thereby avoiding the occurrence of singularities. For Cosserat solids, it is preferable to perform the discretization operation first, followed by the variation operation. This approach leads to considerable computation efficiency and simplicity.


1977 ◽  
Vol 99 (3) ◽  
pp. 773-779 ◽  
Author(s):  
N. Orlandea ◽  
M. A. Chace ◽  
D. A. Calahan

The work described herein is an extension of sparse matrix and stiff integrated numerical algorithms used for the simulation of electrical circuits and three-dimensional mechanical dynamic systems. By applying these algorithms big sets of sparse linear equations can be solved efficiently, and the numerical instability associated with widely split eigenvalues can be avoided. The new numerical methods affect even the initial formulation for these problems. In this paper, the equations of motion and constraints (Part 1) and the force function of springs and dampers (Part 2) are set up, and the numerical solutions for static, transient, and linearized types of analysis as well as the modal optimization algorithms are implemented in the ADAMS (automatic dynamic analysis of mechanical systems) computer program for simulation of three-dimensional mechanical systems (Part 2). The paper concludes with two examples: computer simulation of the front suspension of a 1973 Chevrolet Malibu and computer simulation of the landing gear of a Boeing 747 airplane. The efficiency of simulation and comparison with experimental results are given in tabular form.


1977 ◽  
Vol 99 (3) ◽  
pp. 780-784 ◽  
Author(s):  
N. Orlandea ◽  
D. A. Calahan ◽  
M. A. Chace

The work described herein is an extension of sparse matrix and stiff integrated numerical algorithms used for the simulation of electrical circuits and three-dimensional mechanical dynamic systems. By applying these algorithms, big sets of sparse linear equations can be solved efficiently, and the numerical instability associated with widely split eigenvalues can be avoided. The new numerical methods affect even the initial formulation for these problems. In this paper, the equations of motion and constraints (Part 1) and the force function of springs and dampers (Part 2) are set up, and the numerical solutions for static, transient, and linearized types of analysis as well as the model optimization algorithms are implemented in the ADAMS (automatic dynamic analysis of mechanical systems) computer program for simulation of three-dimensional mechanical systems (Part 2). The paper concludes with two examples: computer simulation of the front suspension of a 1973 Chevrolet Malibu and computer simulation of the landing gear of a Boeing 747 airplane. The efficiency of simulation and comparison with experimental results are given in tabular form.


1976 ◽  
Vol 98 (2) ◽  
pp. 143-150 ◽  
Author(s):  
R. W. Doll ◽  
C. D. Mote

The longitudinal, torsional and 2-transverse equations of motion are formulated for the titled problem through application of Hamilton’s Principle. Curvature-torsion conditions under which linear oscillation in a plane can exist are identified. The finite element method with isoparametric elements is used for discretization prior to spectra analysis. Natural frequency calculations over a range of mass transport velocities and cylinder end conditions were carried out for comparison with constant and variable curvature analyses and experiment. These results support the application of the constant curvature, inextensible centerline model for curved cylinder vibration analysis.


2013 ◽  
Vol 2013 ◽  
pp. 1-10
Author(s):  
Michał Hać

Discussion on equations of motion of planar flexible mechanisms is presented in this paper. The finite element method (FEM) is used for obtaining vibrational analysis of links. In derivation of dynamic equations it is commonly assumed that the shape function of elastic motion can represent rigid-body motion. In this paper, in contrast to this assumption, a model of the shape function specifically dedicated to the rigid-body motion is presented, and its influence on elastic motion is included in equations of motion; the inertia matrix related to the rigid-body acceleration vector depends on both shape functions of the elastic and rigid elements. The numerical calculations are conducted in order to determine the influence of the assumed shape function for rigid-body motion on the vibration of links in the case of closed-loop and open-loop mechanisms. The results of numerical simulation show that for transient analysis and for some specific conditions (e.g., starting range, open-loop mechanisms) the influence of assumed shape functions on vibration response can be quite significant.


2020 ◽  
Vol 46 (3) ◽  
pp. 175-181
Author(s):  
Marcelo Bighetti Toniollo ◽  
Mikaelly dos Santos Sá ◽  
Fernanda Pereira Silva ◽  
Giselle Rodrigues Reis ◽  
Ana Paula Macedo ◽  
...  

Rehabilitation with implant prostheses in posterior areas requires the maximum number of possible implants due to the greater masticatory load of the region. However, the necessary minimum requirements are not always present in full. This project analyzed the minimum principal stresses (TMiP, representative of the compressive stress) to the friable structures, specifically the vestibular face of the cortical bone and the vestibular and internal/lingual face of the medullary bone. The experimental groups were as follows: the regular splinted group (GR), with a conventional infrastructure on 3 regular-length Morse taper implants (4 × 11 mm); and the regular pontic group (GP), with a pontic infrastructure on 2 regular-length Morse taper implants (4 × 11 mm). The results showed that the TMiP of the cortical and medullary bones were greater for the GP in regions surrounding the implants (especially in the cervical and apical areas of the same region) but they did not reach bone damage levels, at least under the loads applied in this study. It was concluded that greater stress observed in the GP demonstrates greater fragility with this modality of rehabilitation; this should draw the professional's attention to possible biomechanical implications. Whenever possible, professionals should give preference to use of a greater number of implants in the rehabilitation system, with a focus on preserving the supporting tissue with the generation of less intense stresses.


1990 ◽  
Vol 18 (4) ◽  
pp. 216-235 ◽  
Author(s):  
J. De Eskinazi ◽  
K. Ishihara ◽  
H. Volk ◽  
T. C. Warholic

Abstract The paper describes the intention of the authors to determine whether it is possible to predict relative belt edge endurance for radial passenger car tires using the finite element method. Three groups of tires with different belt edge configurations were tested on a fleet test in an attempt to validate predictions from the finite element results. A two-dimensional, axisymmetric finite element analysis was first used to determine if the results from such an analysis, with emphasis on the shear deformations between the belts, could be used to predict a relative ranking for belt edge endurance. It is shown that such an analysis can lead to erroneous conclusions. A three-dimensional analysis in which tires are modeled under free rotation and static vertical loading was performed next. This approach resulted in an improvement in the quality of the correlations. The differences in the predicted values of various stress analysis parameters for the three belt edge configurations are studied and their implication on predicting belt edge endurance is discussed.


2018 ◽  
Vol 18 (02) ◽  
pp. 1850017 ◽  
Author(s):  
Iwona Adamiec-Wójcik ◽  
Łukasz Drąg ◽  
Stanisław Wojciech

The static and dynamic analysis of slender systems, which in this paper comprise lines and flexible links of manipulators, requires large deformations to be taken into consideration. This paper presents a modification of the rigid finite element method which enables modeling of such systems to include bending, torsional and longitudinal flexibility. In the formulation used, the elements into which the link is divided have seven DOFs. These describe the position of a chosen point, the extension of the element, and its orientation by means of the Euler angles Z[Formula: see text]Y[Formula: see text]X[Formula: see text]. Elements are connected by means of geometrical constraint equations. A compact algorithm for formulating and integrating the equations of motion is given. Models and programs are verified by comparing the results to those obtained by analytical solution and those from the finite element method. Finally, they are used to solve a benchmark problem encountered in nonlinear dynamic analysis of multibody systems.


Sign in / Sign up

Export Citation Format

Share Document