23.12.2014 Views

OCTOBER 19-20, 2012 - YMCA University of Science & Technology

OCTOBER 19-20, 2012 - YMCA University of Science & Technology

OCTOBER 19-20, 2012 - YMCA University of Science & Technology

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Proceedings <strong>of</strong> the National Conference on<br />

Trends and Advances in Mechanical Engineering,<br />

<strong>YMCA</strong> <strong>University</strong> <strong>of</strong> <strong>Science</strong> & <strong>Technology</strong>, Faridabad, Haryana, Oct <strong>19</strong>-<strong>20</strong>, <strong>20</strong>12<br />

(FEA) on other hand can be used for solving such problems by using available computer programs such as<br />

described in [15-16].<br />

2. Methodology<br />

The Finite Element method is used to solve physical problems in engineering analysis and design. The<br />

idealization <strong>of</strong> physical problem to mathematical model requires certain assumptions that together lead to<br />

differential equation governing the mathematical model. The FEA solves this mathematical model. Here the<br />

mass matrix formulation is carried out. A two-element cantilever beam is used in order to develop the consistent<br />

mass matrix. The six degrees <strong>of</strong> freedom lumped mass is used for constructing the lumped mass matrix The<br />

global mass matrix is built up as an assemblage <strong>of</strong> element mass matrices. A method analogous to static<br />

condensation, Guyan reduction, is used to reduce the size <strong>of</strong> the two-element cantilever problem. The model is<br />

then solved for its eigenvalues using Guyan reduction.<br />

Element Stiffness Matrix<br />

The element stiffness matrix is developed by using basic strength <strong>of</strong> materials techniques to analyse the forces<br />

required to displace each degree <strong>of</strong> freedom a unit value in the positive direction:<br />

Beam Node Definitions<br />

The two-beam elements are made identical, with the same E, I and length; the global stiffness matrix can then be<br />

rewritten as:<br />

⎡ 24 –12 6 ⎤<br />

0<br />

k<br />

g<br />

⎢ 3<br />

1<br />

⎢<br />

⎢ 0<br />

= ⎢<br />

⎢ –12<br />

⎢ 3<br />

1<br />

⎢ 6<br />

⎢ 2<br />

⎣ 1<br />

8<br />

1<br />

– 6<br />

2<br />

1<br />

2<br />

1<br />

3<br />

1<br />

– 6<br />

2<br />

1<br />

12<br />

3<br />

1<br />

– 6<br />

2<br />

1<br />

2<br />

1 ⎥<br />

2 ⎥<br />

⎥<br />

1 ⎥<br />

– 6⎥<br />

2<br />

1 ⎥<br />

4 ⎥<br />

⎥<br />

1 ⎦<br />

For solving the dynamics <strong>of</strong> the cantilever beam, a mass matrix is developed to complete the equations <strong>of</strong><br />

motion. For a beam finite element, there are a number <strong>of</strong> different mass matrix formulations as discussed in [13-<br />

15] out <strong>of</strong> which Consistent mass – distributed mass is used.<br />

391

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!