Abstract
The purpose of this present work is to investigate the geometrical nonlinearity in free and forced vibration of a shallow arch elastically restrained at the ends. The nonlinear governing equilibrium equation of the shallow arch is obtained after the Euler Bernoulli theory and the Von Karman geometrical nonlinearity assumptions. After applying the ends conditions, the eigenvalues problem of the generalized trancendant equation have been determined iteratively using the NewtonRaphson algorithm. The kinetic and total strain energy have been discretized into a series of a finite spatial functions which are a combination of linear modes and basic function contribution coefficients. Using Hamilton’s principle energy and spectral analysis, the problem is reduced into a set of nonlinear algebraic equations that solved numerically using an approximate explicit method developed previously the socalled second formulation. Considering a multimode approach, the effect of initial rise and concentrated force on nonlinear behaviour of system has been illustrated in the backbone curves giving the nonlinear amplitudefrequency dependence. The corresponding nonlinear deflections and curvatures have been plotted for various vibration amplitudes.
1. Introduction
Investigating the geometrical nonlinearity is one of the major consideration on the design process of the arches. The study of the geometrical nonlinearity of a beams, plates ,shells and circular arches were investigated by the authors of [18]. The authors classified the arches following their shallowness ratio into two classes; shallow arches and deep or noshallow arches. The present study will be mainly focused on the nonlinear free and forced vibration of a sinusoidal shallow arch.
Many researchers have studied the static and dynamic of shallow arch, experimentally and theoretically. Raymond H. Plaut et al. [9] determined the form of the arch which maximizes the fundamental vibration frequency.
Nayfeh et al. [10] obtained an exact solution of the eigenvalue problems governing linear undumped free vibrations of buckled beam, W. Lacarbona and G. Rega [11] obtained the general conditions for orthogonality of the nonlinear modes in case of twotoone, twotothree, and onetoone resonance. Using the multiple scales method, Yi Z., Stanciulescu [12] constructed the NNMs by considering the quadratic and cubic nonlinearities for a shallow arch with one end elastically restrained in vertical and rotational directions in the case of twotoone internal resonances.
George C. Tsiatas and Nick G. Babouskos [13] investigated the linear and geometrically nonlinear response of nonuniform shallow arches under a central concentrated force and they assume that the shallow arches have increased overall stiffness contrary to straight beams. The harmonic balance method and pseudo arclength method was applied by SiQin Ye et al. [14] to study nonlinear vibration of curved structure with nonlinear boundary conditions. The author of [15] used a technique that combining harmonic balance coupled with continuation scheme to determine the solutions for an entire range of externally applied loading. The effects of a low frequency parametric excitation for a double hinged shallow arch with lateral and sinusoidal chargement of arch have been investigated by F. Lakrad, W. Schiehlen [16]. E. Ozkaya et al were analyzed the effect of different elastic springs coefficients and locations on nonlinear transverse vibration of curved beam using the multiple scale method [17]. All these authors confirmed that linear and nonlinear frequencies and modes shapes are sensitive to the linear initial curvature.
In this present paper, the analyze of the geometrical nonlinearity in free and forced vibration of a sinusoidal shallow arch elastically restrained at the ends, will be presented following the analogous of [18, 19]. The effects of the parameters such as the initial curvature and concentrated force were examined. The theoretical model is based on the EulerBernoulli beam theory and the von Kármán geometrical nonlinearity assumptions. Harmonic motion is assumed and the transverse displacement is expended into a series of finites spatial functions. Using Hamilton’s principle energy, the problem is reduced into a set of nonlinear algebraic equations solved numerically using an approximate method the socalled second formulation. Based on multimode approach, the nonlinear behavior of system can be clearly observed in the backbone curves giving the nonlinear amplitudefrequency dependence of the structure. The corresponding first and second nonlinear deflections and curvatures of the arch are presented.
2. Mathematical review
2.1. The governing equation
In the preset study, an EulerBernoulli shallow arch with initial curvature ${w}_{0}\left(x\right)$, initial rise $q$ and elastically restrained at the ends by vertical and rotational springs at left and right with stiffness coefficients ${K}_{{t}_{1}}$, ${K}_{{\theta}_{1}}$, ${K}_{{t}_{2}}$ and ${K}_{{\theta}_{2}}$ respectively. The expressions of the kinetic energy $T$, the axial strain energy ${V}_{a}$ and the bending strain energy ${V}_{b}$ of the system represented in Fig. 1 may be written as:
The total strain energy of a shallow arch $V$ can be expressed as the sum of the axial strain energy ${V}_{a}$ and the bending strain energy ${V}_{b}$ as:
where $L$ is the projected length of the beam, A crosssection, $\rho $ is the density per unit length, $I$ the moment of inertia, $E$ the modulus of elasticity, ${w}_{0}\left(x\right)$ is the initial configuration function which is assumed to be:
Fig. 1Schematic of vibration of a shallow arch supported by translational and rotational springs at both ends
The kinetic energy Eq. (1) and total strain energy Eq. (4) due to the bending and axial strain are used in Hamilton’s principle as defined in Eq. (5) to obtain the equilibrium equation of motion of a shallow arch, the longitudinal displacement is eliminated from the equations:
The corresponding equation of motion can be written as follows:
For convenience, the following nondimensional variables are used:
where $r=\sqrt{I/A}$ is the radius of gyration of the cross section, ${x}^{*}$ is the dimensionless length, ${w}^{*}\left(x\right)$ is the dimensionless deflection from the initial configuration, ${{w}_{0}}^{*}\left(x\right)$ the dimensionless initial curvature of the arch., ${q}^{*}$ is the dimensionless initial rise of the arch. The dimensionless equation of motion can be written as follows:
The stars of the dimensionless parameter of equations are vanished to simplify their rewriting. The boundary conditions of system at left rand at right and are respectively.
The left end condition:
The right end condition:
2.2. Linear problem
The linear mode shapes and frequencies satisfy the following differential system:
The motion is assumed to be harmonic $w(x,t)=w\left(x\right){e}^{i\omega t}$ where $i={\left(1\right)}^{\frac{1}{2}}$, $\omega $ is the linear frequency and $w\left(x\right)$ is the mode shape function. The linear differential equation of motion becomes such bellow:
The solution of the differential Eq. (13), may be written as a sum of homogenous a particular solution as [10]:
The eigenvalues system 5×5 has been determined iteratively using NewtonRaphson algorithm.
The first five frequencies dependent to the initial rise of the arch have been plotted in Fig. 2 by assuming the vertical and rotational stiffness ${K}_{{t}_{1}}={K}_{{\theta}_{1}}={K}_{{t}_{2}}={K}_{{\theta}_{2}}=\infty $ and considering the initial rise of the $q$ varies in interval [0, 20].
In Fig. 3 and 4 the first eight symmetrical and antisymmetrical modes shapes are plotted respectively for a value of initial rise $q=$4.
Fig. 2Variation of the natural frequencies of the five modes of clampedclamped shallow arch with nondimensional rise, dash line ( ) for symmetrical modes, continuous line is for antisymmetrical modes
The Fig. 2 shows the variation of the frequencies with the initial dimensionless rise of the arch. We noticed that only the symmetrical mode influenced by the initial rise dimensionless of the arch. The crossing phenomena between two approaching (first order and second order modes) exists at $q=$ 9.168 and (thirdorder and fourthorder modes) at $q$$=$16.713. The results are compared with those given by W. Lacarbonara and G. Rega in [11] and have a very good agreement as presented in Table 1. The first four linear symmetrical and antisymmetric mode shapes of the shallow arch are plotted respectively in Fig. 3 and 4 for an initial value $q=$4.
Table 1The ratios frequencies of clampedclamped shallow arch
Ratios  Frequencies  Present $q$  Ref. [3] $q$ 
1:1  ${\omega}_{1}={\omega}_{2}$  9.168  9.15 
1:1  ${\omega}_{3}={\omega}_{4}$  16.713  16.72 
2:1  ${\omega}_{2}=2{\omega}_{1}$  2.788  2.79 
2:1  ${\omega}_{3}=2{\omega}_{2}$  3.01  3.01 
3:1  ${\omega}_{3}=3{\omega}_{1}$  4.975  4.98 
3:1  ${\omega}_{4}=3{\omega}_{1}$  10.782  10.78 
3:1  ${\omega}_{3}=3{\omega}_{2}$  14.712  14.70 
Fig. 3Symmetrical linear displacement wi for i= 1, 3, 5 and 7 for an initial rise q= 4
Fig. 4Anti symmetrical linear displacement wi for i= 2, 4, 6 and 8 for an initial rise q= 4
2.3. Formulation of the nonlinear problem
To develop the nonlinear theory, the transverse displacement function is expanded as a series of $N$ basic spatial functions:
Using a generalized parameterization and the usual summation convention defined in [4, 18], the kinetic energy, and the total strain energy due to axial and bending strain energy of the shallow arch can be expressed as:
where ${m}_{ij}$, ${k}_{ij}$, ${C}_{ijk}$ and ${b}_{ijkl}$ are respectively, the mass tensor, linear, cubic nonlinear and quadratic nonlinear rigidity of the shallow arch in which their expressions are defined as:
The coefficients ${a}_{i}$ are unknowns as well as the frequency ${\omega}_{i}$. The dynamic behavior of the beam is examined under a centered point force ${F}_{c}$ applied at point ${x}_{f}$ and a distributed harmonic force ${F}_{d}$ defined by:
$\delta $ is the Dirac function, and ${F}_{i}^{c}\left(t\right)$ and ${F}_{i}^{d}\left(t\right)$ are the corresponding concentrated force and generalized forces that expressed as:
${F}_{i}^{d}\left(t\right)={F}^{d}\mathrm{sin}\left(\omega t\right){\int}_{0}^{1}{w}_{i}\left(x\right)dx={f}_{i}^{d}\mathrm{sin}\left(\omega t\right).$
Using Hamilton’s principle defined in Eq. (5) and the matrix notation by taking into account the forcing term, the process leads to nonlinear algebraic equations for the unknown coefficients ${a}_{i}$ that can be written in a matrix form as follows:
where $\left\{A\right\}$ is a column vector of the basic functions contribution coefficients ${a}_{i}$. To evaluate the nondimensional parameter, ones put:
And the dimensionless generalized forces ${f}_{i}^{c*}$ and ${f}_{i}^{d*}$ are given by:
Which leads to:
where ${m}_{ij}^{*}$, ${k}_{ij}^{*}$, ${C}_{ijk}^{*}$ and ${b}_{ijkl}^{*}$ are the nondimensional generalized parameters given by:
For a uniform rectangular crosssection, $\alpha =$ 3, since in this case, $\left(\frac{{h}^{2}A}{I}=12\right)$. The functions $\left\{{w}_{i}^{*},i=1,\dots ,n\right\}$ have been normalized in order to obtain the mass matrix equals the identity matrix. The corresponding integrals, defined in Eqs. (29)(32), have been calculated numerically, using the Simpson’s rule. Eq. (25) can be rewritten in nondimensional form as:
Using the tensorial notation, ones put:
Eq. (33) represents a set of nonlinear algebraic equations in which $\left\{{F}^{*}\right\}$ is a column vector of generalized forces. This equation appears as a generalization to the nonlinear case of the classical linear forced response equation, very well known in linear modal analysis theory, to which the term $\frac{3}{2}\left[{B}^{*}\left(\left\{A\right\}\right)\right]\left\{A\right\}$corresponding to the nonlinear geometrical rigidity is added.
The parameters ${m}_{ij}^{*}$, ${k}_{ij}^{*}$ and ${b}_{ijkl}^{*}$ have been solved numerically using Simpson’s method in range [0, 1] and the nonlinear algebraic Eq. (34) has been solved using the formulation developed in [10] that called the second formulation, it is an approximation which consists to separate the nonlinear term of Eq. (34) into a terms proportional to ${a}_{1}^{3}$ and terms proportional to ${a}_{1}^{2}{\epsilon}_{i}$ and neglecting terms proportional to ${a}_{1}{\epsilon}_{i}{\epsilon}_{j}$ and terms proportional to ${\epsilon}_{i}{\epsilon}_{j}{\epsilon}_{k}$ one can write:
After substituting and rearranging, Eq. (35) can be written in matrix form as:
In which ${\left[{{K}_{r}}^{*}\right]}_{R}$ = ${\left[{{k}_{ij}}^{*}\right]}_{R}$ and ${\left[{{M}_{r}}^{*}\right]}_{R}$ = ${\left[{{m}_{ij}}^{*}\right]}_{R}$ are reduced rigidity and mass matrices associated with the first nonlinear mode, obtained by varying $i$ and $j$, ${\left[{{\alpha}_{r}}^{*}\right]}_{R}$ is a square matrix, depending on ${a}_{1}$, whose general term ${{\alpha}_{ij}}^{*}$ is equal to ${{a}_{1}}^{2}{b}_{11ir}$, and ${a}_{1}^{3}{b}_{111r}$ is a column vector representing the right side of the linear system Eq. (32) in which the reduced unknown vector is ${\left\{{A}_{r}\right\}}_{R}={\left\{{\epsilon}_{2},{\epsilon}_{3},...,{\epsilon}_{r}\right\}}^{T}$.
The modal contributions ${\epsilon}_{i}$ can be obtained easily by solving the nonlinear system Eq. (35) of $n$equations and $n$ unknowns. Higher nonlinear mode shapes may be obtained in a similar manner, using appropriate reduced matrix in each case.
Where the frequency ratio ${\omega}^{*2}$ is given from Eq. (34):
where ${\omega}_{l}^{*}$ is the linear natural frequency parameter that can be written as:
If free vibration case is considered the dependence of the nonlinear frequency ratio $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$ on the dimensionaless vibration amplitude ${w}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$ is obtained from Eq. (37) by putting ${f}_{i}^{*c}=$ 0. The normalized first nonlinear mode shape for a given assigned value ${a}_{1}$ is then obtained as a series involving the arch modal parameters depending on the first six symmetric modes given by [18]:
In which the predominant term, proportional to the first linear mode shape, is ${a}_{1}{w}_{1}^{*}\left(x\right)$ and other terms, proportional to the higher modes ${w}_{3}^{*}\left(x\right),\dots ,{w}_{11}^{*}\left(x\right)$ are the corrections due to the nonlinearity. For the second amplitudedependent nonlinear mode shape ${w}_{nl2}^{*}\left(x,{a}_{2}\right)$ for a given assigned value ${a}_{2}$ of the second function contribution, the predominant term is ${a}_{2}{w}_{2}^{*}\left(x\right)$ and other terms, proportional to the higher modes ${w}_{4}^{*}\left(x\right),\dots ,{w}_{12}^{*}\left(x\right)$ are the corrections due to the nonlinearity.
3. Results and discussion
In this section, numerical results obtained for various values of the initial rise and levels of the excitation force will be presented via solution of the set of nonlinear algebraic equations Eq. (36), based on the second formulation for the value of the stiffness of the translational and rotational springs at both ends ${K}_{{t}_{1}}={K}_{{t}_{2}}={K}_{{\theta}_{1}}={K}_{{\theta}_{2}}=\infty $.
The basic functions contributions ${\epsilon}_{i}$ for $i$th mode shapes, the corresponding maximum nondimensional vibration amplitudes, and the curvatures at $x=$ 0 are given in Table 2 and 3 for each given value of ${a}_{i}$. It can be seen in this tables, that the contribution ${a}_{i}$ remains predominant for the whole range of vibration amplitudes considered, compared to the contributions ${\epsilon}_{i}$ of the other functions, it can be noticed also that the contribution coefficients ${\epsilon}_{i}$ increased according to the maximum nondimensional amplitude ${w}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$.
The numerical solutions obtained, have been substituted in Eq. (36) in case of $\left\{{F}_{r}^{*}\right\}=\left\{0\right\}$ and calculation has been made of the resulting vector {Residual}.
Fig. 5Backbone curves of clampedclamped shallow arch, in the vicinity of the first mode, for various values of initial rises (q= 2, 4 and 8)
Fig. 6The first normalized nonlinear mode shape for q= 4
Fig. 7The second normalized nonlinear mode shape for q= 4
In Fig. 5, the backbone curves giving the nonlinear amplitudefrequency dependence, based on the multimode approach of a C–C shallow arch are plotted in the neighborhood of the first nonlinear mode shape for various initial rise $q=$2, 4 and 8. It may be noticed that increasing the initial rise the nonlinearity decreases, these results are due to the fact that when the initial rise is increased, the flexural rigidity of the arch increases. We can also notice from the Fig. 2 that the linear frequencies of symmetrical modes increase by the increasing of the initial rise, consequently the frequency ratio decreases since the linear frequency presents the denominator of the frequency ratio. However, in the second vibration mode, the nondimensional frequency ratios exhibit the same behavior as the straight beam. These results are compatible with the results of linear regime.
The normalized nondimensional vibration amplitude corresponding to the first and second mode shape are presented respectively in Fig. 6 and 7.
In Fig. 8 and 9, the curvature curves of clampedclamped shallow arch, associated to the first and second nonlinear deflection are plotted for various levels of vibration amplitudes by taking into account an initial rise $q=$ 4.
Table 2Frequency ratios of free vibration in the first nonlinear mode, curvatures and modal participation of a CC shallow arch at various amplitudes obtained by the second formulation
$\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${w}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{d}^{2}w}{d{x}^{2}}\left(0\right)$  ${a}_{1}$  ${\epsilon}_{3}$  ${\epsilon}_{5}$  ${\epsilon}_{7}$  ${\epsilon}_{9}$  ${\epsilon}_{11}$  Residual 
1.000198  0.074687  2.608599  0.05  1,11E06  5,22E07  1,17E07  3,60E08  1,72E08  8.43E35 
1.019601  0.745852  26.66276  0.5  0,00103  0,000506  0,000117  3,63E05  1,73E05  3.79E29 
1.076265  1.48718  56.22632  1  0,006837  0,003686  0,00092  0,000295  0,000139  7.47E27 
1.164595  2.223823  89.63892  1.5  0,018144  0,010791  0,002977  0,001  0,000467  1.49E25 
Table 3Frequency ratios of free vibration in the second nonlinear mode, curvatures and modal participation of a CC shallow arch at various amplitudes obtained by the second formulation
$\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${w}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{d}^{2}w}{d{x}^{2}}\left(0\right)$  ${a}_{2}$  ${\epsilon}_{2}$  ${\epsilon}_{4}$  ${\epsilon}_{6}$  ${\epsilon}_{8}$  ${\epsilon}_{12}$  Residual 
1.002089  0.0754  6.176212  0.05  2.24E05  4.59E06  1.23E06  5.46E07  3.65E08  4.62E32 
1.190862  0.7470  68.37282  0.5  0.016765  0.004105  0.00122  0.000563  4.26E05  4.85E26 
1.634812  1.4742  154.2536  1  0.077353  0.02398  0.008552  0.004223  0.00041  2.10E24 
2.182516  2.1932  247.0735  1.5  0.155322  0.055208  0.022587  0.011836  0.0015  4.24E23 
Table 4Comparison of the first nonlinear frequency ratio ωnl*ωl* and maximal nonlinear amplitude with those obtained in reference [4]
Ref. [4].  Present with $q=$ 2  Present with $q=$ 4  Present with $q=$ 8  
${a}_{1}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$ 
0.05  1.001  0.079  1.00059  0.07824  1.0001  0.0746  1.00006  0.0636 
0.5  1.106  0.789  1.0581  0.7792  1.0196  0.7462  1.0059  0.6369 
1  1.361  1.554  1.2159  1.5432  1.0762  1.489  1.0237  1.2760 
1.5  1.684  2.297  1.441  2.2905  1.1645  2.228  1.0527  1.9175 
Table 5Comparison of the second nonlinear frequency ratio ωnl*ωl* and maximal nonlinear amplitude with those obtained in reference [4]
Ref. [4].  Present with $q=$ 2  Present with $q=$ 4  Present with $q=$ 8  
${a}_{2}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$  $\frac{{\omega}_{nl}^{*}}{{\omega}_{l}^{*}}$  ${W}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}$ 
0.05  1.002  0.075  1.002  0.0754  1.002  0.0754  1.002  0.0754 
0.5  1.187  0.750  1.190  0.7470  1.190  0.7470  1.190  0.7493 
1  1.603  1.488  1.634  1.4741  1.634  1.4742  1.634  1.4809 
1.5  2.101  2.217  2.182  2.1946  2.182  2.1932  2.182  2.1992 
Fig. 8The linear and nonlinear curvature corresponding to the first mode shape for q= 4
Fig. 9The linear and nonlinear curvatures corresponding to the second mode shape for q= 4
Fig. 10Curvature distributions, associated to the first nonlinear deflection, of a C–C shallow arch, for various values of initial rise (q = 2, 4 and 8) and a maximum nondimensional amplitudes wmax*=1.5
The Table 4 and 5, present the values of the nondimensional nonlinear frequency ratios and the nondimensional vibration amplitude of clampedclamped shallow arch which are calculated using the second formulation for various values of the initial rise 𝑞. A comparison is made with a clampedclamped beam straight beam [4].
The Table 5 shows that the frequency ratio remains constant along the variation of initial rise of the arch. In this case the same behavior of straight beam is may be observed. It was found, that for amplitude equal to 1 time the shallow arch thickness, the difference induced does not exceed 1.8 % for the nondimensional frequency ratio and 0.48 % for the nondimensional vibration in comparison with the results given by [4].
In Fig. 12 and 13, the nonlinear frequency response functions of the shallow arch examined, are plotted in the neighborhood of the first nonlinear mode shape, for various values of the concentrated force and initial rise ($q=$ 2, 4 and 8).
In Figs. 10, the curvature distributions are plotted, associated to the first linear and nonlinear deflection, of a C–C shallow arch in free vibration case, for various initial rise (2, 4, 8) and a maximum nondimensional amplitudes $\left({w}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{*}=1.5\right)$. The corresponding values are summarized in Table 6, in which the percentage correction at the clamps and at the middles are given for various values of initial rise. An increase of curvature with the initial rise, especially at the clamped end of the shallow arch may be clearly observed.
Table 6 presents the calculation of the percentage correction between the curvatures via the nonlinear theory compared to the linear theory. it can be seen that the percentage correction at the middles of the arch reaches 17.74 % for $q$ = 2, 2.91 % for $q$ = 4, and 3.29 % for $q$ = 8. In addition, it can be noticed that the maximum value of the linear and nonlinear curvature of the arch in case of $q$ = 2, 4 and 8 is located near the clamped ends.
Table 6Effect of rise on the percentage correction between the curvatures estimated via the linear and nonlinear theories for various values of initial rise. Maximum nondimensional amplitudes wmax*=1.5
${q}^{*}$  $\frac{{d}^{2}{w}^{*}}{d{x}^{*2}}\left(0\right)$  Percentage correction (%)  $\frac{{d}^{2}{w}^{*}}{d{x}^{*2}}\left(\frac{1}{2}\right)$  Percentage correction (%)  
Linear  Nonlinear  Linear  Nonlinear  
2  46.11  58.06  20.58  –22.36  –18.99  17.74 
4  52.15  57.30  8.98  –13.00  –13.39  2.91 
8  73.04  75.59  3.37  5.95  –2.59  3.29 
Table. 7 presents the contributions of the six symmetric modes. The calculation of the generalized forces led to the conclusion that the first mode is predominant.
The hardening behavior type of geometrical nonlinearity can be clearly observed in Fig. 1113. The discontinuities shown in Fig. 1115 are due to the fact of the existence of three amplitudes for one range of frequency, this phenomenon is characteristic on nonlinear systems, and is called a jump phenomenon. Another property of nonlinear system can be clearly observed; the frequencyresponse functions are not proportional to the level of the concentrated forces applied.
By assuming a multimode approach used previously in [18], Eq. (37) is plotted in the neighborhood of the first nonlinear mode shape in case of various values of concentrated force ${F}^{*}$ in Fig. 12 and various values of the initial rise $q$ in Fig. 13. Figure 11 shows a comparison between the linear and nonlinear frequency response functions, where the shallow arch is excited by a concentrated force ${F}^{*}=100$.
Table 7Percentage of generalized forces exciting the first six symmetric modes of a CC shallow arch excited by uniformly distributed force over an interval [0, 1].
Modes  1  3  5  7  9  11 
${\int}_{0}^{1}{w}_{i}^{*}\left({x}^{*}\right)$  0.6522  0.3572  0.2313  0.1697  0.1340  0.1095 
$\frac{{\int}_{0}^{1}{w}_{i}^{*}\left({x}^{*}\right)}{{\sum}_{i=1}^{12}\left{\int}_{0}^{1}{w}_{i}^{*}\left({x}^{*}\right)\right}$  0.3943  0.2159  0.1398  0.1026  0.0810  0.0662 
In Fig. 12 and 13 the nonlinear frequency response functions of the shallow arch are plotted in the neighborhood of the first nonlinear mode shape for various levels of the concentrated force and initial rise $q$. The hardening type of geometrical nonlinearity of the clampedclamped shallow arch excited by a concentrated force is decreased in a remarkable way with increasing the initial rise. This is due to the fact that when the initial rise is increased, the flexural rigidity of the arch decreases.
The response curves have been illustrated for two cases: the first, the CC shallow arch was excited by a concentric harmonic excitation in Fig. 12 and 13, in the second, by a uniformly distributed force in Fig. 14 and 15. It may be observed from these figures, the hardening type of geometrical nonlinearity also the jump phenomenon where for the same frequency range, several results are possible, as well as the nonproportional evolution of the frequency response to the intensity of excitation.
Fig. 11Linear and nonlinear frequency response functions comparison for F*= 100 and q= 2
Fig. 12Nonlinear frequency response functions, based on the multimode approach, of a CC shallow arch and various levels of the concentrated forces
Fig. 13Nonlinear frequency response functions, based on the multimode approach, of a CC shallow arch for a concentrated forces F*= 100 and various levels of q
Fig. 14Nonlinear frequency response functions. based on the multimode approach. of a CC shallow arch and various levels of excitation: in case uniformly distributed harmonic force.
Fig. 15Nonlinear frequency response functions. based on the multimode approach. of a CC shallow arch for a harmonic force uniformly distributed Fd=100 and various levels of q.
4. Conclusions
A theoretical model to analyze the large vibration amplitude of shallow arch was presented to examine the geometrical nonlinearity in transverse vibration of a shallow arch. The linear frequencies were obtained iteratively using the NewtonRaphson algorithm. The mathematical formulation of the problem is based on the EulerBernoulli theory and the Von Karman assumptions. The kinetic and total strain energy were derived and discretized into a series of a finites spatial functions by applying Hamilton’s principle energy and spectral analysis. The problem was reduced into a set of nonlinear algebraic equations solved numerically using an approximate method developed previously the socalled second formulation. The effect of the initial rise and the concentrated force on nonlinear behavior of the arch was analyzed. It is found that only the symmetrical modes shapes are influenced by the initial rise, these results are due to the fact that when the initial rise is increased, the flexural rigidity of the arch increases. The passage from straight beam to curved beam is very sensitive and depends on the functionality of the future structure, but it is known that curved beam or shallow arch are very performed in case of structure which will support a higher load, because it cans resist more and avoid the earlier deflection of structure.
References

