Clculo del Factor de Intensidad de Tensiones en una viga

sometida a flexin en tres puntos con el programa ABAQUS 6.9

Introduccin
ABAQUS es un programa destinado a resolver problemas de ciencias e ingeniera y
est basado en el mtodo de los elementos finitos. El programa puede resolver casi
todo tipo de problemas, desde un simple anlisis lineal hasta simulaciones complejas
no lineales. Abaqus posee una extensa librera de elementos que permite modelar
virtualmente cualquier geometra, as como una extensa lista de modelos que simulan
el comportamiento de una gran cantidad de materiales, permitiendo su aplicabilidad en
distintas reas de ingeniera.
El Mtodo de los Elementos Finitos
Con excepcin de las estructuras de barras, la mayor parte de las estructuras en
ingeniera son de naturaleza continua y, por lo tanto, su comportamiento no puede
expresarse en forma precisa en funcin de un nmero pequeo de variables discretas.
Un anlisis riguroso de dichas estructuras precisa la integracin de las ecuaciones
diferenciales que expresan el equilibrio de un elemento diferencial genrico de las
mismas. Ejemplos de estas estructuras continuas son comunes en las ingenieras
civil, mecnica, aeronutica y naval, y entre las ms usuales podemos citar las placas,
depsitos, cubiertas, puentes, presas, carroceras de vehculos, fuselajes de aviones,
cascos de barcos, etc.
Aunque las estructuras continuas son inherentemente tridimensionales en algunos
casos su comportamiento puede describirse adecuadamente por modelos
matemticos uni o bidimensionales. As ocurre, por ejemplo, con los problemas de
flexin en placas, en los que el anlisis se limita al estudio de la deformacin del plano
medio de la placa, y con todas las estructuras en las que puede hacerse uso de las
hiptesis simplificativas de la elasticidad bidimensional o de revolucin (presas,
tneles, depsitos, etc).
El Mtodo de los Elementos Finitos es hoy da el procedimiento ms potente para el
anlisis de estructuras de carcter uni, bi o tridimensional sometidas a las acciones
exteriores ms diversas. La gran analoga existente entre los conceptos del anlisis
matricial de estructuras de barras y los del mtodo de los elementos finitos facilitan en
gran manera el estudio de ste.
Para visualizar la implementacin del mtodo de los elementos finitos plantearemos la
aplicacin del mismo a un ejemplo concreto como es el clculo de un puente (Figura
1).

Figura 1

La resolucin del problema consta de las siguientes etapas:

Etapa I: A partir de la realidad fsica del puente, sus apoyos y tipos de cargas que
sobre l acten, es necesario primeramente seleccionar un modelo matemtico
apropiado para describir el comportamiento de la estructura. Por ejemplo podra
utilizarse la teora de lminas planas, lminas curvas, o de la elasticidad
tridimensional. Tambin hay que definir con detalle las propiedades mecnicas de los
materiales del puente y el carcter de la deformacin del mismo (pequeos o grandes
movimientos, anlisis esttico o dinmico, etc). As mismo, para el planteamiento de
las ecuaciones de equilibrio haremos uso del Principio de los Trabajos Virtuales (PTV).
Etapa II: Una vez seleccionado el modelo matemtico se procede a discretizar la
estructura en porciones no intersectantes entre s, denominadas elementos finitos,
dentro de los cuales se interpolan las variables principales en funcin de sus valores
en una serie de puntos discretos del elemento denominados nodos. Los elementos
se conectan entre s por los nodos situados en sus contornos. La malla de elementos
finitos puede, por ejemplo, estar constituda, por elementos de diferente geometra, por
ejemplo elementos bidimensionales rectangulares y triangulares. La etapa de
discretizacin constituye una parte esencial de la fase de preproceso que suele incluir
tambin la representacin grfica de la malla de elementos finitos.
Etapa III: A partir de la aplicacin del PTV se obtienen las matrices de rigidez K(e) y el
vector de cargas f(e) para cada elemento. El clculo de K(e) y f(e) es ms complejo que
en estructuras de barras al intervenir integrales sobre el dominio uni, bi o
tridimensional del elemento.
Etapa IV: Se procede al ensamblaje de las matrices de rigidez y el vector de carga
elementales en la matriz de rigidez global de toda la malla de elementos finitos K y el
vector de cargas sobre los nodos f, respectivamente.

