Home Analysis and Numerical Simulation of Time-Fractional Derivative Contact Problem with Friction in Thermo-Viscoelasticity
Article
Licensed
Unlicensed Requires Authentication

Analysis and Numerical Simulation of Time-Fractional Derivative Contact Problem with Friction in Thermo-Viscoelasticity

  • Mustapha Bouallala ORCID logo EMAIL logo , EL-Hassan Essoufi and Youssef Ouafik
Published/Copyright: March 26, 2024

Abstract

The objective of this study is to analyze a quasistatic frictional contact problem involving the interaction between a thermo-viscoelastic body and a thermally conductive foundation. The constitutive relation in our investigation is constructed using a fractional Kelvin–Voigt model to describe displacement behavior. Additionally, the heat conduction aspect is governed by a time-fractional derivative parameter that is associated with temperature. The contact is modeled using the Signorini condition, which is a version of Coulomb’s law for dry friction. We develop a variational formulation for the problem and establish the existence of its weak solution using a combination of techniques, including the theory of monotone operators, Caputo derivative, Galerkin method, and the Banach fixed point theorem. To demonstrate the effectiveness of our approach, we include several numerical simulations that showcase the performance of the method.

A Appendix

In this section, we recall some results about the fractional calculus and nonlinear analysis that can find it as soon as the following references [17, 18, 12, 23, 16]. From a numerical perspective, the time discretization of Caputo derivative can be found in [35].

Definition A.1 (Riemann–Liouville Fractional Integral).

Let X be a Banach space and let ( 0 , T ) be a finite time interval. The Riemann–Liouville fractional integral of order α > 0 for a given function f L 1 ( 0 , T ; X ) is defined by

I t α 0 f ( t ) = 1 Γ ( α ) 0 t ( t - s ) α - 1 f ( s ) 𝑑 s for all  t ( 0 , T ) ,

where Γ ( ) stands for the Gamma function defined by

Γ ( α ) = 0 t α - 1 e - t 𝑑 t .

To complement the definition, we set I t 0 0 = I , where I is the identity operator, which means that I t 0 0 f ( t ) = f ( t ) for a.e. t ( 0 , T ) .

Definition A.2 (Caputo Derivative of Order 0 < α 1 ).

Let X be a Banach space, let 0 < α 1 and let ( 0 , T ) be a finite time interval. For a given function f A C ( 0 , T ; W ) , the Caputo fractional derivative of f is defined by

D t α 0 C f ( t ) = 0 I t 1 - α f ( t ) = 1 Γ ( 1 - α ) 0 t ( t - s ) - α f ( s ) 𝑑 s for all  t ( 0 , T ) .

The notation A C ( 0 , T ; X ) refers to the space of all absolutely continuous functions from ( 0 , T ) into X.

It is obvious that if α = 1 , the Caputo derivative reduces to the classical first-order derivative, that is, we have

D t 1 0 C f ( t ) = I f ( t ) = f ( t ) for a.e.  t ( 0 , T ) .

Proposition A.1.

Let X be a Banach space and α , β > 0 . Then the following statements hold:

  1. For y L 1 ( 0 , T ; X ) , we have

    I t α 0 I t β 0 y ( t ) = 0 I t α + β y ( t ) for a.e.  t ( 0 , T ) .

  2. For y A C ( 0 , T ; X ) and α ( 0 , α ] , we have

    I t α 0 D t α 0 C y ( t ) = y ( t ) - y ( 0 ) for a.e.  t ( 0 , T ) ,

  3. For y L 1 ( 0 , T ; X ) , we have D t α 0 C I t α 0 y ( t ) = y ( t ) for a.e. t ( 0 , T ) .

Theorem A.1.

