Trebuchet Range Simulation and Optimization

This submission provides two main files: treb.m and plot_treb.m.

treb.m predicts the projectile range of a simple trebuchet given its fulcrum position and release angle. Other parameters are held fixed (but can be adjusted by the user). ode45 is used to simulate the dynamic system, including event detection for projectile release and landing.

plot_treb.m samples treb.m to gather range data and then plots it to visualize range as a function of fulcrum position and release angle.

treb.m is intended to be used with an optimization algorithm (e.g., fmincon - bound constraints should be used). The range plot is useful for visualizing the nonsmoothness in objective function that results from discrete elements of the system (such as limits on how high the trebuchet mass can be raised). This nonsmoothness can cause difficulty with gradient-based optimization. Choosing a good starting point that avoids this nonsmoothness helps, but it is also possible to impose the geometric limits as optimization constraints instead of satisfying them implicitly in the simulation, eliminating nonsmoothness from the objective function.

A handwritten derivation of the model is also provided (TrebuchetDerivation.pdf) that includes both the trebuchet model and a projectile model using parameters appropriate for a golf ball.

This model corresponds to an example video created for ME 149 (Engineering System Design Optimization), a graduate course taught in the Mechanical Engineering department at Tufts University. The video can be viewed at:

http://www.youtube.com/watch?v=3QUJNEzloDo