Next Article in Journal
Threat Hunting System for Protecting Critical Infrastructures Using a Machine Learning Approach
Next Article in Special Issue
Some New Estimates of Hermite–Hadamard Inequalities for Harmonical cr-h-Convex Functions via Generalized Fractional Integral Operator on Set-Valued Mappings
Previous Article in Journal
Short-Term Prediction of Time-Varying Passenger Flow for Intercity High-Speed Railways: A Neural Network Model Based on Multi-Source Data
Previous Article in Special Issue
Study of Log Convex Mappings in Fuzzy Aunnam Calculus via Fuzzy Inclusion Relation over Fuzzy-Number Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variational Approach to Modeling of Curvilinear Thin Inclusions with Rough Boundaries in Elastic Bodies: Case of a Rod-Type Inclusion

Lavrentyev Institute of Hydrodynamics of SB RAS, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(16), 3447; https://doi.org/10.3390/math11163447
Submission received: 6 July 2023 / Revised: 3 August 2023 / Accepted: 7 August 2023 / Published: 8 August 2023
(This article belongs to the Special Issue Variational Problems and Applications, 2nd Edition)

Abstract

:
In the framework of 2D-elasticity, an equilibrium problem for an inhomogeneous body with a curvilinear inclusion located strictly inside the body is considered. The elastic properties of the inclusion are assumed to depend on a small positive parameter δ characterizing its width and are assumed to be proportional to δ 1 . Moreover, it is supposed that the inclusion has a curvilinear rough boundary. Relying on the variational formulation of the equilibrium problem, we perform the asymptotic analysis, as δ tends to zero. As a result, a variational model of an elastic body containing a thin curvilinear rod is constructed. Numerical calculations give a relative error between the initial and limit problems depending on δ .

1. Introduction

The present paper deals with an equilibrium problem for an inhomogeneous body with a curvilinear inclusion located strictly inside the elastic body. The elastic properties of the inclusion are assumed to depend on a small positive parameter δ characterizing its width. More precisely, the Lamé coefficients are proportional to δ 1 . Furthermore, it is supposed that the inclusion has a perturbed curvilinear boundary. This means that the boundary of the inclusion is rough. The main purpose of this paper is to pass to the limit as δ tends to zero. As a result, a model of an elastic body containing a thin elastic curvilinear inclusion inside is constructed (in the considered case, this is a rod-type inclusion). To derive the limit problem, we use an approach based on a variational formulation of the equilibrium problem (see, e.g., [1,2,3]). It is worth noting that in our consideration, we take into account the three key features, which, in our opinion, were taken into account simultaneously earlier only in [1], in the case of a soft inclusion. Namely, these three features are that the inclusion is located strictly inside the body, the inclusion has a curvilinear boundary, and the boundary of the inclusion is rough.
Imperfect interface problems and problems of the coupling of different models of elasticity are of the greatest interest in the description of composite materials. In the present time, various types of thin inclusion problems were considered. Papers [4,5] investigated thin elastic inclusions such as the Euler–Bernoulli beam and the Timoshenko beam. Thin rigid inclusions were considered in [6,7,8,9,10]. Problems of thin elastic junctions between elastic bodies were studied in [11,12,13]. In all the mentioned works, as a rule, the main attention was paid to the study of the well-posedness of variational (weak) formulations of corresponding boundary value problems for equilibrium models but not the derivation of such models.
Currently, there is a huge number of papers devoted to the asymptotic analysis of linear and non-linear models of elasticity with the aim to derive different types of imperfect interfaces between elastic bodies. For example, in [14], a higher-order imperfect interface was studied by using complex variables and asymptotic analysis. In [15,16,17], interface models were derived for coupled thermoelasticity (see also [18,19,20,21] for modeling soft and hard interfaces in 2D- and 3D-elasticity).
One can notice that there are not so many papers devoted to studies of curvilinear interfaces with rough boundaries. In this line, we can only mention that the curvilinear interface problems were considered in [1,18,22], while in [1,23,24,25,26], the roughness of the inclusions was taken into account.
Moreover, we mention refs. [25,27,28,29], in which multi-scale asymptotic analysis and computation of static and dynamics problems in curvilinear coordinates was proposed.
The rest of the article is organized as follows. In Section 2, the equilibrium problem for an elastic body with a narrow inclusion depending on a small parameter δ > 0 is formulated. Then, in Section 3, we introduce the special coordinate transformations, which allow us further in Section 4 to reformulate the initial problem and decompose it in such a way that the new variational problem is defined in a fixed domain that is independent of δ . Section 5 is devoted to the justification of the asymptotic procedure, as a result of which the limit problem is derived. The numerical examples of Section 6 confirm the accuracy of the passage to the limit, justified in Section 5.

2. Statement of the Initial Problem

Let Ω R 2 be a convex bounded domain with a Lipschitz boundary Ω ; Γ N and Γ D be the parts of Ω such that Γ ¯ N Γ ¯ D = Ω and meas Γ D > 0 ; and I, I 1 , and I 2 be the three intervals lying on the abscissa axis O y 1 such that I ¯ 1 I 2 , I ¯ 2 I . Moreover, assume I is the intersection of the domain Ω with O y 1 -axis.
We introduce the functions φ C 1 , 1 ( I ¯ ) , ψ , ψ + C 0 , 1 ( I ¯ ) such that
  • ψ + ψ > 0 on I 1 ;
  • ψ + ψ = 0 on I I ¯ 1 ;
  • φ = 0 on I I ¯ 2 ;
  • The graph of the function φ is located strictly inside the domain Ω .
