Error al intentar ejecutar una simulación con el UMAT de Huang usando un solo elemento.
Estoy tratando de usar el UMAT de Huang (http://www.columbia.edu/~jk2079/KysarResearchLaboratory/SingleCrystalUMAT.html) en un elemento cubo en Abaqus que construí siguiendo este tutorial: https://www.youtube.com/watch?v=9GQkNKI09ow.
Cada vez que intento ejecutar la simulación, Abaqus/Standard Analysis se cierra con un error. El archivo de mensaje no existe y en la pestaña de advertencia en el monitor puedo leer lo siguiente: “El parámetro hourglass = enhanced en la opción de control de secciones es relevante para estos elementos: c3d8r, cax4r, cgax4r, cpeg4r, cpe4r, cps4r, m3d4r, s4r, sc8r y sus homólogos híbridos, térmicos y de presión donde corresponda. También es relevante para todos los tipos de elementos triangulares y tetraédricos modificados. Esta advertencia se puede ignorar si la función se aplica solo a estos tipos de elementos”.
También apareció el siguiente mensaje la primera vez que intenté ejecutar las simulaciones: “Se han definido 1 elementos con rigidez cero en sus componentes hour glass. Puede utilizar *hourglass stiffness o cambiar el tipo de elemento. Los elementos han sido identificados en el conjunto de elementos ErrElemZeroHourGlassStiffness”.
Parece que el problema está relacionado con el parámetro hourglass, pero aún no he podido encontrar una solución. ¿Podría alguien ayudarme a entender cuál es el problema? Mi objetivo es ejecutar esta simulación.
“`
Encabezado
* Nombre del trabajo: Job-1 Nombre del modelo: Job-1
* Generado por: Abaqus/CAE 2019
Preimpresión, eco=NO, modelo=NO, historial=NO, contacto=NO
*
PARTES
Parte, nombre=CUBE-1
1, 1., 1., 1.
2, 1., 0., 1.
3, 1., 1., 0.
4, 1., 0., 0.
5, 0., 1., 1.
6, 0., 0., 1.
7, 0., 1., 0.
8, 0., 0., 0.
nobr>Elemento, tipo=C3D8R
1, 5, 6, 8, 7, 1, 2, 4, 3
Nset, nset=SET-1, generar
1, 8, 1
*Elset, elset=SET-1
1,
* Sección: Sección-1-SET-1
,
*Fin de la parte
*
ENSAMBLE
Ensamble, nombre=Ensamble
*
Instancia, nombre=CUBE-1, parte=CUBE-1
*Fin de instancia
*
5, 8, 1
*Elset, elset=SET-1, instancia=CUBE-1
1,
*Nset, nset=SET-2, instancia=CUBE-1
3, 4, 7, 8
*Elset, elset=SET-2, instancia=CUBE-1
1,
*Nset, nset=SET-3, instancia=CUBE-1, generar
2, 8, 2
*Elset, elset=SET-3, instancia=CUBE-1
1,
*Nset, nset=SET-4, instancia=CUBE-1, generar
1, 4, 1
*Elset, elset=SET-4, instancia=CUBE-1
1,
*Fin de ensamble
*
CONTROLES DE ELEMENTOS
1., 1., 1.
*
MATERIALES
APUNTES DE DEPVAR
* =====================
ITRATN , ITRMAX , GAMERR ,
— , — , — ,
ITRATN: parámetro que determina si se utiliza un método iterativo para
resolver incrementos de tensiones y variables de estado en términos de
incrementos de deformación
ITRATN=0. — sin iteración
de lo contrario — con iteración
ITRMAX: número máximo de iteraciones
GAMERR: error absoluto de las deformaciones de corte en sistemas de deslizamiento
*Material, nombre=MATERIAL-1
*Depvar
125,
*Material de usuario, constantes=160, unsymm
168400.,121400., 75400., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
1., 0., 0., 0., 0., 0., 0., 0.
1., 1., 1., 1., 1., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
-1., 0., 1., 0., 0., 1., 0., 0.
0., 1., 0., 0., 1., 0., 0., 0.
10., 0.001, 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
541.5, 109.5, 60.8, 0., 0., 0., 0., 0.
1., 1., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0., 0., 0., 0., 0., 0., 0., 0.
0.5, 1
davy.ai
El mensaje de error sugiere que el problema podría estar relacionado con el parámetro de reloj de arena. En particular, la advertencia indica que el parámetro es relevante para ciertos tipos de elementos y sus contrapartes modificadas. Sin embargo, parece que el elemento de cubo que construyó sigue un tipo diferente, C3D8R, que podría no estar cubierto por la advertencia.
Para resolver el problema, podría intentar establecer el parámetro de reloj de arena en cero, o deshabilitarlo por completo eliminando la opción
hourglass=enhancedde la entradaSection Controls. Esto podría ayudar a evitar el mensaje de error sobre la rigidez del reloj de arena cero y permitir que la simulación se ejecute correctamente.Otra opción sería modificar el tipo de elemento de su cubo a uno de los tipos mencionados en la advertencia, como
c3d8r,s4rosc8r, entre otros, que podrían funcionar mejor con el UMAT de Huang que está utilizando. Sin embargo, esto requeriría trabajo adicional y puede que no sea práctico dependiendo de sus necesidades específicas.En general, vale la pena experimentar con diferentes configuraciones y tipos de elementos para encontrar la mejor solución para su simulación. Además, verificar los datos de entrada y asegurarse de que sean consistentes y precisos también puede ayudar a diagnosticar y resolver cualquier problema que pueda surgir.