       Re: PSE Example 8.6 via Reduce

• To: mathgroup at smc.vnet.net
• Subject: [mg112882] Re: PSE Example 8.6 via Reduce
• From: Eduardo Cavazos <wayo.cavazos at gmail.com>
• Date: Tue, 5 Oct 2010 05:31:45 -0400 (EDT)
• References: <i7mq41\$ovh\$1@smc.vnet.net> <i8c8t0\$g8f\$1@smc.vnet.net>

```On Oct 4, 5:05 am, Eduardo  Cavazos <wayo.cava... at gmail.com> wrote:
> On Sep 26, 1:44 am, Eduardo  Cavazos <wayo.cava... at gmail.com> wrote:
>
>
>
> > Here's example 8.6 from PSE :
>
> >    http://i.imgur.com/HJLqs.png
>
> > Here I'll be working with the "exercise" portion of the example; i.e.
> > the case where there is friction between point A and point B.
>
> > Use Reduce on the set of equations which describe the system between
> > point A and point B:
>
> >    http://i.imgur.com/PLWza.png
>
> > The result shows that vB is about 12.9.
>
> > Now use Reduce on the set of equations which describe the system
> > between point B and point C. Also, plug in our value of 'vB' manually:
>
> >    http://i.imgur.com/rujUz.png
>
> > We get the correct answer for dBC; i.e. about 40.3.
>
> > Above I solved the problem in "piecemeal" fashion; first I used Reduce
> > on the system from A to B and used the result from that to work with
> > the system from B to C, also using Reduce.
>
> > I don't have alot of experience with Reduce, but it seems like, in
> > theory, I should be able to throw the whole set of equations at Reduce
> > in one step, instead of in two steps. However, when I try to do so,
> > Mathematica seems to have trouble with it (at least on this system;
> > Pentium M with about 1.5GB RAM). So something along the lines of:
>
> >    http://i.imgur.com/2swfE.png
>
> > So my questions are... Is this too much for Reduce to handle? Would
> > you recommend the piecemeal method instead?
>
> In particular, the trouble is that the computation ran over 3 hours
> without completing. When I split up the problem, each piece takes a
> few seconds.

Here's the code for the non piecemeal approach:

Timing[Reduce[{
Subscript[\[CapitalEpsilon], A] ==
Subscript[K, A] + Subscript[U, A],
Subscript[K, A] == 1/2 m*Subscript[v, A]^2,
Subscript[U, A] == m*g*Subscript[y, A],
Subscript[\[CapitalEpsilon], B] ==
Subscript[K, B] + Subscript[U, B],
Subscript[K, B] == 1/2 m*Subscript[v, B]^2,
Subscript[U, B] == m*g*Subscript[y, B],
Subscript[\[CapitalDelta]E, AB] ==
Subscript[\[CapitalEpsilon], B] - Subscript[\[CapitalEpsilon],
A], Subscript[\[CapitalDelta]E,
AB] == -Subscript[\[CapitalEpsilon], frictionAB],
Subscript[\[CapitalEpsilon], frictionAB] ==
Subscript[f, kAB]*Subscript[d, AB],
Subscript[f, kAB] == Subscript[\[Mu], k]*Subscript[n, AB],
Subscript[n, AB] == m*g*Cos[\[Theta]],
Sin[\[Theta]] == Subscript[y, A]/Subscript[d, AB],
Subscript[d, AB] > 0,
Subscript[\[CapitalEpsilon], frictionAB] > 0,
Subscript[v, B] > 0,
Subscript[\[CapitalEpsilon], B] ==
Subscript[K, B] + Subscript[U, B],
Subscript[K, B] == 1/2 m*Subscript[v, B]^2,
Subscript[U, B] == m*g*Subscript[y, B],
Subscript[\[CapitalEpsilon], C] ==
Subscript[K, C] + Subscript[U, C],
Subscript[K, C] == 1/2 m*Subscript[v, C]^2,
Subscript[U, C] == m*g*Subscript[y, C],
Subscript[\[CapitalDelta]E, BC] ==
Subscript[\[CapitalEpsilon], C] - Subscript[\[CapitalEpsilon],
B], Subscript[\[CapitalDelta]E,
BC] == -Subscript[\[CapitalEpsilon], frictionBC],
Subscript[\[CapitalEpsilon], frictionBC] ==
Subscript[f, kBC]*Subscript[d, BC],
Subscript[f, kBC] == Subscript[\[Mu], k]*Subscript[n, BC],
Subscript[n, BC] == m*g,
Subscript[d, BC] > 0,
Subscript[\[CapitalEpsilon], frictionBC] > 0
} /. {
Subscript[y, A] -> 20, Subscript[v, A] -> 0,
Subscript[y, B] -> 0,
Subscript[y, C] -> 0, Subscript[v, C] -> 0,
Subscript[\[Mu], k] -> 0.210, g -> 9.8, \[Theta] -> 20 Degree
}, Reals]]

Ed

```

• Prev by Date: Re: What assumptions to use to check for orthogonality
• Next by Date: Re: Mean Values
• Previous by thread: Re: PSE Example 8.6 via Reduce
• Next by thread: double loop do