Register to reply

Stiffness matrix of bending plate

by Ronankeating
Tags: bending, matrix, plate, stiffness
Share this thread:
Ronankeating
#1
Dec20-11, 12:42 PM
P: 64
hi,

I'm studying the book for elasticity theory and I stuck at one of the equation in that book, the stiffness matrix for elements(thin or thick plate) in bending is given in that form

ke= t∫∫BT*Db*B*det(J)dζdη+∫∫BT*Ds*B*det(J)dζdη
  • B is the strain matrix(5x12) where the inner products include the shape function and derivation of shape function regarding the ζ and η.
  • Db is the bending coefficient matrix (5x5) (scalar)
  • Ds is the shear coefficient matrix (5x5) (scalar)
  • det(J) is the determinant of Jacobian matrix (2x2)
and book suggests that the equation can be solved by carrying out the Gauss integration, with 2 Gauss points for 1st term and 1 Gauss point for 2nd term, in order to avoid the shear locking phenomenon.

The question is:
Since I know that shape function(Ni = 1/4(1 + ξ*ξi )(1 + η*ηi)) is bilinear of ζ and η. B matrix is 1st order bilinear, so BT*Db*B yieleds 2nd order. Form of det(J) is also ζ and η dependant. Multiplication of all of this terms will result in at least 3 order form of equation , is that really can be solved with Gauss integration?

Regards,
Phys.Org News Partner Engineering news on Phys.org
Tricorder XPRIZE: 10 teams advance in global competition to develop consumer-focused diagnostic device
Study shows local seismic isolation and damping methods provide optimal protection for essential computing equipment
New filter technology uses inert gas to bore holes in high-quality steel
AlephZero
#2
Dec20-11, 03:10 PM
Engineering
Sci Advisor
HW Helper
Thanks
P: 7,169
Gaussian integration with n points (in one dimension) is exact for polynomials of order 2n-1.

In any case, there are often reasons NOT to want to do "exact" integration. (For example you mentioned one reason, to avoid shear locking). Gaussian integration is similar to fitting an approximate polynomial to the "exact" function by least squares, and integrating the approximate polynomial.
Ronankeating
#3
Dec20-11, 03:39 PM
P: 64
Thanks in advance,

So as a result Gauss integration can be performed?

AlephZero
#4
Dec20-11, 08:38 PM
Engineering
Sci Advisor
HW Helper
Thanks
P: 7,169
Stiffness matrix of bending plate

Well, obviously it CAN always be performed. The question is whether it is a good or bad way to evaluate the integrals. Usually, it's a good way, or at least "good enough", considering the output from the FE model is only an approximation to the "exact" solution of a continuum mechanics problem, and the continuum mechanics problem is only an approximation to the real-world situation.

There are other numerical integration methods which you might find out about later in your course, so don't get the idea that it's the ONLY way to do numerical integration for finite elements.
Ronankeating
#5
Dec21-11, 12:32 AM
P: 64
Thanks a lot


Register to reply

Related Discussions
Bending of a 304 SST plate Mechanical Engineering 13
Stiffness matrix fea plate element Mechanical Engineering 2
Stiffness matrix Linear & Abstract Algebra 1
The stiffness matrix in matrix analysis of strutures Classical Physics 0
Stiffness Influence Coefficient Fea For Plate Engineering Systems & Design 1