Is it reasonable to average von Mises stress plots at discontinuities?
Is it reasonable to average von Mises stress plots at discontinuities?
(OP)
Suppose I have a FEM with two bodies. One body represents the washer underneath a bolt head; it is circular and small relative to the other body, which is a large, thick, flat plate. The washer is bonded to the plate in a circular region with compatible meshes. The edge of the bonded contact is a discontinuity, and of course the stresses peak at this discontinuity.
In one example, the peak stress observed at the discontinuity exceeds the material allowable. Looking at the stress plot, it appears to me that the peak stress is artificial, due to the discontinuity caused by the boundary condition. I would also argue that if this stress were real, localized yielding would occur and would redistribute the stress until the system was in equilibrium. I don't think "localized yielding will occur" is going to fly with my customer, regardless of whether that is accurate.
One method proposed to resolve this is to use a more realistic contact, such as a spring element with friction between the bodies and no penetration. This of course results in a non-linear model, which takes far longer to solve. Suppose for this discussion that non-linear analyses are prohibitively expensive to run.
A simpler method to resolve the peak stress that has been proposed in our group is to average the stress values in a small area around the peak stress. For example, if the observed peak is at a node, we take the stress values from the surrounding nodes at a distance of 3 element edges away. I would argue this is fairly well supported by the Saint-Venant principle, where stresses are equivalent at sufficient distances from the contact, as long as the load is statically equivalent. This raises the question how far away is sufficient? How far away is too far? Or am I completely misunderstanding the Saint-Venant principle and applying it inappropriately?
Now my core questions for you all are:
1) is it reasonable to average von Mises stresses at nodes to reduce the artificially high stresses observed at boundary condition discontinuities?
2) if the above is reasonable, is an arithmetic mean the correct measure? or would RMS, geometric mean, or something else be more appropriate?
3) are there alternative methods for overcoming artificially high stresses at boundary condition discontinuities (aside from running non-linear studies)?
Thank you.
In one example, the peak stress observed at the discontinuity exceeds the material allowable. Looking at the stress plot, it appears to me that the peak stress is artificial, due to the discontinuity caused by the boundary condition. I would also argue that if this stress were real, localized yielding would occur and would redistribute the stress until the system was in equilibrium. I don't think "localized yielding will occur" is going to fly with my customer, regardless of whether that is accurate.
One method proposed to resolve this is to use a more realistic contact, such as a spring element with friction between the bodies and no penetration. This of course results in a non-linear model, which takes far longer to solve. Suppose for this discussion that non-linear analyses are prohibitively expensive to run.
A simpler method to resolve the peak stress that has been proposed in our group is to average the stress values in a small area around the peak stress. For example, if the observed peak is at a node, we take the stress values from the surrounding nodes at a distance of 3 element edges away. I would argue this is fairly well supported by the Saint-Venant principle, where stresses are equivalent at sufficient distances from the contact, as long as the load is statically equivalent. This raises the question how far away is sufficient? How far away is too far? Or am I completely misunderstanding the Saint-Venant principle and applying it inappropriately?
Now my core questions for you all are:
1) is it reasonable to average von Mises stresses at nodes to reduce the artificially high stresses observed at boundary condition discontinuities?
2) if the above is reasonable, is an arithmetic mean the correct measure? or would RMS, geometric mean, or something else be more appropriate?
3) are there alternative methods for overcoming artificially high stresses at boundary condition discontinuities (aside from running non-linear studies)?
Thank you.
RE: Is it reasonable to average von Mises stress plots at discontinuities?
I know that you exclude nonlinear analyses from your considerations but I’ll just say that nowadays this type of simulation became a standard and the computational cost is usually low when compared to how valuable results you can get.
RE: Is it reasonable to average von Mises stress plots at discontinuities?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
You're looking at element results ? or nodal results ?
And the nodal results from different elements have different stresses ?
This is physically "wrong" ... stress is a continuous function.
Averaging stresses at nodes is reasonable, IMHO.
But "surely" you're plotting the results in a post processor (rather than looking at reams of numbers) ?
doesn't your post-processor do this "automatically" ?
Normally you plot either element centroid results or averaged nodal results.
another day in paradise, or is paradise one day closer ?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
Just to be clear, in most traditional finite element methods the solution is approximated with piecewise continuous basis functions (i.e. C0 basis functions), the exception being that some shell and beam elements use smooth (C1) basis functions (e.g. Cubic Hermites). This means that while the solution variable (e.g. displacement) is continuous across element boundaries, derived values such as strain/stress are discontinuous (C-1) across element boundaries and thus it is completely expected that the stress evaluated at a node will be different depending on which element you're evaluating. If you have a node shared by 8 elements, there will be 8 equally-valid stress values at the node. Physically "wrong" but variationally "correct" and satisfies the best-approximation property for the given solution space. Codes like Abaqus automatically use volume-weighted averaging for reporting stress (and other derived values) at nodes and for contour plots.
RE: Is it reasonable to average von Mises stress plots at discontinuities?
another day in paradise, or is paradise one day closer ?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
In the physical world, elastic displacements (on the interior) are not merely continuous, though this is the assumption in the weak-form / FEM, and aren't even "just" smooth (which would allow stress to be continuous). Rather they are infinitely smooth (on the interior). Yet we don't dismiss the displacements computed by FEM as being "physically wrong" -- even though they technically are physically wrong.
RE: Is it reasonable to average von Mises stress plots at discontinuities?
Little from FEMs (other than ridiculously simple FEMs) is truth (or physically right).
It is just a good enough approximation.
And nodal averaged stress is similarly a good enough approximation, although there are times when stress is discontinuous (or, in the real world, very nearly so).
another day in paradise, or is paradise one day closer ?
RE: Is it reasonable to average von Mises stress plots at discontinuities?
Rick Fischer
Principal Engineer
Argonne National Laboratory
RE: Is it reasonable to average von Mises stress plots at discontinuities?
I think that whats said above by rickfischer51 is correct. Definetly use Neuber/Glinka's method to sort out the real stress (you can resort to that in order to keep the simulation linear if you cannot run a non linear one).
RE: Is it reasonable to average von Mises stress plots at discontinuities?
Sounds like you're wanting to use it to say "Stress at a distance is a good substitute for stress at this point that's badly modelled" when what it really says is "Stress at a distance is correct, regardless of the stress at this point that's badly modelled".
The stress at the distance will be right, according to St Venant's principle, but if you care about the stress at the badly modelled point, it doesn't matter any less. And it sounds like your client cares about that point (since apparently they won't accept local yielding).
You have to either;
1) Run a more realistic analysis and redesign/submit as appropriate per it's results.
2) Have enough experience to say "That's a stress singularity/geometric discontinuity and is not real" (and have your client believe you).
3) Have enough experience to say "That will yield and you don't care" (and have your client believe you).
If you don't have the experience/clout for items (2) or (3) to work for you, then your only option is (1), as it should be, because clearly, the reality is - you might be wrong. But not all is lost - treat it as a learning experience. Running the extra analysis and seeing how things change, is how you eventually get the experience/clout/number of gray hairs required to make the judgement call one way or the other, and convince people of it.