Stiffness Influence Coefficient Fea For Plate

AI Thread Summary
The discussion focuses on the calculation of the stiffness influence coefficient (K(IJ)) for rectangular elements in finite element analysis (FEA). It explains that while the stiffness for rod elements can be calculated using the formula AE/L, the derivation for solid and plane elements involves more complex methods such as the principle of virtual work or the Galerkin method. The stiffness matrix for a bilinear 2D solid element is derived through a series of equations that incorporate virtual displacements and interpolation functions. Additionally, there is a query about the implications of dividing a 10mm line into multiple rod or truss elements, with an emphasis on maintaining accuracy in finite element modeling. Understanding these calculations is crucial for effective FEA applications.
chandran
Messages
137
Reaction score
1
i have put this question in other forum also.

I have understood the stiffness matrix generation for a single rectangular element as below. I consider each node having 2 translational dofs

each in x and y coord system. There are 4 nodes in this rectangular plate element and hence 8 dofs overall.


K11U1+K12U2+K13U3+K14U4+K15U5+K16U6+K17U7+K18U8= F1
K21U1+K22U2+K23U3+K24U4+K25U5+K26U6+K27U7+K28U8= F2
K31U1+K32U2+K33U3+K34U4+K35U5+K36U6+K37U7+K38U8= F3
K41U1+K42U2+K43U3+K44U4+K45U5+K46U6+K47U7+K48U8= F4
K51U1+K52U2+K53U3+K54U4+K55U5+K56U6+K57U7+K58U8= F5
K61U1+K62U2+K63U3+K64U4+K65U5+K66U6+K67U7+K68U8= F6
K71U1+K72U2+K73U3+K74U4+K75U5+K76U6+K77U7+K78U8= F7
K81U1+K82U2+K83U3+K84U4+K85U5+K86U6+K87U7+K88U8= F8


where K(IJ) is the stiffness influence coefficient saying the force at ith dof due to a unit displacement at jth dof.
U1,U2,U3 ETC is the displacement at ith dof.
F1,F2,F3 ETC is the force at the ith dof.


Now i have a question. In nastran or other fea packages how this K(IJ) is calculated for rectangular elements. I know that for rod element it is

AE/L where A is the rod area,E youngs modulus and
L is the length of the rod. But for solid elements and plane elements how this is calculated.

One more question regarding rod or truss elements. There is a line of 10mm length horizontal to x axis. Why shouldn't i divide the
line into 10 rod or truss elements. What will be the error? Can i divide the line with 10 beam elements?
 
Engineering news on Phys.org
chandran said:
One more question regarding rod or truss elements. There is a line of 10mm length horizontal to x axis. Why shouldn't i divide the
line into 10 rod or truss elements. What will be the error? Can i divide the line with 10 beam elements?

I'm not sure I'm getting this 'easy' part ... can't see why you couldn't, and basic finite elements such as rods, beams etc. typically produce in practise identical results with analytical theories (as long as we're considering typical bending theories and leave nonlinearities etc. out) == extremely accurate, work well with pretty crude meshes.

chandran said:
Now i have a question. In nastran or other fea packages how this K(IJ) is calculated for rectangular elements. I know that for rod element it is

AE/L where A is the rod area,E youngs modulus and
L is the length of the rod. But for solid elements and plane elements how this is calculated.

This goes into the "heart" of FEM ... you can work with simple element formulations using rod, beam and so forth analogies, but in general it's preferable to derive the balance equations using "valid" method of derivation, which in case of FEM means either principle of virtual work, minimum of potential energy, variational derivation or the "best" of all (personally) using weak forms, such as the Galerkin method.

Considering the bilinear 2D solid element, its stiffness matrix is typically derived as follows (the notation I'm using is similar to links I've added to the end):

*Assuming infinitesimal strain - displacement relations and linear - elastic constitutive relations (plane stress or strain), write the Cauchy's equation of motion (equilibrium equation) substituting the kinematic and constitutive relationship to it.
*Using virtual displacements, write for plane elastic body :

<br /> \int_\textrm{Ve}(\sigma_{ij}\delta\epsilon_{ij}+\rho\ddot{u_{i}}\delta u_{i})dV-\int_\textrm{Ve}f_{i}\delta u_{i}dV-\int_\textrm{Se}t_{i}\delta u_{i}dS=0<br />

... so write the virtual displacement energy balance law element - wise

*Introduce the FE approximation functions / interpolations, along the lines
<br /> u=\sum_{i=1}^n u_{i}^e \psi_{i}^e(x,y)<br />

*substitute them to above and you'll end up with a typical FEA matrix equation :

<br /> [M^e]({\ddot{\Delta^e}})+[K^e](\Delta^e)-(f^e)-(q^e)=0<br />

*where the stiffness matrix you're after is:

<br /> [K^e]=h_{e} \int_\textrm{Ve} [B^e]^T [C^e] [B^e] dxdy<br />

... get the specifics of the derivation for example here :

http://www.vector-space.com/TourPE.pdf
http://www.mathsoft.cse.clrc.ac.uk/felib/Docs/html/Intro/intro-node10.html

although looking at common FEA software manuals or any book is probably the best way to go.
 
Last edited by a moderator:
Thread 'I need a concave mirror with a focal length length of 150 feet?'
I need to cut down a 3 year old dead tree from top down so tree causes no damage with small pieces falling. I need a mirror with a focal length of 150 ft. 12" diameter to 36" diameter will work good but I can't think of any easy way to build it. Nothing like this for sale on Ebay. I have a 30" Fresnel lens that I use to burn stumps it works great. Tree service wants $2000.
Hi all, i have some questions about the tesla turbine: is a tesla turbine more efficient than a steam engine or a stirling engine ? about the discs of the tesla turbine warping because of the high speed rotations; does running the engine on a lower speed solve that or will the discs warp anyway after time ? what is the difference in efficiency between the tesla turbine running at high speed and running it at a lower speed ( as fast as possible but low enough to not warp de discs) and: i...
Thread 'Where is my curb stop?'
My water meter is submerged under water for about 95% of the year. Today I took a photograph of the inside of my water meter box because today is one of the rare days that my water meter is not submerged in water. Here is the photograph that I took of my water meter with the cover on: Here is a photograph I took of my water meter with the cover off: I edited the photograph to draw a red circle around a knob on my water meter. Is that knob that I drew a red circle around my meter...

Similar threads

Back
Top