Let us fix a small parameter δ > 0 and put
Ω ± = { ( x 1 , x 2 ) Ω | ± x 2 > ± φ ( x 1 ) , x 1 I } ,
Ω m δ = { ( y 1 , y 2 ) Ω | φ ( y 1 ) + δ ψ ( y 1 ) < y 2 < φ ( y 1 ) + δ ψ + ( y 1 ) , y 1 I 1 } ,
S ± δ = { ( y 1 , y 2 ) Ω | y 2 = φ ( y 1 ) + δ ψ ± ( y 1 ) , y 1 I 1 } ,
Ω e δ = Ω Ω m δ ¯ , Ω ± δ = Ω e δ Ω ± .
Note that, for all small-enough δ > 0 , the domain Ω m δ lies strictly inside Ω .
In our consideration, the domain Ω is an elastic inhomogeneous body, containing a narrow curvilinear strip Ω m δ with width of order δ and with a rough boundary S ¯ + δ S ¯ δ , where roughness is characterized by functions ψ ± (see Figure 1).
By C δ , C 0 , we denote the tensors characterizing the material constants of the inclusion Ω m δ and the elastic matrix Ω e δ , respectively, with usual symmetrical and elliptical properties (see, e.g., [30]). We prescribe homogeneous Dirichlet’s conditions on Γ D and Neumann’s conditions on the remaining part Γ N of the external boundary Ω . This means that the body Ω is clamped on Γ D , whereas the force g L 2 ( Γ N ) is applied to Γ N .
We formulate the equilibrium problem for the body Ω with the inclusion Ω m δ in the framework of the two-dimensional elasticity. Let u = ( u 1 , u 2 ) be a vector of displacements of the composite body and σ ( u ) = ( σ i j ( u ) ) i , j = 1 , 2 and ε ( u ) = ( ε i j ( u ) ) i , j = 1 , 2 be the stress and the strain tensors, respectively, where
ε i j ( u ) = 1 2 ( u i , j + u j , i ) , i , j = 1 , 2 ,
σ ( u ) = C 0 ε ( u ) in Ω e δ , C δ ε ( u ) = def δ 1 C 1 ε ( u ) = δ 1 λ m I t r ε ( u ) + 2 μ m ε ( u ) in Ω m δ .
Here and further, lower indices after the comma denote the operation of differentiation with respect to the corresponding coordinate, the summation over repeated indices is performed, λ m and λ m are the Lamé coefficients such that λ m + 2 μ m > 0 and μ m > 0 ; hence, the tensor C 1 is positive definite. In general, the elastic matrix characterized by the tensor C 0 is anisotropic and homogeneous with the elastic coefficients belonging to L l o c ( R 2 ) .
The equilibrium problem under study is the following boundary value problem:
σ i j , j ( u δ ) = 0 in Ω , i = 1 , 2 ,
u δ = 0 on Γ D ,
σ i j ( u δ ) n j = g i on Γ N , i = 1 , 2 ,
where n = ( n 1 , n 2 ) is a normal unit vector to Ω , g = ( g 1 , g 2 ) L 2 ( Γ N ) 2 is a given traction.
In the sequel, the following variational (weak) formulation of Problems (1)–(3) is used: find a function u δ H Γ D ( Ω ) satisfying the variational equality
Ω σ ( u δ ) : ε ( v ) d y = Γ N g v d s
for all kinematically admissible displacement test functions v H Γ D 1 ( Ω ) , where
H Γ D ( Ω ) = { v H 1 ( Ω ) 2 | v = 0 a . e . on Γ D } .

3. Transformations of Coordinates

Introduce into consideration the following sets:
Ω m = { ( z 1 , z 2 ) R 2 | ψ ( z 1 ) < z 2 < ψ + ( z 1 ) , z 1 I 1 } ,
S ± = { ( z 1 , z 2 ) R 2 | z 2 = ψ ± ( z 1 ) , z 1 I 1 } ,
S = { ( x 1 , x 2 ) R 2 | x 2 = φ ( x 1 ) , x 1 I 1 } ,
which are independent of δ , and define coordinate transformations mapping domains Ω ± δ and Ω m δ onto domains Ω ± and Ω m , respectively. For this purpose, take two convex open domains D 1 and D 2 such that I 2 = D 2 O y 1 , I ¯ 1 D 1 , D ¯ 1 D 2 , D ¯ 2 Ω , and consider the following domains:
D i φ = { ( z 1 , z 2 ) | z 1 = y 1 , z 2 = φ ( y 1 ) + y 2 , ( y 1 , y 2 ) D i } , i = 1 , 2 .
The inclusions S ¯ D 1 φ and D ¯ 1 φ D 2 φ hold. Additionally, assume that domain D 2 φ lies strictly inside Ω .
Further, consider a smooth cut-off function θ C 1 ( Ω ¯ ) such that
θ = 1 in D 1 φ ; 0 < θ < 1 in D 2 φ ; θ = 0 in Ω D ¯ 2 φ .
Finally, define the coordinate transformations x = Λ ± ( y ) and z = Λ m ( y ) by the formulas
x 1 = y 1 , x 2 = y 2 δ ψ ± ( y 1 ) θ ( y 1 , y 2 ) , ( y 1 , y 2 ) Ω ± δ , ( x 1 , x 2 ) Ω ± ,
z 1 = y 1 , z 2 = y 2 φ ( y 1 ) δ , ( y 1 , y 2 ) Ω m δ , ( z 1 , z 2 ) Ω m ,
respectively. It is easy to see that Transformations (5) and (6) are one-to-one. Moreover, these coordinate transformations generate an isomorphism between H 1 ( Ω ± ) , H 1 ( Ω m ) and H 1 ( Ω ± δ ) , H 1 ( Ω m δ ) , respectively (see, e.g., Chapter 2, Lemma 3.4 in [31]). By y = Λ ± 1 ( x ) and y = Λ m 1 ( z ) , we further denote the inverse transformations to (5) and (6), respectively.

4. Decomposition and Transformation of the Problem