Etapa V: El sistema de ecuaciones resultante Ka = f se resuelve para calcular las

variables incgnitas (movimiento de todos los nodos de la malla) a, utilizando uno
cualquiera de los mtodos conocidos para la solucin de ecuaciones algebraicas
simultneas lineales.
Etapa VI: Una vez calculados los movimientos nodales a se pueden calcular las
deformaciones y, seguidamente, las tensiones en cada elemento, as como las
reacciones en los nodos con movimiento prescritos.
Etapa VII: Para obtener la solucin en las etapas III-VI es necesario utilizar un software
que emplee el mtodo de los elementos finitos, lo cual puede llevarse a acabo a partir
de un programa comercial o bien de uno desarrollado al respecto.
Etapa VIII: Obtenidos los resultados numricos, la etapa siguiente en la interpretacin
y presentacin de los mismos, Para ello suele hacerse uso de tcnicas grficas que
facilitan dicha labor (postproceso).
Etapa IX: Una vez estudiados los resultados, el tcnico analista puede plantearse
efectuar varias modificaciones en cualquiera de las etapas anteriores. As, por
ejemplo, puede encontrar que la teora de clculo de estructuras inicialmente adoptada
es inapropiada y consiguientemente debe modificarse. Por otro lado, la malla de
elementos finitos utilizada en el anlisis puede ser demasiado grosera para reproducir
la distribucin de desplazamientos o tensiones

Desarrollo del ejemplo

En la Figura 2 se presenta la geometra y el estado de cargas de una viga sometida a
flexin en tres puntos.
F= 1N

Figura 2
Materiales:
E= 21000 MPa
n= 0.33
En el modelado se emplear condicin de simetra
1)Crear un archivo nuevo
Para resolver el presente problema trabajaremos en la plataforma CAE
Al ingresar a ABAQUS/CAE elegir la opcin CREATE MODEL DATA BASE

Figura 3

Luego de elegir esa opcin se abrir el entorno de trabajo de ABAQUS, como se

muestra en la Figura 4

Figura 4
2) Crear una parte
Para poder cargar la geometra de la viga lo primero que debemos hacer es crear una
parte, para ello seleccionamos PART en el men desplegable MODULE, como se
muestra en la Figura 5

Figura 5

Una vez que nos encontramos en el men Part hacemos clic en el cono CREATE
PART

y se abrir la siguiente ventana

Figura 6
Una vez all le asignaremos un nombre a la parte (1), en este caso Viga,
seleccionaremos la opcin (2) 2D PLANAR y luego hacemos click en CONTINUE. En
la pantalla aparecer la grilla de trabajo y utilizando las herramientas de dibujo que se
encuentran en la parte izquierda de la misma para dibujar la mitad de la geometra de
la viga, ya que emplearemos condicin de simetra. Una vez terminado el dibujo
hacemos clic en la opcin DONE en la parte inferior izquierda de la pantalla

Herramientas de dibujo

Figura 7
As obtenemos la siguiente geometra

Figura 8

3) Mallado
Un aspecto muy importante en todo programa de elementos finitos es el mallado, ya
que ste tiene una gran influencia en la precisin de los resultados obtenidos con el
modelo.
Para realizar el mallado de la viga debemos desplegar el men MODULE y seleccionar
la opcin MESH, como se observa en la Figura 9

