# American Institute of Mathematical Sciences

December  2017, 10(6): 1539-1561. doi: 10.3934/dcdss.2017079

## An energy based formulation of a quasi-static interface damage model with a multilinear cohesive law

 1 Technical University of Košice, Civil Engineering Faculty, Vysokoškolská 9,042 00 Košice, Slovakia 2 University of Seville, School of Engineering, Camino de los Descubrimientos s/n, 41092 Seville, Spain

Dedicated to Tomáš Roubíček on the occasion of his 60th birthday

Received  October 2016 Revised  February 2017 Published  June 2017

A new quasi-static and energy based formulation of an interface damage model which provides interface traction-relative displacement laws like in traditional trilinear (with bilinear softening) or generally multilinear cohesive zone models frequently used by engineers is presented. This cohesive type response of the interface may represent the behaviour of a thin adhesive layer. The level of interface adhesion or damage is defined by several scalar variables suitably defined to obtain the required traction-relative displacement laws. The weak solution of the problem is sought numerically by a semi-implicit time-stepping procedure which uses recursive double minimization in displacements and damage variables separately. The symmetric Galerkin boundary-element method is applied for the spatial discretization. Sequential quadratic programming is implemented to resolve each partial minimization in the recursive scheme applied to compute the time-space discretized solutions. Sample 2D numerical examples demonstrate applicability of the proposed model.

Citation: Roman VodiČka, Vladislav MantiČ. An energy based formulation of a quasi-static interface damage model with a multilinear cohesive law. Discrete & Continuous Dynamical Systems - S, 2017, 10 (6) : 1539-1561. doi: 10.3934/dcdss.2017079
##### References:

show all references

