How to recreate NULL material model with dynamic viscosity term in LS-Dyna using a UMAT?
I apologize for missing this topic. Seems like an interesting venture to undertake. I'd be interested if you figured out what the issue was.
I still haven't solved this issue. I've narrowed it down a bit in that it looks like the strain rates that I am calculating (using eps/dt1) do not match the strain rates reported from LS-Dyna when using the built in Null material model.
the results match fine when the dynamic viscosity is 0, but I can't get the results to match at all when I assign a value to the dynamic viscosity.
What are you trying to match to when you assign a dynamic viscosity? I would not think it would match the built-in *Mat_null material because that was the point of re-creating the null material but with dynamic viscosity. Maybe I am just not sure what your goal of this is.
@negativevolume, My goal is actually to modify this null material model eventually, but to start I wanted to make sure that I could replicate the built in one first. The built in null model has a dynamic viscosity term already. So the answers should be very similar if not exactly the same. They are good when you ignore the dynamic viscosity term, but this just means that the hydrostatic portion is good, but the deviatoric portion is not right. The deviatoric portion is controlled by the value of dynamic viscosity (which I make the same for both) and by the strain rate (which is what the difference is from my investigation). I just can't figure out how I should be calculating the strain rate. The LS-Dyna theory shows it as simply the strain increment divided by the timestep, which is what I am doing. Perhaps it is actually calculated in a different manner?
@kennedy Okay, I understand now. This is probably obvious, but what type of element are you using to test the materials? I'm pretty sure shell and beam elements are not affected using *Mat_null (except for contacts) so you should be using solid elements.
Ah, also make sure you know how *Mat_null is calculating strain (true vs engineering). I know LS-Dyna outputs strain as True Strain. This might be where you are getting your difference between the material models.