Figura 9
Una vez all debemos particionar la geometra (definir regiones), seleccionar el tipo de
elemento a utilizar, la tcnica de mallado y finalmente realizar el mallado propiamente
dicho.
La particin de la geometra debe llevarse a cabo para poder definir los contornos de
integracin que sern necesarios para poder aplicar la opcin CONTOUR INTEGRAL,
la cual se utiliza para el clculo del factor de intensidad de tensiones (KI), esto se ver
en detalle en el apartado 8)
Para realizar las particiones hacemos click en el cono
el cual nos permite
editar la geometra, por lo que utilizando nuevamente las herramientas de dibujo
definiremos una serie de semicrculos alrededor de la punta de la fisura, como se
observa en la Figura 10. A su vez realizaremos otras particiones adicionales (rectas
verticales y horizontales) que nos permitirn mejorar el aspecto de la malla, a definir
en los pasos siguientes.

Herramientas de dibujo

Figura 10

Luego de particionada la geometra podemos pasar a asignar a la misma el tipo de

elemento a utilizar, para lo cual hacemos click en el cono
y nos pedir que
seleccionemos las regiones a las que queremos asignarle ELEMENT TYPES, para lo
cual hacemos una ventana con el mouse y seleccionamos todas las regiones que
aparecen en pantalla, luego hacemos click en DONE, como se observa en la Figura 11

Figura 11
En segundo lugar se abrir la ventana que se observa en la Figura 12, all
seleccionamos PLANE STRAIN (deformacin plana) y hacemos click en OK

Figura 12

10

El siguiente paso consiste en seleccionar la tcnica de mallado, haciendo click en el

cono
nos pedir que seleccionemos las regiones a las que queremos
asignarles MESH CONTROLS, seleccionamos todas las regiones que tenemos en
pantalla y hacemos click en DONE, de la misma forma que hicimos en la Figura 11.
Luego de realizar esto aparecer la ventana que se observa en la Figura 13, en la cual
seleccionamos la opcin QUAD-DOMINATED en la parte superior y luego
STRUCTURED. Finalmente hacemos click en OK

Figura 13
La geometra aparecer de color verde como se observa en la Figura 14, hacemos
click en DONE

Figura 14

11

Una vez definida la tcnica de mallado procedemos a asignar la cantidad de divisiones

que queramos para cada una de las lneas definidas en la geometra, para ello
mantenemos apretado el cono

y se desplegarn las siguientes opciones

Seleccionamos la opcin marcada con el recuadro haciendo click sobre la misma, sta
nos permitir asignarle a cada lnea la cantidad de divisiones que queramos, para ello
seleccionaos la lnea a la cual vamos a asignarle las divisiones con un click (Figura
15). Si queremos seleccionar ms de una lnea debemos mantener presionada la tecla
SHIFT mientras clickeamos todas las lneas
Seleccionamos lnea

Click en DONE
Figura 15
Luego en la parte inferior de la pantalla nos pide que indiquemos la cantidad de
divisiones, por ejemplo 20, apretamos ENTER y en la pantalla veremos las 20
divisiones solicitadas (Figura 16)

12

Divisiones

Figura 16
Una vez que terminamos de dividir todas las lneas hacemos click en DONE en la
parte inferior de la pantalla (Figura 17).

Figura 17

13

NOTA: Si tenemos que borrar las divisiones en alguna lnea debemos hacer click en el

cono
Una vez definidas las divisiones debemos realizar la malla, para ello mantenemos

presionado el cono

y se desplegarn las siguientes opciones

Seleccionamos la opcin indicada con el recuadro, la cual nos pedir que

seleccionemos la regin que deseamos mallar y luego la confirmemos haciendo click
en DONE en la parte inferior izquierda de la pantalla (Figuras 18a y b)

Seleccionar
Regin a
mallar

Figura 18a

14

Malla

Figura 18b
Repetimos la opresin anterior para todas las regiones en que se dividi la geometra
y obtenemos el mallado que se observa en la Figura 19

