KratosMultiphysics / Kratos

Kratos Multiphysics (A.K.A Kratos) is a framework for building parallel multi-disciplinary simulation software. Modularity, extensibility and HPC are the main objectives. Kratos has BSD license and is written in C++ with extensive Python interface.
https://kratosmultiphysics.github.io/Kratos/
Other
991 stars 244 forks source link

Temperature dependent material properties #7641

Open mashsquad777 opened 3 years ago

mashsquad777 commented 3 years ago

Description Query regarding the decision of involving KRATOS in my thermo-mechanical coupling project.

Scope Which areas of Kratos are involved

I am trying to solve the thermo-mechanical problem using heat equation as well as elastoplastic material model. I have just come to know about the KRATOS and want to choose it for part of my project. As in my case of application, mechanical material properties are temperature dependent, so can or does KRATOS SmallStrainIsotropicPlasticity3DVonMisesVonMises take care of the temperature dependence of material properties (i.e. Young's Modulus, Yield strength, hardening coefficient etc.) ? Thank you!

philbucher commented 3 years ago

@KratosMultiphysics/structural-mechanics

@mashsquad777 as far as I know this is not yet possible, we had a large discussion about this some time ago, but we couldn't come up with a good solution unfortunately

In theory it is not difficult to do if you need it only for temperature, I think you could change the code of the Constitutive Law itself.

OR sth that I do myself but is not memory efficient is to create one property per element that already has the interpolated values

also is your coupling one or two way?

tteschemacher commented 3 years ago

@mashsquad777 we recently did something which could be similar to your problem. We simulated the behavior of epoxy in the curing process. To do this we have one convection diffusion problem. Here we compute a lot of post processing variables, as curing rate, specific heat flux, ... One of those variables is the strain, which happens due to the curing, is computed accordingly. We then map this strain to structural mechanics problem and apply it there as a prestrain. Gauss point by gauss point.

Is this something similar to your problem?

mashsquad777 commented 3 years ago

@tteschemacher thanks for the response. In my case I am applying the thermal expansion and contraction on the structure. I want to use rate independent elastoplastic material model e.g. J2 plasticity with Von Mises criteria to take of the thermal loading. As temperature ranges are high so I want to use the material model which also take care of the material properties changing locally with local temperature. The thermo-mechanical coupling is supposed to be one way.