Averaging results across materials
Averaging results across materials
(OP)
Hi all,
I’m trying to convince myself that I can average von mises results across material boundaries for a this particular case:
Imagine a model with an orthopaedic implant (material 1) and bone (material 2) side by side (no contact elements, just different partitions) in a linear static analysis. Now bone has this amazing ability to remodel itself and as such can grow into the orthopaedic implant creating amazing fixation. Here is the bit I’m trying to justify. Now material 1 is fully ‘bonded’ to material 2 and as such isn’t a clear boundary anymore. Would averaging the results across the materials give a fair representation of the stress in the bone of a fully ‘bonded’ implant and conversely ‘unbonded’ if I don’t average across the material or partition.
Apologies for the non-technical description of the problem but it seemed like the clearest method. Any thoughts that could help would be gratefully received.
Ideas version 9
Martin
I’m trying to convince myself that I can average von mises results across material boundaries for a this particular case:
Imagine a model with an orthopaedic implant (material 1) and bone (material 2) side by side (no contact elements, just different partitions) in a linear static analysis. Now bone has this amazing ability to remodel itself and as such can grow into the orthopaedic implant creating amazing fixation. Here is the bit I’m trying to justify. Now material 1 is fully ‘bonded’ to material 2 and as such isn’t a clear boundary anymore. Would averaging the results across the materials give a fair representation of the stress in the bone of a fully ‘bonded’ implant and conversely ‘unbonded’ if I don’t average across the material or partition.
Apologies for the non-technical description of the problem but it seemed like the clearest method. Any thoughts that could help would be gratefully received.
Ideas version 9
Martin





RE: Averaging results across materials
I hope it can help you.
RE: Averaging results across materials
Let element 1 of material 1 has a common boundary with element 2 of material 2 and let the boundary nodes are x and y. When you calculate stresses in gauss points of E1, you use M1 and extrapolate to nodes x,y. Similarly when you calculate stresses in gauss points in E2, you use M2 and extrapolate to nodes x,y. Now you have two values contributed to nodes x,y which you can average over. By this way you get average contribution of E1 and E2. Most softwares display nodal stresses by making average (ABAQUS for instance). So you just pick up nodal stresses along the boundary. Make sure your software does this.
To me this can be the only explanation (accoeding to my understanding of your problem).
RE: Averaging results across materials
RE: Averaging results across materials
I think Bhat165 clarification summed up the problem. I think I need to average across the material otherwise I do not get the contribution from the other material when viewing the results at the interface of the two materials. In an ideal world it would be nice to model the area of fusion as described by Corus but I don’t think the data is readily available to do this. Also in my problem averaging across the materials actually increases the stress in the area I am interested, not decreases as suggested.
Many thanks for your help; I have a clear argument in my head now.
Martin
RE: Averaging results across materials
you can't average von Mises stresses across material boundaries for two reasons.
1) At a fully bonded boundary it's the strains which are equal not the stresses.
2) You can't average von Mises stresses anyway, you have to average the individual stress components (x, y, z, xy, yz, xz) and then recalculate the von Mises stress.
Hope this helps
Marc
RE: Averaging results across materials
About the secdond point made by Marc: Its an important observation. You really cant average Von Misses stress directly. Its obvious.
Sudip
RE: Averaging results across materials
This will manifest itself visually as a stress discontinuity at the nodes; however this will show a correct representation of the stresses in each material as they approach the nodes. (Note that this visualized stress discontinuity is in fact correct; a strain field for the problem posed will show continuity at the interface).
What code(s) are you using for processing and post-processing?
Brad
RE: Averaging results across materials
cheers
Martin
RE: Averaging results across materials
Never, ever average results across material boundaries! The result is pure rubbish. This is one of the many great dangers of automated contour plotting codes.
Steve
RE: Averaging results across materials
Can you elaborate on your comment? Thanks.
RE: Averaging results across materials
For the same (simple uniaxial tensile) strain of 0.001, the aluminum will have a stress of 70 MPa and the steel will have a stress of 200 MPa. For similar element size, and decent material quality, these stresses when properly evaluated (with the stress discontinuity at the boundary) will be accurate for the respective materials. However, with averaging the resulting nodal stress would be (200+70)/2 = 135. This is understated for steel, and overstated for aluminum. It's a useless answer on both counts. Hence my earlier suggestion to have your averaging criteria set such that no nodal averaging occurs. I'm sorry that I don't know how to do this in IDEAS.
Brad