Figura 19
Una vez finalizada la operacin de mallado hacemos click en DONE
Un aspecto importante a tener en cuenta con respecto al mallado es el orden en que
se malla, para obtener un mallado lo ms regular posible; para ello en este problema
conviene comenzar a mallar desde los semicrculos ms grandes hacia los ms
pequeos y luego el resto de la geometra, como se indica en la Figura 20

15

6
6

Figura 20

16

NOTA: Si necesitamos borrar alguna regin que hayamos mallado utilizamos el cono

y si queremos borrar todo el mallado en un solo paso utilizaremos

Otro opcin interesante que nos ofrece ABAQUS es verificar el mallado que
realizamos, para determinar si existen errores o warnings con respecto al mismo, para

ello seleccionamos el cono

4) Materiales, secciones, asignacin de propiedades
Para crear el material y la seccin que vamos a utilizar y as poder asignarlos a la
geometra generada debemos desplegar el men MODULE y seleccionar la opcin
PROPERTY, como se observa en la Figura 21

Figura 21
Una vez all lo primero que haremos ser crear un material, para ello hacemos click en
el cono
, el cual se encuentra en la parte superior izquierda del rea de trabajo.
Aparecer la siguiente ventana

17

Indicamos un nombre

Figura 22

Escribimos un nombre para el material, por ejemplo Hormign y luego hacemos click
en la opcin MECHANICAL
ELASTICITY
ELASTIC, a partir de lo
cual el programa nos pedir que indiquemos el mdulo de elasticidad del material, en
este caso 21000 MPa, y el coeficiente de Poisson, en este caso 0.33, como se
observa en la Figura 23. Una vez realizado esto hacemos click en OK

18

Indicamos valores

Figura 23

Ahora debemos crear una seccin para la geometra, por lo que hacemos click en el

cono

y aparecer la ventana que se observa en la Figura 24

Indicamos nombre

Figura 24

19

Indicamos un nombre para la seccin, por ejemplo VIGA, en Category dejamos la

opcin por defecto SOLID, al igual que en Type HOMOGENEUS, hacemos click en
CONTINUE y aparecer la ventana que se observa en la Figura 25

Seleccin del material

Figura 25
En MATERIAL aparecer el material creado previamente, si tuviramos ms de un
material tenemos la posibilidad de seleccionar el que corresponda haciendo click en la
flecha de la izquierda y luego hacemos click en OK para finalizar.
Finalmente el ltimo paso consiste en asignar a la geometra la seccin creada, para

ello hacemos click en el cono

, el cual nos pedir que seleccionemos las
regiones a las que queremos asignarle seccin y material, seleccionamos la
geometra, hacemos click en DONE y aparecer la siguiente ventana

Seleccin de seccin

Figura 26
Elegimos la seccin creada anteriormente y hacemos click en OK y finalmente
clickeamos en DONE en la parte inferior izquierda del rea de trabajo para finalizar.

4) Ensamble
En modelos conformados por varias partes el mdulo ASSEMBLY es utilizado por
ABAQUS para vincular dichas partes. Como en el presente ejemplo slo contamos con
una parte lo nico que debemos hacer es llevarla al mdulo ASSEMBLY, para ello

20

desplegamos el MODULE y seleccionamos ASSEMBLY, como se observa en la Figura

27

Figura 27
Una vez all seleccionamos el cono

Aparecer la siguiente ventana

nica parte del modelo

Figura 28
Aparecer seleccionada la nica parte que existe en el modelo y luego hacemos click
en OK

21

5) Asignacin de condiciones de borde

Las condiciones de borde del presente problema consisten en un apoyo simple en el
extremo izquierdo y una condicin de simetra en el extremo derecho.
Para asignar estas condiciones de borde a la geometra debemos seleccionar el men
LOAD dentro del Module

Figura 29

Una vez all seleccionamos el cono

