A NEW STABILIZED ONE POINT INTEGRATED SHEAR ELASTIC PLATE ELEMENT
|
|
- Warren Walker
- 5 years ago
- Views:
Transcription
1 European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) R. Owen and M. Mikkola (assoc. eds.) Jyväskylä, July 2004 A NEW STABILIZED ONE POINT INTEGRATED SHEAR ELASTIC PLATE ELEMENT Friedrich Gruttmann, Werner Wagner Institut für Werkstoffe und Mechanik im Bauwesen Alexanderstr. 7, D Technical University of Darmstadt, Germany gruttmann@iwmb.tu-darmstadt.de, web page: Institut für Baustatik Kaiserstr. 2, D-763 University of Karlsruhe(TH), Germany ww@bs.uni-karlsruhe.de, web page: Key words: Hellinger Reissner variational principle, Mindlin Reissner plate theory, quadrilateral element, effective one point integration and stabilization matrix, bending patch test. Abstract. A new and efficient formulation for a quadrilateral Mindlin Reissner plate element is presented. A Hellinger Reissner functional with independent displacements, rotations and stress resultants is used in case of a linear isotropic elastic material. Within the mixed formulation the stress resultants are interpolated using five parameters for the bending moments and four parameters for the shear forces. The hybrid element stiffness matrix resulting from the stationary condition can be integrated analytically and thus is obtained by a one point integration and a stabilization matrix. The element possesses a correct rank, does not show shear locking and is applicable for the evaluation of displacements and stress resultants within the whole range of thin and thick plates.
2 INTRODUCTION In the past, considerable research efforts have been directed towards the development of efficient and reliable finite plate elements and numerous publications can be found in the literature see e.g. the textbook []. Here we mention among may others the so called DKT and DKQ elements, where the Kirchhoff constraints are only fulfilled at discrete points of a thin plate, see e.g. [2]. Most of the work has been focussed on the Mindlin Reissner model, e.g. [3,4], which by passes the difficulties caused by C requirements of the classical Kirchhoff theory, e.g. [,5]. However for the standard bilinear interpolation for the transverse displacements and rotations severe shear locking occur for thin plates. A simple method to avoid this locking behavior is the application of reduced or selective reduced integration, see e.g. [6,7], which leads then to a rank deficiency of the element stiffness matrix. Hence several authors have developed stabilization techniques to regain the correct rank, e.g. [8, 9]. These techniques have been extended and refined for different boundary value problems in [0], where stabilization matrices on basis of the enhanced strain method have been derived. A further method uses substitute shear strain fields [], subsequently extended and reformulated in [2,3] and [4 6]. For mixed hybrid models the choice of assumed internal stress fields is particularly crucial, e.g. [7]. The paper is organized as follows. The variational formulation for a linear plate accounting for transverse shear strains is based on a Hellinger Reissner functional. The interpolation functions for the displacements, strains and stress resultants are specified and explicit expressions for the element matrices are derived. The analytical integration leads to the element stiffness matrix which is obtained by one-point integration and a stabilization matrix. Several examples demonstrate the efficiency of the developed finite plate element. 2 BASIC EQUATIONS 2. Variational Formulation For a the description of a plate we introduce the domain Ω, the boundary Γ, the thickness h, a transverse load p =[p, 0, 0] T in Ω and boundary loads t =[ p, m x, m y ] T on Γ σ. The variational formulation is based on a Hellinger Reissner functional. Π HR (u, σ) = (Ω) (ε T σ 2 σt C σ) da (Ω) u T p da (Γ σ) u T t ds stat. () Here, the displacement field is u =[w, β x,β y ] T, with the transverse deflection w and the rotations β x and β y, see Fig.. Furthermore, we introduce the vector of stress resultants σ =[m x,m y,m xy,q x,q y ] T with the bending moments m x,m y,m xy and the shear forces 2
3 q x,q y. Within a Mindlin Reissner theory the strains read κ x β x, x κ y β y, y ε = 2κ xy = β x, y +β y, x. (2) γ x β x + w, x γ y β y + w, y Furthermore, the constitutive matrix for linear isotropic elasticity is introduced as [ ] ν 0 [ ] C b 0 C = 0 C s, C b = D ν ν, Cs = κgh 0 2 with the bending rigidity D = Eh3, Young s modulus E, shear modulus G, Poisson s 2 ( ν 2 ) ratio ν and shear correction factor κ = 5. The stationary condition yields 6 δπ HR (u, σ,δu,δσ) = [δε T σ + δσ T (ε C σ) δu T p] da δu T t ds =0 (4) (Ω) with virtual displacements and stresses δu =[δw, δβ x,δβ y ] T, δσ =[δm x,δm y,δm xy,δq x,δq y ] T. (Γ σ) (3) h A 4 p z D 3 w 3 y3 x3 midsurface (z=0) C z y q y m y q x m yx m x m xy z,w y x B 2 x x y Figure : Quadrilateral plate element 2.2 Finite Element Equations For a quadrilateral element we exploit the isoparametric concept with coordinates ξ, η defined in the unit square {ξ, η} [, ], see Fig., and interpolate the transverse displacements and rotations as well as the virtual quantities using bilinear functions w = N T w, β x = N T β x, β y = N T β y. (5) 3
4 Here, w, β x, β y denote the nodal displacements and rotations and N the vector of the bilinear shape functions with N = [N,N 2,N 3,N 4 ] T = a 0 + ξ a + η a 2 + ξη h a 0 = 4 a = 4 a 2 = 4 h = 4 In order to fulfill the bending patch test, see e.g. Taylor et al. [8], we approximate the shear strains with independent interpolation functions proposed in [4] as follows. (6) [ γx γ y ] [ ] = J γξ γ η where γ ξ = 2 [( η)γb ξ +(+η)γ D ξ ] γ η = 2 [( ξ)γa η +(+ξ)γ C η ] (7) with the Jacobian matrix computed with the nodal coordinates x =[x,x 2,x 3,x 4 ] T and y =[y,y 2,y 3,y 4 ] T x, ξ y, ξ x T (a + ηh) y T (a + ηh) J = = x, η y, η x T (a 2 + ξh) y T (8) (a 2 + ξh). Thus, with eq. (7) the covariant components of the shear strains are transformed to the cartesian coordinate system. The determinant of J yields det J = j 0 + ξj + ηj 2 with j 0 = (x T a )(y T a 2 ) (x T a 2 )(y T a ) j = (x T a )(y T h) (y T a )(x T h) j 2 = (y T a 2 )(x T h) (x T a 2 )(y T h). The strains at the midside nodes A, B, C, D, see Fig., are specified as follows γξ M = [x, ξ β x + y, ξ β y + w, ξ ] M M = B,D γη L = [x, η β x + y, η β y + w, η ] L L = A, C (0) (9) where the following quantities are given with the bilinear interpolation (5) βα A = (β 2 α4 + β α ), α = x, y w A, η = (w 2 4 w ) βα B = (β 2 α + β α2 ) w B, ξ = (w 2 2 w ) βα C = (β 2 α2 + β α3 ) w C, η = (w 2 3 w 2 ) βα D = (β 2 α3 + β α4 ) w D, ξ = (w 2 3 w 4 ) [ ] x r A, η = (r 2 4 r ), r = r y, η = (r 2 3 r 2 ) r B, ξ = (r 2 2 r ) r D, ξ = (r 2 3 r 4 ) () 4
5 Considering (2) and (5) - () the approximation of the strains is now obtained by ε h = Bv, B =[B, B 2, B 3, B 4 ], v =[v, v 2, v 3, v 4 ] T, (2) where v I =[w I,β xi,β yi ] T and the submatrices for bending and shear are 0 N B b I, x 0 I B I = B s, B b I = 0 0 N I, y, I Bs I = J 0 N I, y N I, x with b I = ξ I x M, ξ b 2 I = ξ I y M, ξ N I, ξ b I N I, ξ b 2 I N I, ξ N I, η b 2 I N I, η b 22 I N I, η (3) b 2 I = η I x L, η b 22 I = η I y L, η. The coordinates of the unit square are ξ I {,,, }, η I {,,, } and the allocation of the midside nodes to the corner nodes is given by (I,M,L) {(,B,A); (2,B,C); (3,D,C); (4,D,A)}. The stress field σ is interpolated as follows σ h = S β S =[ (5 5), S] β = [ β 0, β ] T (4) S = J 0 J 0 (η η) J 0 2J 0 2(ξ ξ) 0 0 J 0 2J 0 2(η η) J 0 22J 0 22(ξ ξ) 0 0 J 0 J 0 2(η η) J 0 2J 0 22(ξ ξ) J 0 (η η) J 0 2(ξ ξ) 0 0 J 0 2(η η) J 0 22(ξ ξ), (5) where the vectors β 0 and β contain 5 and 4 parameters, respectively. The transformation coefficients Jαβ 0 in (5) denote the components of the Jacobian matrix (8) evaluated at the element center (ξ = 0,η = 0) and transform the contravariant components of the stress resultant tensors to the cartesian basis system. The coefficients have to be constant in order to fulfill the patch test. The constants ξ and η which are introduced to obtain decoupled matrices denote the coordinates of the center of gravity of the element ξ = A e (Ω e) ξda= j 3 j 0 η = A e (Ω e) ηda= j 2. (6) 3 j 0 The element area is given by A e =4j 0. Inserting (2) and (5) and the corresponding equations for the virtual stresses and virtual strains into the stationary condition (4) yields δπ h HR = numel e= [ δβ δv ] T e {[ H G G T 0 5 ][ β v ] [ 0 f ]} e =0, (7)
6 where numel denotes the total number of plate elements to discretize the problem and the virtual element vectors δβ and δv, respectively. The element load vector f =[f, f 2, f 3, f 4 ] T which follows from the external virtual work is identical with a pure displacement formulation. For a constant load p one obtains f I =[f w I, 0, 0] T with fi w = A e p (a 0I + j a I + 3 j 0 3 j 2 j 0 a 2I ), (8) where a 0I,a I,a 2I are the components of the vectors defined in (6). Furthermore the matrices H and G are introduced H := S T C S da, G := S T B da. (9) (Ω e) Since all integrants in (9) involve only polynomials of the coordinates ξ and η the integration for the element matrices can be carried out analytically. The components of the matrix H are given by and H = A e C 0 with h = hb 0 0 h 0 h s (4 4) h b = 4A ef Eh 3 (J 02 + J 02 2 ) 2 h b 22 = 4A ef 22 Eh 3 (J J ) 2 h b 2 = h b 2 = 4A ef 2 Eh 3 [ (J 0 J J 0 22J 0 2) 2 ν (J 0 J 0 22 J 0 2J 0 2) 2] h s = A ef 02 (J + J2 02 ) h s 22 = A ef (J2 + J22 02 ) 3κGh 3κGh h s 2 = h s 2 = A ef 2 3κGh (J 0 J2 0 + J22J 0 2) 0 (Ω e) (20) (2) f = ( ) 2 j2 f 22 = ( ) 2 j f 2 = j j 2. 3 j 0 3 j 0 3 j 0 j 0 Furthermore the matrix G is obtained by analytical integration as follows 0 J G =[G, G 2, G 3, G 4 ] G I = A 0 J 0 2 e B 0 I g g I = I 3 A 0 J2 0 J 0 22 eγ I γ I γi 2 γi 2 γi 22 (22) 6
7 with B 0 I = B I (ξ =0,η =0)and γ I = h I j 2 a I j a 2I j 0 j 0 γi = (b I h I b I a I b 2 I a 2I )/γ I j 0 j 0 γi 2 = (b 2 I h I b 2 I a I b 22 I a 2I )/γ I j 0 j 0 γi 2 = (b 2 I h I b I a I b 2 I a 2I )/γ I j 0 j 0 j 2 j 2 j 2 j 2 γi 22 = (b 22 I h I b 2 I a I b 22 I a 2I )/γ I. j 0 j 0 j j j j (23) The parameters h I,a I,a 2I are the components of the nodal vectors defined in (6), whereas the b αβ I are defined in (4). Since the interpolation of the stress resultants are discontinuous at the element boundaries, the stress parameters are eliminated on element level β = H Gv. (24) Thus considering (20) and (22) one obtains the element stiffness matrix k e = G T H G = k 0 + k stab k IK = A e B 0T I CB 0 K + g T I h g K. (25) Here, k 0 denotes the stiffness of a one point integrated Mindlin Reissner plate element with substitute shear strains and k stab the stabilization matrix. The matrix h according to (20) consists of two submatrices of order two and thus can easily be inverted. The element possesses with three zero eigenvalues the correct rank. 3 EXAMPLES The derived element formulation has been implemented in an extended version of the general purpose finite element program FEAP, see Zienkiewicz and Taylor []. 3. Constant bending patch test First we investigate the element behaviour within a constant bending patch test as is depicted in Fig. 2. A rectangular plate of length a and width b supported at three corners is loaded by a concentrated load at the fourth corner and by bending moments at the corners. The geometrical and material data and the loading parameters are given. The solution of the problem can be computed analytically. The vertical displacement of node is w =2.48 and the bending moments m x = m y = m xy =.0are constant throughout the plate. The present element fulfills the patch test as Fig. 3 shows. 7
8 4 8 y 5 x a b 2 Node F z m x m y a = 40 E = 0 6 b = 20 ν = 0.3 h = 0. Figure 2: Rectangular plate, patch of 5 elements 9.999E-0 min.000e+00 max -.248E+0 min -.53E E E E E E E E E E E E E E-0 max Figure 3: Bending moment m x and displacement w for the present element 3.2 Clamped square plate subjected to a concentrated load The problem with geometrical and material data is defined in Fig. 4. The mesh consists of 2 2 elements over a quarter of the plate, where fourfold symmetry has been used. Here the influence of element distorsion is tested, where one inner node is moved by 0 <s<0 in x- and y-direction. An analytical Kirchhoff solution for the center deflection yields w = Fa 2 /D =, see e.g. [5]. The sensitivity of different element formulations with respect to the distorsion parameter s is depicted in Fig. 5. The DKQ element [2] behaves relatively insensitive with respect to the mesh distorsion and yields for the present coarse mesh a solution which is too weak. The results computed with the new element are slightly better than with the Bathe/Dvorkin element [5]. The clamped plate allows a calculation without stabilization matrix, since the hourglass modes are suppressed by the boundary conditions. Thus for the present example the best results are obtained with an one point integrated element U. However this is not the case for arbitrary boundary conditions. Results for the Belytschko/Tsay [8] element are similar to the element U in the recommended range for the stabilization parameter 0.02 r w A contour plot of w with a distorsion parameter s = 0 is given in Fig. 5 for the present element. 8
9 a = 00 h = F = E = 0000 ν = 0.3 a/2 s y F, w x s a/2 displacement w Figure 4: Clamped square plate subjected to a concentrated load,20,00 0,80 0,60 DKQ Present 0,40 U_ B/D 0,20 0,00 0,00 2,00 4,00 6,00 8,00 0,00 distorion s 0.000E+00 min 5.562E-02.2E-0.669E E E E E E E E-0 6.8E E E E-0 max Figure 5: Influence of mesh distorsion on center deflection and displacements for a mesh distorsion s=0 3.3 Corner supported square plate subjected to uniform load A corner supported plate with edge length 2a subjected to uniform load is discussed. Considering symmetry the mesh consists of 8 8elementsforaquarteroftheplate, see Fig. 6. The geometrical and material data are also given. An approximate ansatz according to [20] reads w(x, y) =c + c 2 x 2 + c 3 y 2 + c 4 x 4 + c 5 x 2 y 2 + c 6 y 4, (26) where the origin of the co-ordinate system lies in the center of the plate. The boundary condition of vanishing bending moments at the edges can only be fulfilled in an integral sense. The other boundary conditions and the partial differential equation can be fulfilled exactly. The constants are determined and thus for y = 0 the approximate Kirchhoff solution reads w(x, y =0)= qa4 2Eh [ 6ν 3 ν2 +( 5+4ν + ν 2 )( x a )2 +(+ ν 2 ν2 2 )(x a )4 ]. (27) The deflections w(x, y = 0) obtained with different elements are plotted in Fig. 7. The Belytschko/Tsay element [8] leads to hourglass modes for parameters r w < 0.02, optimal 9
10 a = 2 h = q = E = ν = 0.38 a y x a Figure 6: Corner supported plate results for 0.02 r w 0.05 and locking for r w > 0.05, see also [8]. The parameter r β =0.02 has been chosen constant in all cases. Furthermore the deformed mesh using the present element amplified by a factor 0 is depicted in Fig. 7. It can be seen that no hourglassing effects occur. 0,4 0,2 0,0 displacement w 0,08 0,06 0,04 0,02 B/T r_w = 0.02 B/T r_w = 0.00 B/T r_w = 0 Present Analytical (approx.) 0,00 0,00 2,00 4,00 6,00 8,00 0,00 2,00 coordinate x Figure 7: Deflection w(x, y = 0), comparison of different elements and deformed mesh For a convergence study of the center displacement the shear correction factor is increased for the present element and the Bathe/Dvorkin element to approximate the Kirchhoff solution. The results according to Table show nearly the same convergence behaviour for all three compared elements against the same value, which differs from the approximate analytical solution. 8*8 6*6 24*24 48*48 96*96 92*92 DKQ B/D (κ = 000) Present (κ = 000) analytical (approx.) Table : Corner supported plate, convergence study 0
11 3.4 Concrete plate structure subjected to uniform load A typical concrete plate structure of a building subjected to uniform load is presented in the final example, see Fig. 8. The plate of constant thickness h =0.25 m is loaded with q = 0 kn/m 2 and simply supported at the walls. A linear elastic material behaviour with E = kn/m 2 and ν = 0.2 is chosen. The plate is discretized with 3628 elements and 3790 nodes. The associated displacements w in cm are depicted in Fig. 8, where the minimum occurs around the hole. Similar results are obtained with the DKQ element, [2]. Here, a comparison of CPU-time shows that in our implementation the DKQ-element needs about 52% more time than the present element for calculating the global stiffness matrix r= min max 3.9 [m] Figure 8: Concrete plate structure: system and displacement w in [cm] 4 CONCLUSIONS The formulation of a quadrilateral plate element with three displacement degrees of freedom (transverse deflection, two rotations) at each node has been presented. The element possesses a correct rank, does not show shear locking and is applicable for the evaluation of displacements and stress resultants within the whole range of thin and thick plates. No parameters have to be adjusted to avoid shear locking or to prevent zero energy modes. The computed results are very satisfactory. The convergence behaviour for the displacements and stresses is slightly better than comparable quadrilateral assumed strain elements. However the essential advantage is the fast stiffness computation due to the analytically derived stiffness.
12 REFERENCES [] O.C. Zienkiewicz and R.L. Taylor. The Finite Element Method, Vol. 3, 5. ed., Butterworth-Heinemann, Oxford, [2] J. L. Batoz and M. B. Tahar. Evaluation of a New Quadrilateral Thin Plate Bending Element, Int. J. Num. Meth. Engng. 2, , 982. [3] E. Reissner. The effect of transverse shear deformation on the bending of elastic plates, J. Appl. Mech. 2, 69 76, 945. [4] R. D. Mindlin. Influence of rotatory inertia and shear flexural motions of isotropic elastic plates, J. Appl. Mech. 8, 3 38, 95. [5] S. P. Timoshenko and S. Woinowsky Krieger. Theory of Plates and Shells, McGraw Hill 2 nd ed, 970. [6] O. C. Zienkiewicz, J. Too and R.L. Taylor. Reduced integration techniques in general analysis of plates and shells, Int. J. Num. Meth. Engng., 3, , 97. [7] T. R. J. Hughes, R. L. Taylor and W. Kanoknukulchai. A simple and efficient finite element for plate bending, Int. J. Num. Meth. Engng.,, , 977. [8] T. Belytschko and S.-H. Tsay. A stabilization procedure for the quadrilateral plate element with one point quadrature, Int. J. Num. Meth. Engng., 9, , 983. [9] T. Belytschko and W. Bachrach. Efficient implementation of quadrilaterals with high coarse mesh accuracy, Comp. Meth. Appl. Mech. Engng., 54, , 986. [0] S. Reese and P. Wriggers. A stabilization technique to avoid hourglassing in finite elasticity, Int. J. Num. Meth. Engng., 48, 79 09, [] R. H. MacNeal. A simple quadrilateral shell element, Comput. Struct., 8, 75 83, 978. [2] T. J. R. Hughes and T. E. Tezduyar. Finite elements based upon Mindlin plate theory, with particular reference to the 4 node bilinear isoparametric element, J. Appl. Mech., 48, , 98. [3] R. H. MacNeal. Derivation of element stiffness matrices by assumed strain distribution, Nuclear Engineering Design, 70, 3 2, 982. [4] E. Dvorkin and K.-J. Bathe. A Continuum Mechanics Based Four Node Shell Element for General Nonlinear Analysis, Engineering Computations,, 77 88,
13 [5] K.-J. Bathe and E. Dvorkin. A 4 Node Plate Bending Element based on Mindlin/Reissner Theory and a Mixed Interpolation, Int. J. Num. Meth. Engng., 2, , 985. [6] E. Hinton and H. C. Huang. A family of quadrilateral Mindlin plate elements with substitute shear strain fields. Comp. Struct., 23, , 986. [7] T. H. H. Pian, D. Kang and C. Wang. Hybrid plate elements based on balanced stresses and displacements, Finite Element Methods for Plate and Shell Structures, Vol., Pineridge Press, Swansea, , 986. [8] R. L. Taylor, J. C. Simo, O. C. Zienkiewicz and A. C. H. Chan. The patch test - A condition for assessing FEM convergence, Int. J. Num. Meth. Engng., 22, 39 62, 986. [9] T. H. H. Pian and K. Sumihara. Rational approach for assumed stress finite elements. Int. J. Num. Meth. Eng., 20, , 984. [20] S. L. Lee and P. Ballesteros. Uniformly loaded rectangular plate supported at corners, Int. J. Mech. Sci., 2, 206 2,
Universität Karlsruhe (TH) Institut für Baustatik. A robust nonlinear mixed hybrid quadrilateral shell element
Universität Karlsruhe (TH) Institut für Baustatik A robust nonlinear mixed hybrid quadrilateral shell element W. Wagner, F. Gruttmann Mitteilung 3(2004) BAUSTATIK Universität Karlsruhe (TH) Institut für
More informationUncorrected Proof Available online at
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering00 (2017) 000 000 www.elsevier.com/locate/procedia Sustainable Civil Engineering Structures and Construction Materials, SCESCM
More informationMIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING
144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes
More informationFinite Element Analysis of Saint-Venant Torsion Problem with Exact Integration of the Elastic-Plastic Constitutive
Finite Element Analysis of Saint-Venant Torsion Problem with Exact Integration of the Elastic-Plastic Constitutive Equations W. Wagner Institut für Baustatik Universität Karlsruhe (TH) Kaiserstraße 12
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationUNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING
UNIVERSITY OF HAWAII COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ACKNOWLEDGMENTS This report consists of the dissertation by Ms. Yan Jane Liu, submitted in partial fulfillment
More informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationA FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS
TRENDS IN COMPUTATIONAL STRUCTURAL MECHANICS W.A. Wall, K.-U. Bletzinger and K. Schweizerhof (Eds.) c CIMNE, Barcelona, Spain 2001 A FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationDevelopment DKMQ Shell Element with Five Degrees of Freedom per Nodal
International Journal of Mechanical Engineering and Robotics Research Vol. 6, No., May 7 Development DKMQ Shell Element with Five Degrees of Freedom per Nodal Herry Irpanni, Irwan Katili, and Imam J. Maknun
More informationPlate analysis using classical or Reissner- Mindlin theories
Plate analysis using classical or Reissner- Mindlin theories L. Palermo Jr. Faculty of Civil Engineering, State Universiv at Campinas, Brazil Abstract Plates can be solved with classical or Reissner-Mindlin
More informationAdvances in Engineering Software
Advances in Engineering Software 41 (21) 712 728 Contents lists available at ScienceDirect Advances in Engineering Software journal homepage: www.elsevier.com/locate/advengsoft The quadratic MITC plate
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationINTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS
INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3
More informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationFinite element modelling of structural mechanics problems
1 Finite element modelling of structural mechanics problems Kjell Magne Mathisen Department of Structural Engineering Norwegian University of Science and Technology Lecture 10: Geilo Winter School - January,
More informationINCOMPATIBLE BUBBLES: A GENERALIZED VARIATIONAL FRAMEWORK INCLUDING INCOMPATIBLE MODES AND ENHANCED ASSUMED STRAINS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationThe problem of isotropic rectangular plate with four clamped edges
Sādhanā Vol. 32, Part 3, June 2007, pp. 181 186. Printed in India The problem of isotropic rectangular plate with four clamped edges C ERDEM İMRAK and ISMAIL GERDEMELI Istanbul Technical University, Faculty
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More informationA Simple Finite Rotation Formulation for Composite Shell Elements W. Wagner & F. Gruttmann
A Simple Finite Rotation Formulation for Composite Shell Elements W. Wagner & F. Gruttmann Institut für Baumechanik und Numerische Mechanik ABSTRACT Universität Hannover Appelstr. 9A, D-3000 Hannover 1,
More informationA smoothed finite element method for plate analysis
A smoothed finite element method for plate analysis H. Nguyen-Xuan T. Rabczuk Stéphane Bordas J.F.Debongnie September 18, 2007 Abstract A quadrilateral element with smoothed curvatures for Mindlin-Reissner
More informationCOMPARISON OF TWO MODERATELY THICK PLATE ELEMENTS BASED ON THE ABSOLUTE NODAL COORDINATE FORMULATION
MULTIBODY DYNAMICS 2009, ECCOMAS Thematic Conference K. Arczewski, J. Frączek, M. Wojtyra (eds.) Warsaw, Poland, 29 June 2 July 2009 COMPARISON OF TWO MODERATELY THICK PLATE ELEMENTS BASED ON THE ABSOLUTE
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationJournal of the Brazilian Society of Mechanical Sciences Print version ISSN
Journal of the Brazilian Society of Mechanical Sciences Print version ISSN 0100-7386 J. Braz. Soc. Mech. Sci. vol.21 no.3 Rio de Janeiro Sept. 1999 http://dx.doi.org/10.1590/s0100-73861999000300007 Geometric
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationUniversität Karlsruhe (TH) Institut für Baustatik. Shear correction factors in Timoshenko s beam theory for arbitrary shaped cross sections
Universität Karlsruhe (TH) Institut für Baustatik Shear correction factors in Timoshenko s beam theory for arbitrary shaped cross sections F. Gruttmann, W. Wagner Mitteilung 6(2001) BAUSTATIK Universität
More informationA simple locking-free discrete shear triangular plate element
Copyright 011 Tech Science Press CMES, vol.77, no.4, pp.1-38, 011 A simple locking-free discrete shear triangular plate element Y.C. Cai 1,,3, L.G. Tian 1 and S.N. Atluri 3 Abstract: A new three node triangular
More informationDYNAMIC ANALYSIS OF MODERATELY THICK DOUBLY CURVED SHELLS VIA EFFICIENT 3D ELEMENTS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More information1 Bending of beams Mindlin theory
1 BENDNG OF BEAMS MNDLN THEORY 1 1 Bending of beams Mindlin theory Cross-section kinematics assumptions Distributed load acts in the xz plane, which is also a plane of symmetry of a body Ω v(x = 0 m Vertical
More informationFolded-plate structures with openings analysis and optimization
Folded-plate structures with openings analysis and optimization Bogdan Szybiński, Andrzej P. Zieliński Institute of Mechanics & Machine Design, Cracow University of Technology Marek Karaś Institute of
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationVerification Examples. FEM-Design. version
FEM-Design 6.0 FEM-Design version. 06 FEM-Design 6.0 StruSoft AB Visit the StruSoft website for company and FEM-Design information at www.strusoft.com Copyright 06 by StruSoft, all rights reserved. Trademarks
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationA study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation
Multibody Syst Dyn (2014) 31:309 338 DOI 10.1007/s11044-013-9383-6 A study of moderately thick quadrilateral plate elements based on the absolute nodal coordinate formulation Marko K. Matikainen Antti
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The M4 Element. Figure 1: Bilinear Mindlin element.
Professor Terje Hakaas University of British Colmbia, ancover www.inrisk.bc.ca The M Element variety of plate elements exist, some being characterized as Kirchhoff elements, i.e., for thin plates, and
More informationCOUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM
VOL. 4, NO. 10, DECEMBER 009 SSN 1819-6608 006-009 Asian Research Publishing Network (ARPN. All rights reserved. COUPLED FNTE-NFNTE ELEMENTS MODELNG OF BULDNG FRAME-SOL NTERACTON SYSTEM Ramakant Agrawal
More informationEffect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test
Effect of Specimen Dimensions on Flexural Modulus in a 3-Point Bending Test M. Praveen Kumar 1 and V. Balakrishna Murthy 2* 1 Mechanical Engineering Department, P.V.P. Siddhartha Institute of Technology,
More informationStudy of a family of hybrid-trefftz folded plate p-version elements
Computer Assisted Mechanics and Engineering Sciences, 4: 453 476, 1997. Copyright c 1997 by Polska Akademia Nauk Study of a family of hybrid-trefftz folded plate p-version elements Jaroslav Jirousek LSC,
More informationA MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES
COMPUTATIONAL MECHANICS New Trends and Applications S. Idelsohn, E. Oñate and E. Dvorkin (Eds.) c CIMNE, Barcelona, Spain 1998 A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationStability of Simply Supported Square Plate with Concentric Cutout
International OPEN ACCESS Journal Of Modern Engineering Research (IJMER) Stability of Simply Supported Square Plate with Concentric Cutout Jayashankarbabu B. S. 1, Dr. Karisiddappa 1 (Civil Engineering
More informationPORTFOLIO. Advanced Written Communication in English for International Students (FLE ) DECEMBER 5, TUHIN ROY (Student ID # )
PORTFOLIO Advanced Written Communication in English for International Students (FLE 402-002) DECEMBER 5, 2016 TUHIN ROY (Student ID # 200157097) Table of Contents Reflective Analysis... 1 Error Log Project
More informationMITOCW MITRES2_002S10linear_lec07_300k-mp4
MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More information2008 by authors and 2008 Springer Science+Business Media
Antti H. Niemi, Harri Hakula, and Juhani Pitkäranta. 28. Point load on a shell. In: Karl Kunisch, Günther Of, and Olaf Steinbach (editors). Numerical Mathematics and Advanced Applications. Proceedings
More informationAN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION
AN INTEGRATED IRCHHOFF PLATE ELEMENT BY GALERIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION Ragesh.P.P, V.Mustafa 2, T.P.Somasundaran 3 (Research Scholar, National Institute of Technology Calicut,
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More informationNumerical Solution of Non-Linear Biharmonic. Equation for Elasto-Plastic Bending Plate
Applied Mathematical Sciences, Vol. 9, 5, no. 6, 769-78 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ams.5.575 Numerical Solution of Non-Linear Biharmonic Equation for Elasto-Plastic Bending Plate
More informationAN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
More informationKirchhoff Plates: Field Equations
20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationQUASI-STATIC AND DYNAMIC SIMULATION OF SHEET METAL FORMING PROCESSES USING LINEAR AND QUADRATIC SOLID- SHELL ELEMENTS
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Azores, 26-30 July 2015 PAPER REF: 5456 QUASI-STATIC AND DYNAMIC
More informationME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?
SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.
More informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationHierarchic Isogeometric Large Rotation Shell Elements Including Linearized Transverse Shear Parametrization
Proceedings of the 7th GACM Colloquium on Computational Mechanics for Young Scientists from Academia and Industry October 11-13, 2017 in Stuttgart, Germany Hierarchic Isogeometric Large Rotation Shell
More informationComputers and Structures
Computers and Structures 87 (2009) 1451 1460 Contents lists available at ScienceDirect Computers and Structures journal homepage: www.elsevier.com/locate/compstruc A triangular six-node shell element Do-Nyun
More informationCIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives
CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationAnalysis of thin plate structures using the absolute nodal coordinate formulation
345 Analysis of thin plate structures using the absolute nodal coordinate formulation K Dufva 1 and A A Shabana 2 1 Department of Mechanical Engineering, Lappeenranta University of echnology, Lappeenranta,
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationA NEW REFINED THEORY OF PLATES WITH TRANSVERSE SHEAR DEFORMATION FOR MODERATELY THICK AND THICK PLATES
A NEW REFINED THEORY OF PLATES WITH TRANSVERSE SHEAR DEFORMATION FOR MODERATELY THICK AND THICK PLATES J.M. MARTÍNEZ VALLE Mechanics Department, EPS; Leonardo da Vinci Building, Rabanales Campus, Cordoba
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationTRIANGULAR, NINE-DEGREES-OF-FREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY*
QUARTERLY OF APPLIED MATHEMATICS 303 OCTOBER, 1973 TRIANGULAR, NINE-DEGREES-OF-FREEDOM, C PLATE BENDING ELEMENT OF QUADRATIC ACCURACY* By ISAAC FRIED and SHOK KENG YANG (Boston University) Introduction.
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationRetrospectives on My Studies of Solid Mechanics (II)
Retrospectives on My Studies of Solid Mechanics (II) - a new energy method based on the stationary principle of total energy By Tadahiko Kawai + and Etsu Kazama ++ ABSTRACT A new energy method is proposed
More informationBending Analysis of Symmetrically Laminated Plates
Leonardo Journal of Sciences ISSN 1583-0233 Issue 16, January-June 2010 p. 105-116 Bending Analysis of Symmetrically Laminated Plates Bouazza MOKHTAR 1, Hammadi FODIL 2 and Khadir MOSTAPHA 2 1 Department
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
An alternative boundary element formulation for plate bending analysis J.B. Paiva, L.O. Neto Structures Department, Sao Carlos Civil Engineering School, f, BrazzY Abstract This work presents an alternative
More informationA Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
Structural Engineering and Mechanics, Vol. 21, No. 3 (2005) 315-332 315 A Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
More informationInvestigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
Print ISSN: 2322-2093; Online ISSN: 2423-6691 DOI: 10.7508/ceij.2016.02.011 Technical Note Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis
More informationSemiloof Curved Thin Shell Elements
Semiloof Curved Thin Shell Elements General Element Name Y,v,θy X,u,θx Z,w,θz Element Group Element Subgroup Element Description Number Of Nodes Freedoms Node Coordinates TSL 1 2 Semiloof 3 QSL8 7 8 1
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationENHANCED STRAIN METHODS FOR ELASTICITY PROBLEMS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationEsben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer
Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT METHOD FOR REISSNER-MINDLIN PLATES
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 1, Pages 48 62 c 2017 Institute for Scientific Computing and Information ANALYSIS OF A MIXED-SHEAR-PROJECTED QUADRILATERAL ELEMENT
More informationMechanics of Materials and Structures
Journal of Mechanics of Materials and Structures IMPROVED HYBRID ELEMENTS FOR STRUCTURAL ANALYSIS C. S. Jog Volume 5, No. 3 March 2010 mathematical sciences publishers JOURNAL OF MECHANICS OF MATERIALS
More informationA *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +"' ( + -"( (' (& -+" % '('%"' +"-2 ( -!"',- % )% -.C>K:GH>IN D; AF69>HH>6,-+
The primary objective is to determine whether the structural efficiency of plates can be improved with variable thickness The large displacement analysis of steel plate with variable thickness at direction
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
More information