En mecánica de fracturas , la tasa de liberación de energía ,, es la velocidad a la que la energía se transforma cuando un material se fractura . Matemáticamente, la tasa de liberación de energía se expresa como la disminución de la energía potencial total por aumento en el área de la superficie de fractura, [1] [2] y, por lo tanto, se expresa en términos de energía por unidad de área. Se pueden construir varios balances de energía relacionando la energía liberada durante la fractura con la energía de la nueva superficie resultante, así como con otros procesos disipativos como la plasticidad.y generación de calor. La tasa de liberación de energía es fundamental para el campo de la mecánica de la fractura al resolver problemas y estimar las propiedades del material relacionadas con la fractura y la fatiga .
Definición
La tasa de liberación de energía se define [3] como la pérdida instantánea de energía potencial total por unidad de área de crecimiento de grietas ,
donde la energía potencial total se escribe en términos de la energía de deformación total , tracción superficial , desplazamiento y fuerza corporal por
La primera integral está sobre la superficie del material, y el segundo sobre su volumen .
La figura de la derecha muestra la gráfica de una fuerza externa. vs el desplazamiento del punto de carga , en el que el área bajo la curva es la energía de deformación. El área blanca entre la curva y el-Eje se conoce como la energía complementaria. En el caso de un material linealmente elástico , es una línea recta y la energía de deformación es igual a la energía complementaria.
Desplazamiento prescrito
En el caso de un desplazamiento prescrito, la energía de deformación se puede expresar en términos del desplazamiento especificado y la superficie de la grieta. , y el cambio en esta energía de deformación solo se ve afectado por el cambio en el área de la superficie de fractura: . En consecuencia, la tasa de liberación de energía en este caso se expresa como [3]
Aquí es donde uno puede referirse con precisión como la tasa de liberación de energía de deformación.
Cargas prescritas
Cuando se prescribe la carga en lugar del desplazamiento, la energía de deformación debe modificarse como . La tasa de liberación de energía se calcula como [3]
Si el material es linealmente elástico, entonces y uno puede en su lugar escribir
G en casos bidimensionales
En los casos de problemas bidimensionales, el cambio en el área de crecimiento de la grieta es simplemente el cambio en la longitud de la grieta multiplicado por el espesor de la muestra. A saber,. Por lo tanto, la ecuación para calcular se puede modificar para el caso 2D:
- Desplazamiento prescrito:
- Carga prescrita:
- Carga prescrita, elástico lineal:
Se puede consultar los cálculos de ejemplo incluidos en la siguiente sección para obtener más información. A veces, la energía de deformación se escribe usando, un espesor de energía por unidad. Esto da
- Desplazamiento prescrito:
- Carga prescrita:
- Carga prescrita, elástico lineal:
Relación con los factores de intensidad del estrés
La tasa de liberación de energía está directamente relacionada con el factor de intensidad de la tensión asociado con un modo de carga bidimensional dado ( Modo I, Modo II o Modo III ) cuando la grieta crece en línea recta. [3] Esto es aplicable a grietas bajo tensión plana , deformación plana , y de cizallamiento antiplano .
Para el Modo I, la tasa de liberación de energía la tasa está relacionada con el factor de intensidad del estrés Modo-I para un material linealmente elástico por
dónde está relacionado con el módulo de Young y la relación de Poisson dependiendo de si el material se encuentra bajo tensión plana o tensión plana:
Para el Modo II, la tasa de liberación de energía se escribe de manera similar como
Para el Modo III (cizallamiento antiplano), la tasa de liberación de energía ahora es una función del módulo de cizallamiento ,
Para una combinación arbitraria de todos los modos de carga, estas soluciones elásticas lineales pueden superponerse como
Representación alternativa de en términos de [4] |
---|
Que puede verse como equivalente a la representación anterior a través de la relación entre el módulo de Young y el módulo de corte: |
Relación con la tenacidad a la fractura
El crecimiento de grietas se inicia cuando la tasa de liberación de energía supera un valor crítico , que es una propiedad material,
Bajo carga Modo-I, la tasa de liberación de energía crítica luego se relaciona con la tenacidad a la fractura Mode-I , otra propiedad material, por
Calculando G
Hay una variedad de métodos disponibles para calcular la tasa de liberación de energía dadas las propiedades del material, la geometría de la muestra y las condiciones de carga. Algunas dependen de que se cumplan ciertos criterios, como que el material sea completamente elástico o incluso linealmente elástico, y / o que la grieta deba crecer en línea recta. El único método presentado que funciona arbitrariamente es el que utiliza la energía potencial total. Si ambos métodos son aplicables, deberían producir tasas de liberación de energía idénticas.
Energía potencial total
El único método para calcular para condiciones arbitrarias es calcular la energía potencial total y diferenciarla con respecto al área de la superficie de la fisura. Por lo general, esto se hace mediante:
- calcular el campo de tensión resultante de la carga,
- calcular la energía de deformación en el material resultante del campo de tensión,
- calcular el trabajo realizado por las cargas externas,
todo en términos de la superficie de la grieta.
Ejemplo de cálculo de utilizando la energía potencial total |
---|
Este problema es bidimensional y tiene una carga fija, por lo que con , Dado que el material es linealmente elástico, y por lo tanto Las tensiones en el DCB se deben a las tensiones de flexión en cada viga en voladizo, dónde es la longitud en la página. Utilizando, uno ahora tiene la energía de tensión donde el factor de 2 en el frente se debe a que hay 2 vigas en voladizo. Resolviendo, luego tomando una derivada con respecto a y dividiendo por , |
Método de cumplimiento
Si el material es linealmente elástico, el cálculo de su tasa de liberación de energía puede simplificarse mucho. En este caso, la curva de carga vs.desplazamiento del punto de carga es lineal con una pendiente positiva, y el desplazamiento por unidad de fuerza aplicada se define como el cumplimiento, [3]
La energía de deformación correspondiente (área bajo la curva) es igual a [3]
Usando el método de cumplimiento, se puede demostrar que la tasa de liberación de energía para ambos casos de carga prescrita y desplazamiento resulta ser [3]
Ejemplo de cálculo de utilizando el método de cumplimiento [3] |
---|
Considere una muestra de viga en voladizo doble (DCB) como se muestra en la figura de la derecha. El desplazamiento de una sola viga es El desplazamiento del punto de carga resultante es, por lo tanto . Sustituya en la ecuación el cumplimiento y simplifique: Ahora, se calcula como Finalmente, la tasa de liberación de energía de esta muestra de DCB se puede expresar como Tenga en cuenta que, alternativamente, la tasa de liberación de energía se puede expresar en términos de y : Indicando que disminuye con la longitud de la grieta para el caso de desplazamiento fijo, y viceversa para el caso de carga fija. |
Múltiples métodos de muestra para materiales no lineales
En el caso de un desplazamiento prescrito, manteniendo fija la longitud de la grieta, la tasa de liberación de energía se puede calcular mediante [3]
mientras que en el caso de carga prescrita, [3]
Como se puede ver que, en ambos casos, la tasa de liberación de energía veces el cambio de superficie devuelve el área entre curvas, que indica la energía disipada para la nueva área de superficie como se ilustra en la figura de la derecha [3]
Cierre de grietas integral
Dado que la tasa de liberación de energía se define como la derivada negativa de la energía potencial total con respecto al crecimiento de la superficie de la fisura, la tasa de liberación de energía se puede escribir como la diferencia entre la energía potencial antes y después de que la fisura crezca. Después de una derivación cuidadosa, esto nos lleva a la integral de cierre de grietas [3]
dónde es la nueva superficie de fractura, son los componentes de la tracción que se liberan en la superficie superior de la fractura a medida que crece la grieta, son los componentes del desplazamiento de la abertura de la grieta (la diferencia en los incrementos de desplazamiento entre las superficies superior e inferior de la grieta), y la integral está sobre la superficie del material .
La integral de cierre de grietas es válida solo para materiales elásticos, pero sigue siendo válida para grietas que crecen en cualquier dirección. Sin embargo, para una grieta bidimensional que de hecho crece en línea recta, la integral de cierre de la grieta se simplifica a [3]
dónde es la nueva longitud de la grieta, y los componentes de desplazamiento se escriben en función de las coordenadas polares y .
Ejemplo de cálculo de utilizando la integral de cierre de fisuras |
---|
Considere la grieta en la muestra de DCB que se muestra en la figura. Los componentes de esfuerzo y desplazamiento distintos de cero vienen dados por [3] como La integral de cierre de la grieta para este material linealmente elástico, asumiendo que la grieta crece en línea recta, es Considere reescalar la integral usando por donde se calcula que la integral más simple es dejando la tasa de liberación de energía como la relación esperada. En este caso, no es sencillo obtenerdirectamente de la carga y la geometría del problema, pero dado que la grieta crece en línea recta y el material es linealmente elástico, la tasa de liberación de energía aquí debería ser la misma que la tasa de liberación de energía calculada usando los otros métodos. Esto permite recuperar indirectamente el factor de intensidad del estrés para este problema como |
J - integral
En determinadas situaciones, la tasa de liberación de energía se puede calcular utilizando la integral J , es decir, usando [3]
dónde es la densidad de energía de deformación elástica, es el componente del vector unitario normal a , la curva utilizada para la integral de línea, son los componentes del vector de tracción , dónde es el tensor de tensión, y son los componentes del vector de desplazamiento.
Esta integral es cero sobre una ruta cerrada simple y es independiente de la ruta , lo que permite utilizar cualquier ruta simple que comience y termine en las caras de la grieta para calcular. Para equiparar la tasa de liberación de energía a la integral J,, deben cumplirse las siguientes condiciones:
- la grieta debe estar creciendo en línea recta, y
- la deformación cerca de la grieta (encerrada por ) debe ser elástico (no plástico).
La integral J puede calcularse con estas condiciones violadas, pero luego . Cuando no se violan, se puede relacionar la tasa de liberación de energía y la integral J con los módulos elásticos y los factores de intensidad de la tensión utilizando [3]
Ejemplo de cálculo de usando la integral J [3] |
---|
Considere la muestra de viga en voladizo doble que se muestra en la figura, donde la grieta se centró en la viga de altura tiene una longitud de y una carga se aplica para abrir la grieta. Suponga que el material es linealmente elástico y que la grieta crece hacia adelante. Considere un camino rectangular que se muestra en la segunda figura: comience en la cara superior de la grieta, (1) suba a la parte superior en, (2) vaya a la derecha pasando la punta de la grieta, (3) vaya hacia abajo en , (4) vaya por la parte inferior hacia la izquierda y (5) vuelva a subir hasta la cara inferior de la grieta. La integral J es cero a lo largo de muchas partes de este camino. El material se descarga efectivamente detrás de la grieta, por lo que tanto la densidad de energía de deformación como la tracción son cero a lo largo de (1) y (2) y, por lo tanto, la integral J. A lo largo de (2) y (4) uno tiene así como (sin tracción en la superficie libre), por lo que la integral J es cero en (2) y (4) también. Esto deja solo (3); asumiendo que uno está lo suficientemente lejos de la grieta en (3), el término de tracción es cero ya que y lejos de la grieta, dejando a lo largo de (3), y se debe a la tensión de flexión de una viga en voladizo dónde es la longitud en la página. Utilizando, uno ahora tiene donde el factor de 2 en el frente se debe a que hay 2 vigas en voladizo. Resolviendo, |
Métodos computacionales en mecánica de fracturas.
Existen varios métodos para calcular con elementos finitos. Aunque es posible un cálculo directo de la integral J (utilizando las deformaciones y tensiones generadas por FEA ), existen enfoques aproximados para algún tipo de crecimiento de grietas y proporcionan una precisión razonable con cálculos sencillos. Esta sección desarrollará algunos métodos relativamente simples para el análisis de fracturas utilizando simulaciones numéricas.
Método de liberación nodal
Si la grieta crece recta, la tasa de liberación de energía se puede descomponer como una suma de 3 términosasociado con la energía en cada 3 modos. Como resultado, el método de liberación nodal (NR) se puede utilizar para determinara partir de los resultados de FEA. La tasa de liberación de energía se calcula en los nodos de la malla de elementos finitos para la grieta en una longitud inicial y se extiende por una pequeña distancia.. Primero, calculamos la variación de desplazamiento en el nodo de interés(antes y después de que se suelte el nodo de la punta de la grieta). En segundo lugar, realizamos un seguimiento de la fuerza nodal.producido por FEA. Finalmente, podemos encontrar cada componente de usando las siguientes fórmulas:
Si la tasa de liberación de energía excede un valor crítico, la grieta crecerá. En este caso, se realiza una nueva simulación FEA (para el siguiente paso de tiempo) donde se libera el nodo en la punta de la grieta. Para un sustrato delimitado, simplemente podemos dejar de aplicar condiciones de límite de Dirichlet fijas en el nodo de la punta de la grieta del paso de tiempo anterior (es decir, los desplazamientos ya no están restringidos). Para una fisura simétrica, necesitaríamos actualizar la geometría del dominio con una abertura de fisura más larga (y por lo tanto generar una nueva malla [5] ).
Cierre de grietas modificado integral
Similar al método de liberación nodal, la integral de cierre de grietas modificada (MCCI) es un método para calcular la tasa de liberación de energía utilizando desplazamientos nodales FEA y fuerzas . [6] [7] Donderepresenta la dirección correspondiente a los vectores base cartesianos con origen en la punta de la grieta, yrepresenta el índice nodal. MCCI es más eficiente computacionalmente que el método de liberación nodal porque solo requiere un análisis por cada incremento de crecimiento de grietas.
Una condición necesaria para el método MCCI es la longitud uniforme del elemento. a lo largo de la cara de la grieta en el dirección. Además, este método requiere una discretización suficiente de modo que en la longitud de un elemento los campos de tensión sean auto-similares . Esto implica quea medida que se propaga la grieta. A continuación se muestran ejemplos del método MCCI con dos tipos de elementos finitos comunes.
Elementos de 4 nodos
Los elementos lineales cuadrados de 4 nodos que se ven en la Figura 2 tienen una distancia entre los nodos y igual a Considere una grieta con su punta ubicada en el nodo Similar al método de liberación nodal, si la grieta se propagara la longitud de un elemento a lo largo de la línea de simetría (paralela a la eje) el desplazamiento de la apertura de la grieta sería el desplazamiento en la punta de la grieta anterior, es decir y la fuerza en la nueva punta de la grieta sería Dado que se supone que el crecimiento de la grieta es auto-similar, el desplazamiento en el nodo después de que la grieta se propaga es igual al desplazamiento en el nodo antes de que se propague la grieta. Este mismo concepto se puede aplicar a las fuerzas en el nodo y Utilizando el mismo método que se muestra en la sección de liberación nodal, recuperamos las siguientes ecuaciones para la tasa de liberación de energía:
Dónde (desplazamiento por encima y por debajo de la cara de la grieta, respectivamente). Debido a que tenemos una línea de simetría paralela a la grieta, podemos suponer
Por lo tanto,
Elementos de 8 nodos
Los elementos rectangulares de 8 nodos que se ven en la Figura 3 tienen funciones de base cuadrática . El proceso para calcular G es el mismo que el de los elementos de 4 nodos con la excepción de que (el crecimiento de la grieta sobre un elemento) es ahora la distancia desde el nodo a Una vez más, asumiendo un crecimiento de grietas rectas auto-similar, la tasa de liberación de energía se puede calcular utilizando las siguientes ecuaciones:
Al igual que con el método de liberación nodal, la precisión de MCCI depende en gran medida del nivel de discretización a lo largo de la punta de la grieta, es decir La precisión también depende de la elección del elemento. Una malla de elementos cuadráticos de 8 nodos puede producir resultados más precisos que una malla de elementos lineales de 4 nodos con el mismo número de grados de libertad [8] en la malla.
Enfoque integral de dominio para J
La integral J se puede calcular directamente usando las funciones de malla y forma de elementos finitos. [9] Consideramos un contador de dominio como se muestra en la figura 4 y elegimos una función suave arbitraria tal que en y en .
Para grietas elásticas lineales que crecen en línea recta, . La tasa de liberación de energía se puede calcular sobre el área delimitada por el contorno utilizando una formulación actualizada:
La fórmula anterior se puede aplicar a cualquier área anular que rodee la punta de la grieta (en particular, se puede usar un conjunto de elementos vecinos). Este método es muy preciso, incluso con una malla gruesa alrededor de la punta de la grieta (se puede elegir un dominio de integración ubicado lejos, con tensiones y desplazamientos menos sensibles al refinamiento de la malla).
Derivación de la integral J para el enfoque integral de dominio El J-intregral se puede expresar en todo el contorno de la siguiente manera:
Con . , en y el trabajo y el estrés se cancelan en y , por lo tanto, mediante la aplicación del teorema de la divergencia, esto conduce a:
Finalmente, al señalar que y usando la ecuación de equilibrio:
Elementos singulares de punta de grieta 2-D
Los métodos antes mencionados para calcular la tasa de liberación de energía se acercan asintóticamente a la solución real con una mayor discretización, pero no logran capturar completamente la singularidad de la punta de la grieta. Se pueden realizar simulaciones más precisas utilizando elementos de un cuarto de punto alrededor de la punta de la grieta. [10] Estos elementos tienen una singularidad incorporada que produce con mayor precisión campos de tensión alrededor de la punta de la grieta. La ventaja del método de un cuarto de punto es que permite mallas de elementos finitos más gruesas y reduce en gran medida el costo computacional. Además, estos elementos se derivan de pequeñas modificaciones a elementos finitos comunes sin requerir programas computacionales especiales para su análisis. Para los propósitos de esta sección se examinarán los materiales elásticos, aunque este método puede extenderse a la mecánica de fractura elástica- plástica . [11] [12] [13] [14] Suponiendo una elasticidad perfecta, los campos de tensión experimentarán un singularidad de la punta de la grieta.
Elemento isoparamétrico de 8 nodos
El elemento cuadrático de 8 nodos se describe en la Figura 5 en ambos espacios principales con coordenadas locales y y por el elemento mapeado en el espacio físico / global por y El elemento padre se asigna desde el espacio local al espacio físico mediante las funciones de forma y las coordenadas del grado de libertad La punta de la grieta se encuentra en o
De manera similar, los desplazamientos (definidos como ) también se pueden mapear.
Una propiedad de las funciones de forma en el método de elementos finitos es el soporte compacto , específicamente la propiedad delta de Kronecker (es decir, en el nodo y cero en todos los demás nodos). Esto da como resultado las siguientes funciones de forma para los elementos cuadráticos de 8 nodos: [8]
Al considerar una línea frente a la grieta que es colineal con la - eje (es decir ) todas las funciones básicas son cero excepto para
Calcular la deformación normal implica usar la regla de la cadena para tomar la derivada del desplazamiento con respecto a
Si los nodos están espaciados uniformemente en el elemento rectangular, la deformación no contendrá la singularidad. Moviendo los nodos 5 y 8 a una cuarta parte de la longitud del elemento más cercano a la punta de la grieta como se ve en la figura 5, el mapeo de se convierte en:
Resolviendo para y tomar la derivada da como resultado:
Conectando este resultado a la ecuación de deformación se obtiene el resultado final:
Al mover los nodos medios a una posición de un cuarto, se obtiene el resultado correcto singularidad de la punta de la grieta.
Otros tipos de elementos
El método del elemento rectangular no permite que los elementos singulares se enreden fácilmente alrededor de la punta de la grieta. Esto impide la capacidad de capturar la dependencia angular de los campos de tensión que es fundamental para determinar la trayectoria de la fisura. Además, excepto a lo largo de los bordes del elemento,La singularidad existe en una región muy pequeña cerca de la punta de la grieta. La Figura 6 muestra otro método de un cuarto de punto para modelar esta singularidad. El elemento rectangular de 8 nodos se puede mapear en un triángulo. [15] Esto se hace colapsando los nodos en la línea. a la ubicación del nodo medio y desplazando los nodos medios en a la ubicación de un cuarto de punto. El rectángulo colapsado puede rodear más fácilmente la punta de la grieta, pero requiere que los bordes del elemento sean rectos o se reducirá la precisión del cálculo del factor de intensidad de la tensión.
Un mejor candidato para el método de un cuarto de punto es el triángulo natural como se ve en la Figura 7. La geometría del elemento permite rodear fácilmente la punta de la grieta y simplificar el mallado. Siguiendo el mismo procedimiento descrito anteriormente, el campo de desplazamiento y deformación de los elementos triangulares son:
Este método reproduce los dos primeros términos de las soluciones de Williams [16] con un término constante y singular.
Una ventaja del método de un cuarto de punto es que se puede generalizar fácilmente a modelos tridimensionales. Esto puede reducir en gran medida el cálculo en comparación con otros métodos tridimensionales, pero puede provocar errores si la punta de la grieta se propaga con un gran grado de curvatura. [17]
Ver también
- Mecánica de fracturas
- Factor de intensidad de estrés
- Tenacidad a la fractura
- J-integral
Referencias
- ^ Li, FZ; Shih, CF; Needleman, A. (1985). "Una comparación de métodos para calcular las tasas de liberación de energía". Ingeniería Mecánica de Fracturas . 21 (2): 405–421. doi : 10.1016 / 0013-7944 (85) 90029-3 . ISSN 0013-7944 .
- ^ Rice, JR; Budiansky, B. (1973). "Leyes de conservación y tasas de liberación de energía". J. Appl. Mech . 40 (1): 201–3. Código bibliográfico : 1973JAM .... 40..201B . doi : 10.1115 / 1.3422926 . S2CID 13910502 .
- ^ a b c d e f g h i j k l m n o p q Alan Zehnder (2012). Mecánica de fracturas . Londres; Nueva York: Springer Science + Business Media. ISBN 9789400725942.
- ^ Soboyejo, WO (2003). "11.6.5 Equivalencia de G y K". Propiedades mecánicas de los materiales de ingeniería. Marcel Dekker. ISBN 0-8247-8900-8 . OCLC 300921090.
- ^ Tradegard, A. (15 de julio de 1998). "Técnica de remallado FEM aplicada a problemas de crecimiento de grietas". Métodos informáticos en mecánica e ingeniería aplicadas . 160 (1-2): 115-131. Código bibliográfico : 1998CMAME.160..115T . doi : 10.1016 / s0045-7825 (97) 00287-9 .
- ^ Rybicki, EF; Kanninen, MF (enero de 1977). "Un cálculo de elementos finitos de los factores de intensidad de tensión mediante una integral de cierre de fisuras modificada". Ingeniería Mecánica de Fracturas . 9 (4): 931–938. doi : 10.1016 / 0013-7944 (77) 90013-3 . ISSN 0013-7944 .
- ^ Sethuraman, R .; Maiti, SK (enero de 1988). "Cálculo basado en elementos finitos de la tasa de liberación de energía de deformación por integral de cierre de fisura modificado". Ingeniería Mecánica de Fracturas . 30 (2): 227–231. doi : 10.1016 / 0013-7944 (88) 90226-3 . ISSN 0013-7944 .
- ^ a b Zehnder, Alan T. (3 de enero de 2012). Mecánica de fracturas . Dordrecht. ISBN 9789400725959. OCLC 773034407 .
- ^ Zehnder, Alan T. (2012). Mecánica de fracturas . Apuntes de clases en Mecánica Aplicada y Computacional. 62 . Dordrecht: Springer Holanda. doi : 10.1007 / 978-94-007-2595-9 . ISBN 9789400725942.
- ^ Henshell, RD; Shaw, KG (1975). "Los elementos finitos de la punta de la grieta son innecesarios". Revista Internacional de Métodos Numéricos en Ingeniería . 9 (3): 495–507. Código bibliográfico : 1975IJNME ... 9..495H . doi : 10.1002 / nme.1620090302 . ISSN 0029-5981 .
- ^ Barsoum, Roshdy S. (1977). "Elementos triangulares de cuarto de punto como elementos de punta de fisura elásticos y perfectamente plásticos". Revista Internacional de Métodos Numéricos en Ingeniería . 11 (1): 85–98. Código bibliográfico : 1977IJNME..11 ... 85B . doi : 10.1002 / nme.1620110109 . ISSN 0029-5981 .
- ^ Sun, CT; Jin, Z.-H. (2012), "Criterios de fractura elástica-plástica", Fracture Mechanics , Elsevier, págs. 171–187, doi : 10.1016 / b978-0-12-385001-0.00007-9 , ISBN 9780123850010
- ^ Stern, Morris (1979). "Familias de elementos conformes consistentes con campos derivados singulares". Revista Internacional de Métodos Numéricos en Ingeniería . 14 (3): 409–421. Código Bibliográfico : 1979IJNME..14..409S . doi : 10.1002 / nme.1620140307 . ISSN 0029-5981 .
- ^ Levy, N .; Marcal, PV; Ostergren, WJ; Rice, JR (junio de 1971). "Producción a pequeña escala cerca de una grieta en la deformación plana: un análisis de elementos finitos". Revista Internacional de Mecánica de Fracturas . 7 (2). doi : 10.1007 / bf00183802 . ISSN 0020-7268 .
- ^ Barsoum, Roshdy S. (1976). "Sobre el uso de elementos finitos isoparamétricos en mecánica de fracturas lineales". Revista Internacional de Métodos Numéricos en Ingeniería . 10 (1): 25–37. Código bibliográfico : 1976IJNME..10 ... 25B . doi : 10.1002 / nme.1620100103 . ISSN 0029-5981 .
- ^ Williams, ML (1959). "Las tensiones en torno a una falla o fisura en medios disímiles" . Boletín de la Sociedad Sismológica de América . 49 : 199-204.
- ^ Peano, A .; Pasini, A. (febrero de 1982). "Una advertencia contra el mal uso de elementos de un cuarto de punto". Revista Internacional de Métodos Numéricos en Ingeniería . 18 (2): 314–320. Código bibliográfico : 1982IJNME..18..314P . doi : 10.1002 / nme.1620180212 . ISSN 0029-5981 .
enlaces externos
- Notas de mecánica de fracturas no lineales por el profesor John Hutchinson (de la Universidad de Harvard)
- Tasa de liberación de energía de deformación de Griffith en www.fracturemechanics.org