Let T > 0 , γ ( 0 , 1 ) and p [ 1 , ) . Let M, B, Y be Banach spaces such that M B compactly and B Y continuously. Suppose W L loc 1 ( ( 0 , T ) ; M ) satisfies the following:

  1. There exist r 1 [ 1 , ) and C 1 > 0 such that for all u W ,

    (A.1) sup t ( 0 , T ) J γ ( u M r 1 ) = sup t ( 0 , T ) 1 Γ ( γ ) 0 t ( t - s ) γ - 1 u M r 1 ( s ) 𝑑 s C 1 ,

    where

    J α u ( t ) = 1 Γ ( α ) 0 t ( t - s ) α - 1 u ( s ) 𝑑 s .

  2. There exists p 1 ( p , ] such that W is bounded in L p 1 ( ( 0 , T ) ; B ) .

  3. There exist r 2 [ 1 , ) and C 2 > 0 such that for all u W , there is an assignment of initial value u 0 so that the weak Caputo derivative satisfies

    (A.2) D c γ u L r 2 ( ( 0 , T ) ; Y ) C 2 .

    Then W is relatively compact in L p ( ( 0 , T ) ; B ) .

Proposition A.2.

Let - α . Assume f : [ 0 , ) × ( α , β ) R is continuous and locally Lipschitz in the second variable. In other words, for all A > 0 and K ( α , β ) compact, there exists L A , K > 0 such that

(A.3) sup 0 t A | f ( t , v 1 ) - f ( t , v 2 ) | L A , K | v 1 - v 2 | for all  v 1 , v 2 K .

Let 0 < γ < 1 and v 0 ( α , β ) . Then the IVP

(A.4) D c γ v ( t ) = f ( t , v ( t ) ) , v ( 0 ) = v 0

has a unique continuous solution v ( ) on [ 0 , T b ) , where

(A.5) T b = sup { h > 0 : the solution  v C [ 0 , h ) , v ( t ) ( α , β ) for all  t [ 0 , h ) }

is the largest time of existence satisfies T b ( 0 , ] . If T b < , then lim inf t T b - v ( t ) = α or lim sup t T b - v ( t ) = β .

For the time discretization of the Caputo derivative in (2.1) and (2.4), letting t n = n k , n = 0 , 1 , , N , and k = T N being the time step size, we have the following formulation of Caputo derivative (see [35]):

D t α 0 C f ( t n + 1 ) = 1 Γ ( 1 - α ) j = 0 n t j t j + 1 ( t n + 1 - s ) - α f ( s ) 𝑑 s
= 1 Γ ( 1 - α ) j = 0 n f ( t j + 1 ) - f ( t j ) k t j t j + 1 ( t n + 1 - s ) - α 𝑑 s + r k n + 1
= 1 Γ ( 1 - α ) j = 0 n f ( t j + 1 ) - f ( t j ) k k 1 - α 1 - α ( ( n + 1 - j ) 1 - α - ( n - j ) 1 - α ) + r k n + 1
= 1 Γ ( 2 - α ) j = 0 n f ( t n + 1 - j ) - f ( t n - j ) k α d j + r k n + 1 ,

where d j = ( j + 1 ) 1 - α - j 1 - α , and r k n + 1 is the truncated error.

Let us define the discrete fractional differential operator L t α by

L t α f ( x , t n + 1 ) = 1 Γ ( 2 - α ) j = 0 n d j f ( x , t n + 1 - j ) - f ( x , t n - j ) k α .

So we have

(A.6) D t α 0 C f ( x , t n + 1 ) = L t α f ( x , t n + 1 ) + r k n + 1 .

References

[1] M. Allen, L. Caffarelli and A. Vasseur, A parabolic problem with a fractional time derivative, Arch. Ration. Mech. Anal. 221 (2016), no. 2, 603–630. 10.1007/s00205-016-0969-zSearch in Google Scholar

[2] M. Barboteu and M. Sofonea, Modeling and analysis of the unilateral contact of a piezoelectric body with a conductive support, J. Math. Anal. Appl. 358 (2009), no. 1, 110–124. 10.1016/j.jmaa.2009.04.030Search in Google Scholar

[3] K. Bartosz, D. Danan and P. Szafraniec, Numerical analysis of a dynamic bilateral thermoviscoelastic contact problem with nonmonotone friction law, Comput. Math. Appl. 73 (2017), no. 5, 727–746. 10.1016/j.camwa.2016.12.026Search in Google Scholar

