Free Vibration Analysis of Arches with Interval-Uncertain Parameters

: The dynamic characteristics of a structure serve as a crucial foundation for structural assessment, fault diagnosis, and ensuring structural safety. Therefore, it is imperative to investigate the impacts of uncertain parameters on the dynamic performance of structures. The dynamic characteristics of arches with uncertain parameters are analyzed in this paper. The uncertain parameters are regarded as non-probabilistic uncertainties and represented as interval variables. A model of an arch considering interval-uncertain parameters is built, and kinematic equations are established. The natural frequencies are obtained using the differential quadrature (DQ) method, and the relationships between natural frequency, radius, and central angle are also analyzed. On this basis, the Chebyshev polynomial surrogate (CPS) model is employed to solve the uncertain dynamic problem, and the nat-ural frequencies are seen to be the objective functions of the CPS model. The accuracy veriﬁcation of the model is achieved by comparison with the Monte Carlo simulation (MCS). Simulations are carried out considering different uncertainties, and the results show that the bounds of natural frequencies are inﬂuenced not only by the types of uncertain parameters, but also by their combinations.


Introduction
Arch structures are often used in civil engineering, such as in the design of long-span bridges, roofs, and tunnels, due to their strong pressure-bearing capacity and unique shapes.The dynamic performance of arch beam structures has always attracted interest.The precise structural dynamic characteristics serve as crucial foundations for operational condition assessment and fault diagnosis, as well as ensuring structural safety [1][2][3][4][5][6][7][8][9].Traditional arch structures primarily comprise masonry arches, which are incapable of withstanding tension and possess an inadequate load-bearing capacity.To address this limitation, lighter and stronger materials such as steel and concrete-filled steel tubes are employed in modern arch structures.Additionally, the cross-section of arch structures is often reduced to achieve larger spans, resulting in arch beam structures.Consequently, this study focuses on investigating steel arch beam structures.The free vibration of linear elastic arch beams is a crucial aspect in the investigation of dynamic performance, typically evaluated using the finiteelement method.In earlier studies, a simplified model neglecting the moment of inertia and shear deformation was commonly employed, resulting in a singular sixth-order partial differential equation that incurred high computational costs and inaccuracies.To achieve enhanced accuracy with reduced computational expense, partial differential equations were established for the free vibration of an arch model by considering the axial deformation, tangential deformation, and moment of inertia, and the differential quadrature method for numerical computation was subsequently employed.The differential quadrature (DQ) method has been demonstrated to be a straightforward and efficient approach for solving differential equations, yielding highly accurate numerical solutions with minimal grid Appl.Sci.2023, 13, 12391 2 of 16 points.It is particularly well-suited for structures exhibiting a curvilinear geometry, such as arch beam structures [10][11][12][13][14][15][16][17][18][19].
Uncertainties are common in the field of engineering; for example, they are frequently encountered in manufacturing tolerances, assembly deviation, dispersion in materials, and so on.Even a slight parametric variation may lead to a drastic change in the natural frequency.Due to the inherent parametric uncertainties associated with arch structures, occasional variations occur in their structural dynamic characteristics which can potentially impact both the robustness and the accuracy of structural diagnoses and assessments.Therefore, it is important to investigate the effect of uncertain parameters on the dynamic characteristics of arch structures.Generally, the probabilistic model [20,21] and non-probabilistic models such as the fuzzy set model [22,23] and the interval model [24][25][26][27][28][29][30][31][32] are employed to characterize the uncertainty inherent in structural systems.The probabilistic model requires the probability distribution function, whereas the fuzzy set model necessitates the membership function.The aforementioned approaches are frequently high-cost and occasionally challenging to accomplish.To address this issue, the interval model has been proposed as a valuable supplement.The interval model only requires the boundaries that define the extreme values of uncertainties, rather than their accuracy distribution.Additionally, there is no need to clarify the evolution in varying ranges, making the model easy to implement.To date, there have been limited studies about the free-vibration characteristics of arches while considering the interval uncertainty.Feng et al. used a probability field model, including interval fields, to represent spatially correlated uncertain parameters and analyze the non-deterministic free vibration of CFST arch structures [33].Wu et al. considered the uncertain material properties and creep effect of CFST arch structures, established a unified interval field model, and determined the limit range of the natural frequency of an arch using the interval finite element method [34].The aforementioned methodologies are considered intrusive and require modification of the original deterministic models, resulting in complex processes.In contrast, nonintrusive interval methods utilize orthogonal polynomial approximation theory combined with interval operations to overcome the limitation of interval solutions having a wide range [27][28][29][30][31][32][35][36][37][38][39][40][41][42][43].
In this paper, the free-vibration differential equations of a linear elastic arch are initially established using the traditional approach, and the natural frequencies are determined by employing the DQ method, which has been shown to be more efficient and accurate.Then, the uncertain parameters are treated as interval variables, and the non-intrusive Chebyshev polynomial surrogate (CPS) model is for the first time employed to address the natural frequencies of arches with uncertainties.The natural frequencies are regarded as the objective functions of the surrogate model, and the CPS model is established based on the original deterministic value of these frequencies.Finally, the model is utilized to obtain the bounds of natural frequencies.In addition, the Monte Carlo simulation (MCS) is employed as a means of validating the accuracy of the CPS model.
The remainder of this paper is organized as follows.The derivation of the freevibration equation and the determination of the natural frequency are elaborated upon in Section 2. The resolution of the natural frequency boundary through the CPS model is addressed in Section 3. Several numerical examples are analyzed in Section 4. In Section 5, the conclusions are summarized.