Following the arguments from [1,2,32] with natural modifications, we now rewrite Problem (4) in the equivalent form. We introduce the set
K δ = { ( v + , v , v m ) H 1 ( Ω + δ ) × H 1 ( Ω δ ) × H 1 ( Ω m δ ) | v ± = 0 a . e . on Γ D Ω ± δ , v ± = v m a . e . on S ± δ , v + = v a . e . on Ω δ Ω + δ } ,
and define a variational problem: find a triplet ( u δ + , u δ , u δ m ) K δ satisfying the equality
Ω + δ C 0 ε ( u δ + ) : ε ( v + ) d y + Ω δ C 0 ε ( u δ ) : ε ( v ) d y + 1 δ Ω m δ C 1 ε ( u δ m ) : ε ( v m ) d y = Γ N Ω + δ g v + d s + Γ N Ω δ g v d s
for all ( v + , v , v m ) K δ .
Clearly, Problem (7) has a unique solution ( u δ + , u δ , u δ m ) for all δ > 0 . Moreover, the relation
u δ = u δ ± a . e . in Ω ± δ , u δ m a . e . in Ω m δ
holds, where u δ is the solution to Problem (4).
Let K H 1 ( Ω ) × H 1 ( Ω + ) × H 1 ( Ω m ) be the image of the set K δ under the coordinate transformations (5) and (6),
K = { ( v + , v , v m ) H 1 ( Ω + ) × H 1 ( Ω ) × H 1 ( Ω m ) | v ± | S = v m | S ± , v ± = 0 a . e . on Γ D Ω ± , v + = v a . e . on ( Ω Ω + ) S ¯ } .
The equality v ± | S = v m | S ± in the definition of K means that
v ± ( x 1 , φ ( x 1 ) ) = v m ( x 1 , ψ ± ( x 1 ) ) for   all x 1 I 1 .
After applying coordinate transformations y = Λ ± 1 ( x ) and y = Λ m 1 ( z ) to the variational equality (7), we conclude that the transformed displacements
u ˜ δ ± ( x ) = u δ ± ( Λ ± 1 ( x ) ) , u ˜ δ m ( z ) = u δ m ( Λ m 1 ( z ) )
satisfy the variational equality
A + δ ( u ˜ δ + , v + ) + A δ ( u ˜ δ , v )       + Ω m C 1 e δ ( u ˜ δ m ) : e δ ( v m ) d z = Ω + Γ N g v + d s + Ω Γ N g v d s
for any triplet ( v + , v , v m ) K , where
A ± δ ( u ± , v ± ) = Ω ± J δ ± C 0 E ( Ψ δ ± ; u ± ) : E ( Ψ δ ± ; v ± ) d x ,
E i j ( Ψ δ ± ; v ) = ( 1 / 2 ) v i , k ( Ψ δ ± ) k j + v j , k ( Ψ δ ± ) k i , i , j = 1 , 2 ;
e δ ( v m ) = v 1 , 1 m φ δ v 1 , 2 m 1 2 1 δ v 1 , 2 m + v 2 , 1 m φ δ v 2 , 2 m 1 2 1 δ v 1 , 2 m + v 2 , 1 m φ δ v 2 , 2 m 1 δ v 2 , 2 m ,
with the Jacobi matrices Ψ δ ± and the Jacobian J δ ± of the transformations y = Λ ± 1 ( x ) .
In the sequel, we use the first-order asymptotic expansions (see, e.g., [33,34])
A ± δ ( u , v ) = A ± ( u , v ) + δ r ± ( δ , u , v ) , A ± ( u , v ) = Ω ± C 0 ε ( u ) : ε ( v ) d x ,
where
| r ± ( δ , u , v ) | c ( u H 1 ( Ω ± ) 2 + v H 1 ( Ω ± ) 2 ) ,
with a constant c independent of δ , u , and v .

5. Justification of the Asymptotic Analysis