##### References:
The used notation for two bonded domains.
Stress-displacement curves for (a) the bilinear and (b) multilinear CZMs with $m_{\rm d}{=}3$.
Stress-displacement relations for multilinear CZMs with $m_{\rm d}{=}3$: (a) Mode 'subsequent', (b) Mode 'at once'.
An example of stress-relative displacement relation for a multilinear CZM with $m_{\rm d}{=}3$, $\sigma_{\max\,{\rm n}}{=}1.25\sigma_{\max\,{\rm t}}$, $G_{\tiny\rm IIc}{=}2G_{\tiny\rm Ic}$.
Simple tension in the two-square example, $a_1{=}200\,$mm: (a) the problem layout, (b) the traction-relative displacement law in the cohesive zone, (c) the loading function $g$ from (24).
The stress-displacement relation at the point $x_2{=}50$mm (the quarter of the interface) and the evolution of the damage parameters at the same point: (a) Mode 'subsequent', (b) Mode 'at once'.
Double cantilever beam: (a) the problem layout: $\ell{=}190\,$mm, $\ell_{\text{ini}}{=}55\,$mm, $w{=}20\,$mm, $h{=}5\,$mm, (b) the traction-relative displacement law in the cohesive zone: $u_0{=}0.014\,$mm, $u_1{=}0.25\,$mm, $u_{\text c}{=}4\,$mm, $\sigma_0{=}62\,$MPa, $\sigma_1{=}0.67\,$MPa.
Deformations of DCB, the damage evolution and normal stress distribution in the partially cracked interface at selected time instants corresponding to prescribed displacement $g$, $\ell_{\rm ini}$ is the initial crack length, cf. Figure 7.
Normal stress-relative displacement graphs at the interface point $x_1{=}\ell_{\text{ini}}{+}4$mm. The damage range is kept the same in the right and left part, only the range for the normal stress $\sigma_{\rm n}$ is changed in the right picture.
The applied forces for DCB calculated at the place where the vertical displacement is imposed.
The mixed mode beam, cf. [40]; (a) the problem layout: $\ell{=}120$mm, $\ell_{\text C}{=}92$mm, $\ell_{\text{ini}}{=}8$mm, $h{=}20$mm, $w{=}2$mm, $s{=}2$mm, (b) used stress-displacement law in the cohesive zone: $u_1{=}2u_0$, $u_2{=}3u_0$, $u_{\text c}{=}4u_0$, $\sigma_2{=}\frac18\sigma_0$ and $\sigma_0{=}7.5\,$MPa, where $u_0{=}0.01\,$mm in normal component; in the tangential component either the same value (no dependence on mode-mixity) or $u_0{=}0.04\,$mm (mixed-mode dependent: $G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$).
The total reaction forces for the mixed-mode beam calculated at the simple support constraint: (a) observing the influences of the fracture mode mixity ($G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$ or $G_{\tiny\rm IIc}{=}G_{\tiny\rm Ic}$) and viscosity (solid lines for no viscosity ${\tau_{\rm r}}{=}0$, dashed lines for ${\tau_{\rm r}}{=}10$ms), (b) changes of the solution for various discretizations.
The energy evolution and fulfillment of the energy balance (7) for various discretizations calculated for the mixed-mode beam, $G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$: (a) no viscosity ${\tau_{\rm r}}{=}0$, (b) viscosity with ${\tau_{\rm r}}{=}10$ms.
Deformations of the beam, the damage evolution and stress distribution in the cracked interface at selected time instants corresponding to the prescribed displacement $g$: comparison of the cases $G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$ or $G_{\tiny\rm IIc}{=}G_{\tiny\rm Ic}$ (referenced respectively by indices 4 and 1) with no viscosity.
Deformations of the beam, the damage evolution and stress distribution in the cracked interface at selected time instants corresponding to the prescribed displacement $g$: comparison of the cases ${\tau_{\rm r}}{=}10$ms and ${\tau_{\rm r}}{=}0$ (referenced respectively by indices 10 and 0) with $G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$.
Stress-relative displacement graphs at selected points of the interface: ${\tau_{\rm r}}{=}0$, $G_{\tiny\rm IIc}{=}4G_{\tiny\rm Ic}$.
Stress-relative displacement graphs at selected points of the interface: ${\tau_{\rm r}}{=}0$, $G_{\tiny\rm IIc}{=}G_{\tiny\rm Ic}$.
 [1] Marita Thomas, Chiara Zanini. Cohesive zone-type delamination in visco-elasticity. Discrete & Continuous Dynamical Systems - S, 2017, 10 (6) : 1487-1517. doi: 10.3934/dcdss.2017077 [2] G. Leugering, Marina Prechtel, Paul Steinmann, Michael Stingl. A cohesive crack propagation model: Mathematical theory and numerical solution. Communications on Pure & Applied Analysis, 2013, 12 (4) : 1705-1729. doi: 10.3934/cpaa.2013.12.1705 [3] Anna Marciniak-Czochra, Andro Mikelić. A nonlinear effective slip interface law for transport phenomena between a fracture flow and a porous medium. Discrete & Continuous Dynamical Systems - S, 2014, 7 (5) : 1065-1077. doi: 10.3934/dcdss.2014.7.1065 [4] Tomáš Roubíček, V. Mantič, C. G. Panagiotopoulos. A quasistatic mixed-mode delamination model. Discrete & Continuous Dynamical Systems - S, 2013, 6 (2) : 591-610. doi: 10.3934/dcdss.2013.6.591 [5] Songqiang Qiu, Zhongwen Chen. An adaptively regularized sequential quadratic programming method for equality constrained optimization. Journal of Industrial & Management Optimization, 2017, 13 (5) : 1-14. doi: 10.3934/jimo.2019075 [6] Tuan Hiep Pham, Jérôme Laverne, Jean-Jacques Marigo. Stress gradient effects on the nucleation and propagation of cohesive cracks. Discrete & Continuous Dynamical Systems - S, 2016, 9 (2) : 557-584. doi: 10.3934/dcdss.2016012 [7] Gianni Dal Maso, Flaviana Iurlano. Fracture models as $\Gamma$-limits of damage models. Communications on Pure & Applied Analysis, 2013, 12 (4) : 1657-1686. doi: 10.3934/cpaa.2013.12.1657 [8] Ye Tian, Cheng Lu. Nonconvex quadratic reformulations and solvable conditions for mixed integer quadratic programming problems. Journal of Industrial & Management Optimization, 2011, 7 (4) : 1027-1039. doi: 10.3934/jimo.2011.7.1027 [9] Tien-Tsan Shieh. From gradient theory of phase transition to a generalized minimal interface problem with a contact energy. Discrete & Continuous Dynamical Systems - A, 2016, 36 (5) : 2729-2755. doi: 10.3934/dcds.2016.36.2729 [10] Christopher J. Larsen. Local minimality and crack prediction in quasi-static Griffith fracture evolution. Discrete & Continuous Dynamical Systems - S, 2013, 6 (1) : 121-129. doi: 10.3934/dcdss.2013.6.121 [11] Tao Lin, Yanping Lin, Weiwei Sun. Error estimation of a class of quadratic immersed finite element methods for elliptic interface problems. Discrete & Continuous Dynamical Systems - B, 2007, 7 (4) : 807-823. doi: 10.3934/dcdsb.2007.7.807 [12] Michael Stiassnie, Raphael Stuhlmeier. Progressive waves on a blunt interface. Discrete & Continuous Dynamical Systems - A, 2014, 34 (8) : 3171-3182. doi: 10.3934/dcds.2014.34.3171 [13] Ben A. Vanderlei, Matthew M. Hopkins, Lisa J. Fauci. Error estimation for immersed interface solutions. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1185-1203. doi: 10.3934/dcdsb.2012.17.1185 [14] Frédéric Lebon, Raffaella Rizzoni. Modeling a hard, thin curvilinear interface. Discrete & Continuous Dynamical Systems - S, 2013, 6 (6) : 1569-1586. doi: 10.3934/dcdss.2013.6.1569 [15] Amadeu Delshams, Marina Gonchenko, Sergey V. Gonchenko, J. Tomás Lázaro. Mixed dynamics of 2-dimensional reversible maps with a symmetric couple of quadratic homoclinic tangencies. Discrete & Continuous Dynamical Systems - A, 2018, 38 (9) : 4483-4507. doi: 10.3934/dcds.2018196 [16] Champike Attanayake, So-Hsiang Chou. An immersed interface method for Pennes bioheat transfer equation. Discrete & Continuous Dynamical Systems - B, 2015, 20 (2) : 323-337. doi: 10.3934/dcdsb.2015.20.323 [17] Zhongyi Huang. Tailored finite point method for the interface problem. Networks & Heterogeneous Media, 2009, 4 (1) : 91-106. doi: 10.3934/nhm.2009.4.91 [18] Jian Hao, Zhilin Li, Sharon R. Lubkin. An augmented immersed interface method for moving structures with mass. Discrete & Continuous Dynamical Systems - B, 2012, 17 (4) : 1175-1184. doi: 10.3934/dcdsb.2012.17.1175 [19] Qiang Du, Manlin Li. On the stochastic immersed boundary method with an implicit interface formulation. Discrete & Continuous Dynamical Systems - B, 2011, 15 (2) : 373-389. doi: 10.3934/dcdsb.2011.15.373 [20] Thomas I. Seidman. Interface conditions for a singular reaction-diffusion system. Discrete & Continuous Dynamical Systems - S, 2009, 2 (3) : 631-643. doi: 10.3934/dcdss.2009.2.631

2018 Impact Factor: 0.545