Establishment of Free-Vibration Equations
The model of the arch is shown in Figure 1.Only the elastic deformation and the in-plane vibration are considered in the model.The term ϑ denotes the central angle formed between the radius of any cross-section and the Y axis, which is positively oriented in the clockwise direction; ϑ∈[−ϑ 0 ,ϑ 0 ], where ϑ 0 is half of the central angle of the whole arch.R is the radius and u(ϑ,t), v(ϑ,t), and ϕ(ϑ,t) are the tangential displacement, radial displacement, and bending slope of the cross-section at a given central angle ϑ and arbitrary time t, respectively.

Establishment of Free-Vibration Equations
The model of the arch is shown in Figure 1.Only the elastic deformation and the inplane vibration are considered in the model.The term ϑ denotes the central angle formed between the radius of any cross-section and the Y axis, which is positively oriented in the clockwise direction; ϑ∈ [−ϑ0,ϑ0], where ϑ0 is half of the central angle of the whole arch.R is the radius and u(ϑ,t), v(ϑ,t), and φ(ϑ,t) are the tangential displacement, radial displacement, and bending slope of the cross-section at a given central angle ϑ and arbitrary time t, respectively.The kinematic equation can be written as follows [6]: , where ρ is the mass density for unit volume; A and I are the area and moment of inertia of the cross-section, respectively.N(ϑ,t), Q(ϑ,t), and M(ϑ,t) represent the axial force, shear force, and bending moment, respectively, and these internal forces can also be expressed by the constitutive relation of displacement as: The kinematic equation can be written as follows [6]: where ρ is the mass density for unit volume; A and I are the area and moment of inertia of the cross-section, respectively.N(ϑ,t), Q(ϑ,t), and M(ϑ,t) represent the axial force, shear force, and bending moment, respectively, and these internal forces can also be expressed by the constitutive relation of displacement as: in which E is Young's modulus and G is shear modulus; κ 0 is the shear factor.Substituting the relations in (2) into Equation (1), and with u(ϑ,t), v(ϑ,t), and ϕ(ϑ,t) indicated as u, v, and ϕ, the equations can be expressed as follows: Equation ( 3) comprises the in-plan free-vibration equations.In general, the separation of variables can be employed to seek solutions that exhibit temporal harmony with frequency ω, so the displacement can be expressed as: U(ϑ), V(ϑ), and Φ(ϑ) are indicated as U, V, and Φ, which represent the tangential, radial, and rotational vibration amplitude values of the structure, respectively.Substituting the relations in (4) into Equation (3), we obtain: The equations describing the boundary conditions are written as follows: Clamped edge boundary U = 0, V = 0, Φ = 0 at ϑ = ±ϑ 0 (6)

