Computational materials science

Which is the correct way of thinking when using the finite element method to calculate the elastic stresses in a component with holes or notches?

  1. ABecause rounding errors accumulate as the number of elements increases, the error grows, so elements should be made as coarse as possible, even where stress concentrates
  2. BIn a linear elastic calculation, even the stress at a sharp crack tip converges to a definite finite value if the elements are fine enough, so that value can be compared with the tensile strength
  3. CIt is an approximate method, so elements are refined where stress concentrates and it is checked that the result hardly changes with element size (convergence)
Show answer and explanation

Answer: C. It is an approximate method, so elements are refined where stress concentrates and it is checked that the result hardly changes with element size (convergence)

In the finite element method, the structure is divided into small elements such as triangles and quadrilaterals, and the displacement within each element is interpolated from the nodal displacements with simple functions. The global stiffness equation [K]{u} = {F} is solved for the nodal displacements, from which strains and stresses are calculated. Because the interpolation is approximate, the basic practice is to refine the elements near holes and notches where the stress changes sharply, and to check that the values stop changing as the mesh is refined. Coarser elements normally give larger errors. On the other hand, in linear elasticity the stress at the tip of a sharp crack is theoretically infinite, so the value keeps growing as the mesh is refined and never converges. Cracks are therefore evaluated not by stress but by the stress intensity factor or the J-integral, with singular elements or similar used at the tip.

Report an error

Practice this topic Back to topic

More Computational materials science questions

  1. In materials informatics, when a machine learning model is built from data on existing materials to predict the properties of new materials, which caution is correct?
  2. Which pairing of a materials computation method with the range of space and time it can usually handle is most appropriate?
  3. Which statement about the ideas behind density functional theory (DFT) is correct?
  4. Which is a correct tendency of the results of standard DFT calculations using LDA or GGA exchange-correlation functionals?
  5. A 108-atom supercell (with periodic boundary conditions) of a metal calculated with DFT has a total energy of −405.00 eV. Removing one atom to create one vacancy and relaxing the structure gives −400.10 eV. Which is closest to the vacancy formation energy calculated as E_v = E(N−1) − [(N−1)/N]E(N)?