Grigoriy Yashin

20 Reputation

3 Badges

5 years, 47 days

MaplePrimes Activity


These are replies submitted by Grigoriy Yashin

@Preben Alsholm, I change several parts of these expressions as constants and I got solution, but I really don't understand why I can't get this solution for full system of six equations.

Simple.mw

 

@Preben Alsholm, you are wrong. I have 8 equations, but I can choose six of them for calculation of 6 variables. At the same time I guess that the remaining two variables are constants.

@Kitonum, unfortunately, it didn't help. The result is the same.

Quasi-static2.mw

@acer, You are right, for me, visually good expression is better, it is E2 (variant of @Kitonum), but I understand that from program view E1 is shorter. Nevertheless, many thanks for the answers, it helps a lot!

@acer, thank you, it works good. But I have misunderstanding, in following example I have not full simplification again:
 

restart; PDEtools[declare](`θ__si`(t), `β__si`(t), psi(t), x(t), z(t)); PDEtools[declare](prime = t); V__1six := diff(x(t), t)-`l__1с`*sin(`θ__si`(t))*cos(psi(t))*(diff(psi(t), t))-sin(psi(t))*`l__1с`*cos(`θ__si`(t))*(diff(`θ__si`(t), t)); V__1siy := -`l__1с`*sin(`θ__si`(t))*(diff(`θ__si`(t), t)); V__1siz := diff(z(t), t)+`l__1с`*sin(`θ__si`(t))*sin(psi(t))*(diff(psi(t), t))-cos(psi(t))*`l__1с`*cos(`θ__si`(t))*(diff(`θ__si`(t), t)); V__1si := simplify(V__1six^2+V__1siy^2+V__1siz^2, size); V__1si := map(simplify, collect(V__1si, [diff, diff(`θ__si`(t), t), `l__1с`, `l__1с`^2]))

`l__1с`^2*(diff(psi(t), t))^2*sin(theta__si(t))^2-2*(diff(psi(t), t))*sin(theta__si(t))*`l__1с`*(cos(psi(t))*(diff(x(t), t))-(diff(z(t), t))*sin(psi(t)))+(diff(x(t), t))^2-2*(diff(x(t), t))*sin(psi(t))*cos(theta__si(t))*(diff(theta__si(t), t))*`l__1с`+(diff(z(t), t))^2-2*cos(psi(t))*cos(theta__si(t))*(diff(theta__si(t), t))*(diff(z(t), t))*`l__1с`+`l__1с`^2*(diff(theta__si(t), t))^2

(1)

NULL


Download trigon.mw

In this example Maple combined part of expression and forgot about this: 

-2*diff(x(t), t)*sin(psi(t))*cos(theta__si(t))*diff(theta__si(t), t)*`l__1с` + diff(z(t), t)^2 - 2*cos(psi(t))*cos(theta__si(t))*diff(theta__si(t), t)*diff(z(t), t)*`l__1с`

I understood, that we can manually optimize our expression, but it's not effective by using freeze for some expression parts. Is it possible to carry out it automatically by set of some commands?

@acer, thank you so much for the full answers. Problem with symbol "c" was due to of Russian language and program for automatically switching keyboard layouts. I apologize for the concern, but there is one more question. How can I simplify my expression more effectively? For example:
 

restart; PDEtools[declare](`θ__l`(t), `β__l`(t), `θ__si`(t), `β__si`(t), psi(t), x(t), z(t)); PDEtools[declare](prime = t); U__l := simplify(g*(m__2*(z-(1/2)*l__b*sin(psi(t))-l__1*sin(`θ__l`(t)+psi(t))-l__2c*sin(180-`β__l`(t)-`θ__l`(t)-psi(t)))+m__f*(z-(1/2)*l__b*sin(psi(t))-l__1*sin(`θ__l`(t)+psi(t))-l__2*sin(180-`β__l`(t)-`θ__l`(t)-psi(t))-l__fc*sin(180-`β__l`(t)-`θ__l`(t)-psi(t))^H)+m__1*(z-(1/2)*l__b*sin(psi(t))-l__1c*sin(`θ__l`(t)+psi(t))))+(1/2)*k__1*`Δx__kl`^2+(1/2)*k__2*`Δx__fl`^2, trig)

-g*m__f*l__fc*(-sin(-180+beta__l(t)+theta__l(t)+psi(t)))^H+g*(l__2*m__f+l__2c*m__2)*sin(-180+beta__l(t)+theta__l(t)+psi(t))-g*(l__1*m__2+l__1*m__f+l__1c*m__1)*sin(theta__l(t)+psi(t))-(1/2)*g*l__b*(m__1+m__2+m__f)*sin(psi(t))+z*(m__1+m__2+m__f)*g+(1/2)*k__2*`Δx__fl`^2+(1/2)*k__1*`Δx__kl`^2

(1)

``


Download PotenEnergy.mw In this expression we can take out the multipliers (m1+m2+mf) and g, also we can replace sin(180-fi) as sin(fi). I try to use factor, but it works like expand.

@Kitonum, I am very grateful to you! Could you comment why previous representation didn't work? Just for my better understanding of Maple logic.

P.S. Simplification wasn't applied, "simplify" just substitutes the expressions.  I used this:

V__2si := t -> simplify((V__2six^2 + V__2siy^2 + V__2siz^2)(t), trig)

, but it's not effective.

I plan to calculate the Lagrange equations in Maple. These equations contain many variables and require working with trigonometric expressions and differentiation. I tried to implement the approach you described, but did not get the result. I attached a file, could you tell me what I'm wrong about?

equations.mw

@tomleslie, is the "trig" in combine function a simplification procedure? Can I use several types of expression simplification?

@acer, Thank you very much for the answer with explanations!

Page 1 of 1