Establishment of Free-Vibration Equations
In this subsection, our primary objective is to derive a solution for the system of ordinary differential equations.Directly integrating the equations incurs significant complexity, thus necessitating the discovery of a concise and efficient approach.The differential quadrature (DQ) method is a highly effective approach for solving partial and ordinary differential equations.Its fundamental principle involves partitioning the domain of independent variables for smooth functions into discrete sample points according to specific rules and approximating the derivative of the function with respect to these independent variables as a weighted sum of their corresponding function values [11][12][13].Therefore, the nth order derivative of function f (ϑ) with respect to ϑ at the ith point ϑ i is approximated as follows: where ξ (n) ij denotes the nth-order derivative weighting coefficients of ϑ at the ith sample point calculated for the jth sample points in the domain, H is the total number of sample points, and f (ϑ i ) is the function value at point ϑ i .The weighting coefficients can be determined by Lagrange polynomial functions, and the first-order derivative weighting coefficients are: Appl.Sci.2023, 13, 12391 5 of 16 The second-and higher-order derivative weighting coefficients are: Each order of the derivative weighting coefficients can be obtained according to the above formulas.It is evident that the computation of derivative weight coefficients beyond the second order necessitates prior acquisition of first-order weight coefficients.The selection of appropriate sample points serves as a fundamental basis for ensuring both accuracy and efficiency in DQ method calculations.The Chebyshev-Gauss-Lobatto (C-G-L) grid distribution [11] is employed in this paper.According to the C-G-L rule, the sample point r i on the standard interval [0,1] is expressed as: and its values can be established by mapping the relationship with ϑ i as follows: The required sample points for this method are now available.After selecting a specific number of sample points based on the calculation requirements, the derivative weighting coefficient can be obtained through relations ( 8) and (9).
Using the DQ method, the arch is discretized according to sample points ϑ i, and the propositions of the original Equation ( 5) can be rewritten as: where U i , V i , and Φ i represent the tangential, radial, and rotational vibration amplitude values, respectively, at θ i , i = 1, 2,. ..H.The set of equations can be reduced to in which [K] and [M] are the global stiffness matrix and mass matrix, respectively, and 6) can be transformed into: Considering the boundary condition, Equation ( 13) can be written as The subscripts 'b' and 'd' in the matrix denote the components related to the boundary and internal domain, respectively, in Equation (15).The preservation of only the internal domain displacement vector is employed to streamline computation and enhance efficiency.Now Equation ( 15) can be expressed as follows: 16) can be written as: The natural frequencies can be obtained by the eigenvalue decomposition of Equation ( 17).

