Taylor McDonnell

The code differs slightly from the paper now. While equilibrium is still enforced in the same manner, compatability equations are enforced for each element instead of at the nodes. For...

Actually, now that I think about it, the simplest solution is probably just to modify PrescribedConditions to allow prescribed displacements to be linked to the displacements of another node.

1. The dropped L term is introduced earlier: https://github.com/byuflowlab/GXBeam.jl/blob/6d9b8ad37ec6fe62173561c9cdb540462c04513b/src/element.jl#L68-L70 2. Yes, I switched to element based compatibility equations as part of a general move from element-based state variables to point-based...

I didn't do a thorough review, but it looks good. Extending it to the other analyses is trivial since it would involve the exact same operations.

You'll probably get better results if you start from a steady state solution and then gradually increase the loads.

One of the hard things about time domain simulations is getting the right initial conditions.

It's totally possible to do a time dependent loading. Just pass in the prescribed conditions as a function of time. See the wind turbine blade example or the joined wing...

To be honest, I'm not sure I have the time to invest in a bake-off. I put this package together because I needed a large deflection beam model in order...

@ChrisRackauckas That's good to know. Does the sensitivity analysis also cover non-constant mass matrices created by using DiffEqArrayOperators? I think trying to combine sensitivity analyses and non-constant mass matrices might...

> You mean backsolve adjoint? That is unlikely to be stable for any DAEs: we know that from tests on the mass matrix case. I think the more important thing...