Next Article in Journal
Analysis of Blockchain in the Healthcare Sector: Application and Issues
Next Article in Special Issue
Positive Solutions for a Fractional Differential Equation with Sequential Derivatives and Nonlocal Boundary Conditions
Previous Article in Journal
Catalytic Asymmetric Synthesis of C-Chiral Phosphonates
Previous Article in Special Issue
Numerical Integration of Stiff Differential Systems Using Non-Fixed Step-Size Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conjugation Conditions for Systems of Differential Equations of Different Orders on a Star Graph

by
Baltabek Kanguzhin
and
Gauhar Auzerkhan
*
Institute of Mathematics and Mathematical Modeling, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(9), 1761; https://doi.org/10.3390/sym14091761
Submission received: 24 July 2022 / Revised: 17 August 2022 / Accepted: 19 August 2022 / Published: 23 August 2022

Abstract

:
In this paper, a one-dimensional mathematical model for investigating the vibrations of structures consisting of elastic and weakly curved rods is proposed. The three-dimensional structure is replaced by a limit graph, on each arc of which a system of three differential equations is written out. The differential equations describe the longitudinal and transverse vibrations of an elastic rod, taking into account the influence of longitudinal and transverse vibrations on each other. Describing conjugation conditions at joints of four or more rods is an important problem. This article assumes new conjugation conditions that guarantee the all-around decidability and symmetry of the resulting boundary value problems for systems of differential equations on a star graph. In addition, the paper proposes a physical interpretation of the conjugation conditions found. Thus, the work presents one more area of knowledge where symmetry phenomena occur. The symmetry here is manifested in the preservation of conjugation conditions when passing to the conjugate operator.

1. Introduction

The work is devoted to the description of well-posed boundary value problems for systems of differential equations on graphs. In contrast to well-known works, on each edge of the graph, there is a system of three differential equations with various orders that do not coincide with each other. At the e n of each arc of star graph G, let there be a system of differential equations at n = 1 , , m + 1 and z n l n 2 , l n 2 .
d 2 d z n 2 μ n z n a n z n d 2 w 1 n z n d z n 2 + b n z n d 2 w 2 n z n d z n 2 d n z n d w 3 n z n d z n = F 1 n z n , d 2 d z n 2 μ n z n b n z n d 2 w 1 n z n d z n 2 + c n z n d 2 w 2 n z n d z n 2 f n z n d w 3 n z n d z n = F 2 n z n , d d z n μ n z n d n z n d 2 w 1 n z n d z n 2 f n z n d 2 w 2 n z n d z n 2 + d w 3 n z n d z n = F 3 n z n .
Assume that the coefficients
μ n ( z n ) , a n ( z n ) , b n ( z n ) , d n ( z n ) , c n ( z n ) , f n ( z n ) , n = 1 , , m + 1 give simple real-valued functions. Furthermore, the physical meaning of the coefficients is explained. The required functions w 1 n ( z n ) , w 2 n ( z n ) , w 3 n ( z n ) , n = 1 , , m + 1 represent the transverse and longitudinal displacements of rods from their respective axes. The rods are considered to be connected at their ends in one knot. The other ends of the rods are hard-fixed. The hard-fixed ends of the rods in the graph correspond to the boundary vertices of the graph, and the connecting node corresponds to the inner vertex of the star graph. The main question that interests us is: what conjugation conditions can be fulfilled at an interior vertex of the star graph so that the corresponding boundary value problem for the system of differential equations is symmetric? In case the boundary value problem is symmetric, then the eigenfrequencies will be real. As far as physical meaning is concerned, these types of boundary tasks are interesting. The system (1) is introduced in [1] in connection with the asymptotic analysis of one-dimensional equations of deformation of thin weakly curved rods. Thus, this paper presents one more area of knowledge where symmetry plays an important role.
Usually, cross-sections of rods [2] are characterized by space, static moments of plane sections, the center of gravity position, moments of inertia, radii of inertia, and moments of resistances. Let z change along the rod. We denote the transverse cross-section of the rod at point z by ω ( z ) . The coefficients of the system (1) are then input by the formulas d ( z ) = ω ( z ) η 1 d η 1 d η 2 ,   f ( z ) = ω ( z ) η 2 d η 1 d η 2 ,   a ( z ) = ω ( z ) η 1 2 d η 1 d η 2 ,   b ( z ) = ω ( z ) η 1 η 2 d η 1 d η 2 ,   c ( z ) = ω ( z ) η 2 2 d η 1 d η 2 .
Then, according to [2], d ( z ) = ω ( z ) η 1 d η 1 d η 2 ,   f ( z ) = ω ( z ) η 2 d η 1 d η 2 are the static moments of cross-section area w ( z ) relative to the axes O η 1 and O η 2 , respectively.
The values a ( z ) = ω ( z ) η 1 2 d η 1 d η 2 ,   c ( z ) = ω ( z ) η 2 2 d η 1 d η 2 represent the axial moments of inertia of section ω ( z ) relative to the axes O η 1 and O η 2 , respectively. The centrifugal moment of the inertia of cross-section ω ( z ) concerning the two co-orthogonal axes O η 1 and O η 2 is equal to b ( z ) = ω ( z ) η 1 η 2 d η 1 d η 2 . If the diameter of the cross-section w ( z ) is considered to be a small of order ε , then
  • The static moments of the cross-sectional area d ( z ) ,   f ( z ) have an order of smallness ε 3 ;
  • The cross-sectional axial moments of inertia a ( z ) ,   c ( z ) have the order of smallness ε 4 ;
  • The centrifugal moment of inertia of the section b ( z ) has an order of smallness ε 4 .
Consider separately the third equation of system (1)
d d z n μ n ( z n ) d n ( z n ) d 2 w 1 n ( z n ) d z n 2 f n ( z n ) d 2 w 2 n ( z n ) d z n 2 + d w 3 n ( z n ) d z n = F 3 n ( z n ) .
As there are coefficients ε + 0 , d n ( z n ) = O ( ε 3 ) ,   f n ( z n ) = O ( ε 3 ) , the first two terms can be neglected in the last equation. As a result, we have (2)
d d z μ ( z ) d w 3 ( z ) d z = F 3 ( z ) .
Similar equations describe the longitudinal vibrations of rods (3). Thus, if the cross-sectional diameter of the rod is ε 0 , then the influence of the transverse vibrations on the longitudinal vibrations can be neglected.
Under natural simplifying assumptions, this is the equation of longitudinal vibrations of the rod [3]. Therefore, the longitudinal offsets w 3 ( z ) can be found first. The transverse offsets can be defined from the first two equations of system (1).
d 2 d z 2 μ ( z ) a ( z ) d 2 w 1 ( z ) d z 2 + b ( z ) d 2 w 2 ( z ) d z 2 d ( z ) d w 3 ( z ) d z = F 1 ( z ) , d 2 d z 2 μ ( z ) b ( z ) d 2 w 1 ( z ) d z 2 + c ( z ) d 2 w 2 ( z ) d z 2 f ( z ) d w 3 ( z ) d z = F 2 ( z ) .
The last system of equations falls into two systems, and each of them can be solved independently from the other. The fourth-order differential equation for transverse displacement along the O y axis:
d 2 d z 2 μ ( z ) ( b 2 ( z ) a ( z ) c ( z ) ) d 2 w 2 ( z ) d z 2 = F ˜ 2 ( z ) ,
where
F ˜ 2 ( z ) = b ( z ) F 1 ( z ) + d 2 d z 2 d ( z ) d w 3 ( z ) d z a ( z ) F 2 ( z ) + d 2 d z 2 f ( z ) d w 3 ( z ) d z .
The fourth-order differential equation for transverse displacement in the 0 x -axis:
d 2 d z 2 μ ( z ) ( a ( z ) c ( z ) b 2 ( z ) ) d 2 w 1 ( z ) d z 2 = F ˜ 1 ( z ) ,
where
F ˜ 1 ( z ) = c ( z ) F 1 ( z ) + d 2 d z 2 d ( z ) d W 3 ( z ) d z b ( z ) F 2 ( z ) + d 2 d z 2 f ( z ) d W 3 ( z ) d z .
Thus, the transverse displacements in the O x and O y axes, in this case, can be computed independently of each other. This confirms Timoshenko’s theory of beam bends, which proposes that the bends are determined from fourth-order differential equations. If the section diameter ω ( z ) is considered to be a small order of ε , then system (1) confirms Timoshenko’s theory of bending beams [4].
In many engineering calculations, it is presumed that movements are separated: transverse vibrations do not affect longitudinal ones and vice versa. However, this division of rod movements is not always justified. Thus, in general, system (1) does not always decompose into equations of (2) and (3) types.
Consider a mechanical system (Figure 1) consisting of n rods in three dimensions. The ends of the rods are connected to each other at X ε , o , which is a node. If the parameter ε tend toward zero, we obtain the limit graph (Figure 2).
The rods are cylinders X ε , e = I e × ε Y e with a cross-sectional radius ε of the Y e variety for each rod e = 1 , 2 , , n . We denote X ε , o by the central manifold. In the limit, the family { X ε , o , X ε , e , e = 1 , 2 , , n } at ε 0 converges to the limit graph. The limit operator represents the limits of the differential operators describing the transverse and longitudinal vibrations of the rods. The limit operator depends on the choice of coupling boundary conditions at the central node where the rods connect to the node. The limit operator depends on the choice of boundary conditions and conjugation conditions at the central node where the rods connect to the node.
Different types of rod–node connections correspond to different conjugation conditions. In the literature, there are different types of connections between a rod and a node (Figure 3 and Figure 4). Figure 3 shows the case where the dimensions of the node tend to zero much more slowly than the dimensions of the cross-sections of the rods. Figure 4 shows the case where the dimensions of the node tend to zero much faster than the cross-sectional dimensions of the rods. In [5,6,7,8,9], limiting conjugation conditions corresponding to the two mentioned connections of knots with rods have been investigated. The slowly decaying and borderline cases were introduced in [7] (see also [8]) showing the convergence of the spectrum for compact graphs and manifolds. The notion “plumber’s shop” was introduced in the article of Rubinstein and Schatzman [5], where the authors used it for the necessary local estimates of the identification operators from the graph to the graph-like manifold and vice versa. In work [9], the authors extend the analysis here to non-compact spaces and show in particular the quasi-unitary quivalence implying, e.g., the convergence of resolvents and the convergence of the entire spectrum. Other types of boundary conditions for limit operators have been studied in [10,11,12,13]. The authors in the work [14] investigate the initial boundary value problem describing the oscillation process on a geometric graph with hysteresis-type boundary conditions.