General Polynomial Surrogate Model for Interval Natural Frequency Analysis
In this section, the uncertain natural frequency bounds of an arch will be analyzed with Chebyshev interval algorithms.It has been demonstrated that for a continuous function f (x) defined in [a, b], there always exists a polynomial g(x) that approximates it with an error no larger than an arbitrary small positive quantity ε.On this basis, the Chebyshev polynomial surrogate (CPS) model is utilized to deal with uncertain dynamic problems; the frequencies are treated as the proposed function f (x) and can be approximated by orthogonal polynomials with acceptable errors in the following.The original dynamic questions can be modeled as simple mathematical functions by using a few interpolation points, which can save the computational cost of the traditional Monte Carlo simulation (MCS) without losing accuracy.As mentioned previously, the geometric dimensions of the cross-sectional area and physical parameters such as the Young's modulus may be uncertain.Due to the lack of necessary information, it is a complicated task to model the uncertainties probabilistically.However, only their bounds are needed in the interval approach, which can be achieved easily.Therefore, it is more reasonable to regard the uncertainties as interval variables.An uncertain parameter is regarded as an interval parameter q I , where superscript I denotes the interval variable, and can be denoted as: in which q L , q U , and q C denote the lower bound, upper bound, and nominal value of the interval parameter q I , respectively.Deviation degree β denotes the deviation coefficients.R is the set of all real numbers.If there is more than one uncertain parameter, the expression above can be expressed in the form of multi-dimensional interval vectors, as: q I is the uncertain parameters vector, where i = 1, 2,. ..r. r is the total number of uncertain parameters, and the parameters are assumed to be mutually independent in the following.Considering the matrices in the free-vibration equations for the uncertain arches as the function of the uncertain parameter, Equation ( 17) can be rewritten in the interval form The goal of this study is to determine the interval range of natural frequencies, and a surrogate model based on Chebyshev orthogonal polynomial approximation theory is introduced to achieve this.First, the one-dimensional uncertain model is analyzed.The first type of Chebyshev orthogonal polynomial is where n is an arbitrary positive integer.T n (x) is orthogonal on [−1,1] with respect to the weight function ρ(x) = (1 − x 2 ) 1/2 , and its recurrence formulas are as follows: The one-dimensional n-order Chebyshev polynomial approximation of function f (x) is formulated as: in which a i is the expansion coefficient and where f (x) is the deterministic value of f (x) and n is the expansion order.In general, the direct integration of a i poses significant complexities.However, when compared to this approach, employing the Gauss-Chebyshev integration formula proves to be a more efficient and expeditious method for obtaining a i .The Gauss-Chebyshev integration formula is a type of interpolation integral formula expressed as follows: in which m is the number of integral points and should satisfy m ≥ n + 1.
In the Gauss-Chebyshev integration, the interpolation points involved in the integration process are the zeros of m-degree Chebyshev polynomials, which are given by: For Equation (24), the weight function (1 − x 2 ) 1/2 of the Chebyshev orthogonal polynomial is chosen as ρ(x), and in this case, all the integral coefficients A k = π/m [27].Substituting this into Equation ( 24), the expansion coefficient a i can be expressed as Substituting this into Equation ( 22), we obtain Now the polynomial surrogate model function based on Chebyshev orthogonal polynomials can be constructed.It is evident that once the original frequency values at interpolation points, as given by Equation (25), are obtained from the deterministic model in Equation ( 17), the expansion coefficient a i can be calculated immediately, thereby simultaneously obtaining the surrogate-model function.The upper and lower bounds of the function being approximated can be obtained by searching the maximum and minimum values of g n (x) on standard interval [−1,1], which can be easily achieved based on the maximum theorem of a continuous function on closed intervals.Moreover, linear transformation does not change the relative locations of the extreme values.Considering the arbitrary interval distribution of the uncertain parameter x, and the fact that the Chebyshev polynomial surrogate model is applicable only for the standard interval [−1,1], it is necessary to transfer the parameter to a standard interval: The roots of the derivative equation are: The extremal value points x min and x max , which correspond to the uncertain parameter x, are determined by comparing the values of g n (x r ), g n (−1), and g n (1).This implies that the natural frequencies attain their maximum and minimum values at these points.The corresponding locations of x will be marked as x min and x max .
Supposing that multiple uncertain parameters are independent, for r-dimensional uncertain parameter vector q I , the process of searching the extremum should be performed r times to obtain the maximum and minimum frequency value locations for each parameter.Two extremum vectors will be obtained: Accordingly, the two deterministic parameter vectors that account for the bounds of the frequencies will be defined as follows: where q C , q U , and q L represent the nominal values, upper bounds, and lower bounds of the uncertain parameters, respectively, which are all deterministic.Substituting Equation (31) into Equation ( 17), the natural-frequency bounds can be obtained from two differential equations: K(q min ) {δ d (q min , ϑ)} + ω(q min ) 2 M(q min ) {δ d (q min , ϑ)} = 0 The solutions to Equation ( 32) are all deterministic and of the same size as the original equation denoted by Equation (17), and the equations can be solved easily.The whole computational procedure of the natural frequency bounds obtained using the CPS model is illustrated in Figure 2.

