TFG 2210 Hernandez
TFG 2210 Hernandez
Autor:
Carlos Hernández Castelao
Tutor:
Jesús Vázquez Valeo
Profesor titular
iii
Proyecto Fin de Grado: Software Excel: tensiones y deformaciones en fatiga por fretting con contacto
cilíndrico
El tribunal nombrado para juzgar el Proyecto arriba indicado, compuesto por los siguientes miembros:
Presidente:
Vocales:
Secretario:
Sevilla, 2019
v
El Secretario del Tribunal
A mi familia
A mis maestros
vii
Resumen
La fatiga por fretting es un tipo de fatiga que puede provocar el fallo en multitud de componentes industriales.
Este fenómeno de fallo se produce en partes de estos componentes en las que existe una unión mecánica,
unión en la que se encuentran en contacto dos o más piezas. Estos contactos poseen una amplia variedad de
formas geométricas, sin embargo, en muchas situaciones, por su similitud geométrica y de acción de cargas
externas, estos pares de contacto se pueden simplificar a un contaco cilíndrico, es decir, por el contacto entre
un cilindro y un plano.
El objetivo de este proyecto ha sido la creación de un software para calcular las tensiones y deformaciones y la
representación gráfica de las mismas producidas en un contacto cilíndrico. A parte de estos objetivos, este
software se ha realizado de forma que sea lo más sencillo posible para que, una vez finalizado, este se pueda
usar fácil e intuitivamente por cualquier persona. Para conseguirlo, se ha usado el software Excel por las
características que este posee, especialmente por el lenguaje de programación que tiene integrado: VBA. Este
lenguaje se ha utilizado para crear dos macros que calculan las tensiones y deformaciones producidas en un
contacto cilíndrico. La diferencia entre ambas es que una realiza el cálculo tensiones y deformaciones a lo
largo de un segmento y otra muestra los resultados de una matriz de puntos. Los resultados numéricos
obtenidos ejecutando cada una de ellas se muestran en hojas de Excel diferentes.
En este documento se realiza una introducción a los contactos mecánicos y sus posibles mecanismos de fallos.
Posteriormente, se explican las principales características de la fatiga por fretting y las particularidades del
contacto cilíndrico, centrándose en las expresiones analíticas que gobiernan este problema y que se usan para
calcular tanto las tensiones como deformaciones en el software desarrollado. El siguiente apartado describe las
dos hojas del archivo Excel construido, así como los pasos que hay que seguir para conseguir obtener los
resultados numéricos. Por último, se muestran las conclusiones de la realización de este trabajo y se anexan los
códigos de Visual Basic desarrollados para la creación de las macros.
ix
Abstract
Fretting fatigue is a type of fatigue that can provoke failures in multiple industrial components. This failure
phenomenon takes place in parts of those components that have a mechanic union, in which two or more
pieces come in contact. This contacts possess a large variety of geometrical shapes, however in several
situations, due to their geometrical resemblance and the action of external forces these contacts can be
simplified to a cylindrical contact, meaning the contact between a cylinder and a plane.
The objective of this project has been the creation of a software to calculate tensions and deformations and the
graphical representation of them produced in a cylindrical contact. Leaving these objectives behind, this
software has been made as simple as possible so that, once completed, it can be used easily and intuitively by
anyone. Excel Software has been used to accomplish this, for the characteristics that it has, and specially for
the programing language integrated: VBA. This language has been used to create two macros that calculate the
stresses and strains produced in cylindrical contacts. The difference between both is that one calculates the
stresses and strains along a segment and the other shows the results of a dot matrix. The results obtained
executing each are shown in different excel pages.
This document will make an introduction to mechanical contacts and the possible failure mechanism.
Afterwards the main characteristics of fretting fatigue and the particularities of cylindrical contact will be
explained, focusing in analytic expressions that rule this problem and are used to calculate stresses as well as
strains in the developed software. The next section describes the two sheets of the Excel generated, as well as
the steps that are needed follow to obtain the numerical results. Finally, the conclusions obtained during the
fulfillment of this project are explained and the codes of Visual Basic developed for the macro’s creation
annexed.
xi
Índice
Resumen ix
Abstract xi
Índice xiii
Índice de Figuras xv
1 Introducción a Contacos Mecánicos 1
1.1 Introducción al Contacto Mecánico 1
1.2 Clasificación Contactos Mecánicos 1
1.2.1 Contacto sin fricción 1
1.2.2 Contacto con fricción 1
1.2.3 Contacto conforme 1
1.2.4 Contacto no conforme 2
1.2.5 Contacto no hertziano 2
1.2.6 Otros contactos 2
1.3 Mecanismos de fallo asociados al contacto mecánico 2
1.3.1 Desgaste Adhesivo 2
1.3.2 Desgaste Abrasivo 3
1.3.3 Corrosión 3
1.3.4 Fatiga superficial 3
2 Fatiga por Fretting 5
2.1 Descripción fatiga por fretting 5
2.2 Reseña histórica 7
2.3 Ensayos de fretting 8
2.3.1 Geometrías reales 8
2.3.2 Geometrías sencillas 8
3 Contacto Cílindrico 11
3.1 Descripción contacto cilindro-plano 11
3.2 Expresiones analíticas del contacto 15
3.2.1 Tensiones 15
3.2.2 Deformaciones 19
4 Descripción del Software 21
4.1 Hoja Principal 21
4.1.1 Tabla de variables 21
4.1.2 Botón Cilindro_Plano 23
4.1.3 Columnas de resultados 25
4.1.4 Lista de nombres 28
4.1.5 Construcción de una gráfica 29
4.1.6 Ejemplos 31
4.2 Hoja Contour 42
4.2.1 Plantilla 42
4.2.2 Tabla de variables 42
xiii
4.2.3 Botón ContourPlot 44
4.2.4 Matriz de datos 46
4.2.5 Ejemplo 47
5 Conclusiones 53
6 Anexos: Macros de Visual Basic 55
6.1 Macro Cilindro_Plano 55
6.2 Macro ContourPlot 63
7 Bibliografía 75
ÍNDICE DE FIGURAS
xvii
1 INTRODUCCIÓN AL CONTACTO MECÁNICO
1
2 Introducción al Contacto Mecánico
1.3.3 Corrosión
El desgaste corrosivo se produce como resultado de la interacción mecánica de las superficies en contacto por
la influencia de las condiciones del medio en que se produce. Este ambiente corrosivo provoca unas reacciones
químicas o electroquímicas cuyos productos se quedan en una o ambas de las superficies. Estos productos
quedan generalmente poco adheridos a la superficie, por lo que se eliminan fácilmente de esta con el simple
movimiento entre sólidos. Sin embargo, la eliminación de las partículas dejan a la superficie sin protección y el
efecto corrosivo se repite de nuevo. De esta manera, la interacción repetitiva entre la corrosión y un
movimiento de deslizamiento puede dar como resultado altas tasas de desgaste.
Este tipo de desgaste se puede prevenir reduciendo el efecto del ambiente corrosivo o minimizando el
movimiento relativo entre ambos cuerpos. También se puede hacer uso de algún lubricante que sea
conveniente para las condiciones del medio en el que se aplica el contacto.
3
4 Introducción al Contacto Mecánico
2 FATIGA POR FRETTING
Las condiciones para que aparezca el fenómeno fatiga por fretting es necesario que actúe, además de la fuerza
normal N y la fuerza tangencial Q, una tensión externa σ en uno de los cuerpos en contacto. Por el efecto de las
fuerzas normal y tangencial en la zona de contacto, aparecen en el interior de los sólidos una distribución de
tensiones σlocal, mientras que la tensión externa σ crea una distribución de tensiones σglobal tanto en la zona
cercana al contacto como en la zona alejada de él. La diferencia entre fatiga por fretting y fretting se muestra
en la Figura 2-2. Tanto la fuerza Q como la tensión σ pueden variar con el tiempo, ya sea por vibraciones o
porque están sometidos a esfuerzos alternos. Este hecho hace que, tras la iniciación de grietas en la superficie
por la acción de las fuerzas de contacto, dichas grietas tengan la posibilidad de propagarse hacia el interior del
sólido llegando a producir la fractura del componente. Dependiendo del estado de las cargas, el fretting puede
llegar a tener influencia en la iniciación, propagación, o incluso en ambas fases.
5
6 Fatiga por Fretting
Bajo la condición de deslizamiento parcial, al existir una tensión externa, la zona de adhesión se desplaza una
distancia e, tal y como se aprecia en la Figura 2-3. Esto ocurre porque la acción de la tensión externa σ provoca
una diferencia entre las deformaciones de cada cuerpo en la superficie de contacto.
La fatiga por fretting puede aparecer, a priori, en cualquier componente donde haya un contacto mecánico,
siendo fácilmente identificable debido a que tiene características apreciables. Una de estas características es la
formación de partículas de óxido que surgen del desgaste de la superficie y que son de diferentes colores
dependiendo del materiales de los sólidos en contacto, como por ejemplo negro en los aluminios o rojo en los
aceros. Otra característica apreciable es la formación de huellas en las superficies en contacto.
Este fenómeno es muy importante en la ingeniería industrial puesto que existen numerosos componentes
expuestos a él. Algunos ejemplos de componentes industriales en los que puede aparecer el fretting o fatiga
por fretting son uniones atornilladas, roblonadas o remachadas, álabes de turbinas, rodamientos, contactos por
chavetas, uniones eje cubo de una transmisión, diferentes engranajes, cables metálicos trenzados, rotores
generadores o incluso implantes artificiales. En la Figura 2-4 se muestran algunos ejemplos anteriores,
señalándose aquellas zonas propensas a sufrir este tipo de fatiga.
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 7
conocida y puede ser tratado como un problema de dos dimensiones. Una de las principales dificultades que
presenta este contacto es que, a diferencia del contacto esférico, es necesario alinear el eje del indentador y la
dirección de la carga tangencial, siendo necesario que se mantengan perpencidulares entre sí. Otro
inconveniente que tiene es que no se conoce donde se iniciará la grieta, por lo que es necesario realizar un
análisis microscópico de toda la zona de contacto.
En lo que se refiere al contacto con punzón plano, es un contacto con varias dificultades. Entre ellas destaca
que se conoce la distribución de tensiones únicamente en la superficie de contacto. Para evaluar los resultados
habría que hacer uso algún método numérico. Otra desventaja es que este contacto necesita un alineamiento
respecto a la carga de los dos ejes. Otra consecuencia del uso de este contacto a nivel experimental es que, por
consecuencia de la fricción, se crea un par de giro en el punzón. Este tipo de contactos es útil cuando se realiza
el ensayo mediante puentes de fretting fatiga.
El contacto de un punzón recto con las esquinas redondeadas tiene inconvenientes similares al contacto con un
punzón plano, pero posee una geometría más próxima a la realidad, puesto que no existen aristas vivas en la
práctica, sino que los sólidos presentan aristas redondeadas.
10 Fatiga por Fretting
3 CONTACTO CÍLINDRICO
Esta clase de contacto está caracterizado por tener un ancho de superficie de contacto 2a, que es mucho menor
que el radio de curvatura del cilindro R. El semiancho a se determina con la siguiente expresión.
Como se ha mencionado en el segundo apartado de este documento, al actuar una fuerza tangencial Q se crea
una zona adherida en el centro de la superficie de contacto y otras dos de deslizamiento en los extremos de la
misma. A partir de esta fuerza y su relación con la fuerza normal N mediante el coeficiente de rozamiento μ se
obtiene el valor de c que nos proporciona la información del semiancho de la zona adherida.
El desplazamiento de la zona adherida e, que aparece por la presencia de una tensión externa σ, depende del
radio de curvatura, de la fricción existente entre el cilindro y el plano y de las características elásticas de los
cuerpos.
11
12
Contacto Cílindrico
En este contacto y bajo condiciones de cargas que provocan un deslizamiento parcial, el deslizamiento relativo
de ambas zonas puede producirse en el mismo sentido o con ambas zonas deslizando en sentidos opuestos. El
segundo de estos se denomina deslizamiento reverso. Este hecho depende directamente de la relación existente
entre las cargas normal N, tangencial Q y tensión σ aplicadas. La condición necesaria para que no ocurra el
deslizamiento reverso es que se cumpla la relación , es decir:
Si esta condición es cierta, la solución del contacto cilíndrico es analítica y conocida. Si no fuera así, el
problema no tendría solución analítica y habría que resolver el problema mediante ecuaciones integrales que
relacionan las deformaciones en la superficie con la distribución de tensiones tangenciales.
Una vez descritas las expresiones que determinan las longitudes características de la zona de contacto, se van a
exponer las referentes a las presiones que provocan tanto la fuerza normal N como la fuerza tangencial Q en
las direcciones en las que se aplican.
La distribución de tensiones que se produce por la acción de la fuerza normal N en dirección perpendicular a
la superficie de contacto hacia el interior del plano es:
Donde p0 es la máxima presión normal producida en este contacto y depende tanto de N como de la longitud
de la zona de contacto a. Esta presión máxima ocurre en el centro de la zona de contacto tal y como se observa
en la Figura 3-2 y Figura 3-3. Se determina mediante la siguiente expresión.
12
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 13
-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
0.0
-0.2
σNyy(x,0)/p0
-0.4
-0.6
-0.8
-1.0
x/a
Figura 3-2. Distribución de tensiones normales producidas por la acción de la carga N en la dirección de esta adimensionalizadas por
p0 en estado de compresión
-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
0.0
-0.2
σNyy(x,0)/p0
-0.4
-0.6
-0.8
-1.0
x/a
Figura 3-3. Distribución de tensiones normales producidas por la acción de la carga N en la dirección de esta adimensionalizadas por
p0 en estado de tracción
La distribución hace referencia a un estado de deslizamiento global, es decir, si Q = μN. Esta tiene
influencia en toda la zona de contacto de nuestro problema. Por su parte la distribución correctiva
afecta solo a la zona adherida, donde se tiene que cumplir la condición Q < μN. Las expresiones de estas
tensiones son las que se muestran a continuación.
14
Contacto Cílindrico
Se han añadido las Figuras 3-4 y 3-5 para representar esta distribución de tensión tangencial.
0.40
0.35
0.30
σQxy(x,0)/p0
0.25
0.20
0.15
0.10
0.05
0.00
-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
x/a
Figura 3-4. Distribución de tensiones tangenciales producidas por la acción de la carga Q en la dirección de esta adimensionalizadas
por p0 en estado de compresión
14
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 15
0.00
-0.05
-0.10
σQxy(x,0)/p0
-0.15
-0.20
-0.25
-0.30
-0.35
-0.40
-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
x/a
Figura 3-5. Distribución de tensiones tangenciales producidas por la acción de la carga Q y tensión σ adimensionalizadas por p0 en
estado de tracción
3.2.1 Tensiones
En este subapartado se van a exponer todas las expresiones de los campos de tensiones del contacto. La
tensión total será la suma de las tensiones provocadas por cada una de las acciones.
Las tensiones provocadas por la fuerza normal N y la fuerza tangencial Q tienen expresiones analíticas
diferentes en la superficie y en el interior del plano, por lo que se van a describir por separado.
3.2.1.1 Superficie
En cuanto a las tensiones ejercidas por la acción de la fuerza tangencial Q hay que recordar que es una
superposición de distribuciones, por lo que la tensión tangencial total será la suma de la tensión en estado de
deslizamiento global y la tensión correctiva. Las expresiones de ambas son las que se muestran a continuación.
16
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 17
3.2.1.2 Interior
Antes de describir las expresiones de las tensiones en el interior es necesario definir los parámetros s1 y s2, que
son las raíces positivas de la siguiente ecuación cuadrática.
Para obtener el parámetro s1, el valor de la amplitud es equivalente al semiancho de la zona de contacto a
mientras que para s2, la amplitud es igual al límite de la zona de adhesión c. En el cálculo de s2 también hay
que tener en cuenta que el valor de x estará desplazado por la excentricidad e.
Las distribuciones de tensiones en el interior del plano provocadas por la fuerza normal N son las siguientes.
18
Contacto Cílindrico
18
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 19
3.2.2 Deformaciones
Al igual que las tensiones, las distribuciónes de deformaciones totales serán la suma de todas las
deformaciones provocadas por cada acción externa. Las componentes del tensor de deformaciones totales se
calculan a partir de las siguientes expresiones.
Al tratarse del caso de deformación plana sólo tendrán un valor distinto de cero las componentes contenidas en
el plano XY. La deformación que causa cada una de las fuerzas en cada una de las componentes del tensor de
deformaciones se determina de la siguiente manera.
Donde el superíndice k se sustituirá por cada una de las fuerzas para así obtener las distribuciones de
20
Contacto Cílindrico
deformaciones.
20
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 21
Las variables a introducir por el usuario son, de superior a inferior posición en la tabla, las siguientes:
E: módulo de elasticidad longitudinal o de Young (Pa)
ν: coeficiente de Poisson
R: radio de curvatura del cilindro (m)
P: fuerza normal aplicada al cilindro (N)
Q: fuerza tangencial aplicada al cilindro (N)
σ: tensión externa aplicada al plano (Pa)
μ: coeficiente de fricción
xi: coordenada inicial respecto al eje X (m)
yi: coordenada inicial respecto al eje Y (m)
xf: coordenada final respecto al eje X (m)
yf: coordenada inicial respecto al eje Y (m)
Divisiones segmento: número de divisiones a analizar en la recta formada por la unión de los
puntos inicial y final
Por otro lado, las variables de salida que se van a generar tras rellenar las celdas correspondientes a las
variables de entrada son:
a: longitud característica del semiancho de contacto (m)
22
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 23
Especificación:
Esfuerzo: se refiere al esfuerzo axil al que se someterá al plano, es decir, tracción o compresión.
Tras presionar sobre el botón, Excel mostrará una pantalla emergente en la cual nos pedirá que se introduzca el
esfuerzo al que se quiere someter al plano que está en contacto con el cilindro. El valor uno corresponde con
un esfuerzo de compresión y dos con un esfuerzo de tracción, tal y como se aprecia en la Figura 4-3.
Si se introduce uno de esos valores la macro programada se ejecutará correctamente y los resultados se
mostrarán en sus respectivas columnas de la hoja Principal, rellenándose además la celda correspondiente a la
especificación Esfuerzo de la tabla de variables, como se muestra en la Figura 4-4. Se ha querido añadir esta
información a la tabla para que el usuario sea siempre consciente de las condiciones del problema que se está
resolviendo.
24
Descripción del Software
Si, por el contrario, se presiona el botón de cancelar, se cierra la pantalla Introduzca esfuerzo, no se introduce
ningún valor y se pulsa aceptar, o se introduce un valor distinto de uno o dos, Excel mostrará un mensaje de
error pidiendo que elija una de las opciones, tal como se muestra en la Figura 4-5.
También puede haber otros errores que impidan que se ejecute la macro. Estos se deben a una introducción
errónea de valores en las variables de entrada de la tabla. Si se introduce un valor mayor que cero en la
coordenada inicial o final del eje Y se mostrará el mensaje de la Figura 4-6. En cambio, si el fallo ha sido en
que el valor de la fuerza tangencial excede el límite marcado por el coeficiente de rozamiento y la fuerza
normal se apreciará un mensaje como el de la Figura 4-7. Por último, si se introducen unos valores en las
acciones externas que impliquen unas condiciones de deslizamiento reverso se observará el mensaje de la
Figura 4-8.
Figura 4-6. Ventana emergente de error si se introduce un valor de coordenada yi o yf mayor que cero
24
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 25
El grupo Número de puntos consta de una sola columna, la cual se encuentra en la columna E de la hoja. Esta
columna depende directamente de la variable de entrada Divisiones segmento y determina el número de filas
de datos que se obtendrán en la hoja. Esta columna tiene un color de fondo de celda amarillo tal y como se
aprecia en la Figura 4-9.
26
Descripción del Software
4.1.3.2 Coordenadas
Este grupo proporciona información acerca de las coordenadas de cada división a analizar. Este abarca desde
la columna F hasta la columna I y para diferenciarlas del resto se les ha asignado un color de fondo naranja
como se puede observar en la Figura 4-10.
Todos los valores que muestran las celdas de este grupo están adimensionalizados por la longitud de
semiancho de contacto a para una mejor visión de los datos obtenidos. La primera columna de este grupo se
llama x/a y depende de las variables de entradas xi y xf . Esta columna muestra los valores de la coordenada X
de cada punto discretizado.
La siguiente columna, denominada y/a, desempeña la misma función que la columna x/a pero proporcionando
la información de las coordenadas respecto al eje Y, por lo que depende de las variables yi e yf.
El título de la tercera columna es s1/a y es necesaria para los cálculos de las tensiones. Esta depende tanto de
los valores de las coordenadas X e Y de cada división del segmento como de la longitud del semiancho de
contacto a.
La última columna de este grupo es s2/a y, al igual que s1/a depende de las coordenadas X e Y de cada punto
intermedio y de la longitud de semiancho de la zona adherida c y de la excentricidad de esta e.
26
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 27
4.1.3.3 Tensiones
El grupo Tensiones se reconocen por el color de fondo verde que posee en la hoja Excel, tal y como se aprecia
en las Figura 4-11 y Figura 4-12. Los valores de estas columnas están adimensionalizados respecto a la presión
nominal p0 o por μp0 en función de la fuerza que provoque las tensiones. Por orden de izquierda a derecha, se
muestran las distribuciones referidas a la fuerza normal P, a la fuerza tangencial Q , tanto deslizamiento global
como correctiva, y las tensiones suma de todas las acciones externas incluyendo, a parte de las tres
distribuciones mencionadas, la tensión externa σ. Por último, cada grupo de distribución de tensiones posee
cuatro columnas que representan las componentes del tensor de tensiones σxx, σyy, σyx, σzz.
Figura 4-12. Columnas de tensiones provocadas por Q correctiva y por la suma de todas las acciones
4.1.3.4 Deformaciones
El grupo de los valores de las deformaciones se identifica por el color azul de fondo de celda que poseen, tal y
como se aprecia en la Figura 4-13. Al igual que las tensiones, los valores que aparecen en la hoja están
adimensionalizados por la presión nominal p0 o por μp0 dependiendo de la fuerza que provoque la distribución.
Este grupo se subdivide de la misma manera que el anterior en lo que a fuerzas externas se refiere pero se
diferencia en que cada una de las acciones solo tiene componentes εxx, εyy y εxy al tratarse el problema en
estudio con hipótesis de deformación plana.
28
Descripción del Software
28
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 29
Tras esto, hay que pulsar el botón derecho del ratón sobre el cuadro blanco donde se va a representar la gráfica
y pinchar en la opción Seleccionar datos, como se muestra en la Figura 4-16.
Llegados a este punto, interesa tener los datos de las tensiones o deformaciones en el eje Y de la gráfica,
mientras que en el eje X los valores de la columna s1. Para conseguir esto e introducir los valores deseados en
30
Descripción del Software
el eje Y hay que pinchar en el botón Agregar que aparece en la pantalla Seleccionar origen de datos, tal y
como se aprecia en la Figura 4-17.
30
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 31
4.1.6 Ejemplos
Mediante la programación de las distintas ecuaciones descritas anteriormente de las tensiones y deformaciones
en Visual Basic se consiguen desarrollar las diferentes gráficas buscadas con el desarrollo del software Excel
construido. A continuación se van a mostrar una serie de gráficas que caracterizan el problema de fatiga por
fretting en el contacto entre un cilindro y un plano. Para diferenciar el caso de compresión del caso de tracción
con mayor facilidad, las gráficas que se muestran a continuación se han representado en color rojo si se trata de
compresión y en color azul si se trata de tracción.
4.1.6.1 Superficie
En primer lugar, se van a representar las tensiones y deformaciones sobre la superficie del plano, es decir, se
van a elegir solamente puntos con coordenada y igual a cero. Para esta ejemplificación se van a considerar los
siguientes datos de entrada:
E = 2.10E11 Pa
ν = 0.3
R = 0.1 m
P = 1.00E6 N
Q = 1.00E5 N
σ = 1.10E8 Pa
μ = 0.7
xi = -1.5757E-3 m
yi = 0 m
xf = 1.5757E-3 m
yf = 0 m
Divisiones segmento = 250
El resultado de las variables de salida para las anteriores variables de entrada es:
a = 1.0505E-3 m
c = 9.7254E-4 m
e = 6.8095E-5 m
p0 = 6.0604E8 Pa
Se han elegido unos valores de fuerzas que no dieran lugar al caso de deslizamiento inverso, es decir, que la
longitud a sea mayor que las longitudes c y e sumadas. También se ha seleccionado un rango en el eje X que
abarcara desde -1.5a hasta 1.5a para conocer que ocurre en ciertas tensiones y deformaciones fuera de la zona
de contacto. Las gráficas de tensiones y deformaciones para el estado de compresión abarcan desde la Figura
4-20 hasta la Figura 4-26 mientras que las gráficas referentes al estado de tracción de la Figura 4-27 hasta la
Figura 4-33.
32
Descripción del Software
0.0
-0.2
-0.4
σxx/p0
-0.6
-0.8
-1.0
-1.2
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-20. Tensión σxx en la superficie para el estado de compresión adimensionalizada por p0
0.0
-0.2
-0.4
σyy/p0
-0.6
-0.8
-1.0
-1.2
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-21. Tensión σyy en la superficie para el estado de compresión adimensionalizada por p0
0.40
0.35
0.30
0.25
σxy/p0
0.20
0.15
0.10
0.05
0.00
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-22. Tensión σxy en la superficie para el estado de compresión adimensionalizada por p0
32
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 33
0.00
-0.10
-0.20
-0.30
σzz/p0
-0.40
-0.50
-0.60
-0.70
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-23. Tensión σzz en la superficie para el estado de compresión adimensionalizada por p0
0.0E+00
-5.0E-13
-1.0E-12
-1.5E-12
-2.0E-12
εxx/p0
-2.5E-12
-3.0E-12
-3.5E-12
-4.0E-12
-4.5E-12
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-24. Deformación εxx en la superficie para el estado de compresión adimensionalizada por p0
1.5E-12
1.0E-12
5.0E-13
0.0E+00
εyy/p0
-5.0E-13
-1.0E-12
-1.5E-12
-2.0E-12
-2.5E-12
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-25. Deformación εyy en la superficie para el estado de compresión adimensionalizada por p0
34
Descripción del Software
2.5E-12
2.0E-12
1.5E-12
εxy/p0
1.0E-12
5.0E-13
0.0E+00
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-26. Deformación εxy en la superficie para el estado de compresión adimensionalizada por p0
1.0
0.8
0.6
0.4
0.2
σxx/p0
0.0
-0.2
-0.4
-0.6
-0.8
-1.0
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-27. Tensión σxx en la superficie para el estado de tracción adimensionalizada por p0
0.0
-0.2
-0.4
σyy/p0
-0.6
-0.8
-1.0
-1.2
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-28. Tensión σyy en la superficie para el estado de tracción adimensionalizada por p0
34
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 35
0.00
-0.05
-0.10
-0.15
σxy/p0 -0.20
-0.25
-0.30
-0.35
-0.40
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-29. Tensión σxy en la superficie para el estado de tracción adimensionalizada por p0
0.30
0.20
0.10
0.00
-0.10
σzz/p0
-0.20
-0.30
-0.40
-0.50
-0.60
-0.70
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-30. Tensión σzz en la superficie para el estado de tracción adimensionalizada por p0
4.0E-12
3.5E-12
3.0E-12
2.5E-12
2.0E-12
1.5E-12
εxx/p0
1.0E-12
5.0E-13
0.0E+00
-5.0E-13
-1.0E-12
-1.5E-12
-2.0E-12
-2.5E-12
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-31. Deformación εxx en la superficie para el estado de tracción adimensionalizada por p0
36
Descripción del Software
0.0E+00
-5.0E-13
-1.0E-12
εyy/p0
-1.5E-12
-2.0E-12
-2.5E-12
-3.0E-12
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-32. Deformación εyy en la superficie para el estado de tracción adimensionalizada por p0
0.0E+00
-5.0E-13
-1.0E-12
εxy/p0
-1.5E-12
-2.0E-12
-2.5E-12
1.12 0.66 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.66 1.12
s1/a
Figura 4-33. Deformación εxy en la superficie para el estado de tracción adimensionalizada por p0
Ahora va a desarrollarse la opción genérica que se puede conseguir con el programa diseñado. Para ello se va a
seleccionar un segmento diagonal que vaya desde un punto de la superficie hasta un punto que se encuentre en
el interior con coordenada y menor que cero. Primero se muestran todas las tensiones y deformaciones para el
estado de compresión del plano, desde la Figura 4-34 a la Figura 4-40, y siguiendo a estas, se muestran las
mismas gráficas pero para el estado de tracción, abarcando desde la Figura 4-41 hasta la Figura 4-47. Los
valores de entrada seleccionados para este caso son los siguientes:
E = 2.10E11 Pa
ν = 0.3
R = 0.1 m
P = 1.00E6 N
Q = 1.00E5 N
σ = 1.10E8 Pa
μ = 0.7
36
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 37
xi = -1.5757E-3 m
yi = 0 m
xf = 1.5757E-3 m
yf = -5.2523E-4 m
Divisiones segmento = 250
El resultado de las variables de salida para las anteriores variables de entrada es:
a = 1.0505E-3 m
c = 9.7254E-4 m
e = 6.8095E-5 m
p0 = 6.0604E8 Pa
0.0
-0.1
-0.2
-0.3
σxx/p0
-0.4
-0.5
-0.6
-0.7
-0.8
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-34. Tensión σxx en el interior para el estado de compresión adimensionalizada por p0
0.0
-0.2
-0.4
σyy/p0
-0.6
-0.8
-1.0
-1.2
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-35. Tensión σyy en el interior para el estado de compresión adimensionalizada por p0
38
Descripción del Software
0.30
0.25
0.20
0.15
0.10
σxy/p0
0.05
0.00
-0.05
-0.10
-0.15
-0.20
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-36. Tensión σxy en el interior para el estado de compresión adimensionalizada por p0
0.00
-0.10
-0.20
σzz/p0
-0.30
-0.40
-0.50
-0.60
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-37. Tensión σzz en el interior para el estado de compresión adimensionalizada por p0
0.0E+00
-5.0E-13
-1.0E-12
εxx/p0
-1.5E-12
-2.0E-12
-2.5E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-38. Deformación εxx en el interior para el estado de compresión adimensionalizada por p0
38
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 39
1.0E-12
5.0E-13
0.0E+00
-5.0E-13
-1.0E-12
εyy/p0 -1.5E-12
-2.0E-12
-2.5E-12
-3.0E-12
-3.5E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-39. Deformación εyy en el interior para el estado de compresión adimensionalizada por p0
2.0E-12
1.5E-12
1.0E-12
5.0E-13
εxy/p0
0.0E+00
-5.0E-13
-1.0E-12
-1.5E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-40. Deformación εxy en el interior para el estado de compresión adimensionalizada por p0
0.2
0.1
0.0
-0.1
σxx/p0
-0.2
-0.3
-0.4
-0.5
-0.6
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-41. Tensión σxx en el interior para el estado de tracción adimensionalizada por p0
40
Descripción del Software
0.2
0.0
-0.2
σyy/p0 -0.4
-0.6
-0.8
-1.0
-1.2
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-42. Tensión σyy en el interior para el estado de tracción adimensionalizada por p0
0.25
0.20
0.15
0.10
σxy/p0
0.05
0.00
-0.05
-0.10
-0.15
-0.20
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-43. Tensión σxy en el interior para el estado de tracción adimensionalizada por p0
0.10
0.00
-0.10
σzz/p0
-0.20
-0.30
-0.40
-0.50
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-44. Tensión σzz en el interior para el estado de tracción adimensionalizada por p0
40
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 41
1.0E-12
5.0E-13
εxx/p0 0.0E+00
-5.0E-13
-1.0E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-45. Deformación εxx en el interior para el estado de tracción adimensionalizada por p0
5.0E-13
0.0E+00
-5.0E-13
-1.0E-12
εyy/p0
-1.5E-12
-2.0E-12
-2.5E-12
-3.0E-12
-3.5E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-46. Deformación εyy en el interior para el estado de tracción adimensionalizada por p0
1.5E-12
1.0E-12
5.0E-13
εxy/p0
0.0E+00
-5.0E-13
-1.0E-12
-1.5E-12
1.12 0.67 0.21 0.19 0.21 0.25 0.31 0.42 0.62 0.93 1.29
s1/a
Figura 4-47. Deformación εxy en el interior para el estado de tracción adimensionalizada por p0
42
Descripción del Software
4.2.1 Plantilla
Antes de la descripción del proceso que hay que realizar dentro de la hoja Contour para obtener el gráfico de
superficie hay que ubicar el archivo Plantilla.crtx. Este debe ubicarse en la carpeta de plantillas de Excel en el
equipo donde se desarrolle el trabajo. Esta plantilla gráfica de Excel es un gráfico de superficie que consta de
13 bandas que representan las distintas áreas de contorno. Estas bandas poseen una gama de colores idéntica a
la que se usa en el Método de los Elementos Finitos, correspondiéndose el valor mínimo con el azul más
oscuro y el valor máximo con el rojo más oscuro.
42
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 43
En esta hoja las variables de entrada son, de superior a inferior posición en la tabla, las siguientes:
E: módulo de elasticidad longitudinal o de Young (Pa)
ν: coeficiente de Poisson
R: radio de curvatura del cilindro (m)
P: fuerza normal aplicada al cilindro (N)
Q: fuerza tangencial aplicada al cilindro (N)
σ: tensión externa aplicada al plano (Pa)
μ: coeficiente de fricción
xi: coordenada inicial respecto al eje X (m)
yi: coordenada inicial respecto al eje Y (m)
xf: coordenada final respecto al eje X (m)
yf: coordenada inicial respecto al eje Y (m)
Divisiones eje X: número de divisiones entre las coordenadas inicial y final en el eje X
Divisiones eje Y: número de divisiones entre las coordenadas inicial y final en el eje Y
En este caso las variables de salida son las mismas que en la hoja Principal.
a: longitud característica del semiancho de contacto (m)
c: longitud característica del semiancho de la zona adherida (m)
44
Descripción del Software
Esfuerzo: se refiere al esfuerzo axil al que se someterá al plano, es decir, tracción o compresión.
Variable: se muestra si se quiere representar una distribución de tensiones o de deformaciones.
Fuerza: hace referencia a la fuerza externa que provoca la distribución de tensiones o
deformaciones que se quiere representar en el gráfico. Se elegirá entre normal, tangencial bajo
estado de deslizamiento global, tangencial distribución correctiva, tangencial total o la suma de
todas las acciones externas, incluyendo la tensión externa.
Componente: muestra que componente del tensor de tensiones o deformaciones se representa. Si
se grafica una tensión habrá que elegir entre las componentes xx, yy, xy o zz. En el caso de que se
represente una deformación se elegirá entre las componentes xx, yy o xy
De la descripción de esta tabla se observa cómo posee las mismas variables de entrada y salida pero se ha
eliminado la variable Divisiones segmento y se ha sustituido por las dos variables, Divisiones eje X y
Divisiones eje Y, para la determinación de los distintos puntos a analizar. Respecto a las especificaciones,
conserva la especificación Esfuerzo y se añaden tres más, Variable, Fuerza y Componente.
Después de pulsar este botón, para obtener la matriz de resultados, se tiene que responder una serie de
preguntas que aparecerán de forma consecutiva en varias pantallas emergentes, las cuales se pueden observar
en las figuras que abarcan desde la Figuera 4-50 a la Figura 4-54. Estas preguntas corresponden a las
especificaciones descritas en el apartado de la explicación de la tabla de variables. Para la consecución de los
resultados hay que elegir, en cada una de las preguntas planteadas, una de las opciones expuestas en cada
pantalla emergente. Cada vez que se elija una de las opciones, la celda correspondiente a la especificación se
rellena de forma automática.
44
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 45
En esta hoja se mostrarán los mismos mensajes de error que en la hoja Principal si no se elige una de las
opciones requeridas en cada una de las especificaciones o se introduce un valor inadmisible en las variables de
entrada.
A diferencia de la hoja Principal, la creación del gráfico de superficie está programado en la macro ejecutada,
por lo que una vez aparezcan los valores en la hoja Contour, la gráfica deseada se generará automáticamente
46
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 47
4.2.5 Ejemplo
Al igual que en el apartado anterior, se van a añadir al documento una serie de gráficas que muestran los
resultados numéricos obtenidos a partir de las ecuaciones analíticas del problema de fretting para el caso de un
contacto entre un cilindro y un plano. Para ello se han elegido los datos de entrada que aparecen a
continuación:
E = 2.10E11 Pa
ν = 0.3
R = 0.1 m
P = 1.00E6 N
Q = 1.00E5 N
σ = 1.10E8 Pa
μ = 0.7
xi = -2.1009E-3 m
yi = 0 m
xf = 2.1009E-3 m
yf = -1.0505E-3 m
Divisiones eje X = 30
Divisiones eje Y = 30
El resultado de las variables de salida para las anteriores variables de entrada es:
a = 1.0505E-3 m
c = 9.7254E-4 m
e = 6.8095E-5 m
p0 = 6.0604E8 Pa
Las gráficas de tensiones y deformaciones para el estado de compresión abarcan desde la Figura 4-56 hasta la
Figura 4-62, mientras que las gráficas referentes al estado de tracción se muestran de la Figura 4-63 hasta la
Figura 4-69.
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
-0.153--0.074
-0.231--0.153
-0.167
-0.309--0.231
-0.387--0.309
-0.333 -0.465--0.387
-0.544--0.465
y/a
Figura 4-56. Gráfico de superficie de la tensión σxx para el estado de compresión adimensionalizada por p0
48
Descripción del Software
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
-0.071-0.006
-0.149--0.071
-0.167
-0.226--0.149
-0.304--0.226
-0.333 -0.381--0.304
-0.459--0.381
y/a
Figura 4-57. Gráfico de superficie de la tensión σyy para el estado de compresión adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
0.234-0.272
0.196-0.234
-0.167
0.157-0.196
0.119-0.157
-0.333 0.081-0.119
0.043-0.081
y/a
Figura 4-58. Gráfico de superficie de la tensión σxy para el estado de compresión adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
-0.069--0.022
-0.115--0.069
-0.167
-0.162--0.115
-0.208--0.162
-0.333 -0.255--0.208
-0.302--0.255
y/a
Figura 4-59. Gráfico de superficie de la tensión σzz para el estado de compresión adimensionalizada por p0
48
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 49
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 -2.10E-13-6.75E-14
-4.88E-13--2.10E-13
-0.167 -7.66E-13--4.88E-13
-1.04E-12--7.66E-13
-0.333 -1.32E-12--1.04E-12
-1.60E-12--1.32E-12
y/a
-0.667 -2.43E-12--2.16E-12
-2.71E-12--2.43E-12
-2.99E-12--2.71E-12
-0.833
-3.27E-12--2.99E-12
-3.55E-12--3.27E-12
-1.000
Figura 4-60. Gráfico de superficie de la deformación εxx para el estado de compresión adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 8.31E-13-1.15E-12
5.08E-13-8.31E-13
-0.167 1.86E-13-5.08E-13
-1.37E-13-1.86E-13
-0.333 -4.59E-13--1.37E-13
-7.82E-13--4.59E-13
y/a
-0.667 -1.75E-12--1.43E-12
-2.07E-12--1.75E-12
-2.40E-12--2.07E-12
-0.833
-2.72E-12--2.40E-12
-3.04E-12--2.72E-12
-1.000
Figura 4-61. Gráfico de superficie de la deformación εyy para el estado de compresión adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 1.45E-12-1.68E-12
1.21E-12-1.45E-12
-0.167 9.74E-13-1.21E-12
7.37E-13-9.74E-13
-0.333 5.01E-13-7.37E-13
2.64E-13-5.01E-13
y/a
-0.667 -4.46E-13--2.09E-13
-6.83E-13--4.46E-13
-9.19E-13--6.83E-13
-0.833
-1.16E-12--9.19E-13
-1.39E-12--1.16E-12
-1.000
Figura 4-62. Gráfico de superficie de la deformación εxy para el estado de compresión adimensionalizada por p0
50
Descripción del Software
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
0.375-0.482
0.268-0.375
-0.167
0.161-0.268
0.054-0.161
-0.333 -0.053-0.054
-0.160--0.053
y/a
Figura 4-63. Gráfico de superficie de la tensión σxx para el estado de tracción adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
-0.077-0.000
-0.154--0.077
-0.167
-0.231--0.154
-0.308--0.231
-0.333 -0.385--0.308
-0.462--0.385
y/a
Figura 4-64. Gráfico de superficie de la tensión σyy para el estado de tracción adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
0.189-0.227
0.150-0.189
-0.167
0.112-0.150
0.074-0.112
-0.333 0.036-0.074
-0.003-0.036
y/a
Figura 4-65. Gráfico de superficie de la tensión σxy para el estado de tracción adimensionalizada por p0
50
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 51
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000
0.089-0.145
0.034-0.089
-0.167
-0.021-0.034
-0.076--0.021
-0.333 -0.131--0.076
-0.187--0.131
y/a
Figura 4-66. Gráfico de superficie de la tensión σzz para el estado de tracción adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 1.77E-12-2.09E-12
1.45E-12-1.77E-12
-0.167 1.13E-12-1.45E-12
8.04E-13-1.13E-12
-0.333 4.84E-13-8.04E-13
1.63E-13-4.84E-13
y/a
-0.667 -8.00E-13--4.79E-13
-1.12E-12--8.00E-13
-1.44E-12--1.12E-12
-0.833
-1.76E-12--1.44E-12
-2.08E-12--1.76E-12
-1.000
Figura 4-67. Gráfico de superficie de la deformación εxx para el estado de tracción adimensionalizada por p0
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 -1.57E-13-1.18E-13
-4.31E-13--1.57E-13
-0.167 -7.05E-13--4.31E-13
-9.80E-13--7.05E-13
-0.333 -1.25E-12--9.80E-13
-1.53E-12--1.25E-12
y/a
-0.667 -2.35E-12--2.08E-12
-2.63E-12--2.35E-12
-2.90E-12--2.63E-12
-0.833
-3.17E-12--2.90E-12
-3.45E-12--3.17E-12
-1.000
Figura 4-68. Gráfico de superficie de la deformación εyy para el estado de tracción adimensionalizada por p0
52
Descripción del Software
x/a
-2.00 -1.33 -0.67 0.00 0.67 1.33 2.00
0.000 1.17E-12-1.40E-12
9.31E-13-1.17E-12
-0.167 6.94E-13-9.31E-13
4.57E-13-6.94E-13
-0.333 2.21E-13-4.57E-13
-1.60E-14-2.21E-13
y/a
-0.667 -7.26E-13--4.89E-13
-9.63E-13--7.26E-13
-1.20E-12--9.63E-13
-0.833
-1.44E-12--1.20E-12
-1.67E-12--1.44E-12
-1.000
Figura 4-69. Gráfico de superficie de la deformación εxy para el estado de tracción adimensionalizada por p0
52
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 53
5 CONCLUSIONES
El objetivo principal de este trabajo ha sido crear y diseñar una herramienta sencilla para la obtención de las
tensiones y deformaciones del contacto entre un cilindro y un plano bajo las condiciones de carga que
provocan la fatiga por fretting. Para ello se ha elegido el software Excel por poseer unas características que
hacen que sea fácil de utilizar y tener una interfaz intuitiva. El software diseñado se ha creado combinando el
manejo de hojas de Excel y programación de macros usando el lenguaje Visual Basic. Aunque se haya
utilizado este lenguaje para el desarrollo del trabajo, el usuario que vaya a utilizar el programa no necesita
tener ningún conocimiento acerca de programación en Visual Basic. Esto ocurre gracias a la posibilidad de
ejecutar las macros programadas en cada hoja pulsando un botón de control de formulario.
El problema de cálculo de tensiones y deformaciones se podría haber resuelto usando softwares como Matlab
o Maple. Estos poseen una potencia de cálculo numérico mucho mayor que Excel. Se ha comprobado que,
para realizar el mismo cálculo de tensiones para la hoja Principal, utilizando el programa Mapple se obtenía
un tiempo de ejecución de resultados ocho veces menor que el necesitado por Excel. Este hecho se ha
estimado insuficiente en la valoración de usar un software distinto de Excel, puesto que, con ningún otro se
habría podido conseguir una herramienta tan fácil de usar como la construida en este trabajo. A parte de
realizar un cálculo numérico, también se han diseñado una serie de gráficas en las que representar dichos
cálculos. En este ámbito, Excel también se presenta como una mejor opción frente a otros programas al poseer
un amplio abanico gráficos a los que se puede acceder de manera intuitiva. En cambio, en los programas
anteriormente mencionados, para poder representar los valores de las tensiones y deformaciones, es necesario
conocer los comandos necesarios para la realización de gráficas. Por esto, si el usuario desconoce el lenguaje
de programación de Maple o Matlab, le será imposible construir las gráficas.
Otra de las características que hacen que el software desarrollado sea sencillo es la actualización automática de
los valores que se representan en las gráficas de la hoja Principal. Esto permite cambiar cualquier dato de las
variables de entrada y solamente tener que volver a pulsar el botón Cilindro_Plano para observar el efecto que
este cambio tiene en las gráficas. Hay que señalar que se puede tener graficadas todas las columnas de
tensiones y deformaciones de esta hoja de forma simultánea. Para visualizar una gráfica nueva hay que seguir
una serie de pasos muy sencillos descritos en el apartado 4.1.5. Respecto a la hoja Contour, hay que decir que
solo se pueden obtener los valores para representar una sola gráfica. Esto es así porque se ha programado para
que en la hoja aparezca solamente una matriz de datos. Se ha tomado esta decisión ya que, si se quisiera tener
las gráficas de todas las tensiones y deformaciones de forma simultánea, el tiempo de ejecución de la macro
sería extremadamente elevado.
A parte de las anteriores, una de las claves para obtener un software fácilmente manejable, ha sido la forma de
determinar el estado de aplicación de cargas que se desea resolver en cada una de las hojas. Esto se refiere al
estado de tracción o compresión y a lo que implica esta elección en los sentidos de la fuerza tangencial
aplicada al cilindro y tensión aplicada al plano. Este método de introducción de requisitos se ha seguido
también en la hoja Contour para seleccionar la tensión o deformación que se quiere representar en el gráfico de
superficie. Para ello el usuario tiene que responder a una serie de preguntas que le aparecen de forma
automática. Sus respuestas dictaminan qué componente del tensor de tensiones o deformaciones se muestra en
el gráfico, eligiendo también si dicha componente es provocada por la fuerza normal, tangencial o por la suma
de todas las acciones externas. Hay que destacar que las opciones para responder cada una de las preguntas
aparecen en la misma pantalla emergente, por lo que simplemente hay que apreciar cual de ellas se desea e
introducir el valor correspondiente.
El apartado cuarto ha sido redactado como una guía para aprender a usar el software construido. Con esto se
54
Conclusiones
pretende conseguir que, aquel que lo utilice, obtenga los resultados numéricos de las tensiones y
deformaciones del contacto propuesto, independientemente de que se conozca o no la teoría relacionada con el
problema del fretting para interpretar los resultados. Aunque se desconozca dicha teoría, algo que si es
realmente importante para poder obtener resultados es el significado de cada una de las variables de entrada.
Es necesario introducir valores adecuados en función de la variable que se trate, por lo que hay que tener un
conocimiento básico de ellas para saber el rango de valores posible de cada una de ellas.
Uno de los objetivos por los que se ha construido este software ha sido para poseer una herramienta que realice
el cálculo numérico de las expresiones analíticas y poder comparar los resultados obtenidos teóricamente con
los resultados obtenidos en experimentos con contactos cilíndricos. Los conocimientos que se han averiguado
sobre la fatiga por fretting desde los comienzos de su estudio, han sido gracias a la realización de estudios
teóricos y experimentales paralelamente. La comparación de ambos proporciona un mejor entendimiento de
este tipo de fatiga.
Otro de los posibles usos que se puede conseguir del desarrollo de este software, es que, al haber otros tipos
de contactos de los que se conocen las expresiones analíticas de las distribuciones de tensiones, sustituyendo
las ecuaciones que gobiernan las distribuciones del contacto cilíndrico por aquellas que determinan las
distribuciones en otro tipo de contacto, se puede conseguir la obtención de otro software para otra clase de
contacto. Las ecuaciones de las mismas para el contacto entre una esfera y un plano se pueden encontrar en los
trabajos realizados por Hamilton y Sackfield. Por tanto, partiendo del software diseñado en este trabajo, se
podría modificar el código de Visual Basic para obtener los resultados de las tensiones y deformaciones para el
caso esférico. Sería necesario estudiar este nuevo tipo de contacto a fondo para tener en cuenta todas las
particularidades que este presente y difieran del contacto cilíndrico.
Por último, la creación de este archivo Excel se considera el primer paso del desarrollo de un software mucho
más complejo. El fin último que se pretende conseguir es obtener un software que permita realizar un conjunto
de cálculos relacionados con la fatiga, como son el cálculo de la vida del componente a fatiga o las tensiones
que se crean en el borde de la grieta. En definitiva, realizar un estudio en profundidad de la iniciación y
propagación de la grieta para determinar cuándo, cómo y dónde falla el componente en estudio, pero en un
software sencillo.
54
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 55
Sub Cilindro_Plano()
''VARIABLES DE LA TABLA''
If Q > mu * P Then
MsgBox "El valor de Q excede el límite, es decir, Q>mu*P" , , "Error"
Exit Sub
End If
'Especificación
esfuerzo = InputBox("¿A qué esfuerzo quiere someter el plano? Compresión=1 o Tracción=2", "Introduzca esfuerzo")
56
Anexos: Macros de Visual Basic
If esfuerzo = 1 Then
sigma = -sigma
Cells(27, 3).Value = "Compresión" ' Q>0 , sigma<0
End If
If esfuerzo = 2 Then
Q = -Q
Cells(27, 3) = "Tracción" ' Q<0 , sigma>0
End If
Else
MsgBox "Tiene que elegir una de las opciones" , , "Error"
Exit Sub
End If
j=2
Do Until j = n + 3
Cells(j, 5).Value = j – 1
56
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 57
Loop
''COLUMNAS DE COORDENADAS''
Length = n + 2
ReDim x(Length)
ReDim y(Length)
ReDim s1(Length)
ReDim s2(Length)
For j = 2 To Length
i=6
Do Until i = 10
With Cells(j, i).Interior
.Pattern = xlSolid
.PatternColorIndex = xlAutomatic
.ThemeColor = xlThemeColorAccent6
.TintAndShade = 0.399975585192419
.PatternTintAndShade = 0
End With
Cells(j, i).NumberFormat = "0.0000"
With Cells(j, i).Borders(xlEdgeLeft)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i).Borders(xlEdgeRight)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i)
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlBottom
.WrapText = False
.Orientation = 0
.AddIndent = False
.IndentLevel = 0
.ShrinkToFit = False
.ReadingOrder = xlContext
.MergeCells = False
End With
i=i+1
Loop
Next j
i=6
Do Until i = 10
j = Length + 1
Do Until IsEmpty(Cells(j, i))
Cells(j, i).Clear
j=j+1
Loop
i=i+1
Loop
''COLUMNAS DE TENSIONES''
ReDim sigmaP_xx(Length)
ReDim sigmaP_yy(Length)
ReDim sigmaP_xy(Length)
ReDim sigmaP_zz(Length)
ReDim sigmaQ_xx_d_total(Length)
ReDim sigmaQ_yy_d_total(Length)
ReDim sigmaQ_xy_d_total(Length)
ReDim sigmaQ_zz_d_total(Length)
ReDim sigmaQ_xx_d_parcial(Length)
ReDim sigmaQ_yy_d_parcial(Length)
ReDim sigmaQ_xy_d_parcial(Length)
ReDim sigmaQ_zz_d_parcial(Length)
ReDim sigmaQ_xx(Length)
ReDim sigmaQ_yy(Length)
ReDim sigmaQ_xy(Length)
ReDim sigmaQ_zz(Length)
ReDim sigma_total_xx(Length)
ReDim sigma_total_yy(Length)
ReDim sigma_total_xy(Length)
ReDim sigma_total_zz(Length)
For j = 2 To Length
If y(j) = 0 Then
Else
sigmaP_xx(j) = -y(j) / A * (2 - (s1(j) / A) / Sqr(1 + (s1(j) / A) ^ 2) - Sqr(1 + (s1(j) / A) ^ 2) / (s1(j) / A) - (x(j) / A) ^ 2 * (s1(j) / A) ^ 3 / (1 +
(s1(j) / A) ^ 2) ^ (3 / 2) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2))
sigmaP_yy(j) = (y(j) / A) ^ 3 * Sqr(1 + (s1(j) / A) ^ 2) / (s1(j) / A) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2)
sigmaP_xy(j) = (x(j) / A) * (y(j) / A) ^ 2 * (s1(j) / A) / Sqr(1 + (s1(j) / A) ^ 2) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2)
End If
58
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 59
If y(j) = 0 Then
Else
sigmaQ_xx_d_total(j) = Sgn(Q) * (-2 * (x(j) / A) * (1 - (s1(j) / A) / Sqr(1 + (s1(j) / A) ^ 2)) + (x(j) / A) * (y(j) / A) ^ 2 * (s1(j) / A) / Sqr(1 +
(s1(j) / A) ^ 2) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2))
sigmaQ_yy_d_total(j) = Sgn(Q) * (x(j) / A) * (y(j) / A) ^ 2 * (s1(j) / A) / Sqr(1 + (s1(j) / A) ^ 2) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2)
sigmaQ_xy_d_total(j) = Sgn(Q) * (y(j) / A) * (2 - (s1(j) / A) / Sqr(1 + (s1(j) / A) ^ 2) - Sqr(1 + (s1(j) / A) ^ 2) / (s1(j) / A) - (x(j) / A) ^ 2 *
(s1(j) / A) ^ 3 / (1 + (s1(j) / A) ^ 2) ^ (3 / 2) / ((s1(j) / A) ^ 4 + (y(j) / A) ^ 2))
End If
If y(j) = 0 Then
Else
sigmaQ_xx_d_parcial(j) = -Sgn(Q) * c / A * (-2 * (x(j) + exc) / c * (1 - (s2(j) / c) / Sqr(1 + (s2(j) / c) ^ 2)) + ((x(j) + exc) / c) * (y(j) / c) ^ 2
* (s2(j) / c) / Sqr(1 + (s2(j) / c) ^ 2) / ((s2(j) / c) ^ 4 + (y(j) / c) ^ 2))
sigmaQ_yy_d_parcial(j) = -Sgn(Q) * c / A * ((x(j) + exc) / c) * (y(j) / c) ^ 2 * (s2(j) / c) / Sqr(1 + (s2(j) / c) ^ 2) / ((s2(j) / c) ^ 4 + (y(j) / c)
^ 2)
60
Anexos: Macros de Visual Basic
sigmaQ_xy_d_parcial(j) = -Sgn(Q) * c / A * (y(j) / c * (2 - (s2(j) / c) / Sqr(1 + (s2(j) / c) ^ 2) - Sqr(1 + (s2(j) / c) ^ 2) / (s2(j) / c) - ((x(j) +
exc) / c) ^ 2 * (s2(j) / c) ^ 3 / (1 + (s2(j) / c) ^ 2) ^ (3 / 2) / ((s2(j) / c) ^ 4 + (y(j) / c) ^ 2)))
End If
i = 10
Do Until i = 26
With Cells(j, i).Interior
.Pattern = xlSolid
.PatternColorIndex = xlAutomatic
.Color = 10092390
.TintAndShade = 0
.PatternTintAndShade = 0
End With
Cells(j, i).NumberFormat = "0.00000"
With Cells(j, i).Borders(xlEdgeLeft)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i).Borders(xlEdgeRight)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i)
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlBottom
.WrapText = False
.Orientation = 0
.AddIndent = False
.IndentLevel = 0
.ShrinkToFit = False
.ReadingOrder = xlContext
.MergeCells = False
60
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 61
End With
i=i+1
Loop
Next j
i = 10
Do Until i = 26
j=n+3
Do Until IsEmpty(Cells(j, i))
Cells(j, i).Clear
j=j+1
Loop
i=i+1
Loop
''COLUMNAS DE DEFORMACIONES''
ReDim defP_xx(Length)
ReDim defP_yy(Length)
ReDim defP_xy(Length)
ReDim defQ_xx_d_total(Length)
ReDim defQ_yy_d_total(Length)
ReDim defQ_xy_d_total(Length)
ReDim defQ_xx_d_parcial(Length)
ReDim defQ_yy_d_parcial(Length)
ReDim defQ_xy_d_parcial(Length)
ReDim defQ_xx(Length)
ReDim defQ_yy(Length)
ReDim defQ_xy(Length)
ReDim def_total_xx(Length)
ReDim def_total_yy(Length)
ReDim def_total_xy(Length)
For j = 2 To Length
i = 26
Do Until i = 38
With Cells(j, i).Interior
.Pattern = xlSolid
.PatternColorIndex = xlAutomatic
.Color = 16756655
.TintAndShade = 0
.PatternTintAndShade = 0
End With
Cells(j, i).NumberFormat = "0.000E+00"
With Cells(j, i).Borders(xlEdgeLeft)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i).Borders(xlEdgeRight)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(j, i)
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlBottom
.WrapText = False
.Orientation = 0
.AddIndent = False
.IndentLevel = 0
.ShrinkToFit = False
.ReadingOrder = xlContext
.MergeCells = False
End With
i=i+1
Loop
Next j
62
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 63
i = 26
Do Until i = 38
j=n+3
Do Until IsEmpty(Cells(j, i))
Cells(j, i).Clear
j=j+1
Loop
i=i+1
Loop
End Sub
Sub ContourPlot()
''VARIABLES DE LA TABLA''
If Q > mu * P Then
MsgBox "El valor de Q excede el límite, es decir, Q>mu*P" , , "Error"
Exit Sub
End If
If yi < yf Then
yi = Cells(17, 3).Value
yf = Cells(15, 3).Value
End If
'Especificación
esfuerzo = InputBox("¿A qué esfuerzo quiere someter el plano? Compresión=1 o Tracción=2 ", "Introduzca esfuerzo")
If esfuerzo = 1 Then
sigma = -sigma
Cells(28, 3).Value = "Compresión" ' Q>0 , sigma<0
End If
If esfuerzo = 2 Then
Q = -Q
Cells(28, 3) = "Tracción" ' Q<0 , sigma>0
End If
Else
MsgBox "Tiene que elegir una de las opciones" , , "Error"
Exit Sub
End If
''MATRIZ''
fila_x = 1
col_x = 7
fila_y = fila_x + 1
col_y = col_x - 1
length_x = nx + col_x
length_y = ny + fila_y
ReDim x(length_x)
ReDim y(length_y)
i = col_x
Do Until i = nx + col_x + 1
64
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 65
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(fila_x, i).Borders(xlBottom)
.LineStyle = xlContinuous
.ColorIndex = xlAutomatic
.TintAndShade = 0
.Weight = xlThin
End With
With Cells(fila_x, i)
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlBottom
.WrapText = False
.Orientation = 0
.AddIndent = False
.IndentLevel = 0
.ShrinkToFit = False
.ReadingOrder = xlContext
.MergeCells = False
End With
i=i+1
Loop
'Borrado de datos si el número de puntos en el eje X es menor que el anterior uso de la macro
j = fila_y
Do Until j = ny + fila_y + 1
.HorizontalAlignment = xlCenter
.VerticalAlignment = xlBottom
.WrapText = False
.Orientation = 0
.AddIndent = False
.IndentLevel = 0
.ShrinkToFit = False
.ReadingOrder = xlContext
.MergeCells = False
End With
j=j+1
Loop
'Borrado de datos si el número de puntos en el eje Y es menor que el anterior uso de la macro
''CÁLCULO DE TENSIONES''
s1(i, j) = Sqr(0.5 * (-(A ^ 2 - x(i) ^ 2 - y(j) ^ 2) + Sqr((A ^ 2 - x(i) ^ 2 - y(j) ^ 2) ^ 2 + 4 * A ^ 2 * y(j) ^ 2)))
s2(i, j) = Sqr(0.5 * (-(c ^ 2 - (x(i) + exc) ^ 2 - y(j) ^ 2) + Sqr((c ^ 2 - (x(i) + exc) ^ 2 - y(j) ^ 2) ^ 2 + 4 * c ^ 2 * y(j) ^ 2)))
If y(j) = 0 Then
66
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 67
sigmaP_xy(i, j) = 0
Else
sigmaP_xx(i, j) = -Sqr(1 - (x(i) / A) ^ 2)
sigmaP_yy(i, j) = -Sqr(1 - (x(i) / A) ^ 2)
sigmaP_xy(i, j) = 0
End If
Else
sigmaP_xx(i, j) = -y(j) / A * (2 - (s1(i, j) / A) / Sqr(1 + (s1(i, j) / A) ^ 2) - Sqr(1 + (s1(i, j) / A) ^ 2) / (s1(i, j) / A) - (x(i) / A) ^ 2 * (s1(i, j) /
A) ^ 3 / (1 + (s1(i, j) / A) ^ 2) ^ (3 / 2) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2))
sigmaP_yy(i, j) = (y(j) / A) ^ 3 * Sqr(1 + (s1(i, j) / A) ^ 2) / (s1(i, j) / A) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2)
sigmaP_xy(i, j) = (x(i) / A) * (y(j) / A) ^ 2 * (s1(i, j) / A) / Sqr(1 + (s1(i, j) / A) ^ 2) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2)
End If
If y(j) = 0 Then
Else
sigmaQ_xx_d_total(i, j) = Sgn(Q) * (-2 * (x(i) / A) * (1 - (s1(i, j) / A) / Sqr(1 + (s1(i, j) / A) ^ 2)) + (x(i) / A) * (y(j) / A) ^ 2 * (s1(i, j) / A) /
Sqr(1 + (s1(i, j) / A) ^ 2) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2))
sigmaQ_yy_d_total(i, j) = Sgn(Q) * (x(i) / A) * (y(j) / A) ^ 2 * (s1(i, j) / A) / Sqr(1 + (s1(i, j) / A) ^ 2) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2)
sigmaQ_xy_d_total(i, j) = Sgn(Q) * (y(j) / A * (2 - (s1(i, j) / A) / Sqr(1 + (s1(i, j) / A) ^ 2) - Sqr(1 + (s1(i, j) / A) ^ 2) / (s1(i, j) / A) - (x(i) /
A) ^ 2 * (s1(i, j) / A) ^ 3 / (1 + (s1(i, j) / A) ^ 2) ^ (3 / 2) / ((s1(i, j) / A) ^ 4 + (y(j) / A) ^ 2)))
End If
If y(j) = 0 Then
End If
Else
sigmaQ_xx_d_parcial(i, j) = -Sgn(Q) * c / A * (-2 * (x(i) + exc) / c * (1 - (s2(i, j) / c) / Sqr(1 + (s2(i, j) / c) ^ 2)) + ((x(i) + exc) / c) * (y(j) /
c) ^ 2 * (s2(i, j) / c) / Sqr(1 + (s2(i, j) / c) ^ 2) / ((s2(i, j) / c) ^ 4 + (y(j) / c) ^ 2))
sigmaQ_yy_d_parcial(i, j) = -Sgn(Q) * c / A * ((x(i) + exc) / c) * (y(j) / c) ^ 2 * (s2(i, j) / c) / Sqr(1 + (s2(i, j) / c) ^ 2) / ((s2(i, j) / c) ^ 4 +
(y(j) / c) ^ 2)
sigmaQ_xy_d_parcial(i, j) = -Sgn(Q) * c / A * (y(j) / c * (2 - (s2(i, j) / c) / Sqr(1 + (s2(i, j) / c) ^ 2) - Sqr(1 + (s2(i, j) / c) ^ 2) / (s2(i, j) / c)
- ((x(i) + exc) / c) ^ 2 * (s2(i, j) / c) ^ 3 / (1 + (s2(i, j) / c) ^ 2) ^ (3 / 2) / ((s2(i, j) / c) ^ 4 + (y(j) / c) ^ 2)))
End If
''CÁLCULO DE DEFORMACIONES''
68
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 69
Next i
Next j
'Borrado de datos si el número de puntos en el eje X o eje Y es menor que el anterior uso de la macro
j = fila_y
i = col_x
Do Until IsEmpty(Cells(j, col_x))
i = col_x
Do Until IsEmpty(Cells(j, i))
If j > length_y Or i > length_x Then
Cells(j, i).Clear
i=i+1
Else
i=i+1
End If
Loop
j=j+1
Loop
'Especificaciones
If fuerza = 1 Then
Cells(30, 3).Value = "Normal"
End If
If fuerza = 2 Then
tangencial = InputBox("¿A qué distribución está asociado el campo deseado? Tensión
tangencial deslizamiento global=1 Tensión tangencial correctiva=2 Tensión tangencial total=3", "Introduzca
tangencial")
If tangencial = 1 Or tangencial = 2 Or tangencial = 3 Then
If tangencial = 1 Then
Cells(30, 3).Value = "Tangencial deslizamiento global"
End If
If tangencial = 2 Then
Cells(30, 3).Value = "Tangencial correctiva"
End If
If tangencial = 3 Then
Cells(30, 3).Value = "Tangencial deslizamiento parcial"
End If
Else
MsgBox "Tiene que elegir una opción" , , "Error"
Exit Sub
End If
End If
If fuerza = 3 Then
Cells(30, 3).Value = "Ambas"
End If
Else
MsgBox "Tiene que elegir una opción" , , "Error"
Exit Sub
End If
If variable = 1 Then
If fuerza = 1 Then
If componente = 1 Then
Cells(j, i).Value = sigmaP_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = sigmaP_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = sigmaP_xy(i, j)
70
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 71
End If
If componente = 4 Then
Cells(j, i).Value = sigmaP_zz(i, j)
End If
End If
If fuerza = 2 Then
If tangencial = 1 Then
If componente = 1 Then
Cells(j, i).Value = sigmaQ_xx_d_total(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = sigmaQ_yy_d_total(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = sigmaQ_xy_d_total(i, j)
End If
If componente = 4 Then
Cells(j, i).Value = sigmaQ_zz_d_total(i, j)
End If
End If
If tangencial = 2 Then
If componente = 1 Then
Cells(j, i).Value = sigmaQ_xx_d_parcial(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = sigmaQ_yy_d_parcial(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = sigmaQ_xy_d_parcial(i, j)
End If
If componente = 4 Then
Cells(j, i).Value = sigmaQ_zz_d_parcial(i, j)
End If
End If
If tangencial = 3 Then
If componente = 1 Then
Cells(j, i).Value = sigmaQ_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = sigmaQ_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = sigmaQ_xy(i, j)
End If
If componente = 4 Then
Cells(j, i).Value = sigmaQ_zz(i, j)
End If
End If
End If
If fuerza = 3 Then
If componente = 1 Then
Cells(j, i).Value = sigma_total_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = sigma_total_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = sigma_total_xy(i, j)
End If
If componente = 4 Then
Cells(j, i).Value = sigma_total_zz(i, j)
End If
End If
End If
72
Anexos: Macros de Visual Basic
If variable = 2 Then
If fuerza = 1 Then
If componente = 1 Then
Cells(j, i).Value = defP_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = defP_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = defP_xy(i, j)
End If
End If
If fuerza = 2 Then
If tangencial = 1 Then
If componente = 1 Then
Cells(j, i).Value = defQ_xx_d_total(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = defQ_yy_d_total(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = defQ_xy_d_total(i, j)
End If
End If
If tangencial = 2 Then
If componente = 1 Then
Cells(j, i).Value = defQ_xx_d_parcial(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = defQ_yy_d_parcial(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = defQ_xy_d_parcial(i, j)
End If
End If
If tangencial = 3 Then
If componente = 1 Then
Cells(j, i).Value = defQ_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = defQ_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = defQ_xy(i, j)
End If
End If
End If
If fuerza = 3 Then
If componente = 1 Then
Cells(j, i).Value = def_total_xx(i, j)
End If
If componente = 2 Then
Cells(j, i).Value = def_total_yy(i, j)
End If
If componente = 3 Then
Cells(j, i).Value = def_total_xy(i, j)
End If
End If
End If
If variable = 1 Then
72
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 73
End If
End Sub
74
Software Excel: tensiones y deformaciones en fatiga por fretting con contacto cilíndrico 75
7 BIBLIOGRAFÍA