y aparecer la ventana que se observa en
la Figura 30, all le asignaremos un nombre a la condicin de contorno, luego en
CATEGORY seleccionamos MECHANICAL y en TYPE FOR SELECTED STEP
seleccionamos DISPLACEMENT/ROTATION. Es importante tener en cuenta que el
STEP en que nos encontramos es el INITIAL.
Cuando terminamos hacemos click en CONTINUE

22

Asignar nombre

Figura 30
El siguiente paso consiste en seleccionar el punto, la lnea o superficie a la que
queremos asignarle condiciones de borde. En este ejemplo seleccionamos el nodo
que se indica en la Figura 31 y luego hacemos click en DONE

Figura 31
Luego aparecer la ventana que se observa en la Figura 32 en la que
seleccionaremos el desplazamiento en la direccin U2, indicando as que el

23

desplazamiento vertical se encuentra restringido en dicho punto. Finalmente hacemos

click en OK.

Figura 32
Para asignar la condicin de simetra el procedimiento es el mismo que el descrito
anteriormente, la diferencia se presenta en que en la ventana CREATE BOUNDARY
CONDITION
debemos
seleccionar
la
opcin
SYMMETRY/ANTISYMMETRY/ENCASTRE en la categora TYPE FOR SELECTED
STEP, como se observa en la Figura 33

Asignar nombre

Figura 33
Luego debemos seleccionar la lnea en la que debe aplicarse la condicin de simetra,
como se observa en la Figura 34 y hacemos click en DONE

24

Figura 34
Aparecer la ventana que se observa en la Figura 35 en la que seleccionamos la
opcin XSYMM (U1=UR2=UR3=0), ya que tenemos simetra en la direccin x.
Finalmente hacemos click en OK

Figura 35
El aspecto de las condiciones de borde es el que se observa en la Figura 36

25

Figura 36
6) Definicin del step
Abaqus divide al anlisis en steps, a travs de stos se introducen cambios en las
condiciones de contorno, en las cargas o en el tipo de anlisis. Es importante destacar
que no existe lmite en el nmero de steps dentro de un anlisis.
En el ejemplo que estamos analizando slo definiremos un step esttico para poder
asignar la nica carga actuante en el modelo, para ello en MODULE seleccionamos el
men STEP

Figura 37

Una vez all hacemos click en el cono

y aparecer la siguiente ventana

26

Asignar nombre

Figura 38
En esta ventana asignamos un nombre al step, luego en PROCEDURE TYPE
seleccionamos GENERAL, y dentro de esta categora STATIC, GENERAL. Finalmente
hacemos click en CONTINUE, lo que dar lugar a la siguiente ventana

27

Figura 39
En la solapa BASIC podemos hacer una descripcin del step, establecer el perodo de
tiempo en el cual se calcula, indicar si se desea trabajar o no con no linealidad
geomtrica y utilizar AUTOMATIC STABILIZATION en el caso de que haya problemas
de convergencia.
En la solapa INCREMENTATION (Figura 40) podemos establecer la cantidad mxima
de incrementos que deseamos que el programa utilice y la longitud de los mismos

28

Figura 40
Para el caso del presente ejemplo dejaremos los parmetros establecidos por defecto
por ABAQUS, por lo que hacemos click en OK.
7) Asignacin de cargas
Una vez definido el step podemos asignar al modelo el estado de cargas
correspondiente, para ello volvemos al men LOAD dentro de la opcin Module

Figura 41

29

En el men Load hacemos click en el cono

observa en la Figura 42

y aparecer la ventana que se

Asignar nombre

Figura 42
All asignamos un nombre a la carga y seleccionamos el STEP GENERAL STATIC
que acabamos de crear. En CATEGORY seleccionamos MECHANICAL y en TYPE
FOR SELECTED STEP seleccionamos CONCETRATE FORCE. Finalmente hacemos
click en CONTINUE
Ahora debemos seleccionar el punto en el que deseamos aplicar la carga, en este
caso el que se encuentra en la esquina superior derecha (ver Figura 43) y clickeamos
en DONE