Fakhreddine H., Adri A., Chajdi M., Rifai S., Benamar R. A multimode approach to geometrically nonlinear forced vibration of beams carrying point masses. Diagnostyka, Vol. 21, Issue 4, 2020, p. 2333.

Chajdi M., Ahmed A., El Bikri K., Benamar R. Analysis of the associated stress distributions to the nonlinear forced vibrations of functionally graded multicracked beams. Diagnostyka, Vol. 22, Issue 1, 2021, p. 101112.

Ahmed A., Rhali B. Linear and geometrically nonlinear frequencies and mode shapes of beams carrying a point mass at various locations. an analytical approach and a parametric study. Diagnostyka, Vol. 18, Issue 2, 2017, p. 1321.

Benamar R., Bennouna M. M. K., White R. G. The effects of large vibration amplitudes on the mode shapes and natural frequencies of thin elastic structures part I: simply supported and clampedclamped beams. Journal of Sound and Vibration, Vol. 149, Issue 2, 1991, p. 179195.

Moussaoui F., Benamar R., White R. G. The effects of large vibration amplitudes on the mode shapes and natural frequencies of thin elastic shells, part I: coupled transversecircumferential mode shapes of isotropic circular cylindrical shells of infinite length. Journal of Sound and Vibration, Vol. 232, Issue 5, 2000, p. 917943.

