Articulo de referencia

Peridinámica

Modelo computacional del estrechamiento de una varilla de aluminio bajo tensión. Los colores indican el aumento de temperatura debido al calentamiento plástico. El cálculo se re...

Modelo computacional del estrechamiento de una varilla de aluminio bajo tensión. Los colores indican el aumento de temperatura debido al calentamiento plástico. El cálculo se realizó con el código informático Emu utilizando un marco de trabajo peridinámico basado en estados.

La peridinámica es una formulación no local de la mecánica de medios continuos que está orientada hacia deformaciones con discontinuidades, especialmente fracturas . Originalmente, se introdujo la peridinámica basada en enlaces , [ 1 ] en la que las fuerzas de interacción internas entre un punto material y todos los demás con los que puede interactuar se modelan como un campo de fuerza central . [ 2 ] Este tipo de campo de fuerza puede imaginarse como una malla de enlaces que conecta cada punto del cuerpo con todos los demás puntos interactuantes dentro de una cierta distancia que depende de una propiedad del material, llamada horizonte peridinámico . Posteriormente, para superar las limitaciones del marco basado en enlaces para la relación de Poisson del material [ 3 ] [ 4 ] (1/3{\displaystyle 1/3}para tensión plana y1/4{\displaystyle 1/4}para deformación plana en configuraciones bidimensionales;1/4{\displaystyle 1/4}Para los tridimensionales), se ha formulado la peridinámica basada en estados . [ 5 ] Su característica principal es que la fuerza intercambiada entre un punto y otro está influenciada por el estado de deformación de todos los demás enlaces en relación con su zona de interacción. [ 1 ]

La característica distintiva de la peridinámica, que la diferencia de la mecánica local clásica, es la presencia de enlaces de alcance finito entre dos puntos cualesquiera del cuerpo material: es una característica que se aproxima a formulaciones como teorías discretas de la materia a mesoescala. [ 1 ]

Etimología

El término peridinámico , como adjetivo, fue propuesto en el año 2000 y proviene del prefijo peri- , que significa alrededor , cerca o que rodea ; y de la raíz dyna , que significa fuerza o potencia . El término peridinámica , como sustantivo, es una forma abreviada de la frase modelo peridinámico de mecánica de sólidos. [ 1 ]

Objetivo

Una fractura es una singularidad matemática a la que no se pueden aplicar directamente las ecuaciones clásicas de la mecánica del continuo . La teoría peridinámica se propuso con el propósito de modelar matemáticamente la formación y la dinámica de fracturas en materiales elásticos. [ 1 ] Se basa en ecuaciones integrales , a diferencia de la mecánica del continuo clásica, que se basa en ecuaciones diferenciales parciales . Dado que no existen derivadas parciales en las superficies de las grietas [ 1 ] ni en otras singularidades geométricas , las ecuaciones clásicas de la mecánica del continuo no se pueden aplicar directamente cuando tales características están presentes en una deformación . Las ecuaciones integrales de la teoría peridinámica también son válidas en las singularidades y se pueden aplicar directamente, porque no requieren derivadas parciales. La capacidad de aplicar las mismas ecuaciones directamente en todos los puntos de un modelo matemático de una estructura deformable ayuda al enfoque peridinámico a evitar la necesidad de técnicas especiales de mecánica de fractura como xFEM . [ 6 ] Por ejemplo, en peridinámica, no hay necesidad de una ley de crecimiento de grietas separada basada en un factor de intensidad de tensión . [ 7 ]

Definición y terminología básica

(a) Cinemática del cuerpo materialΩt{\displaystyle \Omega _{t}}dentro de la teoría peridinámica. (b) Representación del horizonte peridinámico deincógnita{\displaystyle {\bf {x}}}.

En el contexto de la teoría peridinámica, los cuerpos físicos se tratan como constituidos por una malla continua de puntos que pueden intercambiar fuerzas de interacción mutua de largo alcance, dentro de una distancia máxima y bien establecida.δ>0{\displaystyle \delta >0}: el radio del horizonte peridinámico . Esta perspectiva se acerca mucho más a la dinámica molecular que a los cuerpos macroscópicos y, en consecuencia, no se basa en el concepto de tensor de tensión (que es un concepto local) y tiende hacia la noción de fuerza por pares que un punto materialincógnita{\displaystyle {\bf {x}}}intercambios dentro de su horizonte peridinámico. Desde un punto de vista lagrangiano , adecuado para pequeños desplazamientos, el horizonte peridinámico se considera fijo en la configuración de referencia y, por lo tanto, se deforma con el cuerpo. [ 3 ] Consideremos un cuerpo material representado porΩRnorte{\displaystyle \Omega \subset \mathbb {R} ^{n}}, dóndenorte{\displaystyle n}Puede ser 1, 2 o 3. El cuerpo tiene una densidad positiva.ρ{\displaystyle \rho }Su configuración de referencia en el instante inicial se denota porΩ0Rnorte{\displaystyle \Omega _{0}\subset \mathbb {R} ^{n}}. La configuración de referencia puede ser la configuración libre de estrés o una configuración específica del cuerpo elegida como referencia. En el contexto de la peridinámica, cada punto enΩ{\displaystyle \Omega }interactúa con todos los puntosincógnita{\displaystyle {\bf {x}}'}dentro de un determinado vecindario definido pord(incógnita,incógnita)δ{\displaystyle d({\bf {x}},{\bf {x}}')\leq \delta }, dóndeδ>0{\displaystyle \delta >0}yd(,){\displaystyle d(\cdot ,\cdot )}representa una función de distancia adecuada enΩ0{\displaystyle \Omega _{0}}Este barrio suele ser conocido comoBδ(incógnita){\displaystyle B_{\delta }({\bf {x}})}en la literatura. Se conoce comúnmente como el horizonte [ 7 ] [ 8 ] o la familia deincógnita{\displaystyle {\bf {x}}}. [ 3 ] [ 9 ]

