Hello,

I'm trying to solve the effects (deflictions, tensions, etc) of a load on a timoshenko beam. It uses two partial diffential equations, wherein q is the load:

> PDE1 := kappa*G*A*(diff(y(x, t), x, x)-(diff(theta(x, t), x)))-rho*A*(diff(y(x, t), t, t))-q;

> PDE2 := E*J*(diff(theta(x, t), x, x))+kappa*G*A*(diff(y(x, t), x)-theta(x, t))-rho*J*(diff(theta(y, x), t, t));

The boundary conditions are that in the corners the moments (derative of theta) en the deflections (y) are zero, therefore:

> BCS := {y(0, t) = 0, y(L, t) = 0, (D[1](theta))(0, t) = 0, (D[1](theta))(L, t) = 0};

The initial conditions are that there is no movement in the beginning, so y has to be zero everywhere, the derative of y (the velocity). Also the rotation of the normal starts of as zero, so theta at t=0 and the derative of theta also has to be zero.

> ICS := {theta(x, 0) = 0, y(x, 0) = 0, (D[2](theta))(x, 0) = 0, (D[2](y))(x, 0) = 0};

Solving this PDE system numerically for y and theta, ie:

> sol := pdsolve({PDE1, PDE2}, {BCS, ICS}, {theta(x, t), y(x, t)}, numeric, time = t, range = 0 .. L);

Gives me this error:

Error, (in pdsolve/numeric/process_IBCs) invalid initial/boundary condition: {theta(x, 0) = 0, y(x, 0) = 0, (D[2](theta))(x, 0) = 0, (D[2](y))(x, 0) = 0}

And whatever I change, I cannot get it right. I've tried numorous things, but I keep getting this error... If I loose some boundary or initial conditions (I believe you need 8), doesn't make a difference, same error...

What am I doing wrong?

Thank you!

PS: every variable which should be a number is already defined in a previous part in the maple sheet.

shear_beam_pin_pin_t.mw