Unformatted text preview: 16.225  Computational Mechanics of Materials Homework assignment # 5 Handed out: 11/24/03 Due: 12/8/03 Endow your version of sumMIT with dynamics capability. To this end, implement Newmark's algorithm, both in implicit ( = 0) and explicit ( = 0) form. In all cases use a lumped mass matrix computed by the rowsum method. As an example of application, consider the square sheet problem of HA04. Take as initial conditions the final deformation obtained by static stretching, and set the initial velocities to zero. Estimate a stable time step for explicit dynamics calculations as ts = xmin /c, where xmin is one half of the minimum element side among all elements in the mesh (you may estimate this value by creating your mesh and determining the lengths of the sides by inspection) and c = k/0 . Keep the mesh and the time step constant throughout the calculations. Take the calculations up to time T = 5L/c. Consider the following combination of parameters: 1. = 1/2, = 1, t = 10ts and t = 100ts 2. = 1/4, = 1/2, t = 10ts and t = 100ts . 3. = 0, = 1/2, t = 0.2ts , t = ts and t = 10ts . Plot deformed meshes at various stages of deformation. Plot the x1 coordinate of the center point of the right edge of the plate against time. Use for comparison an accurate solution obtained, e. g., by taking = 1/4, = 1/2 and t = 0.01ts . Comment on the conclusions you draw on the accuracy, stability and dissipation properties of the specific algorithm resulting from each set of algorithm parameters. 1 ...
View
Full
Document
This note was uploaded on 02/03/2012 for the course AERO 16.225 taught by Professor Raulradovitzky during the Fall '03 term at MIT.
 Fall '03
 RaulRadovitzky
 Dynamics

Click to edit the document details