IMPORTANCIA DEL ABLANDAMIENTO TÉRMICO EN SIMULACIONES DE IMPACTOS EN PLACAS DE ACERO A ALTAS TEMPERATURAS
B. Erice, F. Gálvez, M. J. Pérez, D. A. Cendón y V. Sánchez-Gálvez
Departamento de Ciencia de Materiales, CISDEM, Universidad Politécnica de Madrid,
Profesor Aranguren s/n, 28040, Madrid E-mail: [email protected]
RESUMEN
En el presente trabajo se presenta un estudio experimental y numérico de impactos balísticos sobre placas de acero inoxidable martensítico a altas temperaturas (400ºC y 700ºC), que pone de manifiesto la importancia del ablandamiento térmico en simulaciones de impactos a altas temperaturas. Mediante un estudio metalográfico de la zona de impacto, se ha observado la aparición de bandas adiabáticas de cortante formadas por el aumento brusco de la temperatura debido a la acumulación del trabajo plástico en el interior del material. La correcta predicción en la formación de estas bandas durante el proceso de penetración es crítica a la hora de obtener resultados representativos de los experimentos realizados. Basándose en datos experimentales de ensayos previamente realizados, se ha calibrado un modelo de material de Johnson-Cook (JC) para su uso con simulaciones numéricas en el código no lineal de elementos finitos LS-DYNA. Mediante estas simulaciones numéricas se demuestra la importancia del ablandamiento térmico en el proceso de perforación de placas, al igual que la incapacidad que un modelo tipo JC tiene para representar el dicho ablandamiento para material estudiado. Esta investigación presenta, finalmente, una modificación a un modelo tipo JC programado como subrutina de material de usuario para LS-DYNA que permite simular correctamente estos procesos de impacto a alta temperatura.
ABSTRACT
An experimental and numerical study of plate impacts at high temperatures (400ºC and 700ºC) is presented in this article. Thermal softening behaviour of a material has been found crucial when simulating dynamic phenomena like plate impacts. A metallographic study on some of the impacted specimens has shown adiabatic shear band formation. The shear bands formation is caused by the adiabatic temperature rise due to the plastic work accumulated inside the material. The correct prediction of the shear band formation in the penetration process has been found critical to obtain results representative from experimental data. Based on previous experiments, a Johnson-Cook (JC) material model has been calibrated. Non-linear finite element code LS-DYNA has been used to carry out the numerical simulations of plate impacts. These simulations show the importance of the thermal softening in plate penetration problems. JC type material model has been proven ineffective to reproduce experimental data. Therefore, a modification of the model is proposed. In order to modify the JC a user material subroutine has been programmed in LS-DYNA. PALABRAS CLAVE: Simulaciones Numéricas; Bandas adiabáticas de cortante; Impacto.
1. INTRODUCCIÓN
La motivación de esta investigación, se basa en los eventos de “blade-off”, que estudian la contención ante un desprendimiento de un alabe en el interior de un de la turbina de un motor de aviación. El material objeto de este estudio es el acero inoxidable martensítico FV535, uno de los materiales más usados para carcasas de turbinas de baja presión de motores de aviación. En estas, las temperaturas de utilización van desde unos 400ºC, hasta unos 800ºC.
El desprendimiento y posterior impacto de un alabe sobre una de estas carcasas, si es que sucediese, es un fenómeno dentro del ámbito de la física o mecánica del impacto. Por lo tanto requiere una caracterización mecánica específica a alta velocidad de deformación y alta temperatura, así como un estudio balístico a altas temperaturas. Una correcta caracterización mecánica del material estudiado, puede permitir simular
numéricamente fenómenos balísticos y reproducir con precisión los datos experimentales.
Una extensa campaña de ensayos mecánicos sobre el acero FV535 [1, 2] ha permitido la calibración de un modelo de material de Johnson-Cook (JC) [3, 4] (ver Tabla 1). Con el objetivo de validar mediante simulaciones numéricas este modelo de material, se han llevado a cabo una serie de impactos balísticos sobre placas de acero FV535 [5]. Dadas las temperaturas de trabajo de las carcasas de turbina de baja presión, se ha optado por realizar esto ensayos balísticos a dos temperaturas, 400ºC y 700ºC. Esto, ha requerido el diseño y construcción de un horno balístico [5] que permite calentar y mantener a la temperatura deseada las placas durante el ensayo, con una temperatura máxima de aproximadamente 900ºC.
La aparición de bandas adiabáticas de cortante en impactos balísticos sobre placas, condiciona el resultado
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
517
0
200
400
600
800
1000
1200
200 400 600 800 1000 1200 1400 1600 1800
Lím
ite e
lást
ico
(MPa
)
Temperatura (ºC)
m = 4.5 ; 10-4
m = 4.5 ; 10+3
Tm
= 870ºC Tmext
= 1500ºC
Figura 1. Evolución del límite elástico con la temperatura. En líneas continuas, ajuste de los datos experimentales con JC. Línea discontinua, el modelo propuesto.
de las simulaciones numéricas [6]. La correcta predicción de la formación de estas bandas, es crítica a la hora de predecir los resultados experimentales. Se ha comprobado que sobre el material objeto de estudio, aparecen bandas de cortante, tal y como se ve en la Figura 2.
Un modelo tipo JC para este material en concreto, no es capaz de reproducir los resultados experimentales, tal y como demostró el presente autor mediante un estudio numérico [7] que simulaba impactos balísticos sobre placas de FV535. La conclusión de dicho estudio, por razones que se explicarán más adelante, es que este tipo de modelo no es capaz de reproducir el ablandamiento térmico que realmente sufre el este material. Para solventar dicho problema, se propone una modificación del modelo de JC que es capaz de reproducir resultados experimentales. Para poder realizar simulaciones con el modelo que se propone como válido, ha sido necesaria la programación de una subrutina de material de usuario en el código de elementos finitos no lineal LS-DYNA [8].
Mediante el uso de curvas de velocidad residual frente a velocidad inicial de impactos sobre placas a altas temperaturas (400ºC y 700ºC), se han comparado ambos modelos, el JC original y el modificado. Tabla 1. Constantes del modelo de JC, así como para el modelo propuesto, para el acero inoxidable FV535.
E (GPa) (kg/m3) pC (J/kgK) 210 0.28 7850 460 0.9
A (MPa) B (MPa) n C m1035 190 0.3 0.016 4.5
mT (ºC) extT (ºC) mextT (ºC) 0 (s-1) ARc 870 845 1500 5x10-4 20
1D 2D 3D 4D 5D 0.1133 2.11 -1.65 0.0125 0.9768
2. MODELO DE MATERIAL
El modelo de Jonhson-Cook (JC) se compone de una relación constitutiva (6) y un criterio de rotura (11) desacoplados entre si, es decir, que no existe ablandamiento del material debido a la acumulación de daño. Se propone la modificación de la relación constitutiva de JC acorde con la expresión (8), para poder representar correctamente el comportamiento térmico del acero FV535. Todas las constantes necesarias para ambos modelos se resumen el en la Tabla 1.
2.1. Implementación de los modelos
El modelo se ha implementado con un criterio de plastificación clásico de von Mises, donde la función de plastificación viene dada por:
, , 0JC p peq eqT (1)
siendo, JCeq la relación constitutiva de JC (ley de
endurecimiento), p la deformación plástica equivalente o de von Mises, p la velocidad de la deformación plástica equivalente, T la temperatura y
eq la tensión equivalente o tensión de von Mises. Con el objetivo de poder tener una medida escalar de una magnitud tensorial se define la tensión equivalente,
23 3 / 2 :eq J s s (2)
donde s es el tensor desviador de tensiones y 2J es el segundo invariante de dicho tensor, y la deformación plástica equivalente,
2 / 3 :p p p (3)
donde p es el tensor de velocidades de deformación plástica.
La regla de flujo adoptada es asociativa y por tanto, se ha tomado la función de plastificación de von Mises como potencial de flujo.
p (4)
donde es la velocidad del multiplicador de flujo plástico, y es el tensor de tensiones. Introduciendo la expresión (3) en la (4) obtenemos la ley de Prandtl-Reuss:
23
p p
eq
s (5)
El algoritmo de retorno a la superficie de plastificación ha sido resuelto mediante el método de escalado radial [9].
2.2. Relación constitutiva de JC
La relación constitutiva de JC [4] consta de tres términos que se multiplican entre si. El primero de ellos
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
518
Figura 2. Imauna banda atransversal al
es una leydeformación endurecimientercero y úmaterial. Tie
JCeq A B
donde A , B*
0p p
de referen* /rT T T
donde rT etemperatura d
La obtenajustando loobservar en como la temexperimentalSabemos quede fusión pertipo de ajuste
2.3. JC con t
Como stemperatura exponer con aunque el maa temperaturaresistencia hconcepto es línea que estno hay evideclaro que es1500ºC. Entresistencia. extendida co
*ext mextT T T
donde mextT ey extT es la tede la cual se
agen de microsadiabática deimpacto de una
y de endureplástica equi
nto por veloúltimo es el ene por expres
1np Cln
B , n , C y mdonde 0 es u
ncia definid/ m rT T eses la tempede fusión. nción de las cos datos expla Figura 1. P
mperatura de fles, tiene un e esta temperaro el modelo e.
temperatura d
se ha expuestde fusión resta nueva r
aterial haya pas cercanas a
hasta su tempeel que se pue
tá trazada conencias experimste material ntre estas dos Así definimomo:
mmext extT T T
es la temperatemperatura exformula la nu
scopio electróne cortante soba de las placas
ecimiento poivalente p . Locidad de de
ablandamiensión:
* *1mp T
m son constanuna velocidadda por e la tempera
eratura inicia
constantes m perimentales Podemos obsfusión, extrapvalor aproxim
atura no es la de JC no per
de fusión exten
to mT en estreal. La ideaelación consti
perdido casi tolos 870ºC, tod
eratura de reaede ver en la
n línea discontmentales de lo no funde a lotemperaturas
os la temper
tura de fusióntendida (ver F
ueva relación c
nico de barridobre la superfanalizadas.
otencial conLa segunda eeformación ynto térmico
(
ntes del materd de deformacel usuario atura homóloal y mT es
y mT [2] se hque se pueervar clarameolando los damado de 870temperatura r
rmite ningún o
ndida (JCT)
te caso no esa que se quiitutiva, es la
oda su resistendavía retiene e
al de fusión. Ea Figura 1 continua. Nótese expuesto pero
os 870ºC sinmantiene cie
ratura homol
(
n extendida o rFigura 3) a paconstitutiva:
o de ficie
la s el
y el del
(6)
rial, ción
y oga,
la
hizo eden ente atos 0ºC. real otro
s la iere que ncia esta Este n la que o es
no a erta loga
(7)
real artir
Figtérm
Jeq
modotsimperNódelimpcaldef
T
donden
2.4
de rotparse pláparausdelsi lmopar
D
donrotconprihuey tdefsig
gura 3. Repremico que sigue
*
JCeqJCT
q JCeq extT
El cálculo deodelos JC y JCtar al materia
mulaciones nrmite dar a lótese que estal modelo MAplementado elculado sumafinidos como:
JC peq
pC
nde es el cnsidad del mate
4. Criterio de
No se han erotura neces
tura de JC [3ra los dos modbasa en la a
ástica en un prámetro, tomsencia de dañol material. Enlos elementosodelo de elemra JC se defin
f
p
nde f es la tura cuya exnstitutiva, forimero, es la coecos con la extercero muestrformación y
guiente expres
sentación esquel modelo JCT
0 ext
ext mext
T T TT T T
e la temperatuCT. Se ha incal de una tempnuméricas, elos elementoscaracterística
AT_015 [8], Jen LS-DYNAando los inc
coeficiente de rial y pC es el c
rotura de JC
encontrado evite modificac], se ha consdelos (JC y JCacumulación parámetro o ina como valoo, y como máx
n el caso de las alcanzan el vmentos finitose como:
deformación xpresión es rmada por treonocida relacixponencial de ran la dependtemperatura ión:
uemática del T.
ura es idénticcorporado la p
mperatura iniciesta temperas la temperaa no está dispJC estándar, qA. La tempe
ncrementos d
Taylor-Quinncalor específico
videncias de qción alguna. Eservado exactCT). Este critede daño por ndicador de dor inicial 0, ximo 1, indicaas simulacionvalor 1D ses. La evoluc
plástica eqhomologa a
es términos sión [10] del c
e la triaxialidadencia con la respectivamen
ablandamiento
(8)
o para ambosposibilidad deial 0T . En lasatura inicialtura deseada.
ponible dentroque ya vieneneratura se hade la misma
(9)
ney, es lao.
que el criterioEl criterio detamente igualerio de rotura,
deformacióndaño D . Esteindicando la
ando la roturaes numéricas,
e eliminan delión del daño
(10)
quivalente dea la relaciónseparados. Elrecimiento ded. El segundovelocidad de
nte. Tiene la
o
s e s l .
o n a a
a
o e l , n e a a , l o
)
e n l e o e a
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
519
2t
Figura 4. Conpara la placa.velocidad inic
2t
Figura 5. Conextendida mexTes de 700rT
0 600v m s .
1 2f D D
donde 1D , material y
1 / 3H tr
de tensiones.
3. MODEL
Todas lasel código exelementos cuadradas deproyectil esfmm de diámmodelo ha sen 2D de cucontrol de ho
Dado el llegar hasta ase ha intentaelemento ad
2.0 s
ntornos de temp La temperaturial del proyecti
2.0 s
ntornos de tem1500ºxt C . La
0ºC y los eleme
*3 1exp D
2D , 3D , 4D* /H eq
es la comp.
LO DE ELEM
s simulacionexplicito no linfinitos simu
e 101mm de lférico de acermetro. Dada ido realizado
uatro nodos yourglassing.
espesor de laproximadameado dentro de ecuado. Se m
t
peratura de las sra inicial es de il es de 0 600v
t
mperatura de laa temperatura eentos se han in
*4 1pD ln
4 y 5D son es la tria
ponente hidros
MENTOS FI
es numéricas sneal LS-DYNAula impactoslado y 2 mm ro de alta resla simetría con elemento
y un punto de
las bandas deente los 25 mlo posible ten
malló la zona
3.4t s
simulaciones nu700ºrT C y l
0m s .
3.4t s
as simulacioneextendida es eT
nicializado con
*5D T (
n constantes axialidad, dostática del ten
NITOS
se realizaronA. El modelos sobre plade espesor desistencia de 5del problemaos axilsimétri
e integración
e cortante pum (ver Figuraner un tamaño
de impacto
uméricas realizlos elementos s
es numéricas re845ºextT C . Co
una temperatu
(11)
del onde nsor
con o de acas e un 5.55 a el icos con
uede a 2), o de con
unaquePardoelede extext
usasolFVmo(MEl mé2D
expaltaelequeprolos
6.4t s
zadas con el mose han inicializa
6.4t s
ealizadas con ontornos solo v
ura de 0 700T
a malla deedándonos asra agilizar el s transiciones
emento según la placa. La p
tremo impidieterior de la miApenas se de
ados en los enlido rígido (M
V535, como sodelizó con
MAT_USER_Dcontacto ent
étodo automátD [8] sin fricci
Uno de los pplícitos de feas deformaci
ementos. En le es calculadoporcional a ls elementos,
0v
odelo estándar dado con una tem
0v
el modelo de Jválidos para la p0ºC . La velocid
e 64 elemesí unos elemecálculo en todde malla, aumnos acercam
placa se simulóendo el movisma. etectó deformnsayos, por loMAT_RIGID se ha explicadun material
DEFINED_MAtre ambas patico de formuón. problemas mánómenos de iones y distlos códigos edo automáticala longitud la
lo que pro
600m s
de JC. Contornmperatura de T
600m s
JC de temperaplaca. La temp
dad inicial del
entos sobre entos de 30 do lo posible mentando así
mos hacía el bó con empotra
vimiento de l
mación sobre loo tanto, se m
[8]). El blando con anterdefinido po
ATERIAL_Martes se simuulación por pe
ás comunes enimpacto balí
torsiones queexplícitos tamamente por ea longitud carovoca en est
nos solo válidos0 700ºT C . La
atura de fusiónperatura inicialproyectil es de
el espesor,m x 30 m.se realizaron
el tamaño delborde exterioramiento en suos nodos del
os proyectilesodelizó comonco de aceroioridad se ha
or el usuarioMODEL [8]).uló usando elenalización en
n los cálculosstico, son lase sufren los
maño de paso,el código, esracterística detos casos la
s a
n l e
, . n l r u l
s o o a o . l n
s s s , s e a
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
520
0
100
200
300
400
500
600
700
200 300 400 500 600 700 800
Experimental
JCT
JC
Velo
cida
d re
sidu
al (m
/s)
Velocidad inicial (m/s)
T = 400ºC
a
0
100
200
300
400
500
600
700
200 300 400 500 600 700 800
Experimental
JCT
JC
Velo
cida
d re
sidu
al (m
/s)
Velocidad inicial (m/s)
T = 700ºC
b
Figura 6. Curvas de velocidad residual frente a velocidad inicial de los datos experimentales y numéricos de impactos de placas a 400ºC (a) y a 700ºC (b). Se puede apreciar que el modelo JCT reproduce bien los resultados experimentales especialmente a 700ºC.
reducción del tamaño de paso hasta valores extremadamente pequeños ( 12~ 10 )La distorsión de los elementos puede ser tal que en muchos casos, dan volúmenes negativos, con lo que el cálculo se para automáticamente teniendo que eliminar dichos elementos manualmente y reinicializando el cálculo.
Para solventar esta clase de problemas numéricos, se implementó un criterio de erosión (eliminación elementos del modelo de elementos finitos) de elementos según que depende de su aspecto de forma. La subrutina del material de usuario programada permite obtener las coordenadas de los nodos de los elementos en todo instante. Definimos aspecto de forma AR como el cociente entre el lado más largo del elemento y el lado más corto. Cuando este aspecto de forma llega a un valor crítico ARc, el elemento se erosiona. Esto permitió disminuir el tiempo de cálculo eliminando elementos problemáticos, que son intrascendentes para el resultado final de las simulaciones numéricas.
4. RESULTADOS Y DISCUSIÓN
Los resultados analizados son en base a las curvas de velocidad residual rv de la masa desprendida por el proyectil del blanco o plug frente a velocidad inicial 0v del proyectil. Para la predicción de las velocidades residuales, se ha usado un modelo analítico modificado de Recht e Ipson (MRI) [11] que tiene por expresión:
22 2
1
1 0
1
kk kbl
r
k v vv
m M (12)
donde 1k y 2k son constantes de ajuste obtenidas por mínimos cuadrados, m es la masa del plug y M es la masa del proyectil. Las constantes para expresión analítica se recogen en la Tabla 2.
Se realizaron simulaciones numéricas con los dos modelos expuestos (JC y JCT) para velocidades iniciales del proyectil parecidas a las lecturas experimentales para las dos temperaturas 400ºC y 700ºC. Los resultados obtenidos pueden verse en la Figura 6.
Para las simulaciones numéricas realizadas a 400ºC (ver Figura 6a) no hay mucha diferencia entre los dos modelo, la predicción de los resultados experimentales es bastante buena para ambos, siendo la del modelo JCT mejor que la del JC. Esto es debido a que el aumento de temperatura del material producido por el impacto del proyectil no es suficientemente grande como para hacer notar mucho la diferencia entre ambos modelos, pero suficiente como para ver que el modelo con la temperatura de fusión modificada (JCT) funciona mejor.
En el caso de las simulaciones numéricas realizadas a 700ºC (ver Figura 6b) las predicciones de las velocidades residuales para el modelo de JC no son buenas, especialmente cuando disminuimos la velocidad inicial del proyectil. En cambio el modelo JCT da unas predicciones excepcionales. La elevada temperatura inicial (700ºC) hace que en el impacto, rápidamente se alcance la temperatura mT con lo que el modelo de JC ya no es capaz de reproducir los resultados experimentales. En cambio el modelo propuesto (JCT) reproduce las velocidades residuales con fidelidad.
En la Figura 4 y la Figura 5 se recogen tres secuencias de las simulaciones numéricas realizadas con los dos modelos JC (Figura 4) y JCT (Figura 5) para una temperatura de 700ºC. Las secuencias se tomaron a para los mismos tiempos, tomando como origen de tiempos el contacto entre proyectil y blanco. Con esta serie de secuencias se quiere recalcar que no solo es importante el predecir las velocidades residuales experimentales. La morfología de la rotura producida por el proyectil en el blanco es fundamental, ya que estamos tratando de reproducir el comportamiento del material y su fractura
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
521
a altas tempecomo las temtemperatura se formen prresiduales qulas experimeun aspecto mdúctil, comoAdemás de eva penetrandaumentan dellegando fina
a
Figura 7. Pluglos modelos den un experim
En la Figcortante se correcta formJCT evita qualcanzan los toda su resistTabla 2. Resimpactos sobr
Experimental
Numérico JC
JCT JC
JCT
5. CONCL
Se ha deválido para modificacióntener cierta rhasta su temmodelo que a
Las simu700ºC han prcapacidad dexperimentallas bandas d
eraturas. Vemmperaturas qu
mT . Esto caurematurament
ue predice estentales. El plug
más parecido o lo es el aesto los elemedo el blanco, e temperaturaalmente a eros
gs obtenidos ende JC (a) y JCT
ento. 0 700ºT
gura 5 puede forman de
mación del plue a pesar de
elementos entencia y se erosultados exper
re placas de ace
l 1k 0.75 0.80
1k 0.76 0.78 0.75 0.76
LUSIONES
emostrado queste material.
n del mismo bresistencia, po
mperatura de fal que se le haulaciones numrobado la valide predicciónles. Se ha vistde cortante po
mos con claridaue alcanzan susa que las bate con lo quee modelo seang formado en al de un mat
acero FV535 entos que ve e
que son los a, pierden todsionarse antes
b
n las simulacioT (b) comparadºC , 700ºrT C
apreciar comtal manera qlug (ver Figurlas elevadas t
n el impacto, osionen. rimentales [5] ero FV535.
2k 1.90 41.90 3
2k 2.50 42.10 42.10 22.20 3
ue un modelo. Se propone basado en queor muy pequefusión real. Aa denominado méricas realizdez de este nun de velocidto que la correor este model
ad en la Figuron cercanas a
andas de cortae las velocidan muy superioel impacto, ti
terial frágil y(ver Figura
el proyectil segque más ráp
da su resistende tiempo.
nes numéricas dos con el recogC y 0 420v m
mo las bandasque permitenra 7). El modtemperaturas estos pierdan
y numéricos
blv T (º440.0 40385.0 70
blv T (420.0 40440.0 40200.0 70380.0 70
o tipo JC noen su lugar
e el material deña que esta s
Así se formulacomo JCT.
zadas a 400ºCuevo modelo ydades residuaecta formaciónlo lleva a que
ra 4 a la ante ades ores iene y no
7). gún
pido ncia
con
gido s .
s de n la delo que
n su
de
ºC) 00 00 (ºC) 00 00 00 00
o es una
debe sea, a el
C y y su ales n de e la
rotun
expEstdetmá
AG
Cietradel
RE[1] Flotria[2] Gámede [3] thretem21 [4] DaandBal[5] Meat Bal[6] FinStr(20[7] andbaltem[8] LivCal[9] Phe199[10meto Phy[11Dy
tura del blancoaspecto parecSe propone
perimental deto requiere reterminar los vás allá de los 8
GRADECIMIE
Los autoresencia e Innovés del progrl proyecto BIA
EFERENCIASErice, B., Gálv
ow and fracturaxialities, strain
Gálvez, F., Clvez, V., High s
etallic materialsPhysique IV, 1Johnson, G.R.
ee metals sumperatures and
(1985) 31-48.Johnson, G.R.
ta for Metals Sd High Temperllistics, The HaErice, B., Gálv
echanical BehavHigh Temperallistics, Beijing
Kane, A., Bønite Element Auck by Blunt P
009) 051302. Erice, B., Gá
d Børvik, T., llistic impacts mperatures, Jour
LS-DYNA KEvermore Softwlifornia, 2007.
Wilkins, Menomena, Sprin99. 0] Hancock, Jechanisms of du
multi-axial strysics of Solids, 1] Recht, R.F.ynamics, Journa
o así como la cido al real.
como trabae la hipótesis ealizar una cavalores del lí850ºC.
ENTOS
s desean agrvación la finrama CONSOA2008-06705-
vez, F., Cendónre behaviour o
n rates and tempCendón, D., Estrain rate and hs for jet engine 34 (2006) 269-, and Cook, W.bjected to vapressures, Eng
and Cook, W.Hubjected to Lar
ratures, in: 7thgue, 1983, pp.
vez, F., Cendónvior of FV535 Satures, in: 25th , China, 2010, p
ørvik, T., HoppAnalysis of PlugProjectiles, Jour
lvez, F., CendóAn experimenon a turbine
rnal of Applied EYWORD USE
ware Technolog
M.L., Computenger-Verlag, B
J.W., and Mauctile failure in ress-states, Jou24 (1976) 147-, and Ipson,
al of Applied M
a morfología d
ajo futuro, lque sigue el ampaña expeímite elástico
radecer al Mnanciación de
OLIDER INGE-C02-01
n, D., and Sánchof FV535 ste
peratures. , enviEnfedaque, A., high temperatu
e turbine contai-274. .H., Fracture charious strains,gineering Fractu
H., A Constitutrge Strains, Hig
h International 541-547.
n, D., and SánchSteel Against B
h Internatonal Spp. 1087-1096.perstad, O., &
ugging Failure rnal of Applied
ón, D., Sánchntal and nume
casing materd Mechanics, (2ER'S MANUALgy Corporatio
er Simulation Berlin Heidelbe
ackenzie, A. high strength s
urnal of the M-169. T. W., Ballis
Mechanics, 30 (1
del plug tenga
la validaciónmodelo JCT.rimental para
o del material
Ministerio deel proyecto aENIO 2010 y
hez-Gálvez, V.,el at differentiado, (2010).
and Sánchez-re behaviour ofinment, Journal
haracteristics of, strain rates,ure Mechanics,
tive Model andgh Strain RatesSymposium on
hez-Gálvez, V.,Ballistic ImpactSymposium on Langseth, M.,
in Steel PlatesMechanics, 76
hez-Gálvez, V.,erical study ofrial at varying010). L Version 971,
on, Livermore,
of Dynamicerg New York,
C. , On thesteels subjectedMechanics and
tic Perforation963) 384-390.
a
n . a l
e a y
, t
-f l
f , ,
d s n
, t n
, s 6
, f g
, ,
c ,
e d d
n
Anales de Mecánica de la Fractura 28, Vol. 2 (2011)
522