# Nonlinear finite elements/Homework 1

## Problem 1: Mathematical Model Development edit

Consider the motion of a pendulum (see Figure 1).

Assume that:

- The rod and the mass at the end are rigid.
- The mass of the rod is negligible compared to the mass of the bob.
- There is no friction at the pivot.

Find the following:

- The equation of motion of the bob (angular displacement as a function of time). State all other assumptions.
- Why is the equation of motion "nonlinear"?
- Assume that is "small". Derive the equation of motion for this situation.
- Why is the small equation of motion "linear"?
- Derive a finite element formulation for the linear equation of motion.
- Use ANSYS or some other tool of your choice (including your own code) to solve the linear equation of motion via finite elements. Compare with the exact solution. You can use any reasonable values for time ( ), mass ( ), and length ( ).
- What happens when you try to solve the nonlinear equation of motion with finite elements? (You will need to formulate the FE formulation to see the difference.)

## Problem 2: Numerical Exercise edit

In the learning resources section, you have read about the process of formulating a FE model for an axial bar under a distributed body load. Ansys 9 shows you how to solve the same problem using ANSYS. Work the example out on your own.

Next, consider a mini centrifuge of the type manufactured by Kisker Biotech (see Figure 2).

This type of mini centrifuge is widely used in the fields of genetic research and microbiology. The centrifuge is rate to run at 15,000 RPM. The allowable temperature range is -20 C to 40 C. The nominal radius is 140 mm. Assume other geometric quantities on the basis of the figure.

We want to extend the temperature range to -80 C - 100 C and the allowable RPM to 25,000.

Assume that the current material is a light metal (any other material of your choice will also do). The only requirement is that it should be cheap and stiff. Ignore bending and the weight of the sample tubes.

Find the thinnest beam of the material of your choice that can stand the axial load that is generated due to rotation at the new speed. Use ANSYS or your own code to solve the problem. State all simplifying assumptions.

We just need a ballpark value. A 1-D simulation will do.