2. Lagrange Formula on a Star Graph in a Conjugate Condition at the Internal Vertex

Let Γ = { V , E } , a graph, be a star (Figure 2), where V is the set of its vertices numbered from 0 to m + 1 , and the set E means its arcs e 1 , , e m . We introduce a vector function [15,16]
Y j ( x j ) = y 1 j ( x j ) y 2 j ( x j ) y 3 j ( x j ) T .
A system of differential equations is given on each arc e j
l 1 j ( Y j ) = d 2 d x j 2 μ j ( x j ) a j ( x j ) d 2 y 1 j ( x j ) d x j 2 + d 2 d x j 2 μ j ( x j ) b j ( x j ) d 2 y 2 j ( x j ) d x j 2
d 2 d x j 2 μ j ( x j ) d j ( x j ) d y 3 j ( x j ) d x j ,
l 2 j ( Y j ) = d 2 d x j 2 μ j ( x j ) b j ( x j ) d 2 y 1 j ( x j ) d x j 2 + d 2 d x j 2 μ j ( x j ) c j ( x j ) d 2 y 2 j ( x j ) d x j 2
d 2 d x j 2 μ j ( x j ) f j ( x j ) d y 3 j ( x j ) d x j ,
l 3 j ( Y j ) = d d x j μ j ( x j ) d j ( x j ) d 2 y 1 j ( x j ) d x j 2 + d d x j μ j ( x j ) f j ( x j ) d 2 y 2 j ( x j ) d x j 2
d d x j μ j ( x j ) d y 3 j ( x j ) d x j .
Furthermore, we assume that the functions μ j ( x j ) ,   a j ( x j ) ,   b j ( x j ) ,   c j ( x j ) ,   d j ( x j ) ,   f j ( x j ) are real continuous on the arc e j . In this case, x j [ 0 , l j ] , and x j belongs to e j . The vertex ( m + 1 ) V is called the internal vertex of the star graph. The vertices 0 , 1 , , m are called the boundary vertices of the star graph. Denote by
L j ( Y j ) = l 1 j ( Y j ) l 2 j ( Y j ) l 3 j ( Y j ) T
at j = 1 , , m + 1 . Introduce formally conjugate differential expressions
l 1 j + ( Y j ) = d 2 d x j 2 μ j ( x j ) a j ( x j ) d 2 y 1 j ( x j ) d x j 2 + d 2 d x j 2 μ j ( x j ) b j ( x j ) d 2 y 2 j ( x j ) d x j 2
+ d d x j μ j ( x j ) d j ( x j ) d 2 y 3 j ( x j ) d x j 2 ,
l 2 j + ( Y j ) = d 2 d x j 2 μ j ( x j ) b j ( x j ) d 2 y 1 j ( x j ) d x j 2 + d 2 d x j 2 μ j ( x j ) c j ( x j ) d 2 y 2 j ( x j ) d x j 2
+ d d x j μ j ( x j ) f j ( x j ) d 2 y 3 j ( x j ) d x j 2 ,
l 3 j + ( Y j ) = d 2 d x j 2 μ j ( x j ) d j ( x j ) d y 1 j ( x j ) d x j d 2 d x j 2 μ j ( x j ) f j ( x j ) d y 2 j ( x j ) d x j
d d x j μ j ( x j ) d y 3 j ( x j ) d x j .
We will also need a designation L j + ( Y j ) = l 1 j + ( Y j ) l 2 j + ( Y j ) l 3 j + ( Y j ) T at j = 1 , , m + 1 . Let us also introduce quasi-derivatives according to the formulas at the point ξ :
D 1 j ( 0 ) ( Y j ; ξ ) = [ y 1 j ] | x j = ξ , D 2 j ( 0 ) ( Y j ; ξ ) = [ y 2 j ] | x j = ξ , D 3 j ( 0 ) ( Y j ; ξ ) = [ y 3 j ] | x j = ξ ,
D 1 j ( 1 ) ( Y j ; ξ ) = [ d y 1 j / d x j ] | x j = ξ , D 2 j ( 1 ) ( Y j ; ξ ) = [ d y 2 j / d x j ] | x j = ξ ,
D 1 j ( 2 ) ( Y j ; ξ ) = μ j a j d 2 y 1 j d x j 2 + μ j b j d 2 y 2 j d x j 2 μ j d j d y 3 j d x j | x j = ξ ,
D 2 j ( 2 ) ( Y j ; ξ ) = μ j b j d 2 y 1 j d x j 2 + μ j c j d 2 y 2 j d x j 2 μ j f j d y 3 j d x j | x j = ξ ,
D 3 j ( 2 ) ( Y j ; ξ ) = μ j d j d 2 y 1 j d x j 2 + μ j f j d 2 y 2 j d x j 2 μ j d y 3 j d x j | x j = ξ ,
D 2 j ( 3 ) ( Y j ; ξ ) = d d x j μ j b j d 2 y 1 j d x j 2 + d d x j μ j c j d 2 y 2 j d x j 2 d d x j μ j f j d y 3 j d x j | x j = ξ ,
D 1 j ( 3 ) ( Y j ; ξ ) = d d x j μ j a j d 2 y 1 j d x j 2 + d d x j μ j b j d 2 y 2 j d x j 2 d d x j μ j d j d y 3 j d x j | x j = ξ .
Sometimes, instead of D k j ( s ) ( Y j ; ξ ) , we will write D k j ( s ) ( Y j ) . Denote by
L ( Y ) = L 1 ( Y 1 ) L 2 ( Y 2 ) L m + 1 ( Y m + 1 ) , Z = Z 1 Z 2 Z m + 1 , Y = Y 1 Y 2 Y m + 1 .
Introduce a scalar product of
< L ( Y ) ; Z > L 2 ( Γ ) = j = 1 m + 1 < L ( Y j ) ; Z J > L 2 ( e j ) = j = 1 m + 1 i = 1 3 < l i j ( Y j ) ; z i j > L 2 ( e j ) .
Lemma 1.
For any two sets of sufficiently smooth functions Y = { Y j ( x j ) , j = 1 , , m + 1 } ,   Z = { Z j ( x j ) , j = 1 , , m + 1 } , the Lagrange identity is valid
< L ( Y ) ; Z > L 2 ( Γ ) < Y , L + ( z ) > L 2 ( Γ ) = j = 1 m + 1 { D 1 j ( 3 ) ( Y j ) D 1 j ( 0 ) ( Z j ¯ ) D 1 j ( 2 ) ( Y j ) D 1 j ( 1 ) ( Z j ¯ ) + D 2 j ( 3 ) ( Y j ) D 2 j ( 0 ) ( Z j ¯ ) D 2 j ( 2 ) ( Y j ) D 2 j ( 1 ) ( Z j ¯ ) + D 3 j ( 2 ) ( Y j ) D 3 j ( 0 ) ( Z j ¯ ) + D 1 j ( 1 ) ( Y j ) D 1 j ( 2 ) ( Z j ¯ ) D 1 j ( 0 ) ( Y j ) D 1 j ( 3 ) ( Z j ¯ ) + D 2 j ( 1 ) ( Y j ) D 2 j ( 2 ) ( Z j ¯ ) D 2 j ( 0 ) ( Y j ) D 2 j ( 3 ) ( Z j ¯ ) D 3 j ( 0 ) ( Y j ) D 3 j ( 2 ) ( Z j ¯ ) } | 0 l j .
Proof of Lemma 1.
Consider a scalar product of
< L ( Y ) ; Z > = j = 1 m + 1 i = 1 3 0 l j l i j ( Y j ) · z ¯ i j ( x j ) d x j .
Using the fractional integration formula, we obtain the following
0 l j l 1 j ( Y j ) · z ¯ 1 j ( x j ) d x j = [ d d t ( μ j ( t ) ( a j ( t ) d 2 y 1 j ( t ) d t 2 + b j ( t ) d 2 y 2 j ( t ) d t 2
d j ( t ) d y 3 j ( t ) d t ) ) z 1 j ( t ) ¯ ] | 0 l j [ μ j ( t ) ( a j ( t ) d 2 y 1 j ( t ) d t 2 + b j ( t ) d 2 y 2 j ( t ) d t 2
d j ( t ) d y 3 j ( t ) d t ) d z 1 j ( t ) ¯ d t ] | 0 l j + [ d y 1 j ( t ) d t μ j ( t ) a j ( t ) d 2 z 1 j ( t ) ¯ d t 2 + d y 2 j ( t ) d t μ j ( t ) b j ( t ) d 2 z 1 j ( t ) ¯ d t 2
y 3 j ( t ) μ j ( t ) d j ( t ) d 2 z 1 j ( t ) ¯ d t 2 ] | 0 l j [ y 1 j ( t ) d d t μ j ( t ) a j ( t ) d 2 z 1 j ( t ) ¯ d t 2
+ y 2 j ( t ) d d t μ j ( t ) b j ( t ) d 2 z 1 j ( t ) ¯ d t 2 ] | 0 l j + 0 l j y 1 j ( x j ) · l 1 j + ( z 1 j ¯ ) d x j .
The integrals are converted in a similar way 0 l j l 2 j ( Y j ) · z 1 j ¯ ( x j ) d x j ,   0 l j l 3 j ( Y j ) · z 1 j ¯ ( x j ) d x j . Using the above relations, we derive identity (4). Lemma 1 is completely proved. From now on, we assume that the boundary vertices of the graph satisfy the rigid anchoring conditions
D 1 j ( 0 ) ( Y j ; ξ j ) = 0 , D 2 j ( 0 ) ( Y j ; ξ j ) = 0 , D 1 j ( 1 ) ( Y j ; ξ j ) = 0 , D 2 j ( 1 ) ( Y j ; ξ j ) = 0 , D 3 j ( 0 ) ( Y j ; ξ j ) = 0 , j = 1 , , m + 1 .
 □
