Interface Models for Multirate Co-Simulation of Nonsmooth Multibody Systems

Author(s):  
Albert Peiret ◽  
József Kövecses ◽  
Francisco González ◽  
Marek Teichmann

Abstract Co-simulation techniques enable the coupling of physically diverse subsystems in an efficient and modular way. Complex engineering applications can be simulated in co-simulation setups, in which each subsystem is solved and integrated using numerical methods tailored to its physical behaviour. Co-simulation implies that the communication between subsystems takes place at discrete-time instants and is limited to a given set of coupling variables, while the internals of each subsystem are generally not accessible to the rest of the simulation environment. In non-iterative co-simulation schemes, this may lead to the instability of the integration. Increasingly demanding requirements in the simulation of machinery have led to the coupling, in real-time co-simulation setups, of multibody models of mechanical systems to computational representations of non-mechanical subsystems, such as hydraulics and electronics. Often, these feature faster dynamics than their mechanical counterparts, which leads to the use of multirate integration in non-iterative co-simulation environments. The stability of the integration in these cases can be enhanced using interface models, i.e., reduced representations of the multibody system, to provide meaningful input values to faster subsystems between communication points. This work describes such interface models that can be used to represent nonsmooth mechanical systems subjected to unilateral contact and friction.

Author(s):  
Albert Peiret ◽  
Francisco González ◽  
József Kövecses ◽  
Marek Teichmann

Abstract Co-simulation techniques enable the coupling of physically diverse subsystems in an efficient and modular way. Communication between subsystems takes place at discrete-time instants and is limited to a given set of coupling variables, while the internals of each subsystem remain undisclosed and are generally not accessible to the rest of the simulation environment. In noniterative co-simulation schemes, commonly used in real-time applications, this may lead to the instability of the numerical integration. The stability of the integration in these cases can be enhanced using interface models, i.e., reduced representations of one or more subsystems that provide physically meaningful input values to the other subsystems between communication points. This work describes such an interface model that can be used to represent nonsmooth mechanical systems subjected to unilateral contact and friction. The dynamics of the system is initially formulated as a mixed linear complementarity problem (MLCP), from which the effective mass and force terms of the interface model are derived. These terms account for contact detachment and stick–slip transitions, and can also include constraint regularization in case of redundancy in the system. The performance of the proposed model is shown in several challenging examples of noniterative multirate co-simulation schemes of a mechanical system with hydraulic components, which feature faster dynamics than the multibody subsystem. Using an interface model improves simulation stability and allows for larger integration step-sizes, thus resulting in a more efficient simulation.


2020 ◽  
Vol 50 (2) ◽  
pp. 143-167
Author(s):  
Jarkko Rahikainen ◽  
Francisco González ◽  
Miguel Ángel Naya ◽  
Jussi Sopanen ◽  
Aki Mikkola

Abstract The simulation of mechanical devices using multibody system dynamics (MBS) algorithms frequently requires the consideration of their interaction with components of a different physical nature, such as electronics, hydraulics, or thermodynamics. An increasingly popular way to perform this task is through co-simulation, that is, assigning a tailored formulation and solver to each subsystem in the application under study and then coupling their integration processes via the discrete-time exchange of coupling variables during runtime. Co-simulation makes it possible to deal with complex engineering applications in a modular and effective way. On the other hand, subsystem coupling can be carried out in a wide variety of ways, which brings about the need to select appropriate coupling schemes and simulation options to ensure that the numerical integration remains stable and accurate. In this work, the co-simulation of hydraulically actuated mechanical systems via noniterative, Jacobi-scheme co-simulation is addressed. The effect of selecting different co-simulation configuration options and parameters on the accuracy and stability of the numerical integration was assessed by means of representative numerical examples.


Author(s):  
Jinzhong Wang ◽  
Zheng-Dong Ma ◽  
Gregory M. Hulbert