Figura 43

30

Ahora aparecer la ventana que se observa en la Figura 44, en la cual en CF2

indicamos -1, esto significa que estamos aplicando una carga de 1N en la direccin 2
(en la terna cartesiana) con sentido negativo. Para finalizar hacemos click en OK.

-0.5

Figura 44
El aspecto de la carga aplicada se observa en la Figura 45

Figura 45
8) Contour integral- Definicin del frente de fisura y de su direccin de
propagacin
ABAQUS/Standard permite realizar la evaluacin de varios parmetros de Mecnica
de Fractura, como la integral J, la integral Ct, el Factor de Intensidad de Tensiones (K),

31

la direccin de propagacin de la fisura y la T-stress, lo cual da lugar al clculo de

diferentes tipos de CONTOUR INTEGRALS.
La evaluacin de varias contour integrals puede llevarse a cabo en cada ubicacin a lo
largo de la fisura. En un modelo de elementos finitos cada evaluacin puede pensarse
como el movimiento virtual de un bloque de material que rodea la punta de la fisura (en
dos dimensiones) o que rodea cada nodo a lo largo de la lnea de la fisura (en tres
dimensiones). Cada tipo de bloque se define por contornos, cada contorno es un anillo
de elementos que rodea completamente la punta de la fisura o un anillo de nodos a lo
largo de la lnea de fisura desde una cara de la fisura a la opuesta. Estos anillos de
elementos se definen de forma recursiva rodeando los contornos previos.
ABAQUS/Standard encuentra automticamente los elementos que forman cada anillo
de las regiones definidas como punta de fisura o lnea de fisura. Cada contorno provee
una evaluacin de la contour integral. El nmero de evaluaciones posibles es el
nmero de estos anillos de elementos, debemos especificar el nmero de contornos a
utilizar en el clculo de las contuor integrals.
El tipo de contour integral que utilizaremos en este ejemplo ser Factor de Intensidad
de Tensiones, para lo cual es necesario que definamos y asignemos la fisura, para ello
debemos seleccionar el men INTERACTON en el MODULE

Figura 46
Una vez all, en la barra de men que se encuentra en la parte superior de la pantalla
seleccionamos SPECIAL
CRACK
CREATE

Figura 47

32

Aparecer la ventana que se observa en la Figura 48, en la que debemos asignarle un

nombre a la fisura y en TYPE seleccionar la opcin CONTOUR INTEGRAL, luego
hacemos click en CONTINUE

Asignar nombre

Figura 48
Luego ABAQUS nos pide que definamos el frente de la fisura, el cual ser el primer
contorno de integracin. Para este problema seleccionaremos la punta de la fisura,
como se observa en la Figura 49 y luego hacemos click en DONE para teminar

Crack front

Figura 49
El siguiente paso consiste en definir la direccin de propagacin de la fisura, la cual
puede ser NORMAL TO CRACK PLANE o Q VECTORS. La primera opcin establece
que la direccin de propagacin de la fisura es normal al plano de la misma, y la
segunda nos permite elegir una direccin de propagacin distinta de la normal.
La opcin NORMAL TO CRACK PLANE debe utilizarse solamente si la direccin es la
misma en todos los puntos a lo largo de la lnea de la fisura, si esto no sucede no

33

podemos seleccionar una nica normal que defina la direccin de propagacin de la

fisura para todos los puntos a lo largo de la misma.
En el presente ejemplo utilizaremos la opcin NORMAL TO CRACK PLANE (Figura
50) seleccionndola en la parte inferior del rea de trabajo

Figura 50
Luego ABAQUS nos pide que definamos la coordenadas de dicho vector normal, para
lo cual ingresaremos en el recuadro que aparece en la parte inferior del rea de
trabajo (0.0,0.0), ENTER y luego (-1.0,0.0) ENTER

