[simrego]

Hi!


I wrote a small total lagrangian solid displacement solver, and I'm trying to create material models for it. The material model is simply returns the Lame's coefficients, BUT. If I loop through the field in my update() function in my material model, the solver gives some weird results.


For example:
Simple bilinear model, I have let's call Ee_, and Ep_, and based on the strain and stress i set the E_ from Ee_ and Ep_. If i loop through the list and set like "E_[i] = Ee_[i];", the results are reallly bad. If i set "E_ = Ee_;" without loop, it's working perfectly.
I wrote out the E_ field for both case and the are exactly the same...


Anyone any ideas? I'm pretty lost now.


Thanks!

[simrego]

E_.correctBoundaryConditions(); did it. :S

[EmadTandis]

Originally Posted by simrego

Hi!


I wrote a small total lagrangian solid displacement solver, and I'm trying to create material models for it. The material model is simply returns the Lame's coefficients, BUT. If I loop through the field in my update() function in my material model, the solver gives some weird results.


For example:
Simple bilinear model, I have let's call Ee_, and Ep_, and based on the strain and stress i set the E_ from Ee_ and Ep_. If i loop through the list and set like "E_[i] = Ee_[i];", the results are reallly bad. If i set "E_ = Ee_;" without loop, it's working perfectly.
I wrote out the E_ field for both case and the are exactly the same...


Anyone any ideas? I'm pretty lost now.


Thanks!

Hi,

I guess the difference is because for the loop case, you don't set boundaryFields but for E_=Ee_, not only cell field but also boundaries are set as well.

[simrego]

Yep you are right. That was the problem.
Thx the feedback.