Dynamic Simulation of Dielectric Elastomer Actuated Multibody Systems

Author(s):  
T. Schlögl ◽  
S. Leyendecker

A three-dimensional electro-mechanically coupled finite element model for dielectric elastomers is used to actuate multibody systems. This setting allows exploring the complex behavior of humanoid robots that are driven by artificial muscles instead of electrical drives. The coupling between the finite element muscle model and the rigid bodies is formulated at configuration level, where Lagrange multipliers account for constraint forces, leading to differential algebraic equations of index-3. A well-chosen set of redundant configuration variables for the multibody system avoids any rotational degrees of freedom and leads to linear coupling constraints. As a result, the coupling between the artificial muscles and the multibody system can be formulated in a very modular way that allows for easy future extension. The applied structure preserving time integration scheme provides excellent long time energy behavior. In addition, the index-3 system is solved directly with numerical accuracy, avoiding index reduction approximations.

Author(s):  
Stefan Reichl ◽  
Wolfgang Steiner

This work presents three different approaches in inverse dynamics for the solution of trajectory tracking problems in underactuated multibody systems. Such systems are characterized by less control inputs than degrees of freedom. The first approach uses an extension of the equations of motion by geometric and control constraints. This results in index-five differential-algebraic equations. A projection method is used to reduce the systems index and the resulting equations are solved numerically. The second method is a flatness-based feedforward control design. Input and state variables can be parameterized by the flat outputs and their time derivatives up to a certain order. The third approach uses an optimal control algorithm which is based on the minimization of a cost functional including system outputs and desired trajectory. It has to be distinguished between direct and indirect methods. These specific methods are applied to an underactuated planar crane and a three-dimensional rotary crane.


Author(s):  
Keisuke Kamiya ◽  
Yusaku Yamashita

The governing equations of multibody systems are, in general, formulated in the form of differential algebraic equations (DAEs) involving the Lagrange multipliers. For efficient and accurate analysis, it is desirable to eliminate the Lagrange multipliers and dependent variables. Methods called null space method and Maggi’s method eliminate the Lagrange multipliers by using the null space matrix for the constraint Jacobian. In previous reports, one of the authors presented methods which use the null space matrix. In the procedure to obtain the null space matrix, the inverse of a matrix whose regularity may not be always guaranteed. In this report, a new method is proposed in which the null space matrix is obtained by solving differential equations that can be always defined by using the QR decomposition, even if the constraints are redundant. Examples of numerical analysis are shown to validate the proposed method.


Author(s):  
Alessandro Fumagalli ◽  
Pierangelo Masarati ◽  
Marco Morandini ◽  
Paolo Mantegazza

This paper discusses the problem of control constraint realization applied to generic under-actuated multibody systems. The conditions for the realization are presented. Focus is placed on the tangent realization of the control constraint. An alternative condition is formulated, based on the practical observation that Differential-Algebraic Equations (DAE) need to be integrated using implicit algorithms, thus naturally leading to the solution of the problem in form of matrix pencil. The analogy with the representation of linear systems in Laplace’s domain is also discussed. The formulation is applied to the solution of simple, yet illustrative problems, related to rigid and deformable bodies. Some implications of considering deformable continua are addressed.


2019 ◽  
Vol 2019 ◽  
pp. 1-11
Author(s):  
Bowen Li ◽  
Jieyu Ding ◽  
Yanan Li

An L-stable method over time intervals for differential-algebraic equations (DAEs) of multibody system dynamics is presented in this paper. The solution format is established based on equidistant nodes and nonequidistant nodes such as Chebyshev nodes and Legendre nodes. Based on Ehle’s theorem and conjecture, the unknown matrix and vector in the L-stable solution formula are obtained by comparison with Pade approximation. Newton iteration method is used during the solution process. Taking the planar two-link manipulator system as an example, the results of L-stable method presented are compared for different number of nodes in the time interval and the step size in the simulation, and also compared with the classic Runge-Kutta method, A-stable method, Radau IA, Radau IIA, and Lobatto IIIC methods. The results show that the method has the advantages of good stability and high precision and is suitable for multibody system dynamics simulation under long-term conditions.


Author(s):  
Sami Holopainen ◽  
Mathias Wallin

The constitutive model for glassy polymers proposed by Arruda and Boyce (BPA model) is reviewed and compared to experimental data for long-term loading. The BPA model has previously been shown to capture monotonic loading accurately, but for unloading and long-term behavior, the response of the BPA model is found to deviate from experimental data. In the present paper, we suggest an efficient extension that significantly improves the predictive capability of the BPA model during unloading and long-term recovery. The new, extended BPA model (EBPA model) is calibrated to experimental data of polycarbonate (PC) in various loading–unloading situations and deformation states. The numerical treatment of the BPA model associated with the finite element analysis is also discussed. As a consequence of the anisotropic hardening, the plastic spin enters the model. In order to handle the plastic spin in a finite element formulation, an algorithmic plastic spin is introduced. In conjunction with the backward Euler integration scheme use of the algorithmic plastic spin leads to a set of algebraic equations that provides the updated state. Numerical examples reveal that the proposed numerical algorithm is robust and well suited for finite element simulations.


Author(s):  
Alfonso Callejo ◽  
Valentin Sonneville ◽  
Olivier A. Bauchau

The gradient-based design optimization of mechanical systems requires robust and efficient sensitivity analysis tools. The adjoint method is regarded as the most efficient semi-analytical method to evaluate sensitivity derivatives for problems involving numerous design parameters and relatively few objective functions. This paper presents a discrete version of the adjoint method based on the generalized-alpha time integration scheme, which is applied to the dynamic simulation of flexible multibody systems. Rather than using an ad hoc backward integration solver, the proposed approach leads to a straightforward algebraic procedure that provides design sensitivities evaluated to machine accuracy. The approach is based on an intrinsic representation of motion that does not require a global parameterization of rotation. Design parameters associated with rigid bodies, kinematic joints, and beam sectional properties are considered. Rigid and flexible mechanical systems are investigated to validate the proposed approach and demonstrate its accuracy, efficiency, and robustness.


Sign in / Sign up

Export Citation Format

Share Document