It is well-known that uniform estimates play an essential role in asymptotic analysis (see, e.g., [20,35,36,37]. Here, we formulate and prove some Korn-type estimates. Primarily, we recall the following Lemmas 1 and 2.
Lemma 1
([2]). For any function v H 1 ( Ω m ) the following inequality holds:
v L 2 ( Ω m ) 2 c v , 2 L 2 ( Ω m ) 2 + v L 2 ( S α ) 2 , α { , + } ,
with a constant c independent of v.
Lemma 2
([1]). For all small-enough δ > 0 and any triplet ( v + , v , v m ) K the following inequality holds:
c v + H 1 ( Ω + ) 2 + v H 1 ( Ω ) 2 + δ v 1 , 1 m φ δ v 1 , 2 m L 2 ( Ω m ) 2 + 1 δ v 1 , 2 m L 2 ( Ω m ) 2 + δ v 2 , 1 m φ δ v 2 , 2 m L 2 ( Ω m ) 2 + 1 δ v 2 , 2 m L 2 ( Ω m ) 2 Ω + C 0 ε ( v + ) : ε ( v + ) d x + Ω C 0 ε ( v ) : ε ( v ) d x     + δ Ω m C 1 e δ ( v m ) : e δ ( v m ) d z
with a constant c independent of δ and ( v + , v , v m ) .
Lemmas 1 and 2 allow us to prove Theorems 1–3 that follow.
Theorems 1 and 2 are preparatory in nature.
Theorem 1.
The solution ( u ˜ δ + , u ˜ δ , u ˜ δ m ) of Problem (9) satisfies the following uniform in δ estimates:
u ˜ δ + H 1 ( Ω + ) 2 c , u ˜ δ H 1 ( Ω ) 2 c ,
δ 1 u ˜ δ 2 , 2 m L 2 ( Ω m ) c ,
u ˜ δ 1 , 1 m φ δ u ˜ δ 1 , 2 m L 2 ( Ω m ) c , u ˜ δ 2 , 1 m + u ˜ δ 1 , 2 m δ L 2 ( Ω m ) c ,
δ 1 2 u ˜ δ 1 , 2 m L 2 ( Ω m ) c .
Proof. 
Let us take ( u ˜ δ + , u ˜ δ , u ˜ δ m ) in (7) as the test functions and apply the coordinate transformations (5) and (6). Taking into account (10), (11), and Korn’s inequality in the domains Ω ± , we obtain (12)–(14). Estimate (15) follows directly from Lemma 2 and (12)–(14). □
Theorem 2.
There exists a sequence, still denoted by δ, and limit functions u ± H 1 ( Ω ± ) and p i L 2 ( Ω m ) , i = 1 , 2 , 3 , such that the following limiting relations
u ˜ δ ± u ± w e a k l y   i n H 1 ( Ω ± ) 2 ,
u ˜ δ m u m w e a k l y   i n L 2 ( Ω m ) 2 ,
u ˜ δ 1 , 1 m φ δ u ˜ δ 1 , 2 m p 1 w e a k l y   i n L 2 ( Ω m ) ,
u ˜ δ 1 , 2 m δ + u ˜ δ 2 , 1 m p 2 w e a k l y   i n L 2 ( Ω m ) ,
u ˜ δ 2 , 2 m δ p 3 w e a k l y   i n L 2 ( Ω m ) ,
u ˜ δ , 2 m u , 2 m = 0 s t r o n g l y   i n L 2 ( Ω m ) 2 .
u ˜ δ 1 , 1 m + φ u ˜ δ 2 , 1 m p 1 + φ p 2 w e a k l y   i n L 2 ( Ω m )
hold true as δ 0 + .
Proof. 
The first six limiting relations follow from Theorem 1 and the weak compactness property of bounded sets in the spaces L 2 ( Ω m ) and H 1 ( Ω ± ) .
Now, let us prove (17). For any function η C 0 ( Ω m ) , we have
Ω m ( p 1 + φ p 2 ) η d z = lim δ 0 Ω m u ˜ δ 1 , 1 m φ δ u ˜ δ 1 , 2 m + φ ( u ˜ δ 1 , 2 m δ + u ˜ δ 2 , 1 m ) η d z
= lim δ 0 Ω m u ˜ δ 1 , 1 m + φ u ˜ δ 2 , 1 m η d z = lim δ 0 Ω m u ˜ δ 1 m η , 1 + u ˜ δ 2 m ( φ η ) , 1 d z
= Ω m ( ( u 1 m + φ u 2 m ) η , 1 + φ η u 2 ) d z ,
which gives (17). Moreover, additionally the following equality holds true:
p 1 + φ p 2 = u 1 , 1 m + φ u 2 , 1 m L 2 ( Ω m ) .
Remark 1.
Note that, in general, functions u 1 , 1 m and u 2 , 1 m do not belong to the space L 2 ( Ω m ) .
Remark 2.
Due to (16), the function u m does not depend on z 2 . From this fact, it follows that the traces of functions u + and u on S coincide with each other.
Theorem 3 is the main result of this paper.
Theorem 3.
Let ( u ˜ δ + , u ˜ δ , u ˜ δ m ) K be a solution of Problem (9) and ( u + , u ) be a solution of the following variational problem: find u = ( u + , u ) satisfying the variational equality
Ω C 0 ε ( u ) : ε ( v ) d x + S 4 μ m ( λ m + μ m ) ( ψ + ψ ) ( λ m + 2 μ m ) ( 1 + φ 2 ) 1 2 ( τ t u τ ) ( τ t v τ ) d s = Γ N g v d s
for all v K l , where u ± = u | Ω ± are the restrictions of u on Ω ± ,
K l = v H Γ D ( Ω ) | ( ψ + ψ ) 1 2 ( τ t v τ ) L 2 ( S ) ,
v = v 1 , 1 v 1 , 2 v 2 , 1 v 2 , 2 ,
τ is a tangent vector to the curve S defined by
τ = ( 1 , φ ) ( 1 + φ 2 ) 1 2 ,
and the superscript t denotes the transposition operator defined for matrices.
Then, the limiting relations
( u ˜ δ + , u ˜ δ ) ( u + , u ) s t r o n g l y   i n H 1 ( Ω + ) × H 1 ( Ω ) ,
u ˜ δ m u m s t r o n g l y   i n L 2 ( Ω m )
hold true as δ 0 + , where u ± are defined above and u m ( z 1 , z 2 ) = u ( z 1 , ψ ( z 1 ) ) is a trace of u on S extended into the domain Ω m .
Proof. 
Multiply (9) by δ and let δ tend to zero. As the result, using Theorem 2, we obtain
φ ( ( λ m + 2 μ m ) p 1 + λ m p 3 ) + μ ( p 2 φ p 3 ) = 0 a . e . in Ω m ,
φ μ m ( p 2 φ p 3 ) + ( λ m + 3 μ m ) p 3 + λ m p 1 = 0 a . e . in Ω m .
Thus, together with (18), we obtain the system of the linear algebraic equations for p 1 , p 2 , p 3 with the non-vanishing discriminant D = μ m ( λ m + 2 μ m ) ( 1 + ( φ ) 2 ) 2 , which is a Lipschitz continuous function on I ¯ 1 .
Solving the system, we find the functions p i L 2 ( Ω m ) , i = 1 , 2 , 3 :
p 1 = λ m + 2 μ m λ m φ 2 ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 ( u m 1 , 1 + φ u m 2 , 1 ) ,
p 2 = φ ( 3 λ m + 4 μ m + ( λ m + 2 μ m ) φ 2 ) ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 ( u m 1 , 1 + φ u m 2 , 1 ) ,
p 3 = φ ( λ m + 2 μ m ) λ m ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 ( u m 1 , 1 + φ u m 2 , 1 ) .
Now, take an arbitrary pair ( v + , v ) ( C 1 ( Ω + ) × C 1 ( Ω ) ) K l and put v m ( z 1 , z 2 ) = v + ( z 1 , φ ( z 1 ) ) for ( z 1 , z 2 ) Ω m . Then, the triple ( v + , v , v m ) belongs to the set K and, therefore, can be substituted into (9) as test functions. After passing to the limit as δ goes to zero, owing to (22)–(24), we obtain
A + ( u + , v + ) + A ( u , v ) + Ω m 4 μ m ( λ m + μ m ) ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 ( u m 1 , 1 + φ u m 2 , 1 ) ( v m 1 , 1 + φ v m 2 , 1 ) d z = Ω + Γ N g v + d s + Ω Γ N g v d s .
Due to the definition of the domain Ω m and the fact that the functions in the third term of (25) do not depend on z 2 , Equation (25) can be rewritten as follows:
A + ( u + , v + ) + A ( u , v ) + I 1 4 μ m ( λ m + μ m ) ( ψ + ψ ) ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 ( u m 1 , 1 + φ u m 2 , 1 ) ( v m 1 , 1 + φ v m 2 , 1 ) d z 1 = Ω + Γ N g v + d s + Ω Γ N g v d s .
In turn, taking into account (8), Equality (26) can be rewritten as follows:
A + ( u + , v + ) + A ( u , v ) + I 1 4 μ m ( λ m + μ m ) ( ψ + ψ ) ( λ m + 2 μ m ) ( 1 + φ 2 ) 2 × ( u 1 , 1 + + φ u 1 , 2 + + φ ( u 2 , 1 + + φ u 2 , 2 + ) ) ( v 1 , 1 + + φ v 1 , 2 + + φ ( v 2 , 1 + + φ v 2 , 2 + ) ) d z 1 = Ω + Γ N g v + d s + Ω Γ N g v d s .
Finally, due to Remark 2 and the density of C 1 ( Ω ) K l in K l , we obtain (19).
The strong limiting relations (20) and (21) follow from the standard arguments (see, e.g., [32]). □
Now, assuming the additional regularity of the solution to Problem (19), we give the equivalent (in the sense of distributions) differential formulation of (19).
First of all, note that
w s = w τ , d τ d s = k ν , d ν d s = k τ ,
where w is an arbitrary rather smooth scalar function, and k is the curvature of the curve S defined by
k = φ ( 1 + ( φ ) 2 ) 3 2 .
The last two formulas in (27) are the Frenet–Serret formulas.
Let w τ and w ν be tangential and normal components of a vector-function w in the ( τ , ν ) -coordinate system, respectively. Due to (27), we have
w s = w τ s τ + k w τ ν + w ν s ν k w ν τ ,
which yields that
τ t w τ = w τ s k w ν .
Then, after the application of Green’s formula to the first integral and integration by parts in the second integral in (19), we arrive at the following boundary value problem:
σ i j , j ( u ) = 0 in Ω S ¯ , i = 1 , 2 ,
u = 0 on Γ D ,
σ i j ( u ) n j = g i on Γ N , i = 1 , 2 ,
s Λ ψ u τ s k u ν = [ σ τ ( u ) ] on S ,
k Λ ψ u τ s k u ν = [ σ ν ( u ) ] on S ,
Λ ψ u τ s k u ν = 0 at S ,
where
Λ ψ = 4 μ m ( λ m + μ m ) λ m + 2 μ m ψ + ψ ( 1 + φ 2 ) 1 2 .
Let us give a mechanical interpretation of the differential equations and boundary conditions (28)–(33). Here, Equation (28) is the standard equilibrium equation of the body Ω containing the thin inclusion S. Condition (29) means that the body is clamped on Γ D , while (30) means that traction g = ( g 1 , g 2 ) is applied to Γ N . Equations (31) and (32) describe deformations of inclusion S taking into account the influence of the elastic matrix surrounding it. It is worth noting that, in the case when k 0 , not only is the jump of the tangential stress σ τ ( u ) is not equal to zero, but the jump of the normal stress σ ν ( u ) on S is not equal to zero as well. Moreover, these jumps are of the Ventcel type (see, e.g., works [15,22,38,39,40,41], in which analogous conditions were used for different models of elasticity). At last, Equation (33) is the boundary condition at the tips of S that are free of loads.

6. Numerical Examples

In this section, by using a standard finite-element method, we compare numerically the solution of the initial problem (4) with the approximated problem (19). In all the examples below, the domain Ω is a square Ω = ( 1 , 1 ) × ( 1 , 1 ) , the body is fixed on the right-hand side Γ D = { 1 } × ( 1 , 1 ) of the square, and the remaining part Γ N = Ω Γ ¯ D of the boundary Ω is divided into three parts such that Γ N = Γ N 1 Γ N 2 Γ N 3 with Γ N 1 = { 1 } × ( 1 , 1 ) , Γ N 2 = ( 1 , 1 ) × { 1 } , and Γ N 3 = ( 1 , 1 ) × { 1 } , which are the left-hand, the upper, and the lower sides of the domain Ω , respectively.
The elastic matrix Ω e δ is assumed to be an isotropic homogeneous material with the Lamé coefficients μ and λ :
μ = E 2 ( 1 + ν ) , λ = 2 ν μ 1 2 ν ,
where E and ν are the Young modulus and the Poisson ratio, respectively. We take the material parameters ν = 0.32 and E = 112 GPa for the elastic matrix Ω e δ . Moreover, we assume that the Lamé coefficients μ m and λ m are defined by the same formulas with the Young modulus E m and the Poisson ratio ν m :
ν = 0.28 , E m = 10 p E ,
where parameter p takes values from the range −2, −1, 0, 1, and 2, characterizing the softness (or hardness) of the inclusion Ω m δ .
The horizontal tensile loading by the following traction forces is imposed: g = ( g 1 , 0 ) with g 1 = 10 3 μ m on Γ N 3 , and g 2 = 0 on Γ N 2 Γ N 1 .
The numerical experiments are implemented using the free software FreeFEM++ (Version 4.11, https://freefem.org/) [42]. The space H 1 ( Ω ) is approximated by the finite-element space consisting of linear functions, namely, of Lagrange’s P1-elements.
Below, we consider several cases of the geometry of the inclusion and get numerical convergence of approximated solutions to the solution of the initial problem with the inclusion of a non-zero width. Let S be a sinusoidal curve, i.e., S be described by the graph of the function y 2 = φ = d sin ( 2 π y 1 ) , where d takes values from 0 (rectilinear line) to 0.5 (curve with non-zero curvature). Moreover, we take functions ψ ± characterizing the roughness of the narrow inclusion Ω m δ as follows:
ψ + ( y 1 ) = 2 y 1 2 0.5 , ψ ( y 1 ) = 16 y 1 4 + 1 .
By ErrL2 and ErrH1, the relative errors of the displacements in the L 2 - and H 1 -spaces, respectively, are denoted.
In Table 1 (see also Figure 2), we investigate the dependence of the relative error on the parameter δ characterizing inclusion’s width (for d = 0.5 , p = 0 ). They show good approximation of the initial problem (with non-zero width) by the thin inclusion problem in both L 2 - and H 1 -norms.
Table 2 shows the dependence of relative error on the parameter p, which characterizes the roughness of the inclusion. It can be seen that the approximation for the soft inclusions is better than the approximation for the hard inclusions. This feature can be explained by the fact that in the case of the hard inclusion, there are large stresses in the vertices of the inclusions and in the places of great curvature. It should be noted that, in [43], the asymptotic behavior of the solution near the tip of the rigid line inclusion was studied in the case of two-dimensional homogeneous isotropic linearized elasticity. In [44], the numerical method of solving bulk and thin rigid inclusion problems was proposed.
Table 3 illustrates the dependence of the relative errors on the parameter d that characterizes the curvature of the inclusion (the case of d = 0 corresponds to the rectilinear inclusion, and the case d = 0.5 corresponds to the sinusoid with an amplitude equal to one).
Let us investigate a mesh sensitivity analysis. Let h . m i n and h . m a x denote the minimal and maximal sizes of mesh triangles of the domain Ω . Table 4 presents numerical calculations in dependence of mesh sizes. Again, it shows a good approximation of the initial problem (with non-zero width) by the thin inclusion problem in both L 2 - and H 1 -norms.
At last, the resulting deformations in the Lagrange coordinates (with the tenfold amplification factor) and the distribution of the von Mises stresses are depicted in Figure 3 for the curvilinear inclusion (case d = 0.5 ) for the applied loading as stated at the beginning of Section 6 and in Figure 4 for the rectilinear inclusion (case d = 0 ) for the non-symmetric tensile loading g = ( g 1 , 0 ) , where g 1 = 10 3 μ m on Γ N 3 { y 2 < 0 } and g 1 = 0 on ( Γ N 3 { y 2 > 0 } ) Γ N 2 Γ N 1 . It can be seen that, in the case of the curvilinear inclusion, there is a singularity in the vicinity of the curved part of the inclusion, while in the case of the rectilinear inclusion singularity, it appears only in the vicinity of its tips. This means that the curvature k plays an essential role in modeling curvilinear inclusions.

7. Conclusions

By using a variational approach, the model of the elastic body containing strictly inside the curvilinear rod-like thin inclusion has been derived. The model approximates the equilibrium problem of the non-homogeneous elastic body with the narrow inclusion, the width and elastic properties of which depend on the small parameter δ (the width is proportional to δ , while elastic properties are proportional to δ 1 ). The approximated model has a variational formulation. This allows the application of the standard finite-element method in order to compare solutions of the initial problem and approximated one. Numerical experiments show a good approximation in both L 2 -norm and H 1 -norm as the parameter δ goes to zero.
In the future, it is planned to perform an asymptotic procedure for other types of dependence of the elastic properties of the inclusion on the small parameter characterizing its width and justify models of thin inclusion of different types (beam-like, rigid, etc.).

Author Contributions

Conceptualization, E.R.; writing—original draft, E.R. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of Science and Higher Education of the Russian Federation, grant number FWGG-2022-0001.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fankina, I.V.; Furtsev, A.I.; Rudoy, E.M.; Sazhenkov, S.A. Asymptotic modeling of curvilinear narrow inclusions with rough boundaries in elastic bodies: Case of a soft inclusion. Sib. Electron. Math. Rep. 2022, 19, 935–948. [Google Scholar]
  2. Furtsev, A.; Rudoy, E. Variational approach to modeling soft and stiff interfaces in the Kirchhoff-Love theory of plates. Int. J. Solids Struct. 2020, 202, 562–574. [Google Scholar] [CrossRef]
  3. Rudoy, E.M. Asymptotic justification of models of plates containing inside hard thin inclusions. Technologies 2017, 8, 59. [Google Scholar] [CrossRef]
  4. Itou, H.; Khludnev, A.M. On delaminated thin Timoshenko inclusions inside elastic bodies. ZAMM Z. Angew. Math. Mech. 2016, 24, 4980–4993. [Google Scholar] [CrossRef]
  5. Khludnev, A.; Negri, M. Crack on the boundary of a thin elastic inclusion inside an elastic body. ZAMM Z. Angew. Math. Mech. 2012, 95, 341–354. [Google Scholar] [CrossRef]
  6. Khludnev, A. On thin inclusions in elastic bodies with defects. Z. Angew. Math. Phys. 2019, 70, 45. [Google Scholar] [CrossRef]
  7. Khludnev, A. T-shape inclusion in elastic body with a damage parameter. J. Comput. Appl. Math. 2021, 393, 113540. [Google Scholar] [CrossRef]
  8. Khludnev, A. Non-coercive problems for Kirchhoff–Love plates with thin rigid inclusion. Z. Angew. Math. Phys. 2022, 73, 54. [Google Scholar] [CrossRef]
  9. Lazarev, N.P. Equilibrium problem for a thermoelastic Kirchhoff-Love plate with a delaminated flat rigid inclusion. Philos. Trans. R. Soc. A 2022, 380, 20210363. [Google Scholar] [CrossRef] [PubMed]
  10. Lazarev, N.P.; Kovtunenko, V.A. Signorini-type problems over non-convex sets for composite bodies contacting by sharp edges of rigid inclusions. Mathematics 2022, 10, 250. [Google Scholar] [CrossRef]
  11. Khludnev, A. Junction problem for thin elastic and volume rigid inclusions in elastic bod. Philos. Trans. R. Soc. A 2022, 380, 20210360. [Google Scholar] [CrossRef]
  12. Khludnev, A. On the Crossing Bridge between Two Kirchhoff–Love Plates. Axioms 2023, 12, 120. [Google Scholar] [CrossRef]
  13. Khludnev, A. Asymptotic analysis of thin elastic junction between two elastic bodies. Appl. Anal. 2023, 102, 800–814. [Google Scholar] [CrossRef]
  14. Baranova, S.; Mogilevskaya, S.G.; Nguyen, T.H.; Schilling, D. Higher-order imperfect interface modeling via complex variables based asymptotic analysis. Int. J. Eng. Sci. 2021, 57, 103399. [Google Scholar]
  15. Fankina, I.V.; Furtsev, A.I.; Rudoy, E.M.; Sazhenkov, S.A. Multiscale analysis of stationary thermoelastic vibrations of a composite material. Philos. Trans. R. Soc. A 2022, 380, 20210354. [Google Scholar]
  16. Mokhtari, H.; Rahmani, L. Asymptotic modeling of the behavior of a reinforced plate governed by a full von Karman thermo-elastic system with nonlinear thermal coupling. Ann. Univ. Ferrara Sez. VII Sci. Mat. 2021, 67, 387–414. [Google Scholar] [CrossRef]
  17. Serpilli, M.; Dumont, S.; Rizzoni, R.; Lebon, F. Interface Models in Coupled Thermoelasticity. Technologies 2021, 9, 17. [Google Scholar] [CrossRef]
  18. Benveniste, Y.; Miloh, T. Imperfect soft and stiff interfaces in two-dimensional elasticity. Mech. Mater. 2001, 33, 309–323. [Google Scholar] [CrossRef]
  19. Bessoud, A.-L.; Krasucki, F.; Serpilli, M. Plate-like and shell-like inclusions with high rigidity. Comptes Rendus Math. 2009, 346, 697–702. [Google Scholar]
  20. Ljulj, M.; Tambaca, J. 3D structure—2D plate interaction model. Math. Mech. Solids 2019, 24, 3354–3377. [Google Scholar] [CrossRef]
  21. Serpilli, M.; Rizzoni, R.; Rodríguez-Ramos, R.; Lebon, F.; Dumont, S. A novel form of imperfect contact laws in flexoelectricity. Compos. Struct. 2022, 300, 116059. [Google Scholar] [CrossRef]
  22. Bessoud, A.-L.; Krasucki, F.; Michaille, G. Multi-materials with strong interface: Variational modelings. Asymptot. Anal. 2012, 61, 1–19. [Google Scholar] [CrossRef] [Green Version]
  23. Dumont, S.; Lebon, F.; Raffa, M.L.; Rizzoni, R. Towards nonlinear imperfect interface models including micro-cracks and smooth roughness. Ann. Solid Struct. Mech. 2017, 9, 13–27. [Google Scholar] [CrossRef] [Green Version]
  24. Fouchal, F.; Lebon, F.; Raffa, M.L.; Vairo, G. An interface model including cracks and roughness applied to masonry. Open Civ. Eng. J. 2014, 8, 263–271. [Google Scholar] [CrossRef]
  25. Kunets, Y.I.; Matus, V.V. Asymptotic approach in the dynamic problems of the theory of elasticity for bodies with thin elastic inclusions. J. Math. Sci. 2023, 270, 87–106. [Google Scholar] [CrossRef]
  26. Mokhtari, H.; Rahmani, L. Asymptotic modeling of a reinforced plate with a thin layer of variable thickness. Meccanica 2022, 57, 2155–2172. [Google Scholar] [CrossRef]
  27. Dong, H.; Zheng, X.; Cui, J.; Nie, Y.; Yang, Z.; Ma, Q. Multi-scale computational method for dynamic thermo-mechanical performance of heterogeneous shell structures with orthogonal periodic configurations. Comput. Methods Appl. Mech. Eng. 2019, 354, 143–180. [Google Scholar] [CrossRef]
  28. Ma, Q.; Li, Z.; Cui, J. Multi-scale asymptotic analysis and computation of the elliptic eigenvalue problems in curvilinear coordinates. Comput. Methods Appl. Mech. Eng. 2018, 340, 340–365. [Google Scholar] [CrossRef]
  29. Ma, Q.; Cui, J.; Li, Z.; Wang, Z. Second-order asymptotic algorithm for heat conduction problems of periodic composite materials in curvilinear coordinates. J. Comput. Appl. Math. 2016, 306, 87–115. [Google Scholar] [CrossRef]
  30. Khludnev, A.M.; Kovtunenko, V.A. Analysis of Cracks in Solids; WIT-Press: Southampton, UK, 2000. [Google Scholar]
  31. Nečas, J. Direct Methods in the Theory of Elliptic Equations; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  32. Rudoy, E.M.; Itou, H.; Lazarev, N.P. Asymptotic Justification of the Models of Thin Inclusions in an Elastic Body in the Antiplane Shear Problem. J. Appl. Ind. Math. 2021, 15, 129–140. [Google Scholar] [CrossRef]
  33. Rudoy, E. First-order and second-order sensitivity analyses for a body with a thin rigid inclusion. Math. Methods Appl. Sci. 2016, 39, 4994–5006. [Google Scholar] [CrossRef]
  34. Rudoy, E.; Shcherbakov, V. First-Order Shape Derivative of the Energy for Elastic Plates with Rigid Inclusions and Interfacial Cracks. Appl. Math. Optim. 2021, 84, 2775–2802. [Google Scholar] [CrossRef]
  35. Freddi, L.; Morassi, A.; Paroni, R. Thin-Walled Beams: The Case of the Rectangular Cross-Section. J. Elast. 2004, 76, 45–66. [Google Scholar] [CrossRef]
  36. Gaudiello, A.; Sili, A. Limit models for thin heterogeneous structures with high contrast. J. Differ. Equ. 2021, 302, 37–63. [Google Scholar] [CrossRef]
  37. Nazarov, S.A. Weighted anisotropic Korn’s inequality for a junction of a plate and a rod. Sb. Math. 2004, 195, 553–583. [Google Scholar] [CrossRef]
  38. Baranova, S.; Mogilevskaya, S.G. On the Bovik–Benveniste methodology and related approaches for modelling thin layers. Philos. Trans. R. Soc. A 2022, 380, 20210420. [Google Scholar] [CrossRef] [PubMed]
  39. Fankina, I.V.; Furtsev, A.I.; Rudoy, E.M.; Sazhenkov, S.A. The homogenized quasi-static model of a thermoelastic composite stitched with reinforcing threads. J. Comput. Appl. Math. 2023, 434, 115346. [Google Scholar] [CrossRef]
  40. Kushch, V.I.; Mogilevskaya, S.G. Higher order imperfect interface models of conductive spherical interphase. Math. Mech. Solids 2022. [Google Scholar] [CrossRef]
  41. Ventcel, A.D. On boundary conditions for multi-dimensional diffusion processes. Theory Probab. Appl. 1959, 42, 164–177. [Google Scholar]
  42. Hecht, F. New development in FreeFem++. J. Numer. Math. 2012, 20, 251–265. [Google Scholar] [CrossRef] [Green Version]
  43. Itou, H.; Khludnev, A.M.; Rudoy, E.M.; Tani, A. Asymptotic behaviour at a tip of a rigid line inclusion in linearized elasticity. ZAMM Z. Angew. Math. Mech. 2012, 92, 716–730. [Google Scholar] [CrossRef]
  44. Rudoy, E.M. On numerical solving a rigid inclusions problem in 2D elasticity. ZAMP Z. Angew. Math. Phys. 2017, 68, 19. [Google Scholar] [CrossRef]
Figure 1. Domain configuration.
Figure 1. Domain configuration.
Mathematics 11 03447 g001
Figure 2. Relative error vs. δ for the case d = 0.5 , p = 0 .
Figure 2. Relative error vs. δ for the case d = 0.5 , p = 0 .
Mathematics 11 03447 g002
Figure 3. Deformed body with the sinusoidal inclusion and distribution of the von Mises stresses.
Figure 3. Deformed body with the sinusoidal inclusion and distribution of the von Mises stresses.
Mathematics 11 03447 g003
Figure 4. Deformed body with the rectilinear inclusion and distribution of the von Mises stresses.
Figure 4. Deformed body with the rectilinear inclusion and distribution of the von Mises stresses.
Mathematics 11 03447 g004
Table 1. Relative error vs. δ for the case d = 0.5 , p = 0 .
Table 1. Relative error vs. δ for the case d = 0.5 , p = 0 .
δ = 0.1 δ = 0.05 δ = 0.025 δ = 0.01
ErrL20.040.0210.0110.004
ErrH10.1480.1070.0770.048
Table 2. Relative error vs. p for the case d = 0.5 , δ = 0.01 .
Table 2. Relative error vs. p for the case d = 0.5 , δ = 0.01 .
p = 2 p = 1 p = 0 p = 1 p = 2
ErrL20.0010.0040.0050.0080.019
ErrH10.0130.0430.0480.0560.086
Table 3. Relative error vs. d for the case δ = 0.01 , p = 0 .
Table 3. Relative error vs. d for the case δ = 0.01 , p = 0 .
00.10.20.30.40.5
ErrL20.0020.0040.0050.0050.0040.005
ErrH10.0540.0510.0520.050.0480.048
Table 4. Relative error vs. mesh sizes h . m i n and h . m a x for the case δ = 0.01 , p = 0 .
Table 4. Relative error vs. mesh sizes h . m i n and h . m a x for the case δ = 0.01 , p = 0 .
h . min h . max ErrL2ErrH1
0.080.280.0130.074
0.030.140.0050.057
0.0130.0760.0050.051
0.0090.0480.0050.049
0.0050.0420.0040.048
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rudoy, E.; Sazhenkov, S. Variational Approach to Modeling of Curvilinear Thin Inclusions with Rough Boundaries in Elastic Bodies: Case of a Rod-Type Inclusion. Mathematics 2023, 11, 3447. https://doi.org/10.3390/math11163447

AMA Style

Rudoy E, Sazhenkov S. Variational Approach to Modeling of Curvilinear Thin Inclusions with Rough Boundaries in Elastic Bodies: Case of a Rod-Type Inclusion. Mathematics. 2023; 11(16):3447. https://doi.org/10.3390/math11163447

Chicago/Turabian Style

Rudoy, Evgeny, and Sergey Sazhenkov. 2023. "Variational Approach to Modeling of Curvilinear Thin Inclusions with Rough Boundaries in Elastic Bodies: Case of a Rod-Type Inclusion" Mathematics 11, no. 16: 3447. https://doi.org/10.3390/math11163447

APA Style

Rudoy, E., & Sazhenkov, S. (2023). Variational Approach to Modeling of Curvilinear Thin Inclusions with Rough Boundaries in Elastic Bodies: Case of a Rod-Type Inclusion. Mathematics, 11(16), 3447. https://doi.org/10.3390/math11163447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop