Solving System of Non-Linear Diff. Equations

  • Context: Graduate 
  • Thread starter Thread starter izotop
  • Start date Start date
  • Tags Tags
    Non-linear System
Click For Summary

Discussion Overview

The discussion revolves around solving a system of non-linear differential equations related to a mechanical engineering problem involving a pendulum and a block. Participants explore methods for linearization, Laplace transforms, and the interdependence of the equations governing the system's dynamics.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested
  • Mathematical reasoning

Main Points Raised

  • One participant seeks to linearize the first equation and apply the Laplace transform to both equations, followed by using Cramer's rule to solve the resulting linear system.
  • Another participant suggests solving the second equation for time, indicating it resembles simple harmonic motion, and proposes rearranging terms in the first equation.
  • Some participants discuss the implications of initial conditions being zero and how they affect the solution, particularly regarding the angle θ and its derivatives.
  • There is a mention of friction acting on the mass M and how it complicates the system's dynamics, with emphasis on the relationship between the motions of the two masses.
  • Participants express uncertainty about the application of torque and angular momentum in the context of the problem, questioning the validity of certain assumptions regarding the motion of point O.

Areas of Agreement / Disagreement

Participants express differing views on the initial conditions and their implications for solving the equations. There is no consensus on the best approach to linearize the equations or the impact of the input force on the system's behavior.

Contextual Notes

Participants note that the initial conditions and the nature of the input function may not be straightforward, leading to potential complications in solving the equations. The discussion reflects a reliance on specific assumptions that may not be universally accepted.

izotop
Messages
5
Reaction score
0
Hi,

I'm new to this forum. It looks like an amazing source of information and people here seem to be helpful and knowledgeable.

Briefly, I'm an undergrad student currently studying Mechanical Engineering at Concordia University, Montreal, Canada.

The first equation seems a bit tricky to me. In fact, I've been working for a while on these two equations without being able to find the right way for solving them.

Here's what I would like to do:
1) Linearize the first equation
2) Take the laplace transform of both equations
3) Use Cramer's rule to solve the system of linear equations
4) Take the inverse laplace of the two equations thus obtained

Is it possible to linearize the first equation?
 

Attachments

  • Sans titre.jpg
    Sans titre.jpg
    3.1 KB · Views: 511
Physics news on Phys.org
welcome to pf!

hi izotop! welcome to pf! :smile:

(have a theta: θ :wink:)

can't you just completely solve the second equation for t (it looks like ordinary shm), then chuck all the θ stuff in the first equation over to the RHS?
 
Regarding the unknown function x(t), the first equation is a linear ODE.
 


tiny-tim said:
hi izotop! welcome to pf! :smile:

(have a theta: θ :wink:)

can't you just completely solve the second equation for t (it looks like ordinary shm), then chuck all the θ stuff in the first equation over to the RHS?


I first thought of solving the second equation but since my initial conditions are all zero and there is no input function in my second equation, theta is simply going to end up being zero.

I attached picture of the actual physical problem I need to solve.
 

Attachments

  • problem Pendulum mass.jpg
    problem Pendulum mass.jpg
    17.1 KB · Views: 536
hi izotop! :wink:
izotop said:
I first thought of solving the second equation but since my initial conditions are all zero and there is no input function in my second equation, theta is simply going to end up being zero.

(hmm … surprisingly neat handwriting! :smile:)

so it's a pendulum hanging from a block free to roll along a horizontal rail in the plane of the swing …

but the initial conditions aren't all zero … when θ' = 0, θ is some maximum θ0 :wink:
 
tiny-tim said:
hi izotop! :wink:


(hmm … surprisingly neat handwriting! :smile:)

so it's a pendulum hanging from a block free to roll along a horizontal rail in the plane of the swing …

but the initial conditions aren't all zero … when θ' = 0, θ is some maximum θ0 :wink:



Regarding the hand-writing, if it's not clear on paper, it won't be clear in my mind!

Now. The problem.
You're right, it's a pendulum "m" free to rotate about point O located on Mass "M". The mass "M" is not quite "free" to move since there is some friction acting on it (Friction = -bx'). Other than that, as you can see, there is the inertial force, the tension of the cable and the input force u(t) acting on it.

I know what you mean: when [tex]\theta[/tex]' is maximum, [tex]\theta[/tex]=0. Similarily, when [tex]\theta[/tex]'=0, [tex]\theta[/tex] is maximum.

But at time=0-, I assume [tex]\theta[/tex] and [tex]\theta[/tex]' to be zero since the system is in static equilibrium.

At time=0+, the input function u(t) which I assume to be the unit step function comes into play. But still, at t=0+, [tex]\theta[/tex] and [tex]\theta[/tex]' are zero.

It seems like I can't solve the two equations independently since the motions of the two masses are definitely linked.

Does that make sense?
 
Oh and I may add that on my drawing, the angle [tex]\theta[/tex] looks like its initial value is NOT zero because I put a certain angle there. I deliberately put an angle to be able to build up my equations.
 
ah, i get it now! :smile:

(i didn't appreciate the input force u :redface:)

yes initially θ = θ' = 0

i think the problem is that you can't do τ = dL/dt about O because O isn't moving parallel to the centre of mass (there'd be an error term of vO x mvm)
 
T = dL/dt? I don't get that part. L is the length of the cable joining the two masses. T is the tension in the cable.
 
  • #10
no, not T, τ (for torque :wink:) …

you have J = ml2 and Jθ'' = -lmgsinθ …

that's τ = dL/dt with angular momentum L = Jθ' …

it doesn't work for moments (torques) about O.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
Replies
3
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 36 ·
2
Replies
36
Views
6K
  • · Replies 7 ·
Replies
7
Views
2K