Hello,

On my multibody model, I meet a mistake during the simulation which is the following :

"Constraint projection failure at 10s".

1) Do you have a idea on this error message so as to correct the settings for the simulation?

2) I would like to know more about the projection method. If, I'm not wrong, there is not a lot of documentation about the projection method used in Maplesim. 

   a) Does someone have some references (scientific papers) which present the theory used for the projection method in Maplesim?

   b) In ths sense, I would be very interesting about learning the theory of the projection method in order to be able to tune correctly the following options linked to the projection method :

- Projection Iterations
- Projection Tolerance
- Event Projection
- Event Iterations
- Event Hysteresis

Thanks a lot for your help.

Hello,

Is it possible to define the orientation of a body with a matrix transformation ?

I have seen this option for the fixed frame but not for the rigid bodies.

I would like to define the rigid bodies thanks to transformation matrix so as to make some checking.

Thank you for your help.

Hello,

In my mechanism, I use absolute coordinates.

For the rotation (namely the initial angles/ angular velocites), I didn't see the difference between Inboard and Euler.

It seems to me that in the two cases, the initial angles are defined with regard to the inertial frame and for both, different kind of rotations can be defined.

Have you some ideas on these two options ?

Thank you for your help.

I am having difficulty with the contourplot3d.  When I hit enter it comes back with a blank plot.  Either Im doing somthing wrong or my machine can't handle it.

Hi,

The system below is successfully solved by using rfk45_dae method

I want to write the Runge Kutta scheme according to rfk45 method so as to have

S'=F(x,U,S,..), U'=G(x,U,S,..), ...

The presence of diff(S,x,x) in the second equation doesn't help me to definie F and G functions.

Thanks for your help 

restart;

with(plots);

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2-M^2*(diff(f(eta), eta))+B(f(eta)*(diff(f(eta), eta, eta))*(diff(f(eta), eta))-f(eta)^2*(diff(f(eta), eta, eta, eta))) = 0;

Eq2 := (diff(theta(eta), eta, eta))/Pr+f(eta)*(diff(theta(eta), eta))-2*(diff(f(eta), eta))*theta(eta) = 0;

Pr := 1

M := 1

S := 0

epsilon := 1

blt := 10

bcs1 := f(0) = S, (D(f))(0) = epsilon, (D(f))(blt) = 0;

bcs2 := theta(0) = 1, theta(blt) = 0;

L := [0, .2, .4, .6, .8, 1.2];

for k to 6 do R := dsolve(eval({Eq1, Eq2, bcs1, bcs2}, B = L[k]), [f(eta), theta(eta)], numeric, output = listprocedure); X1 || k := rhs(R[3]); X2 || k := rhs(R[4]); Y1 || k := rhs(R[5]); Y2 || k := -rhs(R[6]) end do:

print([(X2 || (1 .. 6))(0)])

Consider the equation

eqn1:=5^(2*abs(y-1)+2) = (1/15625)*625^abs(y-1);

and note that

eval( eqn1, y=-3);
eval( eqn1, y=5);

both return valid solutions. So why does

testeq(eqn1);

throw a numeric exception.

Given the manual description of testeq(), I would expect it to return either true or FAIL - but not a numeric exceprion!!!!

Any ideas?

Hi there,

I've created some nice 3D-plots. I want to export them and insert them in LaTex. Exporting to EPS gives me a zero-byte-file but no printable data.

I've already read that I'm not the first to have this issue and there is obviously a work-around by enabling the Classic-Worksheet-mode. Unfortanely this mode appears to be not available in Maple 17 under MAC OSX. 

Is there another way to solve the problem?
Would help me a lot.

Thanks!
 

Hi there,

I used to work with Maple 15. In Maple 15 the input-font-color was a beautiful plain red. Now I have to work with Maple 17, where the font-color was changed to a way darker red. 

I found a way to change the color in a single worksheet with Format-->Character-->Color. But I want to change it globally and permanently, so that every time I create a new worksheet the color is automatically red.

Do you have any ideas how to get this managed?

Thanks a lot!

Hi, is it possible to instantly start the display of plot animations after they are generated in maple? I searched for answers in the documentation and on the internet but couldn't find a solution. I would like to skip the step of manually clicking the plot and the "play animation" button.

Hi,

  Suppose I have a file, 1.mw, which includes a subroutine

aver:=proc(a,b)

  local res;
  res:=(a+b)/2;

  return res;

end proc:

   Now I am writing another file, 2.mw. I would like to load the subroutine, aver, in 1.mw. How shall I write in 2.mw to use that subroutine in 1.mw?

Thank you very much

PS,
PiFast43 is freely available at,

http://numbers.computation.free.fr/Constants/PiProgram/pifast.html

i dont know what`s wrong with it help me guys

Hello everyone,

I want to inseret a text arrow in the graph obtained by plot command to show increasing or decreasing trends of the graph, but I could not find any option to insert the textarrow or textbox in the graph. 

Thanks in advance.