Here, ξ j = l j ,   j = 1 , , m . At the same time, ξ m + 1 = 0 . Let for a set of functions Z = { Z j ( x j ) , j = 1 , , m + 1 } Equation (5) be also fulfilled. Then, the Lagrange identity follows from Lemma 1.
< L ( Y ) ; Z > L 2 ( Γ ) < Y , L + ( z ) > L 2 ( Γ ) = { D 1 m + 1 ( 3 ) ( Y m + 1 ) D 1 m + 1 ( 0 ) ( Z m + 1 ¯ ) D 1 m + 1 ( 2 ) ( Y m + 1 ) D 1 m + 1 ( 1 ) ( Z m + 1 ¯ ) + D 2 m + 1 ( 3 ) ( Y m + 1 ) D 2 m + 1 ( 0 ) ( Z m + 1 ¯ ) D 2 m + 1 ( 2 ) ( Y m + 1 ) D 2 m + 1 ( 1 ) ( Z m + 1 ¯ ) + D 3 m + 1 ( 2 ) ( Y m + 1 ) D 3 m + 1 ( 0 ) ( Z m + 1 ¯ ) + D 1 m + 1 ( 1 ) ( Y m + 1 ) D 1 m + 1 ( 2 ) ( Z m + 1 ¯ ) D 1 m + 1 ( 0 ) ( Y m + 1 ) D 1 m + 1 ( 3 ) ( Z m + 1 ¯ ) + D 2 m + 1 ( 1 ) ( Y m + 1 ) D 2 m + 1 ( 2 ) ( Z m + 1 ¯ ) D 2 m + 1 ( 0 ) ( Y m + 1 ) D 2 m + 1 ( 3 ) ( Z m + 1 ¯ ) D 3 m + 1 ( 0 ) ( Y m + 1 ) D 3 m + 1 ( 2 ) ( Z m + 1 ¯ ) } | l m + 1 j = 1 m { D 1 j ( 3 ) ( Y j ) D 1 j ( 0 ) ( Z j ¯ ) D 1 j ( 2 ) ( Y j ) D 1 j ( 1 ) ( Z j ¯ ) + D 2 j ( 3 ) ( Y j ) D 2 j ( 0 ) ( Z j ¯ ) D 2 j ( 2 ) ( Y j ) D 2 j ( 1 ) ( Z j ¯ ) + D 3 j ( 2 ) ( Y j ) D 3 j ( 0 ) ( Z j ¯ ) + D 1 j ( 1 ) ( Y j ) D 1 j ( 2 ) ( Z j ¯ ) D 1 j ( 0 ) ( Y j ) D 1 j ( 3 ) ( Z j ¯ ) + D 2 j ( 1 ) ( Y j ) D 2 j ( 2 ) ( Z j ¯ ) D 2 j ( 0 ) ( Y j ) D 2 j ( 3 ) ( Z j ¯ ) D 3 j ( 0 ) ( Y j ) D 3 j ( 2 ) ( Z j ¯ ) } | 0
At the inner vertex of the graph, we require a continuity condition.
D s m + 1 ( 0 ) ( Z m + 1 ¯ l m + 1 ) = D s j ( 0 ) ( Z j ¯ , 0 ) , D s m + 1 ( 0 ) ( Y m + 1 l m + 1 ) = D s j ( 0 ) ( Y j , 0 ) j = 1 , , m , s = 1 , 2
Since it is easier to bend a rod than to stretch it, the longitudinal displacements of the rod are influenced by the first derivatives of the transverse displacements. In particular, such an effect has been highlighted by the authors of the paper [17]. From now on, we will assume that the total longitudinal displacement of the rod end is determined by the expression
D 3 j ( 0 ) ( Y j , ξ j ) η 1 j D 1 j ( 1 ) ( Y j , ξ j ) η 2 j D 2 j ( 1 ) ( Y j , ξ j ) .
where ξ j j are the end rod corresponding to the connection node, η 1 j ,   η 2 j are some constants. Let there now be several rods connected at the same node. Then, three rods can be selected as reference rods (independent), and the total longitudinal displacements of the other rods can be considered dependent on the selected three reference rods. A similar effect only in the case of a flat bar connection has been observed in reference work [18]. Let the numbered i , j , k rods be the base rods.
Then, according to our assumption, the total longitudinal displacement of any rod with a number s { 1 , 2 , , m + 1 } { i , j , k } is a linear combination of the total longitudinal displacements of the rods with numbers i , j , k . Thus, in the inner vertex of the graph
D 3 s ( 0 ) ( Y s , ξ s ) η 1 s D 1 s ( 1 ) ( Y s , ξ s ) η 2 s D 2 s ( 1 ) ( Y s , ξ s ) = t = i , j , k α s t D 3 t ( 0 ) ( Y t , ξ t ) η 1 t D 1 t ( 1 ) ( Y t , ξ t ) η 2 t D 2 t ( 1 ) ( Y t , ξ t ) .
where α s t represent certain constants. Consider that ξ m + 1 = l m + 1 ,   ξ i = 0 ,   i = { 1 , , m } .
To record the conjugation conditions at the inner vertex of the star graph, let us introduce the following linear forms:
C j ( Y j , Y m + 1 ) = D 1 m + 1 ( 0 ) ( Y m + 1 , l m + 1 ) D 1 j ( 0 ) ( Y j , 0 ) j = 1 , , m ,
C j + m ( Y j , Y m + 1 ) = D 2 m + 1 ( 0 ) ( Y m + 1 , l m + 1 ) D 2 j ( 0 ) ( Y j , 0 ) j = 1 , , m ,
C s + 2 m ( Y s , Y m 1 , Y m , Y m + 1 ) = D 3 s ( 0 ) ( Y s , ξ s ) η 1 s D 1 s ( 1 ) ( Y s , ξ s ) η 2 s D 2 s ( 1 ) ( Y s , ξ s )
t = m 1 m + 1 α s t D 3 t ( 0 ) ( Y t , ξ t ) η 1 t D 1 t ( 1 ) ( Y t , ξ t ) η 2 t D 2 t ( 1 ) ( Y t , ξ t ) , s = 1 , , m 2 ,
C 3 m 1 ( Y 1 , , Y m , Y m + 1 ) = D 1 m + 1 ( 3 ) ( Y m + 1 , l m + 1 ) j = 1 m D 1 j ( 3 ) ( Y j , 0 ) ,
C 3 m ( Y 1 , , Y m , Y m + 1 ) = D 2 m + 1 ( 3 ) ( Y m + 1 , l m + 1 ) j = 1 m D 2 j ( 3 ) ( Y j , 0 ) ,
C 3 m + 1 ( Y 1 , , Y m , Y m + 1 ) = D 3 m + 1 ( 2 ) ( Y m + 1 , l m + 1 ) j = 1 m 2 α j m + 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 2 ( Y 1 , , Y m , Y m + 1 ) = D 3 m ( 2 ) ( Y m , 0 ) + j = 1 m 2 α j m D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 3 ( Y 1 , , Y m , Y m + 1 ) = D 3 m 1 ( 2 ) ( Y m 1 , 0 ) + j = 1 m 2 α j m 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 4 ( Y 1 , , Y m , Y m + 1 ) = D 1 m 1 ( 2 ) ( Y m 1 , 0 ) + η 1 m 1 j = 1 m 2 α j m 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 5 ( Y 1 , , Y m , Y m + 1 ) = D 1 m ( 2 ) ( Y m , 0 ) + η 1 m j = 1 m 2 α j m D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 6 ( Y 1 , , Y m , Y m + 1 ) = D 1 m + 1 ( 2 ) ( Y m + 1 , l m + 1 ) + η 1 m + 1 j = 1 m 2 α j m + 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 7 ( Y 1 , , Y m , Y m + 1 ) = D 2 m 1 ( 2 ) ( Y m 1 , 0 ) + η 2 m 1 j = 1 m 2 α j m 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 8 ( Y 1 , , Y m , Y m + 1 ) = D 2 m ( 2 ) ( Y m , 0 ) + η 2 m j = 1 m 2 α j m D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 9 ( Y 1 , , Y m , Y m + 1 ) = D 2 m + 1 ( 2 ) ( Y m + 1 , l m + 1 ) + η 2 m + 1 j = 1 m 2 α j m + 1 D 3 j ( 2 ) ( Y j , 0 ) ,
C 3 m + 9 + j ( Y 1 , , Y m , Y m + 1 ) = η 1 j D 3 j ( 2 ) ( Y j , 0 ) D 1 j ( 2 ) ( Y j , 0 ) , j = 1 , , m 2 ,
C 4 m + 7 + j ( Y 1 , , Y m , Y m + 1 ) = η 2 j D 3 j ( 2 ) ( Y j , 0 ) D 2 j ( 2 ) ( Y j , 0 ) , j = 1 , , m 2 .
Then, by substituting relations (7) and (8) in the Lagrange identity (6), we obtain the equality
< L ( Y ) ; Z > L 2 ( Γ ) < Y , L + ( z ) > L 2 ( Γ )
= D 1 m + 1 ( 0 ) ( Y m + 1 , l m + 1 ) C 3 m 1 ( Z 1 , , Z m , Z m + 1 ) D 2 m + 1 ( 0 ) ( Y m + 1 , l m + 1 ) C 3 m ( Z 1 , , Z m , Z m + 1 )
+ D 1 m + 1 ( 0 ) ( Z m + 1 , l m + 1 ) C 3 m 1 ( Y 1 , , Y m , Y m + 1 ) + D 2 m + 1 ( 0 ) ( Z m + 1 , l m + 1 ) C 3 m ( Y 1 , , Y m , Y m + 1 )
+ D 3 m + 1 ( 0 ) ( Z m + 1 , l m + 1 ) C 3 m + 1 ( Y 1 , , Y m , Y m + 1 ) D 3 m + 1 ( 0 ) ( Y m + 1 , l m + 1 ) C 3 m + 1 ( Z 1 , , Z m , Z m + 1 )
D 3 m ( 0 ) ( Z m , 0 ) C 3 m + 2 ( Y 1 , , Y m , Y m + 1 ) + D 3 m ( 0 ) ( Y m , 0 ) C 3 m + 2 ( Z 1 , , Z m , Z m + 1 )
D 3 m 1 ( 0 ) ( Z m 1 , 0 ) C 3 m + 3 ( Y 1 , , Y m , Y m + 1 ) + D 3 m 1 ( 0 ) ( Y m 1 , 0 ) C 3 m + 3 ( Z 1 , , Z m , Z m + 1 )
+ j = 1 m 2 D 1 j ( 1 ) ( Y j , 0 ) C 3 m + 9 + j ( Z 1 , , Z m , Z m + 1 ) j = 1 m 2 D 1 j ( 1 ) ( Z j , 0 ) C 3 m + 9 + j ( Y 1 , , Y m , Y m + 1 )
+ j = 1 m 2 D 2 j ( 1 ) ( Y j , 0 ) C 4 m + 7 + j ( Z 1 , , Z m , Z m + 1 ) j = 1 m 2 D 2 j ( 1 ) ( Z j , 0 ) C 4 m + 7 + j ( Y 1 , , Y m , Y m + 1 )
D 1 m 1 ( 1 ) ( Y m 1 , 0 ) C 3 m + 4 ( Z 1 , , Z m , Z m + 1 ) D 1 m ( 1 ) ( Y m , 0 ) C 3 m + 5 ( Z 1 , , Z m , Z m + 1 )
D 1 m + 1 ( 1 ) ( Y m + 1 , l m + 1 ) C 3 m + 6 ( Z 1 , , Z m , Z m + 1 ) D 1 m 1 ( 1 ) ( Y m 1 , 0 ) C 3 m + 7 ( Z 1 , , Z m , Z m + 1 )
D 1 m ( 1 ) ( Y m , 0 ) C 3 m + 8 ( Z 1 , , Z m , Z m + 1 ) D 1 m + 1 ( 1 ) ( Y m + 1 , l m + 1 ) C 3 m + 9 ( Z 1 , , Z m , Z m + 1 )
+ D 1 m 1 ( 1 ) ( Z m 1 , 0 ) C 3 m + 4 ( Y 1 , , Y m , Y m + 1 ) + D 1 m ( 1 ) ( Z m , 0 ) C 3 m + 5 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m + 1 ( 1 ) ( Z m + 1 , l m + 1 ) C 3 m + 6 ( Y 1 , , Y m , Y m + 1 ) + D 1 m 1 ( 1 ) ( Z m 1 , 0 ) C 3 m + 7 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m ( 1 ) ( Z m , 0 ) C 3 m + 8 ( Y 1 , , Y m , Y m + 1 ) + D 1 m + 1 ( 1 ) ( Z m + 1 , l m + 1 ) C 3 m + 9 ( Y 1 , , Y m , Y m + 1 ) .
Taking into account the last identity, we introduce the main object of study of the article. We also consider the inhomogeneous operator equation L ( Y ) = F with conjugation conditions in the interior vertex C j ( Y ) = 0 ,   j = 1 , , 5 m + 5 and with rigid fixing conditions in the boundary vertices (5) on the star graph Γ . The number of boundary conditions is also 5 m + 5 . The operator corresponding to the above inhomogeneous problem is denoted by B .
Theorem 1.
The domains of definition of operators B and B * coincide, i.e., they are given by the same conjugation conditions and fixation conditions.
The proof of Theorem 1 immediately follows from the Lagrange relationship, which is written in the form (9).
Remark 1.
The meaning of Theorem 1 is that the operators B and B * are self-adjoint in the sense of definitions.
Remark 2.
The mechanical interpretation of a part of conjugation conditions C j ( Y ) = 0 ,   j = 1 , , 3 m 2 is given above when introducing relations (7) and (8). The conjugation conditions C j ( Y ) = 0 ,   j = 3 m 1 , , 5 m + 5 also have a mechanical interpretation, indicating the distribution of forces and moments between the rods in the function. The conjugation conditions given in Theorem 1 seem to be new.
In fact, Theorem 1 is one of the main results of this paper. It states that we have found self-adjoint conjugate conditions in the interior vertex of the graph. The conjugation conditions we have found generalize the well-known Kirchhoff conditions, which are written for second-order differential equations. We have been able to generalize Kirchhoff conditions for systems of differential equations consisting of differential equations of different orders. In [10,11,12,13], different variants of Kirchhoff analogs conditions of the graph’s inner vertex for the second-order differential equations are given. In a sense, Theorem 1 generalizes the results of [10,11,12,13] for systems of differential equations on a star graph.

