1 / 24

Finite Element Analysis Of Composite Layered Structures Connor Kaufmann – B. Sc. ‘14

Finite Element Analysis Of Composite Layered Structures Connor Kaufmann – B. Sc. ‘14 Neola Putnam – M. Eng. ‘14 Ethan Seo – M. Eng. ‘14 Ju Hwan (Jay) Shin – B. Sc. ‘14 Cornell University Sibley School of Mechanical & Aerospace Engineering

madaline-leon
Download Presentation

Finite Element Analysis Of Composite Layered Structures Connor Kaufmann – B. Sc. ‘14

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Finite Element Analysis Of Composite Layered Structures Connor Kaufmann – B. Sc. ‘14 Neola Putnam – M. Eng. ‘14 Ethan Seo – M. Eng. ‘14 Ju Hwan (Jay) Shin – B. Sc. ‘14 Cornell University Sibley School of Mechanical & Aerospace Engineering Spring 2014 – Professor N. Zabaras

  2. Objective • Develop a linear 3D finite element analysis from scratch using MATLAB. • Consider a uniaxial loading of a symmetric laminate. • Verify the results with expected results, namely the state of out-of-plane stresses near the free edges. • Observe the effect of h-refinement (convergence of results through mesh refinements). • Perform a simple sanity check by doing a force reaction balance with the applied traction. • Compare the numerical results with a commercial FE software, or Ansys Composite PrepPost (ACP). C. Kaufmann, N. Putnam, E. Seo, J. Shin

  3. Some stress contours! C. Kaufmann, N. Putnam, E. Seo, J. Shin These results refer to a cross-ply laminate.

  4. Overview of composite materials • Composite materials are commonly used in aerospace structures to minimize mass. • We considered symmetric, unidirectional, fiber-reinforced composites. • Composite lamina (sheets) can be stacked to form high strength laminates. • Laminate stack-ups are characterized by the orientation angles of the fibers, and the materials used. • Anisotropy of the laminates allows one to tailor designs for stiffness and strength in specific directions. C. Kaufmann, N. Putnam, E. Seo, J. Shin i.e.

  5. Complexities of layered structures • Composite materials can give a much more complicated mechanical response than monolithic materials. • Stress equilibrium must be satisfied in the laminate by way of interlaminar stresses. • Special care must be taken to consider the free-edge and free-corner effect in composite samples! • Non-intuitive effects, such as normal-shear coupling, can occur in anisotropic materials. • As a result, finite element analysis is often found useful for predicting the behavior of complicated composites. C. Kaufmann, N. Putnam, E. Seo, J. Shin Normal-shear coupling Delamination

  6. FEM Formulation • Pre-processing: • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Processing: • “Globalize” and assemble the local stiffness matrices and the local load vectors. • Partition and rearrange the global stiffness equation. • Solve for the nodal displacement field! • Post-processing: • Compute the strain field by applying the kinematic equation (displacement ↔ strain). • Compute the stress field by applying the constitutive equation (Hooke’s Law). C. Kaufmann, N. Putnam, E. Seo, J. Shin

  7. Pre-processing (1/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a 8-noded hexahedron, or tri-linear element (3 translational DOF per node). • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Define the coordinate axes, based on fiber orientations. • x-axis: longitudinal direction • y-axis: transverse direction • z-axis: normal direction (or thickness direction) • Specify the size dimensions.  , , and • Specify the material properties.  , , , , , , , , and • Specify the fiber orientations of the off-axis plies. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  8. Pre-processing (2/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a 8-noded hexahedron, or tri-linear element (3 translational DOF per node). • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Store the global nodes as local nodes for each element. • Adhere to the given node-numbering scheme to ensure that the determinant of the Jacobian matrix is positive. • Guarantee invertible mapping to natural coordinate system. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  9. C. Kaufmann, N. Putnam, E. Seo, J. Shin • Pre-processing (3/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Use an eight-noded hexahedron element. • Each node has three translational degrees-of-freedom. • , , and

  10. Pre-processing (4/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Bias Factor: Allows us to have more concentrated mesh density near a particular region of interest. • The spacing between the node becomes a geometrical series. C. Kaufmann, N. Putnam, E. Seo, J. Shin Finer mesh

  11. Pre-processing (5/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • The three-dimensional elasticity (stiffness) matrix is defined by applying the generalized Hooke’s Law. • Take into account the anisotropy, assuming a transversely isotropic layer. • Material nonlinearity (plasticity) is neglected! C. Kaufmann, N. Putnam, E. Seo, J. Shin

  12. Pre-processing (6/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). C. Kaufmann, N. Putnam, E. Seo, J. Shin

  13. Pre-processing (7/7) • Define the required size dimensions, material properties, and the laminate configuration. • Discretize the model into finite elements. • Consider a tri-linear hexahedron element. • Apply any bias factor when discretizing (optional). • Calculate the 3D elasticity matrix, .  i.e. • Compute the elemental matrices necessary in developing the stiffness equation. • Specify the boundary conditions (includes the external load). • Essential Boundary Condition • Natural Boundary Condition • Pressure-based load, , at C. Kaufmann, N. Putnam, E. Seo, J. Shin

  14. Processing (1/3) • “Globalize” and assemble the local stiffness matrices and the local load vectors. • Partition and rearrange the global stiffness equation. • Solve for the nodal displacement field! • The weak form of our finite element formulation is given below. • Use Gauss Quadrature rule to numerically evaluate the local integration. • Apply the transformation rule to and using the connectivity matrix. • Sum individual matrices and vectors for global assembly. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  15. Processing (2/3) • “Globalize” and assemble the local stiffness matrices and the local load vectors. • Partition and rearrange the global stiffness equation. • Solve for the nodal displacement field! • Partition the global stiffness equation into the known and unknown components. • Apply the transformation rule to rearrange them as shown below. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  16. Processing (3/3) • “Globalize” and assemble the local stiffness matrices and the local load vectors. • Partition and rearrange the global stiffness equation. • Solve for the nodal displacement field! • Solve the system of equations efficiently by using the Gaussian elimination method. • In MATLAB, a built-in function, d=K\f can be employed. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  17. Post-processing (1/2) • Compute the strain field by applying the kinematic equation (displacement ↔ strain). • Compute the stress field by applying the constitutive equation (Hooke’s Law). • The elemental strain vector can be computed as shown below. • Weighted average of the strain values, evaluated the Gauss points. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  18. Post-processing (2/2) • Compute the strain field by applying the kinematic equation (displacement ↔ strain). • Compute the stress field by applying the constitutive equation (Hooke’s Law). • The elemental stress vector can be computed as shown below. • Weighted average of the stress values, evaluated the Gauss points. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  19. Pathwise-results! C. Kaufmann, N. Putnam, E. Seo, J. Shin These results refer to a angle-ply laminate. Stress singularity! Ply interface 16800 seconds ≈ 4.5 hours!

  20. Contour Plots! C. Kaufmann, N. Putnam, E. Seo, J. Shin These results refer to a angle-ply laminate.

  21. More contours… C. Kaufmann, N. Putnam, E. Seo, J. Shin These results refer to a angle-ply laminate.

  22. Sanity Check! • A force reaction balance check would indicate that our analysis was modeled correctly! • External load ≟ Reaction load • Reaction load is extracted from nodes (DOF to be more specific), where an essential BC is specified. C. Kaufmann, N. Putnam, E. Seo, J. Shin Nodal Reaction Forces, External Pressure Load, 

  23. Error Analysis • and Energy error norms are considered. • In the below formula, is equal to 1, since a linear element is considered. C. Kaufmann, N. Putnam, E. Seo, J. Shin

  24. Comparison to Ansys? C. Kaufmann, N. Putnam, E. Seo, J. Shin Averaged over an element 5.3539e7 Pa (Ansys ACP) 5.358e7 Pa (MATLAB)

More Related