[4] M. Bouallala and E.-H. Essoufi, A thermo-viscoelastic fractional contact problem with normal compliance and Coulomb’s friction, J. Math. Phys. Anal. Geom. 17 (2021), no. 3, 280–294. 10.15407/mag17.03.280Search in Google Scholar

[5] M. Campo and J. R. Fernández, Numerical analysis of a quasistatic thermoviscoelastic frictional contact problem, Comput. Mech. 35 (2005), no. 6, 459–469. 10.1007/s00466-004-0635-4Search in Google Scholar

[6] A. Curnier and P. Alart, A generalized Newton method for contact problems with friction, J. Mech. Theory. Appl. 7 (1988), no. 1, 67–82. Search in Google Scholar

[7] K. Diethelm, The Analysis of Fractional Differential Equations, Lecture Notes in Math. 2004, Springer, Berlin, 2010. 10.1007/978-3-642-14574-2Search in Google Scholar

[8] L. C. Evans, Partial Differential Equations, 2nd ed., Grad. Stud. Math. 19, American Mathematical Society, Providence, 2010. Search in Google Scholar

[9] J. Han, S. Migórski and H. Zeng, Weak solvability of a fractional viscoelastic frictionless contact problem, Appl. Math. Comput. 303 (2017), 1–18. 10.1016/j.amc.2017.01.009Search in Google Scholar

[10] R. Herrmann, Fractional Calculus: An Introduction for Physicists, World Scientific, Singapore, 2011. 10.1142/9789814340250Search in Google Scholar

[11] B. Jin, Fractional Differential Equations, Springer, Cham, 2021. 10.1007/978-3-030-76043-4Search in Google Scholar

[12] A. A. Kilbas, H. M. Srivastava and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland Math. Stud. 204, Elsevier Science, Amsterdam, 2006. Search in Google Scholar

[13] R. C. Koeller, Applications of fractional calculus to the theory of viscoelasticity, Trans. ASME J. Appl. Mech. 51 (1984), no. 2, 299–307. 10.1115/1.3167616Search in Google Scholar

[14] A. Kubica, K. Ryszewska and M. Yamamoto, Time-Fractional Differential Equations—A Theoretical Introduction, Springer Briefs in Math., Springer, Singapore, 2020. 10.1007/978-981-15-9066-5Search in Google Scholar

[15] G. F. A. L’Hôpital, Analyse des infiniment petits pour l’intelligence des lignes courbes, 2nd ed., François Montalant, Paris, 1715. Search in Google Scholar

[16] C. Li and M. Cai, Theory and Numerical Approximations of Fractional Integrals and Derivatives, Society for Industrial and Applied Mathematics, Philadelphia, 2020. 10.1137/1.9781611975888Search in Google Scholar

[17] L. Li and J.-G. Liu, A generalized definition of Caputo derivatives and its application to fractional ODEs, SIAM J. Math. Anal. 50 (2018), no. 3, 2867–2900. 10.1137/17M1160318Search in Google Scholar

[18] L. Li and J.-G. Liu, Some compactness criteria for weak solutions of time fractional PDEs, SIAM J. Math. Anal. 50 (2018), no. 4, 3963–3995. 10.1137/17M1145549Search in Google Scholar

[19] Q. Ma, X. Lin and H. Li, A comparative teaching of fractional calculus and integer calculus, 2018 International Conference on Social Science and Education Reform (ICSSER 2018), Atlantis Press, Amsterdam (2018), 1–4. 10.2991/icsser-18.2018.22Search in Google Scholar

[20] P. G. Nutting, A general stress-strain time formula, J. Franklin Inst. 285 (1943), 513–524. 10.1016/S0016-0032(43)91483-8Search in Google Scholar

[21] Y. Ouafik, Numerical analysis of a frictional contact problem for thermo-electro-elastic materials, J. Theoret. Appl. Mech. 58 (2020), no. 3, 673–683. 10.15632/jtam-pl/122258Search in Google Scholar