Numerical Simulation and Discussion
In this section, a circular arch with interval uncertainties will be simulated numerically.The first fourth-order deterministic natural frequencies are presented in the first subsection with the deterministic parameters, using the DQ method.Then, the accuracy and efficiency of the proposed CPS model in computing the natural frequencies are validated by comparison with the MCS in the second subsection.In the third subsection, different uncertain parameters such as the cross-sectional height and Young's modulus are considered, and the corresponding bounds of the natural frequencies are shown.The numerical simulations are conducted using the MATLAB R2014a software.

Numerical Simulation and Discussion
In this section, a circular arch with interval uncertainties will be simulated numerically.The first fourth-order deterministic natural frequencies are presented in the first subsection with the deterministic parameters, using the DQ method.Then, the accuracy and efficiency of the proposed CPS model in computing the natural frequencies are validated by comparison with the MCS in the second subsection.In the third subsection, different uncertain parameters such as the cross-sectional height and Young's modulus are considered, and the corresponding bounds of the natural frequencies are shown.The numerical simulations are conducted using the MATLAB R2014a software.

The Natural Frequencies Based on the Traditional Deterministic Method
According to the calculation example provided in reference [19], the arch cross-section is assumed to be rectangular; its physical parameters are presented in Table 1.In this paper, the total number of discretized sample points of arch H is 21.Waterfall plots of the first fourth-order natural frequencies with the central angle varying in the interval [π/6, π] and the radius varying in the interval [1,3] are shown in Figure 3. Firstly, it is shown that, under the conditions of equal radius and central angle, the natural frequency demonstrates a discernible upward trend with increasing orders, which aligns convincingly with empirical evidence.Then, it is also shown that the frequency of each order increases as the radius and central angle decrease, with a more pronounced increase observed in higher-order frequencies compared to lower-order frequencies.Furthermore, the frequency remains numerically continuous within the range of the central angle, thus confirming the appropriateness of the surrogate model utilized in this paper.Specifically, it is essential for the objective function of the model to maintain continuity.

The Natural Frequencies Based on the Traditional Deterministic Method
According to the calculation example provided in reference [19], the arch cross-section is assumed to be rectangular; its physical parameters are presented in Table 1.In this paper, the total number of discretized sample points of arch H is 21.Waterfall plots of the first fourth-order natural frequencies with the central angle varying in the interval [π/6, π] and the radius varying in the interval [1,3] are shown in Figure 3. Firstly, it is shown that, under the conditions of equal radius and central angle, the natural frequency demonstrates a discernible upward trend with increasing orders, which aligns convincingly with empirical evidence.Then, it is also shown that the frequency of each order increases as the radius and central angle decrease, with a more pronounced increase observed in higher-order frequencies compared to lower-order frequencies.Furthermore, the frequency remains numerically continuous within the range of the central angle, thus confirming the appropriateness of the surrogate model utilized in this paper.Specifically, it is essential for the objective function of the model to maintain continuity.