3. The Reversibility of Operator B

The operator B depends on the real parameters η 1 j ,   η 2 j ,   α s t at all valid indexes. For operator B to be reversible, it is necessary to impose constraints on the specified parameters. In this section, we will find out the values of the above parameters that make possible the existence of the inverse operator B 1 . Let ξ m + 1 = 0 ,   ξ j = b j ,   j = 1 , , m . Denote by
p j ( t j ) = μ j ( t j ) a j ( t j ) b ( t j ) d j ( t j ) b ( t j ) c j ( t j ) f j ( t j ) d j ( t j ) f j ( t j ) 1 , q 1 j 1 ( t j ) = c j ( t j ) f j ( t j ) f j ( t j ) 1 ,
q 1 j 2 ( t j ) = b j ( t j ) d j ( t j ) f j ( t j ) 1 , q 1 j 3 ( t j ) = b j ( t j ) d j ( t j ) c j ( t j ) f j ( t j ) ,
q 2 j 1 ( t j ) = b j ( t j ) f j ( t j ) d j ( t j ) 1 , q 2 j 2 ( t j ) = a j ( t j ) d j ( t j ) d j ( t j ) 1 , q 2 j 3 ( t j ) = a j ( t j ) d j ( t j ) b j ( t j ) f j ( t j ) ,
q 3 j 1 ( t j ) = b j ( t j ) c j ( t j ) d j ( t j ) f j ( t j ) , q 3 j 2 ( t j ) = a j ( t j ) b j ( t j ) d j ( t j ) f j ( t j ) , q 3 j 3 ( t j ) = a j ( t j ) d j ( t j ) b j ( t j ) c j ( t j ) .
Then, the solutions of the homogeneous system of differential Equations
l 1 j ( Y j ) = 0 , l 2 j ( Y j ) = 0 , l 3 j ( Y j ) = 0
with rigid anchoring conditions at the boundary vertices
D 1 j ( 0 ) ( Y j ; ξ j ) = 0 , D 2 j ( 0 ) ( Y j ; ξ j ) = 0 , D 1 j ( 1 ) ( Y j ; ξ j ) = 0 , D 2 j ( 1 ) ( Y j ; ξ j ) = 0 , D 3 j ( 0 ) ( Y j ; ξ j ) = 0
at j = 1 , , m + 1 will be
y 1 j ( x j ) = K 1 j ξ j x j ( t j ξ j ) t j q 1 j 1 ( t j ) p j ( t j ) d t j + N 1 j ξ j x j ( t j ξ j ) q 1 j 1 ( t j ) p j ( t j ) d t j
K 2 j ξ j x j ( t j ξ j ) t j q 1 j 2 ( t j ) p j ( t j ) d t j N 2 j ξ j x j ( t j ξ j ) q 1 j 2 ( t j ) p j ( t j ) d t j + K 3 j ξ j x j ( t j ξ j ) q 1 j 3 ( t j ) p j ( t j ) d t j ,
We see solution y 2 j ( x j )
y 2 j ( x j ) = K 1 j ξ j x j ( t j ξ j ) t j q 2 j 1 ( t j ) p j ( t j ) d t j N 1 j ξ j x j ( t j ξ j ) q 2 j 1 ( t j ) p j ( t j ) d t j
+ K 2 j ξ j x j ( t j ξ j ) t j q 2 j 2 ( t j ) p j ( t j ) d t j + N 2 j ξ j x j ( t j ξ j ) q 2 j 2 ( t j ) p j ( t j ) d t j K 3 j ξ j x j ( t j ξ j ) q 2 j 3 ( t j ) p j ( t j ) d t j ,
y 3 j ( x j ) = K 1 j ξ j x j t j q 3 j 1 ( t j ) p j ( t j ) d t j N 1 j ξ j x j q 3 j 1 ( t j ) p j ( t j ) d t j
+ K 2 j ξ j x j t j q 3 j 2 ( t j ) p j ( t j ) d t j + N 2 j ξ j x j q 3 j 2 ( t j ) p j ( t j ) d t j K 3 j ξ j x j q 3 j 3 ( t j ) p j ( t j ) d t j .
where K 1 j , N 1 j , K 2 j , N 2 j , K 3 j are arbitrary constants.
In accordance with the above representation of the solution, let us introduce the following solutions to the homogeneous system of differential Equation (10):
θ s j 1 ( x j ) = ( 1 ) s 1 ξ j x j ( t j ξ j ) t j q s j 1 ( t j ) p j ( t j ) d t j , θ s j 2 ( x j ) = ( 1 ) s 1 ξ j x j ( t j ξ j ) q s j 1 ( t j ) p j ( t j ) d t j ,
θ s j 3 ( x j ) = ( 1 ) s ξ j x j ( t j ξ j ) t j q s j 2 ( t j ) p j ( t j ) d t j , θ s j 4 ( x j ) = ( 1 ) s ξ j x j ( t j ξ j ) q s j 2 ( t j ) p j ( t j ) d t j ,
θ s j 5 ( x j ) = ( 1 ) s 1 ξ j x j ( t j ξ j ) q s j 3 ( t j ) p j ( t j ) d t j , s = 1 , 2 ,
θ s j 1 ( x j ) = ( 1 ) s ξ j x j t j q s j 1 ( t j ) p j ( t j ) d t j , θ s j 2 ( x j ) = ( 1 ) s ξ j x j q s j 1 ( t j ) p j ( t j ) d t j ,
θ s j 3 ( x j ) = ( 1 ) s 1 ξ j x j t j q s j 2 ( t j ) p j ( t j ) d t j , θ s j 4 ( x j ) = ( 1 ) s 1 ξ j x j q s j 2 ( t j ) p j ( t j ) d t j ,
θ s j 5 ( x j ) = ( 1 ) s ξ j x j q s j 3 ( t j ) p j ( t j ) d t j , s = 3 .
Denote by Θ j s ( x j ) the following vector
Θ j s ( x j ) = [ θ 1 j s ( x j ) , θ 2 j s ( x j ) , θ 3 j s ( x j ) ] T .
For further purposes, it is convenient to introduce a matrix M = [ m i k ] with the following elements
m 1 5 ( j 1 ) + s = D 1 j ( 3 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 1 5 ( j 1 ) + s = 0 , j = m 1 , m , m 1 5 m + s = D 1 m + 1 ( 3 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
Hereon in, consider that index j changes only after the index s .
m 2 5 ( j 1 ) + s = D 2 j ( 3 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 2 5 ( j 1 ) + s = 0 , j = m 1 , m , m 1 5 m + s = D 2 m + 1 ( 3 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 3 5 ( j 1 ) + s = α j m + 1 D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 3 5 ( j 1 ) + s = 0 , j = m 1 , m , m 3 5 m + s = D 3 m + 1 ( 3 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 4 5 ( j 1 ) + s = α j m D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 4 5 ( j 1 ) + s = 0 , j = m 1 , m , m 4 5 m + s = D 3 m ( 3 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 5 5 ( j 1 ) + s = α j m 1 D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 5 5 ( j 1 ) + s = 0 , j = m 1 , m , m 5 5 m + s = D 3 m 1 ( 3 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 6 5 ( j 1 ) + s = η 1 m 1 α j m 1 D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 6 5 ( j 1 ) + s = 0 , j = m + 1 , m , m 6 5 ( m 2 ) + s = D 1 m 1 ( 2 ) ( Θ m 1 s , 0 ) , s = 1 , , 5 ,
m 7 5 ( j 1 ) + s = η 1 m α j m D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 7 5 ( j 1 ) + s = 0 , j = m + 1 , m 1 , m 7 5 ( m 1 ) + s = D 1 m ( 2 ) ( Θ m s , 0 ) , s = 1 , , 5 ,
m 8 5 ( j 1 ) + s = η 1 m + 1 α j m + 1 D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 8 5 ( j 1 ) + s = 0 , j = m , m 1 , m 8 5 m + s = D 1 m + 1 ( 2 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 9 5 ( j 1 ) + s = η 2 m α j m D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 9 5 ( j 1 ) + s = 0 , j = m + 1 , m , m 9 5 ( m 2 ) + s = D 2 m 1 ( 2 ) ( Θ m 1 s , 0 ) , s = 1 , , 5 ,
m 10 5 ( j 1 ) + s = η 2 m α j m D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 10 5 ( j 1 ) + s = 0 , j = m + 1 , m 1 , m 10 5 ( m 1 ) + s = D 2 m ( 2 ) ( Θ m s , 0 ) , s = 1 , , 5 ,
m 11 5 ( j 1 ) + s = η 2 m + 1 α j m + 1 D 3 j ( 2 ) ( Θ j s , 0 ) , j = 1 , , m 2 ,
m 11 5 ( j 1 ) + s = 0 , j = m , m 1 , m 11 5 m + s = D 2 m + 1 ( 2 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 ,
m 11 + p 5 ( j 1 ) + s = 0 , j = 1 , , p 1 , p + 1 , , m + 1 ,
m 11 + p 5 ( p 1 ) + s = η 1 p D 3 p ( 2 ) ( Θ p s , 0 ) D 1 p ( 2 ) ( Θ p s , 0 ) , s = 1 , , 5 , p = 1 , , m 2 ,
m 9 + m + p 5 ( j 1 ) + s = 0 , j = 1 , , p 1 , p + 1 , , m + 1 ,
m 9 + m + p 5 ( p 1 ) + s = η 2 p D 3 p ( 2 ) ( Θ p s , 0 ) D 2 p ( 2 ) ( Θ p s , 0 ) , s = 1 , , 5 , p = 1 , , m 2 ,
m 7 + 2 m + p 5 ( j 1 ) + s = 0 , j = 1 , , p 1 , p + 1 , , m , m 7 + 2 m + p 5 ( p 1 ) + s = D 1 p ( 0 ) ( Θ p s , 0 ) ,
m 7 + 2 m + p 5 m + s = D 1 m + 1 ( 0 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 , p = 1 , , m .
m 7 + 3 m + p 5 ( j 1 ) + s = 0 , j = 1 , , p 1 , p + 1 , , m , m 7 + 3 m + p 5 ( p 1 ) + s = D 2 p ( 0 ) ( Θ p s , 0 ) ,
m 7 + 3 m + p 5 m + s = D 2 m + 1 ( 0 ) ( Θ m + 1 s , l m + 1 ) , s = 1 , , 5 , p = 1 , , m .
m 7 + 4 m + p 5 ( j 1 ) + s = 0 , j = 1 , , p 1 , p + 1 , , m 2 ,
m 7 + 4 m + p 5 ( p 1 ) + s = D 3 p ( 0 ) ( Y s , ξ s ) η 1 p D 1 p ( 1 ) ( Θ p s , 0 ) η 2 p D 2 p ( 1 ) ( Θ p s , 0 ) ,
m 7 + 4 m + p 5 ( t 1 ) + s = α p t ( D 3 t ( 0 ) ( Θ t s , 0 ) η 1 t D 1 t ( 1 ) ( Θ t s , 0 ) η 2 t D 2 t ( 1 ) ( Θ t s , 0 ) ) , t = m 1 , m ,
m 7 + 4 m + p 5 m + s = α p m + 1 ( D 3 m + 1 ( 0 ) ( Θ m + 1 s , l m + 1 ) η 1 m + 1 D 1 m + 1 ( 1 ) ( Θ m + 1 s , l m + 1 )
η 2 m + 1 D 2 m + 1 ( 1 ) ( Θ m + 1 s , l m + 1 ) ) , s = 1 , , 5 , p = 1 , , m 2 .
We can now formulate the theorem on the reversibility of operator B .
Theorem 2.
If det M 0 , then a bounded operator B 1 exists in space L 2 ( Γ ) .
Proof of Theorem 2.
Consider the inhomogeneous operator equation L ( Y ) = F with conjugation conditions in the interior vertex C j ( Y ) = 0 ,   j = 1 , , 5 m + 5 and with rigid fixing conditions in the boundary vertices (5) on the star graph Γ . To prove Theorem 2, it suffices to prove the single-valued solvability of this boundary value problem for any right-hand side F L 2 ( Γ ) .
It is known that the following system of equations
D 1 j ( 2 ) ( Y j ; ξ ) = ξ j x j ( t j ξ j ) f 1 j ( t j ) d t j , D 2 j ( 2 ) ( Y j ; ξ ) = ξ j x j ( t j ξ j ) f 2 j ( t j ) d t j , D 3 j ( 2 ) ( Y j ; ξ ) = ξ j x j f 3 j ( t j ) d t j . .
with initial conditions
D 1 j ( 0 ) ( Y j ; ξ j ) = 0 , D 2 j ( 0 ) ( Y j ; ξ j ) = 0 , D 1 j ( 1 ) ( Y j ; ξ j ) = 0 , D 2 j ( 1 ) ( Y j ; ξ j ) = 0 , D 3 j ( 0 ) ( Y j ; ξ j ) = 0 .
is uniquely solvable. Denote the only solution to the above problem by Y j 0 ( x j ) ,   j = 1 , , m + 1 . It is a partial solution of the inhomogeneous operator equation L ( Y ) = F . The solution to the required problem is found at j = 1 , , m + 1 as follows
Y j ( x j ) = Y j 0 ( x j ) + K 1 j Θ j 1 ( x j ) + N 1 j Θ j 2 ( x j ) + K 2 j Θ j 3 ( x j ) + N 2 j Θ j 4 ( x j ) + K 3 j Θ j 5 ( x j ) .
We have to prove that the numbers K 1 j , N 1 j , K 2 j , N 2 j , K 3 j ,   j = 1 , , m + 1 are determined from the conjugation conditions in the inner vertex. If the conditions of Theorem 2 are fulfilled, a vector h is uniquely found from a system of linear algebraic equations M h = g for any g C 5 m + 5 . The elements of vector h can be interpreted as numbers K 1 j , N 1 j , K 2 j , N 2 j , K 3 j ,   j = 1 , , m + 1 . In this case, g is the numeric column that depends on partial solutions Y j 0 ( x j ) . Thus, the existence of a solution to the required problem is proved. The uniqueness of the solution follows from general statements about systems of linear differential equations. Since the coefficients in the system of linear differential equations represent continuous functions, it follows that the corresponding a priori estimates are correct. This means that the boundedness of the inverse operator B 1 . Theorem 2 is fully proved. □
Theorem 2 immediately implies the following corollary.
Corollary 1.
Let det M 0 . Then, for any F L 2 ( Γ ) and arbitrary constants φ i ,   i = 1 , , 5 m + 5 ,   ψ 5 ( j 1 ) + s ,   s = 1 , , 5 ,   j = 1 , , m + 1 the solution of the problem
L ( Y ) = F , C j ( Y ) = φ j , j = 1 , , 5 m + 5 ,
D 1 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 1 , D 2 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 2 ,
D 1 j ( 1 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 3 , D 2 j ( 1 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 4 ,
D 3 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 5 , j = 1 , , m + 1
exists and is unique.

4. Everywhere Solvable Reversible Boundary Value Problems for Systems of Differential Equations on a Star Graph

In this section, we describe all possible universally solvable reversible boundary value problems for the equation L ( Y ) = F . In the previous paragraphs, we wrote out conjugation conditions using linear forms C 1 ( · ) , , C 5 m + 5 ( · ) and fixation conditions using linear forms D 1 j ( 0 ) ( 0 ) ,   D 2 j ( 0 ) ( 0 ) ,   D 1 j ( 1 ) ( 0 ) ,   D 2 j ( 1 ) ( 0 ) ,   D 3 j ( 0 ) ,   j = 1 , , m + 1 . Under the above boundary conditions, the solvability and reversibility of the corresponding problem are proved everywhere.
How do we obtain new conjugation conditions and new anchoring conditions that guarantee all-around solvability and reversibility? This paragraph answers this question.
Let us choose arbitrarily two sets of linear bounded functionals in the function space L 2 ( Γ ) . The first set of linear functionals is denoted by
φ j ( · ) , j = 1 , , 5 m + 5 .
It is convenient to denote the second set of linear continuous functionals by
ψ 5 ( j 1 ) + s ( · ) , s = 1 , , 5 , j = 1 , , m + 1 .
Let us write a new boundary value problem corresponding to the chosen sets of linear continuous functionals.
Consider the inhomogeneous operator equation L ( Y ) = F with conjugation conditions in the interior vertex C j ( Y ) = φ j ( L ( Y ) ) ,   j = 1 , , 5 m + 5 and with fixation conditions in the boundary vertices on the star graph Γ
D 1 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 1 ( L ( Y ) ) , D 2 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 2 ( L ( Y ) ) , D 1 j ( 1 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 3 ( L ( Y ) ) , D 2 j ( 1 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 4 ( L ( Y ) ) , D 3 j ( 0 ) ( Y j ; ξ j ) = ψ 5 ( j 1 ) + 5 ( L ( Y ) ) , j = 1 , , m + 1 .
The written boundary value problem is everywhere solvable in the space L 2 ( Γ ) . The inversion of the problem also follows from Theorem 2. Indeed, if the expression L ( Y ) is replaced by F , in the boundary conditions, we find ourselves in the conditions of Theorem 2. Therefore, a solution to the new problem exists and is unique. The boundedness of the inverse operator in L 2 ( Γ ) space follows from the boundedness of the chosen functionals φ j ( · ) and ψ 5 ( j 1 ) + s ( · ) . Thus, we have proved the conclusion.
Theorem 3.
Let φ j ( · ) ,   j = 1 , , 5 m + 5 and ψ 5 ( j 1 ) + s ( · ) ,   s = 1 , , 5 ,   j = 1 , , m + 1 represent arbitrary sets of linear bounded functionals in space L 2 ( Γ ) . Then, the operator K , given by the differential expression K Y = L ( Y ) on the domain of definition
D ( K ) = { Y L 2 ( Γ ) , L ( Y ) L 2 ( Γ ) : C j ( Y ) = φ j ( L ( Y ) ) ,
j = 1 , , 5 m + 5 , c o n d i t i o n s ( 11 ) } .
has a bounded inverse in L 2 ( Γ ) space.
Remark 3.
The converse (in a sense) to Theorem 3 is true. However, we do not clarify this issue in this paper.

5. Examples

In this paragraph, we will give illustrative examples to show the meaning of Theorem 3. To do this, let us write down the general form of a linear bounded functional in the space L 2 ( Γ ) .
φ k ( F ) = < F , Φ j > L 2 ( Γ ) = s = 1 m + 1 0 l s F ( x s ) Φ k ( x s ) ¯ d x s
We substitute F with L ( Y ) and transform according to the Lagrange formula. In this case, the functions Y are considered to satisfy the rigid fixing conditions (5). Functionals are selected so that the ratios are true
L + ( Φ k ) = 0 , k = 1 , , 5 m + 5 .
We also assume that the functions Φ k ,   k = 1 , , 5 m + 5 also satisfy the rigid anchoring conditions (5). In this case, the value of the functional φ k ( F ) will be the same as the right-hand side of equality (9), where the vector of function Z should be replaced by the vector of function Φ k . According to Corollary 1, for each fixed k , we can arbitrarily choose values of linear forms C j ( Φ k ) ,   j = 1 , , 5 m + 5 . Let C j ( Φ k ) = 0 ,   j , k = 1 , , 5 m + 5 . Thus, the unambiguous choice of functions Φ k ,   k = 1 , , 5 m + 5 is fulfilled. As a result, we have
φ k ( L ( Y ) ) = D 1 m + 1 ( 0 ) ( Φ k m + 1 , l m + 1 ) C 3 m 1 ( Y 1 , , Y m , Y m + 1 )
+ D 2 m + 1 ( 0 ) ( Φ k m + 1 , l m + 1 ) C 3 m ( Y 1 , , Y m , Y m + 1 )
D 3 m ( 0 ) ( Φ k m , 0 ) C 3 m + 2 ( Y 1 , , Y m , Y m + 1 )
D 3 m 1 ( 0 ) ( Φ k m 1 , 0 ) C 3 m + 3 ( Y 1 , , Y m , Y m + 1 )
j = 1 m 2 D 1 j ( 1 ) ( Φ k j , 0 ) C 3 m + 9 + j ( Y 1 , , Y m , Y m + 1 )
j = 1 m 2 D 2 j ( 1 ) ( Φ k j , 0 ) C 4 m + 7 + j ( Y 1 , , Y m , Y m + 1 )
+ D 1 m 1 ( 1 ) ( Φ k m 1 , 0 ) C 3 m + 4 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m ( 1 ) ( Φ k m , 0 ) C 3 m + 5 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m + 1 ( 1 ) ( Φ k m + 1 , l m + 1 ) C 3 m + 6 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m 1 ( 1 ) ( Φ k m 1 , 0 ) C 3 m + 7 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m ( 1 ) ( Φ k m , 0 ) C 3 m + 8 ( Y 1 , , Y m , Y m + 1 )
+ D 1 m + 1 ( 1 ) ( Φ k m + 1 , l m + 1 ) C 3 m + 9 ( Y 1 , , Y m , Y m + 1 ) .
Here, some constants are denoted by D 1 m + 1 ( 0 ) ( Φ k m + 1 , l m + 1 ) + D 2 m + 1 ( 0 ) ( Φ k m + 1 , l m + 1 ) and so on. That is, the value of the functional φ k ( L ( Y ) ) represents a linear combination of the values of C j ( Y ) ,   j = 1 , , 5 m + 5 .
Thus, the conjugation condition C k ( Y ) φ k ( L ( Y ) ) = 0 at a fixed k after making the aforesaid choice of functions Φ k ,   k = 1 , , 5 m + 5 will be
j = 1 5 m + 5 a k j C j ( Y ) = 0 ,
where a k j are some constants.
Given the above reasoning, Theorem 3 will become as follows.
Theorem 4.
Let the set of numbers a k j ,   k , j = 1 , , 5 m + 5 be implemented in the above-mentioned way. Then, on the star graph Γ, the inhomogeneous operator equation L ( Y ) = F with conjugation conditions at the inner vertex
j = 1 5 m + 5 a k j C j ( Y ) = 0 , k = 1 , , 5 m + 5
and with fixing conditions (5) in the boundary vertices is uniquely solvable for any function F from the space L 2 ( Γ ) .
Thus, this example generalizes the conjugation conditions under fixed anchoring conditions. Now, we give another example where the anchoring conditions are generalized.
Let the vector function Y also satisfy the conjugation conditions.
C j ( Y ) = 0 , j = 1 , , 5 m + 5 ,
and otherwise be arbitrary. The set of functions Φ k ,   k = 1 , , 5 m + 5 satisfy the relation L + ( Φ k ) = 0 and the same conjugation conditions (12) and rigid anchoring conditions (5). In this case, the function Φ k is uniquely defined. That is, the function Φ k ,   k = 1 , , 5 m + 5 is the same as in Example 1.
In this case, the value of the functional φ k ( L ( Y ) ) is a linear combination of the form values.
D 1 j ( 0 ) ( Y j ; ξ j ) , D 2 j ( 0 ) ( Y j ; ξ j ) , D 1 j ( 1 ) ( Y j ; ξ j ) , D 2 j ( 1 ) ( Y j ; ξ j ) , D 3 j ( 0 ) ( Y j ; ξ j ) , j = 1 , , m + 1 .
As a result, a statement can be formulated.
Theorem 5.
There exists a set of numbers b k j s that depends on boundary values of functions Φ k ,   k = 1 , , 5 m + 5 . Then, on the star graph Γ, the inhomogeneous operator equation L ( Y ) = F with anchoring conditions in the boundary vertices
j = 1 m + 1 ( b k j 1 D 1 j ( 0 ) ( Y j ; ξ j ) + b k j 2 D 2 j ( 0 ) ( Y j ; ξ j ) + b k j 3 D 1 j ( 1 ) ( Y j ; ξ j )
+ b k j 4 D 2 j ( 1 ) ( Y j ; ξ j ) + b k j 5 D 3 j ( 0 ) ( Y j ; ξ j ) ) = 0 , k = 1 , , 5 m + 5 .
and with conjugation conditions (12) in the interior vertex is uniquely solvable for any function F from the space L 2 ( Γ ) .

6. Conclusions

The paper presents boundary value problems for systems of differential equations on a star graph. A distinctive feature of this paper is that the system consists of differential equations of different orders. In this process, it was necessary to work out a technique for making conjugation conditions at the nodes where several spatial rods are joined.
The result of this work can be used to calculate eigenvalues and eigenforms of spatially connected rod systems. Up to now, only plane eigenforms separately for longitudinal and for transverse vibrations have been used in engineering practice. In the present study, the determination of spatial eigenforms of vibrations of a rod system is proposed.
Furthermore, it is necessary to evaluate how much the eigenfrequencies and spatial eigenforms of the model problem on a graph approximate the original three-dimensional problem for a structure consisting of several rods. Similar asymptotic problems for the Laplace operator have been studied in [7,11]. In our case, instead of a Laplace operator, we consider equations of linear elasticity theory for a system of rods connected in a node. Similar tasks have been investigated in [19,20,21,22,23]. However, the results require further progress.
Note that this paper obtains new classes of conjugacy conditions in the inner vertex of a star graph in the case of systems of differential equations. In this case, the boundary conditions in the boundary vertices of the star graph are considered fixed. In particular, only the case of rigidly fixed boundary vertices is considered. In this case, the boundary conditions in the boundary vertices of the star graph are considered fixed. In particular, only the case of rigidly fixed boundary vertices is considered. We have not investigated other types of boundary fixings. Other types of boundary clauses for second-order differential equations are investigated in [10,11,12,13]. In the future, for systems of differential equations on graphs, we should study all possible kinds of boundary clauses.

Author Contributions

Introduction and canonical problem, B.K. and G.A.; proof theorem, B.K. and G.A.; writing—review and editing, G.A. and B.K.; project administration, B.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Education and Science of the Republic of Kazakhstan (Project AP08855402).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the financial support provided by the Ministry of Education and Science of the Republic of Kazakhstan.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are available from the authors.

References

  1. Nazarov, S.A.; Slutskii, A.S. One-dimensional equations of the deformation of thin slightly curved rods. Asymptotic analysis and justification. Izv. Math. 2000, 643, 531–562. [Google Scholar] [CrossRef]
  2. Andreev, V.I. Technical Mechanics (for Construction Universities and Faculties); Association of Building Education Publishers: Moscow, Russia, 2012; p. 256. [Google Scholar]
  3. Komkov, V. Optimal Control Theory for the Damping of Vibrations of Simple Elastic Systems; Lecture Notes in Mathemetics; Springer: Berlin/Heidelberg, Germany, 1972; Volume 253, Russian Translation: Mir, Moscow, 1975. (Foreword by T. Selezov). [Google Scholar]
  4. Timoshenko, S. On the correction factor for shear of the differential equation for transverse vibrations of bars of uniform cross-section. Philos. Mag. 1921, 41, 744–746. [Google Scholar] [CrossRef]
  5. Rubinstein, J.; Schatzman, M. Variational problems on multiply connected thin strips. I. Basic estimates and convergence of the Laplacian spectrum. Arch. Ration. Mech. Anal. 2001, 160, 271–308. [Google Scholar] [CrossRef]
  6. Kuchment, P.; Zeng, H. Convergence of spectra of mesoscopic systems collapsing onto a graph. J. Math. Anal. Appl. 2001, 258, 671–700. [Google Scholar] [CrossRef]
  7. Kuchment, P.; Zeng, H. Asymptotics of spectra of Neumann Laplacians in thin domains. In Advances in Differential Equations and Mathematical Physics, Contemporary Mathematics; American Mathematical Society: Providence, RI, USA, 2003; Volume 327, pp. 199–213. [Google Scholar]
  8. Exner, P.; Post, O. Convergence of spectra of graph-like thin manifolds. J. Geom. Phys. 2005, 54, 77–115. [Google Scholar] [CrossRef]
  9. Post, O. Spectral Analysis on Graph-like Spaces; Springer: Berlin/Heidelberg, Germany, 2012; pp. 292–366. [Google Scholar]
  10. Post, O. Branched quantum wave guides with Dirichlet boundary conditions:the decoupling case. J. Phys. Math. Gen. 2005, 38, 4917–4931. [Google Scholar] [CrossRef]
  11. Molchanov, S.; Vainberg, B. Scattering solutions in networks of thin fibers: Small diameter asymptotics. Commun. Math. Phys. 2007, 273, 533–559. [Google Scholar] [CrossRef]
  12. Grieser, D. Thin tubes in mathematical physics, global analysis and spectral geometry. arXiv 2008, arXiv:0802.2687v1. [Google Scholar]
  13. Grieser, D. Spectra of graph neighborhoods and scattering. Proc. Lond. Math. Soc. 2008, 9, 718–752. [Google Scholar] [CrossRef]
  14. Kamenskii, M.; Liou, Y.-C.; Wen, C.-F. Zvereva Margarita On a hyperbolic equation on a geometric graph with hysteresis type boundary conditions. J. Math. Program. Oper. Res. 2019, 69, 1029–4945. [Google Scholar]
  15. Kanguzhin, B.; Aimal Rasa, G.H.; Kaiyrbek, Z. Identification of the domain of the sturm–liouville operator on a star graph. Symmetry 2021, 13, 1210. [Google Scholar] [CrossRef]
  16. Nurakhmetov, D.; Jumabayev, S.; Aniyarov, A.; Kussainov, R. Symmetric Properties of Eigenvalues and Eigenfunctions of Uniform Beams. Symmetry 2020, 12, 2097. [Google Scholar] [CrossRef]
  17. Kanguzhin, B.E. Recovering of two-point boundary conditions by finite set of eigenvalues of boundary value problems for higher order differential equations. UFA. Math. J. 2020, 12, 22–29. [Google Scholar] [CrossRef]
  18. Nazarov, S.A. Substantiation of the asymptotic theory of thin rods. Integral and pointwise estimates. Probl. Math. Anal. 1997, 17, 101–152. [Google Scholar]
  19. Nazarov, S.A. Asymptotic Theory of Thin Plates and Rods.Volume 1: Dimensionality Reduction and Integral Estimates; Nauchnaya Kniga: Novosibirsk, Russia, 2001; p. 408. [Google Scholar]
  20. Exner, P.; Post, O. Convergence of resonances on thin branched quantum wave guides. J. Math. Phys. 2007, 48, 092104. [Google Scholar] [CrossRef]
  21. Maz’ya, V.; Nazarov, S.; Plamenevskij, B. Asymptotic Theory of Elliptic Boundary Value Problems in Singularly Perturbed Domains; Birkhauser Verlag: Basel, Switzerland, 2000. [Google Scholar]
  22. Nazarov, S.A. Korn’s inequalities for junctions of spatial bodies and thin rods. Math. Methods Appl. Sci. 1997, 3, 219–243. [Google Scholar] [CrossRef]
  23. Exner, P.; Post, O. Approximation of quantum graph vertex coupling by scaled Schrodinger operators on thin branched manifolds. J. Geom. Phys. 2009, 42, 415305. [Google Scholar] [CrossRef]
Figure 1. A 2-dimensional graph-like manifold with boundary.
Figure 1. A 2-dimensional graph-like manifold with boundary.
Symmetry 14 01761 g001
Figure 2. Limit star graph.
Figure 2. Limit star graph.
Symmetry 14 01761 g002
Figure 3. The case where the value of a node tends to zero much more slowly than the cross-sectional values of the rods.
Figure 3. The case where the value of a node tends to zero much more slowly than the cross-sectional values of the rods.
Symmetry 14 01761 g003
Figure 4. The case where the value of a node tends to zero much faster than the cross-sectional values of the rods.
Figure 4. The case where the value of a node tends to zero much faster than the cross-sectional values of the rods.
Symmetry 14 01761 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kanguzhin, B.; Auzerkhan, G. Conjugation Conditions for Systems of Differential Equations of Different Orders on a Star Graph. Symmetry 2022, 14, 1761. https://doi.org/10.3390/sym14091761

AMA Style

Kanguzhin B, Auzerkhan G. Conjugation Conditions for Systems of Differential Equations of Different Orders on a Star Graph. Symmetry. 2022; 14(9):1761. https://doi.org/10.3390/sym14091761

Chicago/Turabian Style

Kanguzhin, Baltabek, and Gauhar Auzerkhan. 2022. "Conjugation Conditions for Systems of Differential Equations of Different Orders on a Star Graph" Symmetry 14, no. 9: 1761. https://doi.org/10.3390/sym14091761

APA Style

Kanguzhin, B., & Auzerkhan, G. (2022). Conjugation Conditions for Systems of Differential Equations of Different Orders on a Star Graph. Symmetry, 14(9), 1761. https://doi.org/10.3390/sym14091761

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