Boutahar L., El Bikri K., Benamar R. A Homogenization procedure for geometrically nonlinear free vibration analysis of functionally graded annular plates with porosities, resting on elastic foundations. Ain Shams Engineering Journal, Vol. 7, Issue 1, 2016, p. 313333.

Bekkaye Merrimi E., El Bikri K., Benamar R. Geometrically nonlinear steady state periodic forced response of a clampedclamped beam with an edge open crack. Comptes Rendus Mécanique, Vol. 339, Issue 11, 2011, p. 727742.

Outassafte O., Adri A., Yassine El Khouddar, Rifai S., Benamar R. Geometrically nonlinear free inplane vibration of circular arch elastically restrained against rotation at the two ends. International Journal of Engineering Trends and Technology, Vol. 261, 2021, p. 113515.

Raymond Plaut H., Niels Olhoff Optimal forms of shallow arches with respect to vibration and stability. Journal of Structural Mechanics, Vol. 11, Issue 1, 1982, p. 8188.

Nayfeh A. H., Kreider W. Investigation of Natural Frequencies and Mode Shapes of Buckled Beams. AIAA Journal, Vol. 33, Issue 6, 1995, p. 11211126.

Lacarbonara W., Rega G. Resonant nonlinear normal modes. part ii: activation/orthogonality conditions for shallow structural systems. International Journal of NonLinear Mechanics, Vol. 38, Issue 6, 2003, p. 873887.