Comparisons of CPS Model with MCS in Calculating Natural Frequencies
In order to assess the accuracy and efficiency of the CPS model in estimating an uncertain natural frequency, the MCS will be used for comparison.A sufficient number of samples is the guarantor ensuring the accuracy of the simulation in the MCS, and 500 samples are utilized in this MCS.The expansion order n of the CPS model takes the value 3 and the number of the integral point m takes the value represent the upper bound and lower bound of the ith natural frequency, respectively.Variables e U i and e L i represent the upper bound and lower bound of the ith frequency changing rate, respectively, which can be regarded as an indication of the fluctuation in natural frequency and can be formulated as Accounting for inevitable manufacturing tolerances, the cross-sectional height is considered as an uncertainty parameter, with a deviation degree of 5%.The radius takes values of 1 m and 3 m, while the central angle still varies in the interval [π/6, π] for legibility.Figures 4 and 5 show the comparisons of changing rates of the first fourth-order uncertain natural frequencies.It is shown that the results obtained from both models exhibit excellent agreement in terms of the changing rate for each order's frequencies, thereby formalizing the accuracy of the proposed interval model and its applicability in determining frequency intervals.Furthermore, the CPS model demonstrates a computational time that is approximately 85% shorter than that of MCS, thus confirming its high efficiency.Additionally, it is observed that when maintaining an equal radius, there exist fluctuation characteristics and a minimum value in the relationship between the changing rate and central angle.Moreover, as the radius increases, the position of the minimum value shifts towards the left along the horizontal axis representing the central angle.
the upper bound and lower bound of the ith natural frequency, respectively.Variables e U i and e L i represent the upper bound and lower bound of the ith frequency changing rate, respectively, which can be regarded as an indication of the fluctuation in natural frequency and can be formulated as Accounting for inevitable manufacturing tolerances, the cross-sectional height is considered as an uncertainty parameter, with a deviation degree of 5%.The radius takes values of 1 m and 3 m, while the central angle still varies in the interval [π/6, π] for legibility.Figures 4 and 5 show the comparisons of changing rates of the first fourth-order uncertain natural frequencies.It is shown that the results obtained from both models exhibit excellent agreement in terms of the changing rate for each order's frequencies, thereby formalizing the accuracy of the proposed interval model and its applicability in determining frequency intervals.Furthermore, the CPS model demonstrates a computational time that is approximately 85% shorter than that of MCS, thus confirming its high efficiency.Additionally, it is observed that when maintaining an equal radius, there exist fluctuation characteristics and a minimum value in the relationship between the changing rate and central angle.Moreover, as the radius increases, the position of the minimum value shifts towards the left along the horizontal axis representing the central angle.

Single Uncertain Parameter
In the first part of this subsection, a single uncertain parameter is considered.Parts of the numerical simulation results are illustrated regarding the effect of the uncertain cross-sectional height in the comparative analysis previously presented, and here the width of the cross-section is considered with the same deviation degree.Taking into consideration material dispersion, Young's modulus and the density of mass are regarded as uncertain parameters, with respective deviation degrees of 5%.The simulation results are shown in Figure 6.It is found that the bounds of the changing rate are drawn as straight lines in the presence of each of the uncertainties above, and further investigations demonstrate that these phenomena remain unaffected by variations in radius.Furthermore, the natural frequencies are minimally influenced by the uncertain width of the cross-section; however, they exhibit greater sensitivity to variations in the uncertain cross-sectional height compared with Young's modulus and the density of mass.It can be seen that the physical parameters, including Young's modulus and the density of mass, demonstrate a linear relationship with the natural frequencies through eigenvalue decomposition of the free-vibration equation.Additionally, the cross-sectional width can be canceled on both sides of the equation.Moreover, a nonlinear correlation exists between the cross-sectional height and the moment of inertia, subsequently affecting the natural frequency.This observation elucidates the distinct influences of each parameter on the uncertain dynamic characteristics.In the first part of this subsection, a single uncertain parameter is considered.Parts of the numerical simulation results are illustrated regarding the effect of the uncertain cross-sectional height in the comparative analysis previously presented, and here the width of the cross-section is considered with the same deviation degree.Taking into consideration material dispersion, Young's modulus and the density of mass are regarded as uncertain parameters, with respective deviation degrees of 5%.The simulation results are shown in Figure 6.It is found that the bounds of the changing rate are drawn as straight lines in the presence of each of the uncertainties above, and further investigations demonstrate that these phenomena remain unaffected by variations in radius.Furthermore, the natural frequencies are minimally influenced by the uncertain width of the cross-section; however, they exhibit greater sensitivity to variations in the uncertain cross-sectional height compared with Young's modulus and the density of mass.It can be seen that the physical parameters, including Young's modulus and the density of mass, demonstrate a linear relationship with the natural frequencies through eigenvalue decomposition of the free-vibration equation.Additionally, the cross-sectional width can be canceled on both sides of the equation.Moreover, a nonlinear correlation exists between the cross-sectional height and the moment of inertia, subsequently affecting the natural frequency.This observation elucidates the distinct influences of each parameter on the uncertain dynamic characteristics.

