I put together a FEM solver to do some structural analysis. At the moment I am trying to add stress visualization feature to my code. The post-processing is done with Paraview and the results with which I am trying to verify my code are from the following problem:

The structure is two square plates welded on one edge together to form the figure shown in the images below. The structure is clamped at two points as shown by red circles and a uniformely distributed load is applied on the upper edge of the structure as shown below. The simulation was run once using structured mesh and second time using unstructured mesh the results using structured mesh are somehow reasonable even though not completely satisfactory but the results from unstructured mesh are absolutely off base!

At a first glimpse what could be done wrong in transition from structured mesh to unstructured mesh (from coding point of view) that results in such discrepancy?

Stress visualization for an structured mesh using Paraview

Stress visualization for an unstructured mesh using Paraview

  • 1
    $\begingroup$ We can't know -- it may be that you have a bug in your computation of the displacement field, or that you have a bug computing the stresses from the displacements. $\endgroup$ – Wolfgang Bangerth May 5 '20 at 18:37
  • $\begingroup$ I would suggest that you start with the simplest problem, i.e., a problem with uniaxial stress. I agree with Wolfgang, you don't provide enough information so we don't know why the results don't match. $\endgroup$ – nicoguaro May 6 '20 at 17:07
  • $\begingroup$ @nicoguaro I suspected that the information I provided may not be sufficient, I will check the displacement field calculation to see if the problem comes from there. Thanks for the comments. $\endgroup$ – Dude May 12 '20 at 6:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.