Elasticity/Kinematics example 1

From testwiki
Jump to navigation Jump to search

Example 1

Take a unit cube of material. Rotate it 90 degrees in the clockwise direction around the z-axis. Calculate the strains. Discuss your results - their accuracy and the reasons for your conclusions.

Solution

The strains are related to displacements by

ϵxx=ux;ϵyy=vy;ϵzz=wz;γxy=uy+vx;γyz=vz+wy;γzx=wx+uz

Let us consider rotation about the center of the cube. Since the problem concerns a pure rotation, a cylindrical co-ordinate system is appropriate. This problem also provides us a easy way of trying out Maple. Here are the steps that you can follow to find the strains at a point in the cube.


r := sqrt(x^2+y^2);

r:=x2+y2


theta := arctan(y/x);

θ:=arctan(yx)


x1 := r*cos(theta);

x1:=x2+y21+y2x2


y1 := r*sin(theta);

y1:=x2+y2yx1+y2x2


x2 := r*cos(theta+Pi/2);

x2:=x2+y2yx1+y2x2


y2 := r*sin(theta+Pi/2);

y2:=x2+y21+y2x2


u := x2 - x1;

u:=x2+y2yx1+y2x2x2+y21+y2x2


v := y2 - y1;

v:=x2+y21+y2x2x2+y2yx1+y2x2


epsx := simplify(diff(u,x));

epsx:=x2+y2xx2+y2x2


epsy := simplify(diff(v,y));

epsy:=x2+y2xx2+y2x2


gamxy := simplify(diff(u,y) + diff(v,x));

gamxy:=0

From the above Maple calculation, and noting that there is no motion in the z direction, the strains in the cube are

ϵxx=1;ϵyy=1;ϵzz=0;γxy=0;γyz=0;γzx=0

A pure rigid body rotation should not result in any non-zero strains.

Therefore, the measure of strain we have used is not appropriate for large rigid body motions.

Template:Subpage navbar