Figura 51
Aparecer la siguiente ventana, en la que debemos indicar que estamos trabajando
con simetra para que el programa calcule el factor de intensidad de tensiones para la
geometra completa. Finalmente hacemos click en OK para terminar

34

(1,0,0)

Figura 52
9) Output request
Una vez cargado completamente el modelo debemos indicarle al programa los
parmetros que deseamos que calcule, para ello debemos volver al men STEP y

, los cuales corresponden al HISTORY OUTPUT

trabajar con los conos
MANAGER. Al hacer click sobre el cono que se encuentra a la derecha aparecer la
siguiente ventana

Figura 53
Hacemos click en EDIT y aparecer la ventana que se observa en la Figura 54

35

Figura 54
All en DOMAIN seleccionamos CONOTUR INTEGRAL, en FRECUENCY la opcin
LAST INCREMENT, en NUMBER OF CONTOURS indicamos la cantidad contornos de
integracin que deseamos que el programa calcule, por ejemplo 10, y en TYPE
seleccionamos STRESS INTENSITY FACTORS, en CRACK INITIATION CRITERION
dejamos la opcin por defecto MAXIMUN TANGENTIAL STRESS. Finalmente
hacemos click en OK
En la ventana HISTORY OUTPUT REQUESTS MANAGER hacemos click en
DISMISS.
IMPORTANTE: La cantidad de contornos de integracin requeridos viene determinada
por el objetivo de alcanzar una solucin convergente, es decir, que el valor de K
calculado no vare de un contorno al otro, esto se logra a travs de un proceso de
prueba y error, y tiene una gran influencia el grosor del mallado utilizado. Se
recomienda comenzar con 10 contornos de integracin.

36

Una vez terminada la edicin del HISTORY OUTPUT pasamos a observar las

opciones por defecto que se encuentran en el FIELD OUTPUT

haciendo
click en el cono que se encuentra a la derecha. Aparecer una ventana similar a la
correspondiente al HISTORY OUTPUT, como se observa en la Figura 55

Figura 55
All hacemos click en EDIT y aparecer la ventana siguiente

37

Figura 56
En la Figura 56 se observan los parmetros que el programa puede calcular, los que
se encuentran marcados con verde son los seleccionados por defecto, los cuales
pueden aceptarse o deseleccionarse. En nuestro caso dejaremos los parmetros por
defecto, por lo que hacemos click en OK y luego en DISMISS.
10) Generacin del Job
El paso previo antes de que el programa comience con la resolucin del problema
consiste en darle un nombre al modelo que hemos cargado, para ello debemos
acceder al men JOB en el MODULE

38

Figura 57

Una vez que ingresamos al men JOB hacemos click en el cono

la siguiente ventana

y aparecer

Asignar nombre

Figura 58

En esta ventana le asignaremos un nombre al job, por ejemplo Viga-1,

seleccionamos el modelo que deseamos correr, en este caso Viga y hacemos click
en CONTINUE. Aparecer la ventana que se observa en la Figura 59, en la cual
tenemos distintas solapas con opciones referidas a la cantidad de memoria a utilizar,
etc, dejaremos las opciones por defecto y haremos click en OK.

39

Figura 59
Para que ABAQUS resuelva el modelo debemos ir al rbol que se encuentra en la
parte izquierda de la pantalla, all haciendo click en el signo + que se encuentra al lado
de JOB aparecer el job que acabamos de crear, hacemos click derecho sobre el
mismo y seleccionamos la opcin SUBMIT (Figura 60)

40

Figura 60
11) Visualizacin de resultados
Cuando ABAQUS completa un anlisis aparecer la palabra complete entre
parntesis al lado del job que solicitamos calcule. Para visualizar los resultados
obtenidos hacemos click derecho sobre el job, en este caso Viga-1 y seleccionamos
la opcin RESULTS (Figura 61)

