About

In this ME 2016 - Computing Techniques project, I used the finite-difference numerical method to asses torsional vibrations in a wind turbine shaft. The MATLAB script written for this project plots twist angle as a function of position, calculates the maximum twist angle, and calculates shaft lengths for a maximum twist angle of 2 degrees.

Based on the given conditions, plots for twist angle based on shaft position were created. The plot on the left depicts the twist angle for a 2 meter long turbine shaft, and the plot on the right depicts the twist angle for a calculated 1.402 meter long turbine shaft ensuring the maximum twist is less than 2 degrees. The Finite Difference Method results are compared to the built-in BVP4c function in MATLAB in each figure.

This is a sample of the Finite Difference Method script which calculates the node locations and twist angles based on the desired number of nodes.