Multiple Uncertain Parameters
In order to investigate the characteristics of frequency under the influence of multiple uncertain parameters, all the aforementioned uncertainties are considered simultaneously.However, employing the MCS requires a substantial number of samples for each parameter in order to ensure accuracy, resulting in unacceptable computational costs.Utilizing the CPS model can significantly reduce these costs and provide notable efficiency advantages.The numerical simulation results depicting multi-dimensional uncertainty are presented in Figures 7 and 8

Multiple Uncertain Parameters
In order to investigate the characteristics of frequency under the influence of multiple uncertain parameters, all the aforementioned uncertainties are considered simultaneously.However, employing the MCS requires a substantial number of samples for each parameter in order to ensure accuracy, resulting in unacceptable computational costs.Utilizing the CPS model can significantly reduce these costs and provide notable efficiency advantages.The numerical simulation results depicting multi-dimensional uncertainty are presented in Figures 7 and 8

Multiple Uncertain Parameters
In order to investigate the characteristics of frequency under the influence of multiple uncertain parameters, all the aforementioned uncertainties are considered simultaneously.However, employing the MCS requires a substantial number of samples for each parameter in order to ensure accuracy, resulting in unacceptable computational costs.Utilizing the CPS model can significantly reduce these costs and provide notable efficiency advantages.The numerical simulation results depicting multi-dimensional uncertainty are presented in Figures 7 and 8.The changing rates of frequencies surpass those observed when considering individual parameters separately, indicating a linear superposition of the effects of all uncertainties.Furthermore, although resembling single cross-sectional height considerations in terms of image shape, fluctuations exhibit smaller amplitudes.Additionally, as the radius increases, the phenomenon that the position of the minimum value shifts to the left along the horizontal axis representing the central angle does not change.

Conclusions
The linear natural frequencies of arches under interval-uncertain parameters were analyzed in this study.By coupling the DQ method and the Chebyshev approximation theory, a surrogate model was established to solve the uncertainty problem with respect to the dynamic characteristics of arches.The non-intrusive interval method is valid in freevibration analysis, for which probabilistic methods cannot be utilized due to a lack of essential prior information.The accuracy and effectiveness of the proposed CPS model were demonstrated via a comparative study of the natural frequency ranges of the arch obtained using the MCS and the model.It is clear that the CPS model is more efficient, and this will be more obvious in multiple-uncertainty situations, because the sample points in multiple dimensions are enormous, while the interpolation numbers used in the CPS model remain acceptable.Based on the results of several cases exhibited in the foregoing sections, the following conclusions can be made: 1. Some geometrical parameters, such as the cross-sectional height, were found to be sensitive, which will cause large fluctuations in the dynamic properties.2. The influence of most uncertain physical parameters on the natural frequencies exhibits a linear relationship which remains unaffected by variations in radius and central angle.However, different parameters yield distinct rates of change in frequency.3.In the case of multi-dimensional uncertainty, the changing rate of frequencies surpasses that of any individual parameter, as it represents a linear superposition of the combined influence exerted by each parameter.
The analysis method employed in this study is non-intrusive and convenient, and has been demonstrated to be an efficient approach for addressing the free-vibration problems of arches with uncertainties.It is widely accepted that, despite the uncertain parameters, there should be a linear relationship between the natural frequency of the linear elastic arch and the deviation degree of an uncertain parameter.However, this study has

Conclusions
The linear natural frequencies of arches under interval-uncertain parameters were analyzed in this study.By coupling the DQ method and the Chebyshev approximation theory, a surrogate model was established to solve the uncertainty problem with respect to the dynamic characteristics of arches.The non-intrusive interval method is valid in free-vibration analysis, for which probabilistic methods cannot be utilized due to a lack of essential prior information.The accuracy and effectiveness of the proposed CPS model were demonstrated via a comparative study of the natural frequency ranges of the arch obtained using the MCS and the model.It is clear that the CPS model is more efficient, and this will be more obvious in multiple-uncertainty situations, because the sample points in multiple dimensions are enormous, while the interpolation numbers used in the CPS model remain acceptable.Based on the results of several cases exhibited in the foregoing sections, the following conclusions can be made: 1.
Some geometrical parameters, such as the cross-sectional height, were found to be sensitive, which will cause large fluctuations in the dynamic properties.

