Solution for mixed differential and algebraic equations

  • Context: Graduate 
  • Thread starter Thread starter WLamers
  • Start date Start date
  • Tags Tags
    Differential Mixed
Click For Summary
SUMMARY

The discussion centers on solving a mixed differential and algebraic equation (DAE) system related to a mechanical model involving springs and dampers. The equations of motion presented include nonlinear terms and a sign function, complicating the numerical integration process. The user seeks guidance on transforming the DAE into a set of ordinary differential equations (ODEs) without losing essential system characteristics, particularly the spring and friction damper dynamics. The challenge lies in managing the sign function associated with the velocity of the third element.

PREREQUISITES
  • Understanding of Differential Algebraic Equations (DAEs)
  • Familiarity with numerical integration techniques
  • Knowledge of mechanical systems involving springs and dampers
  • Proficiency in using ODE solvers, such as MATLAB's ODE toolbox
NEXT STEPS
  • Research methods for converting DAEs to ODEs in mechanical systems
  • Explore numerical integration techniques for DAEs, focusing on MATLAB or Python libraries
  • Study the implications of nonlinear terms in differential equations
  • Investigate the use of sign functions in dynamic systems and their numerical treatment
USEFUL FOR

Mechanical engineers, applied mathematicians, and researchers working on dynamic systems involving mixed differential and algebraic equations, particularly those interested in numerical methods for solving complex mechanical models.

WLamers
Messages
4
Reaction score
0
I'm confused finding a solution (numerical integration is ok) for the following model:

http://imageshack.com/a/img27/7080/brmb.jpg

Free body diagrams:

http://imageshack.com/a/img22/9523/r0is.jpg And the equations of motion:

\begin{align}
k_{s,1}x_1+k_{s,2}x_1^3+k_m(x_1-x_2)+k_f(x_1-x_3) &= F\\
d_m\dot{x}_2-k_m(x_1-x_2) &= 0\\
F_fsign(\dot{x}_3)-k_f(x_1-x_3) &= 0
\end{align}

How can I solve these equations, preferably using a ODE solver. Problem is the lack of first derivatives. Substitution is possible when the lowest element is not present. But with this system I cannot find a solution. Help appreciated!
 
Last edited by a moderator:
Physics news on Phys.org
Nobody?

I guess this is a DAE problem (Differential Algebraic Equation) and will not be easy to solve with numerical integration. Maybe it will be easier to just modify the problem in a way I end up with just a set of ordinary differential equations. But this requires to seprate the lower spring + friction damper (because of the $$sign(\dot{x}_3)$$, which in turn results in a different than intented behaviour. Maybe someone has a hint to end up with ordinary differential equations without having to sacrifice the sping+friction damper element? Thanks!
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K