[22] Y. Ouafik, Numerical treatment of a thermoelastic contact problem, Eng. Appl. Sci. Lett. 6 (2023), no. 1, 11–18. Search in Google Scholar

[23] I. Podlubny, Fractional Differential Equations, Math. Sci. Eng. 198, Academic Press, San Diego, 1998. Search in Google Scholar

[24] J. Pommier and Y. Renard, Getfem++, an open source generic C++ library for finite element methods, http://getfem.org␣(accessed␣22.08.2023). Search in Google Scholar

[25] Y. Renard, Generalized Newton’s methods for the approximation and resolution of frictional contact problems in elasticity, Comput. Methods Appl. Mech. Engrg. 256 (2013), 38–55. 10.1016/j.cma.2012.12.008Search in Google Scholar

[26] B. Riemann, Versuch einer allgemeinen Auffassung der Integration und Differentiation (1847), Bernard Riemann’s Gesammelte Mathematische Werke und Wissenschaftlicher Nachlass (Cambridge Library Collection—Mathematics), Cambridge University, Cambridge (2014), 331–344. 10.1017/CBO9781139568050.020Search in Google Scholar

[27] P. E. Rouse, The theory of the linear viscoleastic properties of dilute solutions of coiling polymers, J. Chem. Phys. 21 (1953), 1272–1280. 10.1063/1.1699180Search in Google Scholar

[28] A. Soleiman, A. E. Abouelregal, K. M. Khalil and M. E. Nasr, Generalized thermoviscoelastic novel model with different fractional derivatives and multi-phase-lags, European Phys. J. Plus 135 (2020), no. 10, Paper No. 851. 10.1140/epjp/s13360-020-00842-6Search in Google Scholar

[29] I. Ubero-Martínez, L. Rodríguez-Tembleque, J. Cifuentes-Rodríguez and J. Vallepuga-Espinosa, 3D thermoelastic solids under non-linear interface thermal and orthotropic frictional contact conditions, Internat. J. Numer. Methods Engrg. 123 (2022), no. 11, 2631–2659. 10.1002/nme.6953Search in Google Scholar

[30] H. Xuan and X. Cheng, Numerical analysis of a thermal frictional contact problem with long memory, Commun. Pure Appl. Anal. 20 (2021), no. 4, 1521–1543. 10.3934/cpaa.2021031Search in Google Scholar

[31] Q. Yang, I. Turner, F. Liu and M. Ilić, Novel numerical methods for solving the time-space fractional diffusion equation in two dimensions, SIAM J. Sci. Comput. 33 (2011), no. 3, 1159–1180. 10.1137/100800634Search in Google Scholar

[32] F. Zeng, C. Li, F. Liu and I. Turner, Numerical algorithms for time-fractional subdiffusion equation with second-order accuracy, SIAM J. Sci. Comput. 37 (2015), no. 1, 55–78. 10.1137/14096390XSearch in Google Scholar

[33] S. Zeng, Z. Liu and S. Migorski, A class of fractional differential hemivariational inequalities with application to contact problem, Z. Angew. Math. Phys. 69 (2018), no. 2, Paper No. 36. 10.1007/s00033-018-0929-6Search in Google Scholar

[34] S. Zeng and S. Migórski, A class of time-fractional hemivariational inequalities with application to frictional contact problem, Commun. Nonlinear Sci. Numer. Simul. 56 (2018), 34–48. 10.1016/j.cnsns.2017.07.016Search in Google Scholar

[35] T. Zhao and H. Xuan, Error analysis of spectral method for the space and time fractional Fokker–Planck equation, Recent Advances in Applied Nonlinear Dynamics with Numerical Analysis, Interdiscip. Math. Sci. 15, World Scientific, Hackensack (2013), 83–104. 10.1142/9789814436465_0005Search in Google Scholar

Received: 2023-09-02
Revised: 2024-02-21
Accepted: 2024-03-03
Published Online: 2024-03-26
Published in Print: 2025-01-01

© 2024 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 13.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/cmam-2023-0192/html
Scroll to top button