La cinemática deincógnita{\displaystyle {\bf {x}}}se describe en términos de su desplazamiento desde la posición de referencia, denotada como(incógnita,t):Ω0×R+Rnorte{\displaystyle {\bf {u}}({\bf {x}},t):\Omega _{0}\times \mathbb {R} ^{+}\rightarrow \mathbb {R} ^{n}}. En consecuencia, la posición deincógnita{\displaystyle {\bf {x}}}en un momento específicot{\displaystyle t}está determinado pory(incógnita,t):=incógnita+(incógnita,t){\displaystyle {\bf {y}}({\bf {x}},t):={\bf {x}}+{\bf {u}}({\bf {x}},t)}Además, para cada par de puntos interactuantes, se realiza un seguimiento del cambio en la longitud del enlace con respecto a la configuración inicial a lo largo del tiempo mediante la tensión relativa.s(incógnita,incógnita,t){\displaystyle s({\bf {x}},{\bf {x}}',t)}, que puede expresarse como:

s(incógnita,incógnita,t)=|(incógnita,t)(incógnita,t)||incógnitaincógnita|,{\displaystyle s\left({\bf {x}},{\bf {x}}',t\right)={\frac {\left|{\bf {u}}\left({\bf {x}}^{\prime },t\right)-{\bf {u}}({\bf {x}},t)\right|}{\left|{\bf {x}}^{\prime }-{\bf {x}}\right|}},}

dónde||{\displaystyle |\cdot |}denota la norma euclidiana [ 3 ] yincógnitaBδ(incógnita)Ω0{\displaystyle {\bf {x}}'\in B_{\delta }({\bf {x}})\cap \Omega _{0}}.

La interacción entre cualquierincógnita{\displaystyle {\bf {x}}}yincógnita{\displaystyle {\bf {x'}}}se denomina enlace . Estos enlaces por pares tienen longitudes variables en el tiempo en respuesta a la fuerza por unidad de volumen al cuadrado, denotada como [ 3 ].

FF(incógnita,incógnita,(incógnita),(incógnita),t){\displaystyle {\bf {f}}\equiv {\bf {f}}({\bf {x}}',{\bf {x}},{\bf {u}}({\bf {x}}'),{\bf {u}}({\bf {x}}),t)}.

Esta fuerza se conoce comúnmente como la función de fuerza por pares o núcleo peridinámico , y abarca todas las propiedades constitutivas (dependientes del material). Describe cómo las fuerzas internas dependen de la deformación. Vale la pena señalar que la dependencia de{\displaystyle {\bf {u}}}ent{\displaystyle t}se ha omitido aquí por razones de simplicidad en la notación. Además, un término de forzamiento externo,b(incógnita,t){\displaystyle \mathbf {b} ({\bf {x}},t)}, se introduce, lo que da como resultado la siguiente ecuación de movimiento, que representa la ecuación fundamental de la peridinámica: [ 3 ]

ρtt(incógnita,t)=F(incógnita,t).{\displaystyle {\rho {\bf {u}}_{tt}({\bf {x}},t)={\bf {F}}({\bf {x}},t)}\,.}

donde el término integralF(incógnita,t){\displaystyle {\bf {F}}({\bf {x}},t)}es la suma de todas las fuerzas internas y externas por unidad de volumen que actúan sobreincógnita{\displaystyle {\bf {x}}}:

F(incógnita,t):=Ω0Bδ(incógnita)F(incógnita,incógnita,(incógnita),(incógnita))dVincógnita+b(incógnita,t).{\displaystyle {{\bf {F}}({\bf {x}},t):=\int _{\Omega _{0}\cap B_{\delta }({\bf {x}})}{\bf {f}}\left({\bf {x}}',{\bf {x}},{\bf {u}}\left({\bf {x}}'\right),{\bf {u}}({\bf {x}})\right)dV_{{\bf {x}}'}+{\bf {b}}({\bf {x}},t)}\,.}

La función con valores vectorialesF{\displaystyle {\bf {f}}}es la densidad de fuerza queincógnita{\displaystyle {\bf {x'}}}ejerce sobreincógnita{\displaystyle {\bf {x}}}Esta densidad de fuerza depende de los vectores de desplazamiento relativo y posición relativa entreincógnita{\displaystyle {\bf {x'}}}yincógnita{\displaystyle {\bf {x}}}. La dimensión deF{\displaystyle {\bf {f}}}es[norte/metro6]{\displaystyle [N/m^{6}]}. [ 3 ]

Peridinámica basada en enlaces

En esta formulación de la peridinámica, el núcleo está determinado por la naturaleza de las fuerzas internas y las restricciones físicas que rigen la interacción entre solo dos puntos materiales. Para mayor brevedad, se definen las siguientes cantidades. ξ:=incógnitaincógnita{\displaystyle {\bf {\bf {\xi }}}:={\bf {x}}'-{\bf {x}}}yη:=(incógnita)(incógnita){\displaystyle {\bf {\eta }}:={\bf {u}}({\bf {x}}')-{\bf {u}}({\bf {x}})}para que [ 1 ]

F(incógnitaincógnita,(incógnita)(incógnita))F(ξ,η){\displaystyle {\bf {f}}({\bf {x}}'-{\bf {x}},{\bf {u}}({\bf {x}}')-{\bf {u}}({\bf {x}}))\equiv {\bf {{f}({\bf {\xi }},{\bf {\eta }})}}}

Principio de acción y reacción

Para cualquier incógnita{\displaystyle {\bf {x}}} y incógnita{\displaystyle {\bf {x'}}} perteneciente al vecindario Bδ(incógnita){\displaystyle B_{\delta }({\bf {x}})}Se cumple la siguiente relación: F(η,ξ)=F(η,ξ){\displaystyle {\bf {f}}(-\eta ,-\xi )=-{\bf {f}}(\eta ,\xi )} Esta expresión refleja el principio de acción y reacción, comúnmente conocido como la tercera ley de Newton. Garantiza la conservación del momento lineal en un sistema compuesto por partículas que interactúan entre sí. [ 1 ]

Conservación del momento angular

Para cualquierincógnita{\displaystyle {\bf {x}}} yincógnita{\displaystyle {\bf {{x}'}}} perteneciente al vecindarioBδ(incógnita){\displaystyle B_{\delta }({\bf {x}})}Se cumple la siguiente condición:(ξ+η)×F(ξ,η)=0{\displaystyle (\xi +\eta )\times {\bf {f}}(\xi ,\eta )=0}Esta condición surge al considerar el rayo- vector deformado relativo que conectaincógnita{\displaystyle {\bf {x}}} yincógnita{\displaystyle {\bf {{x}'}}}comoξ+η{\displaystyle \xi +\eta }. La condición se satisface si y solo si el vector de densidad de fuerza por pares tiene la misma dirección que el vector de rayo deformado relativo. En otras palabras,F(ξ,η)=F(ξ,η)(ξ+η){\displaystyle {\bf {f}}(\xi ,\eta )=f(\xi ,\eta )(\xi +\eta )} a pesar deξ{\displaystyle \xi } yη{\displaystyle \eta }, dóndeF(ξ,η){\displaystyle f(\xi ,\eta )} es una función escalar. [ 1 ]

Material hiperelástico

Un material hiperelástico es un material con una relación constitutiva tal que: [ 1 ]

ΓF(ξ,η)dη=0, curva cerrada Γ,    ξ0,{\displaystyle \int _{\Gamma }{\bf {f}}({\bf {\xi }},{\bf {\eta }})\cdot d{\bf {\eta }}=0\,,\quad \forall {\text{ closed curve }}\Gamma ,\ \ \ \ \forall {\bf {\xi }}\neq {\bf {{0},}}}

o, equivalentemente, por el teorema de Stokes

η×F(ξ,η)=0{\displaystyle \nabla _{\bf {\eta }}\times {\bf {f}}({\bf {\xi }},{\bf {\eta }})={\bf {{0}\,}}} ,ξ,η{\displaystyle \forall \,{\bf {\xi }},\,{\bf {\eta }}}

y, por lo tanto,

F(ξ,η)=ηΦ(ξ,η)ξ,η.{\displaystyle {\bf {f}}({\bf {\xi }},{\bf {\eta }})=\nabla _{\bf {\eta }}\Phi ({\bf {\xi }},\,{\bf {\eta }})\,\forall {\bf {\xi }},\,{\bf {\eta }}\,.}

En la ecuación anteriorΦ(ξ,η){\displaystyle \Phi ({\bf {\xi }},{\bf {\eta }})}es la función potencial de valor escalar endo2(Rnorte{0}×Rnorte){\displaystyle C^{2}(\mathbb {R} ^{n}\setminus {\bf {{\{0\}}\times \mathbb {R} ^{n})}}}. [ 1 ] Debido a la necesidad de satisfacer la conservación del momento angular , la condición siguiente sobre la función escalarF(ξ,η){\displaystyle f({\bf {\xi }},{\bf {\eta }})}sigue [ 1 ]

F(ξ,η)η=gramo(ξ,η)(ξ+η).{\displaystyle {\frac {\partial f({\bf {\xi }},{\bf {\eta }})}{\partial {\bf {\eta }}}}=g({\bf {\xi }},{\bf {\eta }})({\bf {\xi }}+{\bf {\eta }}).}

dóndegramo(ξ,η){\displaystyle g({\bf {\xi }},{\bf {\eta }})}es una función escalar. Integrando ambos lados de la ecuación, la siguiente condición engramo(ξ,η){\displaystyle g({\bf {\xi }},{\bf {\eta }})}se obtiene [ 1 ]

F(ξ,η)=h(|ξ+η|,ξ)(ξ+η){\displaystyle {\bf {f}}({\bf {\xi }},{\bf {\eta }})=h(|{\bf {\xi }}+{\bf {\eta }}|,{\bf {\xi }})({\bf {\xi }}+{\bf {\eta }})},

parah(|ξ+η|,ξ){\displaystyle h(|{\bf {\xi }}+{\bf {\eta }}|,{\bf {\xi }})}una función de valor escalar. La naturaleza elástica deF{\displaystyle {\bf {f}}}Es evidente: la fuerza de interacción depende únicamente de la posición relativa inicial entre los puntos.incógnita{\displaystyle {\bf {x}}}yincógnita{\displaystyle {\bf {x}}'}y el módulo de su posición relativa,|ξ+η|{\displaystyle |{\bf {\xi }}+{\bf {\eta }}|}, en la configuración deformadaΩt{\displaystyle \Omega _{t}}en ese momentot{\displaystyle t}Aplicando la hipótesis de isotropía , la dependencia del vectorξ{\displaystyle {\bf {\xi }}}puede sustituirse por una dependencia de su módulo|ξ|{\displaystyle |{\bf {\xi }}|}, [ 1 ]

F(ξ,η)=h(|ξ+η|,|ξ|)(ξ+η).{\displaystyle {\bf {f}}({\bf {\xi }},{\bf {\eta }})=h(|{\bf {\xi }}+{\bf {\eta }}|,|{\bf {\xi }}|)({\bf {\xi }}+{\bf {\eta }}).}

Por lo tanto, las fuerzas de enlace pueden considerarse como un modelo de una red de resortes que conecta cada punto. incógnitaΩ0{\displaystyle {\bf {x}}\in \Omega _{0}}por pares conincógnitaBδ(incógnita)Ω0{\displaystyle {\bf {x}}'\in B_{\delta }({\bf {x}})\cap \Omega _{0}}.

Material elástico lineal

Si|η|1{\displaystyle |{\bf {\eta }}|\ll 1}, el núcleo peridinámico puede linealizarse alrededor deη=0{\displaystyle {\bf {\eta }}={\bf {0}}}: [ 1 ]

F(ξ,η)F(ξ,0)+F(ξ,η)η|η=0η;{\displaystyle {\bf {f}}({\bf {\xi }},{\bf {\eta }})\approx {\bf {f}}({\bf {\xi }},{\bf {{0})+\left.{\frac {\partial {\bf {f}}({\bf {\xi }},{\bf {\eta }})}{\partial {\bf {\eta }}}}\right|_{{\bf {\eta }}={\bf {0}}}{\bf {\eta }};}}}

Entonces, se puede definir un tensor de micromódulo de segundo orden como

do(ξ)=F(ξ,η)η|η=0=ξF(ξ,η)η|η=0+F0I{\displaystyle {\bf {C}}({\bf {\xi }})=\left.{\frac {\partial {\bf {f}}({\bf {\xi }},{\bf {\eta }})}{\partial {\bf {\eta }}}}\right|_{{\bf {\eta }}={\bf {0}}}={\bf {\xi }}\otimes \left.{\frac {\partial f({\bf {\xi }},{\bf {\eta }})}{\partial {\bf {\eta }}}}\right|_{{\bf {\eta }}={\bf {0}}}+f_{0}I}

dóndeF0:=F(ξ,0){\displaystyle f_{0}:=f({\bf {\xi }},{\bf {0}})}yI{\displaystyle I}es el tensor identidad. Tras la aplicación del balance de momento lineal, la elasticidad y la condición de isotropía, el tensor de micromódulo puede expresarse de esta forma [ 1 ].

do(ξ)=λ(|ξ|)ξξ+F0I.{\displaystyle {\bf {C}}({\bf {\xi }})=\lambda (|{\bf {\xi }}|){\bf {\xi }}\otimes {\bf {\xi }}+f_{0}I.}

Por lo tanto, para un material hiperelástico linealizado, su núcleo peridinámico tiene la siguiente estructura [ 1 ].

F(ξ,η)F(ξ,0)+(λ(|ξ|)ξξ+F0I)η.{\displaystyle {\bf {f}}({\bf {\xi }},{\bf {\eta }})\approx {\bf {f}}({\bf {\xi }},{\bf {0}})+\left(\lambda (|{\bf {\xi }}|){\bf {\xi }}\otimes {\bf {\xi }}+f_{0}I\right){\bf {\eta }}.}

Expresiones para el núcleo peridinámico

El núcleo peridinámico es una función versátil que caracteriza el comportamiento constitutivo de los materiales dentro del marco de la teoría peridinámica. Una formulación comúnmente empleada del núcleo se utiliza para describir una clase de materiales conocidos como materiales microelásticos frágiles prototipo (PMB). En el caso de materiales PMB isotrópicos, se supone que la fuerza entre pares es linealmente proporcional al estiramiento finito [ 7 ] experimentado por el material, definido como

s:=(|ξ+η||ξ|)/|ξ|{\displaystyle s:=(|{\bf {\xi }}+{\bf {\eta }}|-|{\bf {\xi }}|)/|{\bf {\xi }}|},

de modo que

F(η,ξ)=F(|ξ+η|,|ξ|)norte,{\displaystyle \mathbf {f} ({\bf {\eta }},{\bf {\xi }})=f(|{\bf {\xi }}+{\bf {\eta }}|,|{\bf {\xi }}|){\bf {{n},}}}

dónde

norte:=(ξ+η)/|ξ+η|{\displaystyle {\bf {{n}:=({\bf {\xi }}+{\bf {\eta }})/|{\bf {\xi }}+{\bf {\eta }}|}}}

y donde la función escalarF{\displaystyle f}se define de la siguiente manera [ 7 ]

F=dosμ(s,t)=do|ξ+η||ξ||ξ|μ(s,t),{\displaystyle f=cs\mu (s,t)=c\;{\frac {|{\bf {\xi }}+{\bf {\eta }}|-|{\bf {\xi }}|}{|{\bf {\xi }}|}}\mu (s,t),}con

μ(s,t)={1, si s(t,ξ)<s0,0, de lo contrario,      a pesar de 0tt;{\displaystyle \mu (s,t)=\left\{{\begin{array}{ll}1\,,&{\text{ if }}s\left(t^{\prime },{\bf {\xi }}\right)<s_{0}\,,\\0\,,&{\text{ otherwise, }}\end{array}}\ \ \ \ {\text{ for all }}0\leq t^{\prime }\leq t\right.;}

La constantedo{\displaystyle c}se denomina constante de micromódulo y la funciónμ(s,t){\displaystyle \mu (s,t)}sirve para indicar si, en un momento dadott{\displaystyle t'\leq t}, el estiramiento del enlaces{\displaystyle s}asociado con el par(incógnita,incógnita){\displaystyle ({\bf {x,\,x'}})}ha superado el valor críticos0{\displaystyle s_{0}}. Si se supera el valor crítico, el enlace se considera roto y se asigna una fuerza de pares de cero para todostt{\displaystyle t\geq t'}. [ 1 ]

Después de una comparación entre el valor de la densidad de energía de deformación obtenido bajo extensión isotrópica empleando respectivamente el marco de la teoría de la peridinámica y la teoría del continuo clásica, el valor físico coherente del micromódulodo{\displaystyle c}se puede encontrar [ 7 ]

do=18kπδ4,{\displaystyle c={\frac {18k}{\pi \delta ^{4}}},}

dóndek{\displaystyle k}es el módulo de compresibilidad del material .

Siguiendo el mismo enfoque [ 10 ] la constante del micromódulodo{\displaystyle c}puede extenderse ado(ξ,δ){\displaystyle c({\bf {\xi }},\delta )}, dóndedo{\displaystyle c}Ahora es una función de micromódulo . Esta función proporciona una descripción más detallada de cómo se distribuye la intensidad de las fuerzas entre pares a lo largo del horizonte peridinámico.Bδ(incógnita){\displaystyle B_{\delta }({\bf {x}})}Intuitivamente, la intensidad de las fuerzas disminuye a medida que aumenta la distancia entreincógnita{\displaystyle {\bf {x}}}yincógnitaBδ(incógnita){\displaystyle {\bf {x}}'\in B_{\delta }({\bf {x}})}aumenta, pero la forma específica en que se produce esta disminución puede variar.

La función del micromódulo se expresa como [ 11 ].

do(ξ,δ):=do(0,δ)k(ξ,δ),{\displaystyle c({\bf {\xi }},\delta ):=c({\bf {{0},\delta )k({\bf {\xi }},\delta )\,,}}}

donde la constantedo(0,δ){\displaystyle c({\bf {{0},\delta )}}}se obtiene comparando la densidad de deformación peridinámica con las teorías mecánicas clásicas; [ 12 ]k(ξ,δ){\displaystyle k({\bf {\xi }},\delta )}es una función definida enΩ0{\displaystyle \Omega _{0}}con las siguientes propiedades (dadas las restricciones de conservación del momento e isotropía) [ 11 ]

{k(ξ,δ)=k(ξ,δ),límiteξ0k(ξ,δ)=máximoξ Rnorte{k(ξ,δ)},límiteξδk(ξ,δ)=0,Rnortelímiteδ0k(ξ,δ)dincógnita=RnorteΔ(ξ)dincógnita=1,{\displaystyle \left\{{\begin{array}{l}k({\bf {\xi }},\delta )=k(-{\bf {\xi }},\delta )\,,\\\lim _{{\bf {\xi }}\rightarrow {\bf {0}}}k({\bf {\xi }},\delta )=\max _{{\bf {\xi }}\ \in \mathbb {R} ^{n}}\{k({\bf {\xi }},\delta )\}\,,\\\lim _{{\bf {\xi }}\rightarrow \delta }k({\bf {\xi }},\delta )=0\,,\\\int _{\mathbb {R} ^{n}}\lim _{\delta \rightarrow 0}k({\bf {\xi }},\delta )d{\bf {x}}=\int _{\mathbb {R} ^{n}}\Delta ({\bf {\xi }})d{\bf {x}}=1\,,\end{array}}\right.}

dóndeΔ(ξ){\displaystyle \Delta ({\bf {\xi }})}es la función delta de Dirac .

Micromódulo Peridynamics ampliamente utilizado
Representación gráfica de algunas funciones de micromódulo ampliamente utilizadas.do(ξ,δ)=do(0,δ)k(ξ,δ){\displaystyle c({\bf {\xi }},\delta )=c({\bf {{0},\delta )k({\bf {\xi }},\delta )}}}.

Micromódulo cilíndrico

La expresión más simple para la función de micromódulo es [ 11 ].

do(0,δ)k(ξ,δ)=do1Bδ(incógnita){\displaystyle c({\bf {{0},\delta )k({\bf {\xi }},\delta )=c{\bf {{1}_{B_{\delta }({\bf {x}}')}}}}}},

dónde1A{\displaystyle {\bf {{1}_{A}}}}:incógnitaR{\displaystyle X\rightarrow \mathbb {R} }es la función indicadora del subconjuntoAincógnita{\displaystyle A\subset X}, definido como

1A(incógnita):={1,incógnitaA,0,incógnitaA,;{\displaystyle \mathbf {1} _{A}(x):={\begin{cases}1,&x\in A\,,\\0,&x\notin A\,,\end{cases}}\;\;;}

Micromódulo triangular

Se caracteriza pork(ξ,δ){\displaystyle k({\bf {\xi }},\delta )}ser una función lineal [ 13 ]

k(ξ,δ)=(1|ξ|δ)1Bδ(incógnita).{\displaystyle k({\bf {\xi }},\delta )=\left(1-{\frac {|{\bf {\xi }}|}{\delta }}\right){\bf {{1}_{B_{\delta }({\bf {x}}')}.}}}

micromódulo normal

Si se quiere reflejar el hecho de que la mayoría de los sistemas físicos discretos comunes se caracterizan por una distribución de Maxwell-Boltzmann , para incluir este comportamiento en la peridinámica, la siguiente expresión parak(ξ,δ){\displaystyle k({\bf {\xi }},\delta )}puede utilizarse [ 14 ]

k(ξ,δ)=mi(|ξ|/δ)21Bδ(incógnita);{\displaystyle k({\bf {\xi }},\delta )=e^{-(|{\bf {\xi }}|/\delta )^{2}}{\bf {{1}_{B_{\delta }({\bf {x}}')};}}}

micromódulo cuártico

En la literatura también se puede encontrar la siguiente expresión para lak(ξ,δ){\displaystyle k({\bf {\xi }},\delta )}función [ 11 ]

k(ξ,δ)=(1(ξδ)2)21Bδ(incógnita).{\displaystyle k({\bf {\xi }},\delta )=\left(1-\left({\frac {\xi }{\delta }}\right)^{2}\right)^{2}{\bf {{1}_{B_{\delta }({\bf {x}}')}.}}}

En general, dependiendo de la propiedad específica del material que se desee modelar, existe una amplia gama de expresiones para el micromódulo y, en general, para el núcleo peridinámico. Por lo tanto, la lista anterior no es exhaustiva. [ 11 ]

Daño

Representación de la función de fuerza peridinámica por paresF(ξ,η){\displaystyle {\bf {f}}(\xi ,\eta )}con función de ruptura de enlacesμ(s,t){\displaystyle \mu (s,t)}; después del valor de estiramiento críticos0{\displaystyle s_{0}}Si se supera ese límite, se considera que el vínculo se ha roto y no existe fuerza entre los dos puntos materiales involucrados.

El daño se incorpora a la función de fuerza por pares al permitir que los enlaces se rompan cuando su elongación supera un valor preestablecido. Tras la rotura de un enlace, este deja de soportar fuerza alguna y sus extremos quedan efectivamente desconectados. Cuando un enlace se rompe, la fuerza que soportaba se redistribuye a otros enlaces que aún no se han roto. Este aumento de carga incrementa la probabilidad de que estos otros enlaces se rompan. El proceso de rotura de enlaces y redistribución de la carga, que conduce a nuevas roturas, es el mecanismo por el cual se propagan las grietas en el modelo peridinámico. [ 7 ]

Analíticamente, la ruptura del enlace se especifica dentro de la expresión del núcleo peridinámico, mediante la función [ 7 ].

μ(s,t)={1, si s(t,ξ)<s0,0, de lo contrario,      a pesar de 0tt;{\displaystyle \mu (s,t)=\left\{{\begin{array}{ll}1\,,&{\text{ if }}s\left(t^{\prime },{\bf {\xi }}\right)<s_{0}\,,\\0\,,&{\text{ otherwise, }}\end{array}}\ \ \ \ {\text{ for all }}0\leq t^{\prime }\leq t\right.;}

Si la gráfica deF(s,t){\displaystyle {\bf {f}}(s,t)}versus estiramiento de enlacess{\displaystyle s}Se representa gráficamente la acción de la función de ruptura de enlaces.μ{\displaystyle \mu }en la formación de fracturas es claro. Sin embargo, no solo la fractura abrupta puede modelarse en el marco peridinámico y expresiones más generales paraμ{\displaystyle \mu }puede emplearse. [ 7 ]

Peridinámica basada en estados

La teoría descrita anteriormente supone que cada enlace peridinámico responde independientemente de los demás. Esto es una simplificación excesiva para la mayoría de los materiales y conlleva restricciones en los tipos de materiales que se pueden modelar. En particular, esta suposición implica que cualquier sólido elástico lineal isótropo está restringido a una relación de Poisson de 1/4. [ 3 ]

Para abordar esta falta de generalidad, se introdujo el concepto de estados peridinámicos . Este marco permite que la densidad de fuerza en cada enlace dependa de los estiramientos en todos los enlaces conectados a sus extremos, además de su propio estiramiento. Por ejemplo, la fuerza en un enlace podría depender de los cambios de volumen netos en los extremos. El efecto de este cambio de volumen, en relación con el efecto del estiramiento del enlace, determina el coeficiente de Poisson . Con los estados peridinámicos, cualquier material que pueda modelarse dentro de la teoría estándar de la mecánica de medios continuos puede modelarse como un material peridinámico, conservando las ventajas de la teoría peridinámica para la fractura. [ 5 ]

Matemáticamente, la ecuación del término de fuerza interna y externa es

F(incógnita,t):=Ω0Bδ(incógnita)F(incógnita,incógnita,(incógnita),(incógnita))dVincógnita+b(incógnita,t).{\displaystyle {{\bf {F}}({\bf {x}},t):=\int _{\Omega _{0}\cap B_{\delta }({\bf {x}})}{\bf {f}}\left({\bf {x}}',{\bf {x}},{\bf {u}}\left({\bf {x}}'\right),{\bf {u}}({\bf {x}})\right)dV_{{\bf {x}}'}+{\bf {b}}({\bf {x}},t)}\,.}

El utilizado en las formulaciones basadas en enlaces se sustituye por [ 5 ].F(incógnita,t):=Bδ(incógnita){T_[incógnita,t]incógnitaincógnitaT_[incógnita,t]incógnitaincógnita}dVincógnita+b(incógnita,t),{\displaystyle {\bf {F}}({\bf {x}},t):=\int _{B_{\delta }({\bf {x}})}\left\{{\underline {\mathbf {T} }}[\mathbf {x} ,t]\left\langle \mathbf {x} ^{\prime }-\mathbf {x} \right\rangle -{\underline {\mathbf {T} }}\left[\mathbf {x} ^{\prime },t\right]\left\langle \mathbf {x} -\mathbf {x} ^{\prime }\right\rangle \right\}dV_{\mathbf {x} ^{\prime }}+\mathbf {b} (\mathbf {x} ,t),}

dóndeT_{\displaystyle {\underline {\mathbf {T} }}}es el campo de estado del vector de fuerza.

Un estado general de orden mA_:Bδ(incógnita)Lmetro.{\displaystyle {\underline {\mathbf {A} }}\langle \cdot \rangle :B_{\delta }({\bf {x}})\rightarrow {\mathcal {L}}_{m}.}es un objeto matemático similar a un tensor , con la excepción de que es [ 5 ]

  • En general, no lineal;
  • en general no continuo;
  • no es de dimensión finita.

Los estados vectoriales son estados de orden igual a 2. Para el llamado material simple ,T_{\displaystyle {\underline {\mathbf {T} }}}se define como

T_:=T^_(Y_){\displaystyle {\underline {\mathbf {T} }}:={\underline {\mathbf {\hat {T}} }}({\underline {\mathbf {Y} }})}

dóndeT^_:VV{\displaystyle {\underline {\mathbf {\hat {T}} }}:{\mathcal {V}}\rightarrow {\mathcal {V}}}es una función integrable de Riemann enBδ(incógnita){\displaystyle B_{\delta }({\bf {x}})}, y Y_{\displaystyle {\underline {\mathbf {Y} }}}se denomina campo de estado del vector de deformación y se define mediante la siguiente relación [ 5 ]

Y_[incógnita,t]ξ=y(incógnita+ξ,t)y(incógnita,t)incógnitaΩ0,ξBδ(incógnita),t0{\displaystyle {\underline {\mathbf {Y} }}[\mathbf {x} ,t]\langle {\boldsymbol {\xi }}\rangle =\mathbf {y} (\mathbf {x} +{\boldsymbol {\xi }},t)-\mathbf {y} (\mathbf {x} ,t)\quad \forall \mathbf {x} \in \Omega _{0},\xi \in B_{\delta }({\bf {x}}),t\geq 0}

de este modoY_incógnitaincógnita{\displaystyle {\underline {\mathbf {Y} }}\left\langle \mathbf {x} ^{\prime }-\mathbf {x} \right\rangle }es la imagen del vínculoincógnitaincógnita{\displaystyle \mathbf {x} ^{\prime }-\mathbf {x} }bajo la deformación

de tal manera que

Y_ξ=0 si y solo si ξ=0,{\displaystyle {\underline {\mathbf {Y} }}\langle {\boldsymbol {\xi }}\rangle =\mathbf {0} {\text{ if and only if }}{\boldsymbol {\xi }}=\mathbf {0} ,}

lo que significa que dos partículas distintas nunca ocupan el mismo punto a medida que avanza la deformación. [ 5 ]

Se puede demostrar [ 5 ] que el equilibrio del momento lineal se deriva de la definición deF(incógnita,t){\displaystyle {\bf {F}}({\bf {x,\,t}})}, mientras que, si la relación constitutiva es tal que

Bδ(incógnita)Y_ξ×T_ξdVξ=0Y_V{\displaystyle \int _{B_{\delta }({\bf {x}})}{\underline {\mathbf {Y} }}\langle {\boldsymbol {\xi }}\rangle \times {\underline {\mathbf {T} }}\langle {\boldsymbol {\xi }}\rangle dV_{\boldsymbol {\xi }}=0\quad \forall {\underline {\mathbf {Y} }}\in {\mathcal {V}}}

El campo de estado del vector de fuerza satisface el equilibrio del momento angular. [ 5 ]

Aplicaciones

Fractura dúctil de una aleación de Al-Mg- Si

El creciente interés en la peridinámica [ 6 ] proviene de su capacidad para llenar el vacío entre las teorías atomísticas de la materia y la mecánica clásica de medios continuos locales. Se aplica eficazmente a fenómenos a microescala, como la formación y propagación de grietas , [ 15 ] [ 16 ] [ 17 ] la dispersión de ondas , [ 18 ] [ 19 ] la fractura intragranular. [ 20 ] Estos fenómenos pueden describirse mediante un ajuste apropiado del radio del horizonte peridinámico, que está directamente vinculado a la extensión de las interacciones no locales entre puntos dentro del material. [ 21 ]

Además de los campos de investigación mencionados, el enfoque no local de la peridinámica para las discontinuidades ha encontrado aplicaciones en varias otras áreas. En geomecánica , se ha empleado para estudiar grietas en el suelo inducidas por agua, [ 22 ] [ 23 ] fallas de geomateriales , [ 24 ] fragmentación de rocas, [ 25 ] [ 26 ] y demás. En biología , la peridinámica se ha utilizado para modelar interacciones de largo alcance en tejidos vivos , [ 27 ] rupturas celulares, agrietamiento de biomembranas , [ 28 ] y más. [ 6 ] Además, la peridinámica se ha extendido a la teoría de difusión térmica , [ 29 ] [ 30 ] permitiendo el modelado de la conducción de calor en materiales con discontinuidades, defectos, inhomogeneidades y grietas. También se ha aplicado para estudiar fenómenos de advección-difusión en fluidos multifásicos [ 31 ] y para construir modelos para problemas transitorios de advección-difusión. [ 32 ] Gracias a su versatilidad, la peridinámica se ha utilizado en diversos análisis multifísicos , incluyendo análisis microestructurales , [ 33 ] fatiga y conducción de calor en materiales compuestos, [ 34 ] [ 35 ] corrosión galvánica en metales, [ 36 ] grietas inducidas por electricidad en materiales dieléctricos y más. [ 6 ]

Véase también

Referencias

  1. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 Silling , SA ( enero de 2000 ) . " Reformulación de la teoría de la elasticidad para discontinuidades y fuerzas de largo alcance" . Journal of the Mechanics and Physics of Solids . 48 (1): 175– 209. Bibcode : 2000JMPSo..48..175S . doi : 10.1016/S0022-5096(99)00029-0 . S2CID 122055539 . 
  2. Dimola, Nunzio; Coclite, Alessandro; Fanizza, Giuseppe; Politi, Tiziano (23-10-2022). "Peridinámica basada en enlaces, una revisión que explora teorías no locales de la dinámica de fluidos" . Advances in Continuous and Discrete Models . 2022 (1) 60. arXiv : 2207.06194 . doi : 10.1186/s13662-022-03732-6 . ISSN 2731-4235 . 
  3. 1 2 3 4 5 6 7 8 Madenci, Erdogan; Oterkus, Erkan (2014). Teoría peridinámica y sus aplicaciones . Nueva York, NY: Springer. págs. 19–43 . ISBN  978-1-4614-8464-6.
  4. Macek, Richard W.; Silling, Stewart A. (noviembre de 2007). "Peridinámica mediante análisis de elementos finitos". Finite Elements in Analysis and Design . 43 (15): 1169– 1178. doi : 10.1016/j.finel.2007.08.012 . ISSN 0168-874X . OSTI 1725746 .  
  5. 1 2 3 4 5 6 7 8 Silling, SA; Epton, M.; Weckner, O.; Xu, J.; Askari, E. (2007-08-08). "Estados peridinámicos y modelado constitutivo" . Journal of Elasticity . 88 (2): 151– 184. doi : 10.1007/s10659-007-9125-1 . ISSN 0374-3535 . S2CID 30571789 .  
  6. 1 2 3 4 Javili, Ali; Morasata, Rico; Oterkus, Erkan; Oterkus, Selda (noviembre de 2019). "Revisión de peridinámica" . Matemáticas y mecánica de sólidos . 24 (11): 3714– 3739. doi : 10.1177/1081286518803411 . hdl : 11693/53217 . ISSN 1081-2865 . S2CID 162176799 .  
  7. 1 2 3 4 5 6 7 8 Silling, SA; Askari, E. (junio de 2005). "Un método sin malla basado en el modelo peridinámico de la mecánica de sólidos" . Computers & Structures . 83 ( 17–18 ): 1526–1535 . doi : 10.1016/j.compstruc.2004.11.026 .
  8. Ren, Huilong; Zhuang, Xiaoying ; Cai, Yongchang; Rabczuk, Timon (2016-12-21). "Peridinámica de doble horizonte: Peridinámica de doble horizonte" . International Journal for Numerical Methods in Engineering . 108 (12): 1451– 1476. arXiv : 1506.05146 . doi : 10.1002/nme.5257 . S2CID 117201049 . 
  9. Chen, Ziguang; Bakenhus, Drew; Bobaru, Florin (noviembre de 2016). "Un núcleo peridinámico constructivo para la elasticidad" . Métodos informáticos en mecánica aplicada e ingeniería . 311 : 356–373 . Bibcode : 2016CMAME.311..356C . doi : 10.1016/j.cma.2016.08.012 .
  10. Bobaru, Florin; Duangpanya, Monchai (septiembre de 2010). "La formulación peridinámica para la conducción transitoria de calor" . International Journal of Heat and Mass Transfer . 53 ( 19–20 ): 4047–4059 . Bibcode : 2010IJHMT..53.4047B . doi : 10.1016/j.ijheatmasstransfer.2010.05.024 .
  11. 1 2 3 4 5 Huang, Dan; Lu, Guangda; Wang, Chongwen; Qiao, Pizhong (junio de 2015). "Un enfoque peridinámico extendido para el análisis de deformación y fractura" . Engineering Fracture Mechanics . 141 : 196–211 . doi : 10.1016/j.engfracmech.2015.04.036 .
  12. Chen, Zhiyong; Woody Ju, J.; Su, Guoshao; Huang, Xiaohua; Li, Shuang; Zhai, Lianjun (julio de 2019). "Influencia de las funciones de micromódulo en la simulación peridinámica de la propagación y ramificación de grietas en materiales frágiles" . Engineering Fracture Mechanics . 216 106498. doi : 10.1016/j.engfracmech.2019.106498 . S2CID 197621064 . 
  13. Ha, Youn Doh; Bobaru, Florin (marzo de 2010). "Estudios de propagación dinámica de grietas y ramificación de grietas con peridinámica" . International Journal of Fracture . 162 ( 1–2 ): 229–244 . doi : 10.1007/s10704-010-9442-4 . ISSN 0376-9429 . S2CID 8462707 .  
  14. Kilic, Bahattin (2008). "Teoría peridinámica para la predicción de fallas progresivas en materiales homogéneos y heterogéneos" .{{cite journal}}: Para citar una revista se requiere |journal=( ayuda )
  15. Agwai, Abigail; Guven, Ibrahim; Madenci, Erdogan (septiembre de 2011). "Predicción de la propagación de grietas con peridinámica: un estudio comparativo" . International Journal of Fracture . 171 (1): 65–78 . doi : 10.1007/s10704-011-9628-4 . ISSN 0376-9429 . S2CID 136475045 .  
  16. Lipton, Robert (octubre de 2014). "Fractura frágil dinámica como límite de horizonte pequeño de la peridinámica" . Journal of Elasticity . 117 (1): 21– 50. arXiv : 1305.4531 . doi : 10.1007/s10659-013-9463-0 . ISSN 0374-3535 . S2CID 254462294 .  
  17. Silling, SA; Weckner, O.; Askari, E.; Bobaru, F. (marzo de 2010). "Nucleación de grietas en un sólido peridinámico" . International Journal of Fracture . 162 ( 1–2 ): 219–227 . doi : 10.1007/s10704-010-9447-z . ISSN 0376-9429 . S2CID 209225 .  
  18. ^ Coclita, GM; Dipierro, S.; Fanizza, G.; Magdalena, F.; Romano, M.; Valdinoci, E. (marzo de 2023). "Aspectos cualitativos en la dinámica no local" . Revista de peridinámica y modelado no local . 5 (1): 1– 19. arXiv : 2106.13596 . doi : 10.1007/s42102-021-00064-z . ISSN 2522-896X . S2CID 235652235 .  
  19. Seleson, Pablo; Parks, Michael L.; Gunzburger, Max; Lehoucq, Richard B. (enero de 2009). "Peridynamics as an Upscaling of Molecular Dynamics" . Multiscale Modeling & Simulation . 8 (1): 204–227 . doi : 10.1137/09074807X . ISSN 1540-3459 . OSTI 1678881 .  
  20. Behzadinasab, Masoud; Foster, John T. (abril de 2020). "Un marco de correspondencia constitutiva semilagrangiana para la peridinámica" . Journal of the Mechanics and Physics of Solids . 137 103862. Bibcode : 2020JMPSo.13703862B . doi : 10.1016/j.jmps.2019.103862 . S2CID 212784700 . 
  21. Askari, E; Bobaru, F; Lehoucq, RB; Parks, ML; Silling, SA; Weckner, O (2008-07-01). "Peridynamics for multiscale materials modeling" . Journal of Physics: Conference Series . 125 (1) 012078. Bibcode : 2008JPhCS.125a2078A . doi : 10.1088/1742-6596/125/1/012078 . ISSN 1742-6596 . S2CID 250694017 .  
  22. Ni, Tao; Pesavento, Francesco; Zaccariotto, Mirco; Galvanetto, Ugo; Zhu, Qi-Zhi; Schrefler, Bernhard A. (julio de 2020). "FEM híbrido y simulación peridinámica de propagación de fracturas hidráulicas en medios porosos saturados" . Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 366 113101. arXiv : 2307.10929 . Código Bib : 2020CMAME.366k3101N . doi : 10.1016/j.cma.2020.113101 . S2CID 219519506 . 
  23. Zhou, Xiao-Ping; Wang, Yun-Teng; Shou, Yun-Dong (agosto de 2020). "Modelo peridinámico basado en enlaces hidromecánicos para procesos de fracturación impulsados ​​por fluidos y a presión en rocas porosas fisuradas" . International Journal of Rock Mechanics and Mining Sciences . 132 104383. Bibcode : 2020IJRMM.13204383Z . doi : 10.1016/j.ijrmms.2020.104383 . S2CID 225382857 . 
  24. Song, Xiaoyu; Khalili, Nasser (enero de 2019). "Un modelo peridinámico para el análisis de localización de deformaciones de geomateriales" . Revista Internacional de Métodos Numéricos y Analíticos en Geomecánica . 43 (1): 77– 96. Bibcode : 2019IJNAM..43...77S . doi : 10.1002/nag.2854 . ISSN 0363-9061 . S2CID 125649306 .  
  25. Panchadhara, Rohan; Gordon, Peter A.; Parks, Michael L. (marzo de 2017). "Modelado de la estimulación de un pozo mediante propelente con peridinámica" . International Journal of Rock Mechanics and Mining Sciences . 93 : 330–343 . Bibcode : 2017IJRMM..93..330P . doi : 10.1016/j.ijrmms.2017.02.006 .
  26. Zhou, Xiao-Ping; Wang, Yun-Teng (enero de 2021). "Revisión del estado del arte sobre las características de falla progresiva de los geomateriales en la teoría peridinámica" . Journal of Engineering Mechanics . 147 (1) 03120001. doi : 10.1061/(ASCE)EM.1943-7889.0001876 . ISSN 0733-9399 . S2CID 228906748 .  
  27. Lejeune, Emma; Linder, Christian (agosto de 2017). "Modelado del crecimiento tumoral con peridinámica" . Biomechanics and Modeling in Mechanobiology . 16 (4): 1141– 1157. doi : 10.1007/s10237-017-0876-8 . ISSN 1617-7959 . PMID 28124191. S2CID 254169636 .   
  28. Taylor, Michael; Gözen, Irep; Patel, Samir; Jesorka, Aldo; Bertoldi, Katia (2016-11-09). van Veen, Hendrik W. (ed.). "Modelado peridinámico de rupturas en biomembranas" . PLOS ONE . 11 (11) e0165947. Bibcode : 2016PLoSO..1165947T . doi : 10.1371/journal.pone.0165947 . ISSN 1932-6203 . PMC 5102442. PMID 27829001 .   
  29. Bobaru, Florin; Duangpanya, Monchai (abril de 2012). "Una formulación peridinámica para la conducción de calor transitoria en cuerpos con discontinuidades en evolución" . Journal of Computational Physics . 231 (7): 2764– 2785. Bibcode : 2012JCoPh.231.2764B . doi : 10.1016/j.jcp.2011.12.017 . S2CID 6929467 . 
  30. Oterkus, Selda; Madenci, Erdogan; Agwai, Abigail (mayo de 2014). "Difusión térmica peridinámica" . Journal of Computational Physics . 265 : 71–96 . Bibcode : 2014JCoPh.265...71O . doi : 10.1016/j.jcp.2014.01.027 . S2CID 22835224 . 
  31. Foster, John (2019). "Métodos no locales y de orden fraccional para la turbulencia cerca de la pared, simulación de grandes remolinos e interacción fluido-estructura" . Informe técnico, Universidad de Texas en Austin, Austin, Estados Unidos .
  32. Zhao, Jiangming; Chen, Ziguang; Mehrmashhadi, Javad; Bobaru , Florin (noviembre de 2018). "Construcción de un modelo peridinámico para problemas transitorios de advección-difusión" . International Journal of Heat and Mass Transfer . 126 : 1253–1266 . Bibcode : 2018IJHMT.126.1253Z . doi : 10.1016/j.ijheatmasstransfer.2018.06.075 . S2CID 125321481 . 
  33. Buryachenko, Valeriy A. (octubre de 2020). "Método generalizado de campos efectivos en la micromecánica peridinámica de compuestos de estructura aleatoria" . International Journal of Solids and Structures . 202 : 765–786 . Bibcode : 2020IJSS..202..765B . doi : 10.1016/j.ijsolstr.2020.06.022 . S2CID 225577923 . 
  34. Hu, YL; Madenci, E. (enero de 2017). "Peridinámica para la predicción de la vida a fatiga y la resistencia residual de laminados compuestos" . Composite Structures . 160 : 169–184 . doi : 10.1016/j.compstruct.2016.10.010 .
  35. Oterkus, Erkan; Madenci, Erdogan (28-03-2012). "Análisis peridinámico de materiales compuestos reforzados con fibras" . Journal of Mechanics of Materials and Structures . 7 (1): 45– 84. doi : 10.2140/jomms.2012.7.45 . ISSN 1559-3959 . 
  36. Zhao, Jiangming; Jafarzadeh, Siavash; Rahmani, Mohammad; Chen, Ziguang; Kim, Yong-Rak; Bobaru, Florin (septiembre de 2021). "Un modelo peridinámico para la corrosión galvánica y la fractura" . Electrochimica Acta . 391 138968. doi : 10.1016/j.electacta.2021.138968 .

Lecturas adicionales

  • Bobaru, Florin; Foster, John T.; Geubelle, Philippe H.; Silling, Stewart A., eds. (2016). Manual de modelado peridinámico . Avances en matemáticas aplicadas. Boca Raton, Londres, Nueva York: CRC Press, Taylor & Francis Group, un libro de Chapman & Hall. ISBN 978-1-4822-3044-4.
  • Oterkus, Erkan; Oterkus, Selda; Madenci, Erdogan (24 de abril de 2021). Modelado peridinámico, técnicas numéricas y aplicaciones . Elsevier. ISBN 978-0-12-820441-2.
  • Rabczuk, Timon; Ren, Huilong; Zhuang, Xiaoying (15 de febrero de 2023). Métodos computacionales basados ​​en peridinámica y operadores no locales: teoría y aplicaciones . Springer Nature. ISBN 978-3-031-20906-2.
  • D'Elia, Marta; Li, Xingjie; Seleson, Pablo; Tian, ​​Xiaochuan; Yu, Yue (marzo de 2022). " Una revisión de los métodos de acoplamiento local a no local en difusión no local y mecánica no local" . Journal of Peridynamics and Nonlocal Modeling . 4 (1): 1– 50. arXiv : 1912.06668 . doi : 10.1007/s42102-020-00038-7 . ISSN 2522-896X . OSTI 1833985. S2CID 257114051 .   
  • Bobaru, Florin; Chen, Ziguang; Jafarzadeh, Siavash (1 de diciembre de 2023). Daños por corrosión y fractura asistida por corrosión: modelado y cálculos peridinámicos . Elsevier. ISBN 978-0-12-823174-6.
  • PeriLab es un repositorio de código abierto para peridinámica escrito en Julia.
  • Peridigm, un código de peridinámica computacional de código abierto
  • PeriDoX es un repositorio de código abierto para la peridinámica y su documentación.
  • Laboratorio Sandia - Peridinámica
  • Sitio web sobre peridinámica
  • Implementación de la aproximación de elementos finitos y diferencias finitas de modelos no locales.