# FEM derivation for plates on elastic foundation ?

by Ronankeating
Tags: derivation, elastic, foundation, plates
 P: 64 Hi All, I'm familiar with the theory of how the soil can represented by means of FEM and how soil-structure interaction is being established. Most of books are pointing out at the Winkler -springs implemented mostly for beams and quadrilaterals (4 node) FEM elements. There are obviously other methods e.g. two-parameter modified Vlasov foundation which in my idea is best approximation. If it isn't the best, at least it's better than Winkler. The approach for deriving the stiffness matrices for Vlasov foundaiton are identical to those of Reissner-Mindlin plates, where transverse shear is being taken into account. In order to, represent the vertical and horizontal spring coeff. of soil, corresponding spring coefficients should be found through strain energy equation. The displacements are given as : u(x,y,z), v(x,y,z), w(x,y,z)=w(x,y)*$\phi(z)$ Where $\phi(z)$ is mode shape but not the one obtained from eigen-mode analysis, simply, mode shape that gives the variation of the deflection in the z-direction. For example $\phi(0)$=1.0, $\phi(H)$=0 Obtainig the spring coefficients is: for vertical spring coeffcient k=$\int^{H}_{0}E*(∂\phi/∂z)^2*dz$ for horizontal spring coefficient 2t = $\int^{H}_{0}G*(\phi)^2*dz$ Stiffness matrix for vertical direction is formed by: Minimizing the strain energy by respect to each component of displacement vector. (Uk)e=1/2$\int[w(x,y)^T*k*w(x,y)*dA]$ and [Kek]ij=∂(Uk)e2 / ( ∂2wi*∂2wj) similar procedure is valid horizontal stiffness matrix Having those in hand, I wish to addopt that for triangular FEM element which has the shape function as : Ni =ai +bi*x+ci*y where ai, bi, and ci are known from nodal coordinates. My question is, how can that be addopted to the energy equations and how am I supposed to get derivates for dz where I don't have the any representation for Z-direction in my shape function for triangular element? Your helps will be appreciated! Regards,