Solving System of Non-Linear Diff. Equations

[izotop]

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?

 

[tiny-tim]

welcome to pf!

hi izotop! welcome to pf!

(have a theta: θ )

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?

 

[JJacquelin]

Regarding the unknown function x(t), the first equation is a linear ODE.

 

[izotop]

hi izotop! welcome to pf!

(have a theta: θ )

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.

 

[tiny-tim]

hi izotop!

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! )

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 θ₀

 

[izotop]

hi izotop!


(hmm … surprisingly neat handwriting! )

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 θ₀

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 \theta' is maximum, \theta=0. Similarily, when \theta'=0, \theta is maximum.

But at time=0-, I assume \theta and \theta' 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+, \theta and \theta' 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?

 

[izotop]

Oh and I may add that on my drawing, the angle \theta 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.

 

[tiny-tim]

ah, i get it now!

(i didn't appreciate the input force u )

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 v_O x mv_m)

 

[izotop]

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.

 

[tiny-tim]

no, not T, τ (for torque ) …

you have J = ml² and Jθ'' = -lmgsinθ …

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

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