Yi Z., Stanciulescu I. Nonlinear normal modes of a shallow arch with elastic constraints for twotoone internal resonances. Nonlinear Dynamics, Vol. 83, 2016, p. 15771600.

Tsiatas G. C., Babouskos N. G. Linear and geometrically nonlinear analysis of nonuniform shallow arches under a central concentrated force. International Journal of NonLinear Mechanics, Vol. 92, 2017, p. 92101.

Ye S.Q., Mao X.Y., Ding H., Ji J.C., Chen L.Q. Nonlinear vibrations of a slightly curved beam with nonlinear boundary conditions. International Journal of Mechanical Sciences, Vol. 168, 2020, p. 105294.

Blair K. B., Krousgrill C. M., Farris T. N. Nonlinear dynamic response of shallow arches to harmonic forcing. Journal of Sound and Vibration, Vol. 194, Issue 3, 1996, p. 353367.

Lakrad F., Schiehlen W. Effects of a low frequency parametric excitation. Chaos, Solitons and Fractals, Vol. 22, Issue 5, 2004, p. 11491164.

Ozkaya E., Sarigul M., Boyaci H. Nonlinear transverse vibrations of a slightly curved beam resting on multiple springs. International Journal of Acoustics and Vibration, Vol. 21, Issue 4, 2016, p. 379391.

El Kadiri M., Benamar R., White R. G. Improvement of the semianalytical method, for determining the geometrically nonlinear response of thin straight structures. Part I: application to clampedclamped and simply supportedclamped beams. Journal of Sound and Vibration, Vol. 249, Issue 2, 2002, p. 263305.

Azrar L., Benamar R., White R. G. A semianalytical approach to the nonlinear dynamic response problem of beams at large vibration amplitudes, part ii: multimode approach to the steady state forced periodic response. Journal of Sound and Vibration, Vol. 255, Issue 1, 2002, p. 141.