30.- TRATADO COMPLETO SOBRE CÓMO RESOLVER PROBLEMAS LINEALES DE “ZUNCHADO” CON FEMAP y NX NASTRAN

Hola!,

Os he preparado un vídeo muy completo sobre cómo calcular con FEMAP y NX NASTRAN las tensiones de vonMises de Zunchado (también conocido como “Snap-Fit“, “Press-Fit“, “Interference-Fit“, “Overlapping“, etc..) por contacto lineal “superficie-a-superficie” que se producen en el montaje de piezas eje-agujero en donde el diámetro del agujero es ligeramente menor que el diámetro del eje. El zunchado se puede realizar por presión (“Press-Fit“) o utilizar un calentador por inducción para dilatar suficientemente la pieza hembra y proceder el montaje sin esfuerzo. Finalmente el conjunto se deja enfriar hasta la temperatura ambiente, garantizando una unión perfecta. El Zunchado por presión (“Press-Fit“) únicamente se puede resolver mediante análisis no lineal avanzado (SOL601/701), aquí os voy a enseñar únicamente cómo simular las tensiones de contacto por zunchado térmico mediante un sencillo análisis estático lineral (SOL101).

Método-1: INTERFERENCIA “EXPLÍCITA”

En el ejemplo propuesto se trata de estudiar las tensiones de zunchado entre un eje de Acero y un casquillo de Nylon. El radio del eje es 0.05 mm mayor que el radio del agujero, y la geometría de base para el mallado presenta dicha interferencia de forma “explícita”.

En la siguiente imagen definimos los parámetros MINDi y MAXDi correspondientes a la mínima y máxima distancia de búsqueda en la cual el solver NX NASTRAN determina inicialmente si la distancia entre caras de elementos pertenecientes a las regiones de contacto “source” y “target” están dentro del umbral para crear elementos de contacto. Estos valores se utilizan sólo una vez, y al principio del análisis estático lineal (SOL101), y sirven para determinar dónde se deben crear inicialmente elementos de contacto. La mínima distancia debe ser negativa y ligeramente mayor del valor de la interferencia, tomaremos por ejemplo MINDi=-0.5mm. Por supuesto, en este caso utilizamos el parámetro INIPENE=0, así el contacto se evalúa tal como se ha modelizado …

Os recuerdo cómo funcionan los contactos lineales en NX NASTRAN: el solver proyecta normales a partir de las caras de los elementos pertenecientes a la región origen (“source” o “slave“) y mira a ver si alguna de esas normales intersecta con caras libres de elementos de la región destino (“target” o “master“). Si la normal proyectada intersecta con una cara de un elemento, y la distancia entre las dos caras está dentro del rango definido por los valores MINDi y MAXDi, entonces se crea un elemento de contacto. Recordar también que el análisis SOL101 se usa en problemas lineales con pequeños desplazamientos y deformaciones, por tanto el nº de elementos de contacto creados inicialmente permanece constante para el resto del análisis, no se actualiza con cada iteración de contacto.

La siguiente imagen muestra la deformada (exagerada!!) y las tensiones de vonMises (MPa) de zunchado obtenidas mallando con elementos lineales 3-D sólidos hexaédricos CHEXA de 8-nodos (bajo orden): la calidad de los resultados tanto de tensiones como desplazamientos deja mucho que desear!!.

La siguiente imagen muestra la deformada (exagerada!!) y las tensiones de vonMises (MPa) de zunchado obtenidas mallando con elementos parabólicos 3-D sólidos hexaédricos CHEXA de 20-nodos: la calidad de los resultados es excelente!!.

Y finalmente la siguiente imagen muestra la deformada (de nuevo exagerada!!) y las tensiones de vonMises (MPa) de zunchado obtenidas mallando con elementos parabólicos 3-D sólidos tetraédricos de alto orden CTETRA de 10-nodos: sorprende la calidad de los resultados, prácticamente al nivel de los elementos hexaédricos CHEXA de 20-nodos!!.

Método-2: OFFSET

El parámetro OFFSET es una opción de las regiones de contacto que permite tener en cuenta una capa rígida que exista entre caras de elementos en contacto. Puede tener diversas aplicaciones, por ejemplo, pensemos en un modelo que tenga dos superficies metálicas en contacto y una de ellas tenga un revestimiento cerámico. Si la rigidez del material cerámico no es lo bastante significativa como para incluirla en el análisis, entonces no es necesario discretizarla, pero el espesor que añade la capa cerámica a la cara metálica pude ser importante cuando se resuelva el problema de contacto.

La distancia OFFSET también se utiliza para estudiar interferencias en problemas de zunchado cuando se modelizan como caras coincidentes. El valor del OFFSET es precisamente el zunchado o interferencia de contacto entre cuerpos. La ventaja es que tanto el radio del eje como del anillo es exactamente el mismo, pudiendo estudiar rápidamente el efecto de diferentes valores de zunchado si necesidad de modificar la geometría ni mallar de nuevo, toda una ventaja!!. Por supuesto, en este caso utilizamos el parámetro INIPENE=3, así eliminamos interferencias no deseadas …

La siguiente imagen muestra la deformada (exagerada!!) y las tensiones de vonMises (MPa) de zunchado obtenidas mallando con elementos lineales 3-D sólidos hexaédricos CHEXA de 8-nodos considerando un OFFSET = 0.05 mm: la calidad de los resultados es excelente, y el tiempo de cálculo y los recursos de memoria RAM utilizados mínimos!!.

EN RESUMEN …

En la siguiente tabla tenéis compiladas las pruebas realizadas: a destacar la buena precisión de los resultados obtenidos por el elemento tetraédrico parabólico CTETRA de 10-nodos de NX NASTRAN, a costa claro está del tiempo de cálculo y consumo de recursos de memoria RAM, y por supuesto la grata sorpresa del uso del parámetro OFFSET en regiones de contacto que ofrece grandes posibilidades para “jugar” con diferentes valores de zunchados, y le saca partido al uso de elementos lineales CHEXA de 8-nodos sin apenas penalizar la precisión de los resultados, manteniendo una velocidad de cálculo elevada y reducidos recursos de memoria RAM — una joya!!.

Número
de Nodos y
Elementos
URES
máx.
(mm)
vMises
máx.
(MPa)
Tiempo
de
Cálculo
HEX8 15875 nodos
13213 elem.
(46000 gdl)
0.0845 21.2 0:1:15
HEX20 60708 nodos
13213 elem.
(177500 gdl)
0.0527 11.01 0:4:45
TET10 91654 nodos
61328 elem.
(268950 gdl)
0.0527 10.54 0:7:0
HEX8 +
OFFSET
15875 nodos
13213 elem.
(46000 gdl)
0.0529 10.71 0:0:25

Si quieres repetir este tutorial en tu propio ordenador pídenos los modelos con la geometría de entrada y te lo remitimos por e-mail, es un servicio gratuito y exclusivo para nuestros clientes de IBERISA.

Saludos,
Blas.

2 thoughts on “30.- TRATADO COMPLETO SOBRE CÓMO RESOLVER PROBLEMAS LINEALES DE “ZUNCHADO” CON FEMAP y NX NASTRAN

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s