2.
The influence of most uncertain physical parameters on the natural frequencies exhibits a linear relationship which remains unaffected by variations in radius and central angle.However, different parameters yield distinct rates of change in frequency.

3.
In the case of multi-dimensional uncertainty, the changing rate of frequencies surpasses that of any individual parameter, as it represents a linear superposition of the combined influence exerted by each parameter.
The analysis method employed in this study is non-intrusive and convenient, and has been demonstrated to be an efficient approach for addressing the free-vibration problems of arches with uncertainties.It is widely accepted that, despite the uncertain parameters, there should be a linear relationship between the natural frequency of the linear elastic arch and the deviation degree of an uncertain parameter.However, this study has revealed intriguing variations in the influences of different parameters on the natural frequency, sometimes exhibiting elusive changes.The methodology employed can be extended to calculate both natural frequencies and dynamic responses for nonlinear arches in future studies, which is expected to yield more valuable insights.

Figure 1 .
Figure 1.Model of the arch.

Figure 1 .
Figure 1.Model of the arch.

Figure 2 .
Figure 2. The computational procedure of the natural-frequency bounds via the CPS model.

Figure 2 .
Figure 2. The computational procedure of the natural-frequency bounds via the CPS model.

Figure 3 .Figure 3 .
Figure 3. Waterfall diagrams of the first fourth-order natural frequencies of the arch versus the central angle with different radii: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 4 .
Figure 4.The changing rates of the first fourth-order natural frequencies of the arch with uncertain cross-sectional height and a radius of 1 m using the CPS model and MCS: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 4 .
Figure 4.The changing rates of the first fourth-order natural frequencies of the arch with uncertain cross-sectional height and a radius of 1 m using the CPS model and MCS: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 5 .
Figure 5.The changing rates of the first fourth-order natural frequencies of the arch with uncertain cross-sectional height and a radius of 3 m using the CPS model and MCS: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 5 .
Figure 5.The changing rates of the first fourth-order natural frequencies of the arch with uncertain cross-sectional height and a radius of 3 m using the CPS model and MCS: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 6 .
Figure 6.The changing rates of the first fourth-order natural frequencies of an arch with uncertain cross-sectional width, Young's modulus, and density of mass, respectively: (a) first order; (b) second order; (c) third order; and (d) fourth order.
. The changing rates of frequencies surpass those observed when considering individual parameters separately, indicating a linear superposition of the effects of all uncertainties.Furthermore, although resembling single cross-sectional height considerations in terms of image shape, fluctuations exhibit smaller amplitudes.Additionally, as the radius increases, the phenomenon that the position of the minimum value shifts to the left along the horizontal axis representing the central angle does not change.

Figure 6 .
Figure 6.The changing rates of the first fourth-order natural frequencies of an arch with uncertain cross-sectional width, Young's modulus, and density of mass, respectively: (a) first order; (b) second order; (c) third order; and (d) fourth order.

17 Figure 6 .
Figure 6.The changing rates of the first fourth-order natural frequencies of an arch with uncertain cross-sectional width, Young's modulus, and density of mass, respectively: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 7 .
Figure 7.The changing rates of the first fourth-order natural frequencies of an arch with multiple simultaneously uncertain parameters and a radius of 1 m: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 7 .
Figure 7.The changing rates of the first fourth-order natural frequencies of an arch with multiple simultaneously uncertain parameters and a radius of 1 m: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 8 .
Figure 8.The changing rates of the first fourth-order natural frequencies of an arch with multiple simultaneously uncertain parameters and a radius of 3 m: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Figure 8 .
Figure 8.The changing rates of the first fourth-order natural frequencies of an arch with multiple simultaneously uncertain parameters and a radius of 3 m: (a) first order; (b) second order; (c) third order; and (d) fourth order.

Table 1 .
Geometrical and physical parameters of arch.

Table 1 .
Geometrical and physical parameters of arch.