Figura 61
ABAQUS nos enviar al postprocesador, el cual se denomina ABAQUS VIEWER, all
aparecer una imagen de la viga como se observa en la Figura 62

41

Figura 62
Para observar el campo tensional para la geometra deformada hacemos click en el
que se encuentra en la parte izquierda del rea de trabajo. Aparecen por
cono
defecto las tensiones de Von Mises, como se observa en la Figura 63

Figura 63
Para seleccionar otra direccin de tensiones u otra variable debemos ir en la barra de
men a RESULTS
FIELD OUTPUT y aparecer la ventana que se
observa en la Figura 64

42

Figura 64
All en OUTPUT VARIABLE elegimos la variable que deseamos observar, por ejemplo
TENSIN (S) y en INVARIANT y COMPONENT el tipo o direccin, por ejemplo
MISES. Haciendo click en OK ABAQUS la mostrar en pantalla.

NOTA: Las herramientas

para manipular la pieza modelada, sobre todo en modelo 3D.

son de utilidad

Para observar el valor del factor de intensidad de tensiones debemos acceder a uno
de los archivos que ABAQUS crea cuando analiza un caso. ABAQUS genera archivos
con distintas extensiones en los que se puede observar la forma en cmo opera el
programa, por ejemplo
archivos *.inp: son los archivos de entrada, como el que generamos en la
plataforma CAE
archivos *.msg: son archivos en los que se observa el proceso de clculo y en
los que se indican errores o warnings
archivos *.dat: son archivos en los que se verifican los datos de entrada,
indican errores o warnings vinculados a los mismos.

43

archivos *.sta: son archivos en los que se muestra el progreso del proceso
iterativo
archivos *.odb: son los archivos de resultados

Los valores del Factor de Intensidad de Tensiones se encuentran el al archivo *.dat, en

nuestro caso Viga-1.dat, para acceder al mismo debemos seguir la ruta C:\Temp\Viga1.dat. En la Figura 65 observamos una captura del archivo Viga-1.dat en donde se
indican los valores del factor de intensidad de tensiones calculado para los 10
contornos de integracin

Puede observarse que a partir del contorno 5 los valores comienzan a converger a
2.37 MPa mm1/2.
Para que el valor del Factor de Intensidad de Tensiones calculado con ABAQUS sea
representativo se promedian los valores del mismo obtenidos para los distintos
contornos, descartando el valor del primer contorno, ya que no es representativo.
Por lo tanto para este ejemplo el valor del Factor de Intensidad de Tensiones obtenido
con ABAQUS es 2.341 MPa mm1/2.
Finalmente compararemos este valor con el obtenido en los manuales; para ello
obtenemos, para la geometra y estado de carga analizado los valores de las
relaciones:
L/b= 4
a/b= 0.5
K0= 3LP/2b2 *(pa)1/2 =1.68 MPa mm1/2
Ingresando con estos datos a la Figura A.7 de la pgina 102 del libro Mecnica de la
Fractura Aplicada a Slidos Elsticos Bidimensionales obtenemos:
K/K0= 1.41
Por lo tanto K= 2.369 MPa mm1/2
La diferencia entre el valor obtenido numricamente y el obtenido analticamente es
del 1.18%, por lo tanto podemos concluir que utilizando un programa de elementos

44

finitos para calcular el Factor de Intensidad de Tensiones obtenemos una muy buena
aproximacin.
Referencias bibliogrficas

Eugenio Oate. Clculo de Estructuras por el Mtodo de Elementos Finitos

Anlisis esttico lineal. 1992

ABAQUS. Users Manual Versin 6.9. Dessaults Systmes, 2009

ABAQUS/CAE. Users Manual Versin 6.9. Dessaults Systmes, 2009

Manuel Elices Calafat. Mecnica de la Fractura Aplicada a Slidos Elsticos

Bidimensionales.1998

Mayllet Y. Uzctegui Flores. ABAQUS Programa de Elementos Finitos.

45