The problem to solve

by Nasser M. Abbasi  (oct 2009)

find_B_matrix_1.gif

In this solution, I start directly by solving for the vector field {u,v} and starting from the general degrees of freedom, and from it by matrix inversion, find the shape function matrix N (in terms of nodal degrees of freedom). This involves inversting a 6 by 6 matrix. But Ok, I am using a computer. By hand, I would use the method I showed in the analytical note part of this assignment which involves inverting only a 3 by 3 matrix.

find_B_matrix_2.gif

Start by defining the u and v trial functions (linear polynomials in x and y)

find_B_matrix_3.gif

set up the  u=X a equation

find_B_matrix_4.gif

find_B_matrix_5.gif

Now find the shape functions. Start by expression nodal unknowns in terms of nodal coordinates

find_B_matrix_6.gif

find_B_matrix_7.gif

find_B_matrix_8.gif

find_B_matrix_9.gif

find_B_matrix_10.gif

find_B_matrix_11.gif

find_B_matrix_12.gif

Write the u=A a equation

find_B_matrix_13.gif

find_B_matrix_14.gif

Find find_B_matrix_15.gif from the above by matrix inversion

find_B_matrix_16.gif

Now find the B matrix from the above N matrix by multiplying the N matrix by the following differetial operartors matrix

find_B_matrix_17.gif

find_B_matrix_18.gif

Now find B = oper * N

find_B_matrix_19.gif

find_B_matrix_20.gif

Factor the determinant term from the above to the outside.

find_B_matrix_21.gif

find_B_matrix_22.gif

But area of triangle is

find_B_matrix_23.gif

find_B_matrix_24.gif

Hence B matrix becomes

find_B_matrix_25.gif

find_B_matrix_26.gif

find_B_matrix_27.gif

find_B_matrix_28.gif

Spikey Created with Wolfram Mathematica 8.0