A new gluing algorithm is presented that can be used to couple subsystem models for dynamics simulation of mechanical systems. The gluing algorithm developed relies only on information available at the subsystem interfaces. This strategy not only improves the efficiency of the algorithm, but also engenders model security by limiting model access only to the exposed interface information. These features make the algorithm suitable for a real and practical distributed simulation environment.


Author(s):  
Weitao Chen ◽  
Shenhai Ran ◽  
Canhui Wu ◽  
Bengt Jacobson

AbstractCo-simulation is widely used in the industry for the simulation of multidomain systems. Because the coupling variables cannot be communicated continuously, the co-simulation results can be unstable and inaccurate, especially when an explicit parallel approach is applied. To address this issue, new coupling methods to improve the stability and accuracy have been developed in recent years. However, the assessment of their performance is sometimes not straightforward or is even impossible owing to the case-dependent effect. The selection of the coupling method and its tuning cannot be performed before running the co-simulation, especially with a time-varying system.In this work, the co-simulation system is analyzed in the frequency domain as a sampled-data interconnection. Then a new coupling method based on the H-infinity synthesis is developed. The method intends to reconstruct the coupling variable by adding a compensator and smoother at the interface and to minimize the error from the sample-hold process. A convergence analysis in the frequency domain shows that the coupling error can be reduced in a wide frequency range, which implies good robustness. The new method is verified using two co-simulation cases. The first case is a dual mass–spring–damper system with random parameters and the second case is a co-simulation of a multibody dynamic (MBD) vehicle model and an electric power-assisted steering (EPAS) system model. Experimental results show that the method can improve the stability and accuracy, which enables a larger communication step to speed up the explicit parallel co-simulation.


1991 ◽  
Vol 56 (1) ◽  
pp. 20-41 ◽  
Author(s):  
Dieter Britz ◽  
Merete F. Nielsen

In finite difference simulations of electrochemical transport problems, it is usually tacitly assumed that λ, the stability factor Dδt/δx2, should be set as high as possible. Here, accuracy contours are shown in (nT, λ) space, where nT is he number of finite difference steps per unit (dimensionless) time. Examples are the Cottrell experiment, simple chronopotentiometry and linear sweep voltammetry (LSV) on a reversible system. The simulation techniques examined include the standard explicit (point- and box-) methods as well as Runge-Kutta, Crank-Nicolson, hopscotch and Saul’yev. For the box method, the two-point current approximation appears to be the most appropriate. A rational algorithm for boundary concentrations with explicit LSV simulations is discussed. In general, the practice of choosing as high a λ value when using the explicit techniques, is confirmed; there are practical limits in all cases.


Author(s):  
Swathi Kommamuri ◽  
P. Sureshbabu

Power system stability improvement by a coordinate Design ofThyristor Controlled Series Compensator (TCSC) controller is addressed in this paper.Particle Swarm Optimization (PSO) technique is employed for optimization of the parameterconstrained nonlinear optimization problem implemented in a simulation environment. The proposed controllers are tested on a weakly connected power system. The non-linear simulation results are presented. The eigenvalue analysis and simulation results show the effectiveness and robustness of proposed controllers to improve the stability performance of power system by efficient damping of low frequency oscillations under various disturbances.


Author(s):  
E. Bayo ◽  
J. M. Jimenez

Abstract We investigate in this paper the different approaches that can be derived from the use of the Hamiltonian or canonical equations of motion for constrained mechanical systems with the intention of responding to the question of whether the use of these equations leads to more efficient and stable numerical algorithms than those coming from acceleration based formalisms. In this process, we propose a new penalty based canonical description of the equations of motion of constrained mechanical systems. This technique leads to a reduced set of first order ordinary differential equations in terms of the canonical variables with no Lagrange’s multipliers involved in the equations. This method shows a clear advantage over the previously proposed acceleration based formulation, in terms of numerical efficiency. In addition, we examine the use of the canonical equations based on independent coordinates, and conclude that in this second case the use of the acceleration based formulation is more advantageous than the canonical counterpart.


Sign in / Sign up

Export Citation Format

Share Document