mbs_euler_explicit.c File Reference

This file implements the functions of the Euler Explicit integration method in C (+ Eulaire method). Eulaire is a slightly different method from Euler explicit, see the implementation for more informations. More...

#include "mbs_euler_explicit.h"
#include "mbs_message.h"

Functions
void	initialize_eulerEx (MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Initialize all needed fields in MbsDirdyn for Euler Explicit integrator. More...
int	loop_eulerEx (double t0, double tf, MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Integrate the motion from t0 to tf. More...
void	finish_eulerEx (MbsData *mbs_data, MbsDirdyn *dd)
	Finalize and clear memory allocated for Euler Explicit integrator. More...
void	error_eulerEx (MbsData *mbs_data, MbsDirdyn *dd, int err)
	If an error occurs during time integration, this function prints the error message and call finish_eulerEx. More...
void	euler_explicit (double h, MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Source code for the Euler Explicit integration. More...
void	eulaire (double h, MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Source code for unknown integrator but based partially on euler explicit We call it Eulaire. More...
void	alpha_method (double h, MbsData *mbs_data, MbsDirdyn *mbs_dd)
	Source code for a simplified alpha-method We call it AlphaM. More...

Detailed Description

This file implements the functions of the Euler Explicit integration method in C (+ Eulaire method). Eulaire is a slightly different method from Euler explicit, see the implementation for more informations.

Creation date: April 2018

Author

Sebastien Timmermans

(c) Universite catholique de Louvain

Function Documentation

alpha_method()

void alpha_method	(	double	h,
		MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Source code for a simplified alpha-method We call it AlphaM.

Reference : Chen, Qiong-zhong and Acary, Vincent and Virlez, Geoffrey and Bruls, Olivier. (2013). A nonsmooth generalized-alpha scheme for flexible multibody systems with unilateral constraints. International Journal for Numerical Methods in Engineering. 96. 487-511. 10.1002/nme.4563.

Parameters

h	The integration step size
mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

error_eulerEx()

void error_eulerEx	(	MbsData *	mbs_data,
		MbsDirdyn *	dd,
		int	err
	)

If an error occurs during time integration, this function prints the error message and call finish_eulerEx.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure
err	The error number

eulaire()

void eulaire	(	double	h,
		MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Source code for unknown integrator but based partially on euler explicit We call it Eulaire.

Parameters

h	The integration step size
mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

euler_explicit()

void euler_explicit	(	double	h,
		MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Source code for the Euler Explicit integration.

Parameters

h	The integration step size
mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

finish_eulerEx()

void finish_eulerEx	(	MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Finalize and clear memory allocated for Euler Explicit integrator.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

initialize_eulerEx()

void initialize_eulerEx	(	MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Initialize all needed fields in MbsDirdyn for Euler Explicit integrator.

Parameters

mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

loop_eulerEx()

int loop_eulerEx	(	double	t0,
		double	tf,
		MbsData *	mbs_data,
		MbsDirdyn *	mbs_dd
	)

Integrate the motion from t0 to tf.

Parameters

t0	The initial time
tf	The final time
mbs_data	The computed MBS structure
mbs_dd	The associated MbsDirdyn structure

Returns

Error status, <0 in case of failure.

user loop

user loop