Elmer Discussion Forum

Hi! Yes, that's right. I didn't think of that. Could this be a better way: In ElasticSolver.F90 544 UseUMAT = ListCheckPresentAnyMaterial(Model, 'UMAT Subroutine') 545 IF (UseUMAT) THEN 546 CALL ListAddNewLogical( SolverParams,'UseUMAT',.TRUE.) 547 END IF then in CalculateLoads in SolverUtils.F90 14...

Hello! To follow up from on an old post. UMAT with Mortar BC didn't seem to give proper displacement contact load values at the mortar boundaries. I think it has to do with ElasticSolver is normally (without UMAT) solving Ku = F using the total displacement while with UMAT it is solving for the incr...

Hi!

I posted a pull request. First time for me using GitHub. I hope I did it correctly and as intended. If not, let me know.

Best regards Fredrik

Hi!

Yes sure! I have a simple test case. I just need to clean up the code a little. I will notify later.

Regarding issue 2 and 3 in the first post I think I have found the cause. I will submit some comments for those in a short while as well.

Best regards Fredrik

I think I have found something with the issue with Normal-Tangential Displacement BC with UMAT, (1) in previous post. Here in ElasticSolve, DValues refer to the NT-rotated coords while Displacement refers to cartesian. 1201 IF (UseUMAT) THEN 1202 ! ---------------------------------------------------...

Hello! About ElasticSolver with UMAT and Mortar BC. I have run into some problems 1. Using UMAT without Mortar BC and Normal-Tangential Displacement = True for Dirichlet BC and together with Neumann BC at the load boundary just output a 0 result. After a little debugging I can see in ElasticSolver t...

Hello!

Yes, the updated solver together with exported variable works. Thanks'

Best regards Fredrik

Hi! I'm trying to use SaveMaterials to see material parameters in paraview. It shows the material parameters at nodes. However, the material parameter is (I guess) defined for the element. For nodes that is shared by 2 elements where the 2 elements have different material parameters, it looks like t...

kevinarden, yes thats the way I modified the ElasicSolver also. I'm not exactly sure how to store the previous displacement. | 756 previ = 0 | 757 IF (UseUMAT) THEN |- 758 IF( TransientSimulation ) THEN ||- 759 previ = 3 || 760 ELSE IF( Scanning ) THEN ||- 761 previ = 1 || 762 END IF | 763 END IF In...

Thanks Peter! With the steady state condition keyword I can run StatCurrentSolver and HeatSolver in transient and ElasticSolver with UMAT in steady state. Exactly was I was asking for. However, the displacement is not evolving from time step to time step. Unfortunately I also get the same error as I...