On 11/27/2016 12:06 PM, Mjay21 wrote:
> I am working with a modified Euler Bernoulli beam. The energy
> functional is a function of curvature (theta') and also the slope
> (tan(theta)). I want to obtain the theta profile which minimizes the
> energy. The problem is defined on the domain which spans the length of
> the beam (0 to L). The beam is inextensible so that can be included in
> the energy functional with a Lagrange multiplier. I have been looking up
> a lot of examples and documentation but I can't figure out if MATLAB
> optimization can handle this sort of problem, where the derivatives are
> also included. Any suggestions or additional information would be very
> appreciated. Thanks in advance!
Probably your energy is an integral over the length of the beam. You
will have to figure out a formula for the local energy as a function of
position, and have your objective function be an integral of that. For a
somewhat related example, see
http://www.mathworks.com/help/gads/optimize-an-ode-in-parallel.html
(though you will have to make extensive modifications to that example,
it is simply the only example I know offhand that has an integral as an
objective function).
Good luck,
Alan Weiss
MATLAB mathematical toolbox documentation
> I am working with a modified Euler Bernoulli beam. The energy
> functional is a function of curvature (theta') and also the slope
> (tan(theta)). I want to obtain the theta profile which minimizes the
> energy. The problem is defined on the domain which spans the length of
> the beam (0 to L). The beam is inextensible so that can be included in
> the energy functional with a Lagrange multiplier. I have been looking up
> a lot of examples and documentation but I can't figure out if MATLAB
> optimization can handle this sort of problem, where the derivatives are
> also included. Any suggestions or additional information would be very
> appreciated. Thanks in advance!
Probably your energy is an integral over the length of the beam. You
will have to figure out a formula for the local energy as a function of
position, and have your objective function be an integral of that. For a
somewhat related example, see
http://www.mathworks.com/help/gads/optimize-an-ode-in-parallel.html
(though you will have to make extensive modifications to that example,
it is simply the only example I know offhand that has an integral as an
objective function).
Good luck,
Alan Weiss
MATLAB mathematical toolbox documentation