Implementation of orthotropic thermal effects in tractionBoundaryGradient class
Hello,
I've modified the elasticOrthoSolidFoam solver such that it accounts for thermal residual stresses because I got oscillations in the stress and epsilon fields which seems to be introduced by the BC. The thermal contribution has to be incorporated into the tractionBoundaryGradient class in order to get the BC's right but that hasn't been done for orthotropic boundaries so far. My implementation in tractionBoundaryGradient.C (foam ext 3.1) looks like this: Code:
// Thermal effects I'm running a simulation right now to compare the results but was wondering if someone else has already tried this and can give me some feedback wether this implementation is correct or not. EDIT: Looks like the stress oscillations are gone =) However the cpu time increased by about 11%... Thanks, Alex |
Hi Alex,
As far as I can see, the way you are doing it should be correct. Is it possible for you to share your solver and a simple case? And the tractionBoundaryGradient.C file? Philip |
Hi Philip,
that would be possible but I have to clean things up before that... :D Thanks for checking the code. Alex |
All times are GMT -4. The time now is 11:57. |