Next Article in Journal
On Third Hankel Determinant for Certain Subclass of Bi-Univalent Functions
Previous Article in Journal
Analysis of Mixing Efficiency in a Stirred Reactor Using Computational Fluid Dynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamic Analysis of a Curved Sandwich Beam with a Time-Dependent Viscoelastic Core Using the Generalized Differential Quadrature Method (GDQM)

1
Department of Naval Architecture and Marine Engineering, Yildiz Technical University, TR-34349 Istanbul, Turkey
2
Department of Aeronautical Engineering, Istanbul Technical University, TR-34469 Istanbul, Turkey
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 238; https://doi.org/10.3390/sym16020238
Submission received: 5 December 2023 / Revised: 3 January 2024 / Accepted: 4 January 2024 / Published: 15 February 2024

Abstract

:
This paper focuses on the geometrically nonlinear dynamic analyses of a three-layered curved sandwich beam with isotropic face layers and a time-dependent viscoelastic core. The boundary conditions and equations of motion governing the forced vibration are derived by using Hamilton’s principle. The first-order shear deformation theory is used to obtain kinematic relations. The spatial discretization of the equations is achieved using the generalized differential quadrature method (GDQM), and the Newmark-Beta algorithm is used to solve the time variation of the equations. The Newton–Raphson method is used to transform nonlinear equations into linear equations. The validation of the proposed model and the GDQM solution’s reliability are provided via comparison with the results that already exist in the literature and finite element method (FEM) analyses using ANSYS. Then, a series of parametric studies are carried out for a curved sandwich beam with aluminum face layers and a time-dependent viscoelastic core. The resonance and cancellation phenomena for the nonlinear moving-load problem of curved sandwich beams with a time-dependent viscoelastic core are performed using the GDQM for the first time, to the best of the authors’ knowledge.

1. Introduction

The beams, which are made of conventional isotropic materials, have been studied by many researchers for decades. However, in recent years, with the rapid development of manufacturing and material technologies, studies about new kinds of materials, namely composites and FGM, have been attractive so far. Nonetheless, sandwich structures have been widely preferred in the marine, aerospace, civil, and railroad vehicle industries due to their high strength, low weight, and good vibration damping properties. There are many studies devoted to the static, dynamic, and vibration analyses of structures with different kinds of geometries, namely flat, curved, and variable-thickness beams. Also, both linear and nonlinear analyses have been performed so far. Apart from the exact analytical solutions, many numerical methods for solving partial differential equations are found by scientists.
Beams are used in a wide variety of engineering applications, as they are basic structural elements. They are used in rotor blades, shafts, frames, and robotic structures. Flat beam problems are widely found in literature. Singh et al. [1] presented various formulations about the large-amplitude free vibration of beam structures. Gafsi et al. [2] studied the vibration of nonlinear flexible beams with variable geometric properties. After the GDQM was developed by Richard Bellman and his associates in the early 1970s, many studies have been conducted so far. Bert and Malik [3] investigated freely vibrating cantilever beams by using the GDQM. They presented the effect of different sampling points on the natural frequency values by comparing them with the analytical results. Jang et al. [4] studied the application of the GDQM to the static analyses of structural components, such as beams and rectangular and circular plates. Feng and Bert [5] applied the quadrature method for solving flexural vibration problems, including geometrical nonlinearities. They used only one degree of freedom. Karami and Malekzadeh [6] presented a new differential quadrature methodology for stability, static, and dynamic analyses of beams and also calculated the error regarding the number of sampling points. Zhong and Guo [7] performed nonlinear vibration analyses of thick beams with Timoshenko beam theory using the GDQM and mostly computed the ratio of the linear frequency to the nonlinear frequency. Ghasemi and Mohandes [8] investigated the finite strain effect for nonlinear free vibration analyses of unidirectional composite beams with Timoshenko beam theory by using the GDQM. Mahmoud et al. [9] used GDQM for free vibration analysis of delaminated beam plates. Eftekhari conducted several studies of beams and plates on both linear [10,11,12] and nonlinear [13] dynamic problems, specifically for moving loads using the GDQM. He investigated many effects, such as sampling points and Dirac-delta function coefficients. He found perfect agreement with the analytical solutions.
Also, several studies are found in the literature for curved beams. Yang and Yau [14] investigated the dynamic response of horizontally deep-curved beams under both horizontally and vertically moving loads. Hajianmaleki and Qatu [15] conducted static and vibration analyses for thick and generally laminated deep-curved beams with various boundary conditions. Poojary and Roy [16] used the finite element method in order to study the in-plane vibration of curved beams under moving loads. Moreover, the GDQM is also used to solve problems with curved beams. Kurtaran [17,18] used the GDQM for solving several large displacement transient and static analyses of curved beams. Nikkhoo et al. [19] applied the GDQM in order to study the dynamics of curved beams under moving loads. Kang et al. [20] performed vibration analyses for horizontally curved beams with warping using the GDQM. It is also possible to use the GDQM for solving more complex structural problems than beams, such as plates and shells. Akgün and Kurtaran [21] performed transient analyses of laminated composite super-elliptic shells with geometrical nonlinearities using the GDQM. Malik and Bert [22] showed the implementation of multiple boundary conditions for the GDQM solutions of higher-order partial derivative equations, which govern the free vibration of the plates. Fung [23] solved dynamic problems using the GDQM and examined the stability and accuracy of the method. Most of the sampling point types are used, such as Legendre, Radau, Chebsevey, Chebyshev–Gauss–Lobatto, and uniformly. Bacciocchi et al. [24] also conducted a vibration analysis of variable-thickness shells and plates using the GDQM.
Viscoelastic materials are widely preferred in sandwich structures due to their vibration and energy damping properties. Azevedo and Coben [25] performed a dynamic analysis of a viscoelastic Timoshenko beam numerically using the Kelvin–Voigt material model. Chen and Lin [26] conducted a dynamic analysis of viscoelastic material by incremental FEM. Shafei [27] investigated the nonlinear transient vibration of viscoelastic plates with a NURBS-based iso-geometric higher-order shear deformation theory approach. Tolpekina et al. [28] studied both linear and nonlinear viscoelastic moduli of rubber materials by creep compliance and stress relaxation analyses. Arikoglu [29] made a new fractional derivative linearly viscoelastic material model and parameter identification by genetic algorithms. Demir et al. [30] performed a vibration analysis of a curved composite sandwich beam with a viscoelastic material core layer using the GDQM. Arikoglu and Ozkol [31] and Taskin et al. [32] performed damping and vibration analyses of composite sandwich plates and sandwich cylindrical shells, respectively, in which the core layer is a frequency-dependent viscoelastic material made using the GDQM.
This study aims to predict the nonlinear dynamic characteristics of a three-layered curved sandwich beam with isotropic face layers and a time-dependent viscoelastic core. The boundary conditions and equations of motion governing the forced vibration are derived by using Hamilton’s principle. The first-order shear deformation theory is used to obtain kinematic relations. The spatial discretization of the equations is achieved by using the generalized differential quadrature method (GDQM), and the Newmark-Beta algorithm is used to solve the time variation of the equations. Validation of the proposed model and the GDQ solution’s reliability are provided by comparison with the results in the literature and finite element method (FEM) analyses performed via ANSYS. Then, a series of parametric studies are carried out for a curved sandwich beam with aluminum face layers and a time-dependent viscoelastic core. The dynamic analyses of curved sandwich beams with a time-dependent viscoelastic core are performed using GDQM. Cancellation and resonance phenomena have been studied widely for linear moving load problems; however, it is the first time these phenomena have been considered for a nonlinear problem to the best of the authors’ knowledge.

2. Theory and Formulations

2.1. Kinematics, Constitutive, and Equilibrium Equations

The governing equations are derived from Hamilton’s principle. To obtain the equation of motion of a curved beam, some assumptions are made as follows:
  • The connection between layers is simply perfect, so delamination does not occur.
  • First-order shear deformation theory is used.
  • Face layers are assumed to be incompressible, unlike the core layer.
  • Since the deformations are not small, Von Karman-type nonlinear strain terms are included.
Displacement relations of the layers are obtained for the sandwich beams, which require continuity and are represented in Figure 1 between the face layers for i = 1,3 , where subscript represents the layer number.
u i θ , r , t = u 0 i θ , t + r i φ i θ , t w i θ , r , t = w 0 i θ , t
For the viscoelastic core layer:
u 2 θ , r , t = u 02 θ , t + r 2 φ 2 θ , t w 2 θ , r , t = w 02 θ , t + r 2 φ r θ , t
Regarding the beam theory with Von Karman-type nonlinearity, the above displacement relations are obtained. Here u i and w i represent the circumferential and radial displacements of the beam, respectively. u 0 i , w 0 i , and φ i represent the circumferential displacement, radial displacement, and circumferential rotation of mid-surface of the beam, respectively. φ r and θ are the radial rotation of the core layer and the circumferential coordinate. The displacement relations of the core layer can be expressed in terms of the face layers to eliminate 1/3 of the degrees of freedom as follows:
u 02 θ , r , t = 1 2 ( u 01 θ , r , t + u 03 θ , r , t ) + 1 4 ( h 3 φ 3 θ , r , t h 1 φ 1 θ , r , t ) φ 2 θ , r , t = u 01 θ , r , t u 03 θ , r , t h 2 h 3 φ 3 θ , r , t + h 1 φ 1 θ , r , t 2 h 2 w 02 θ , r , t = 1 2 ( w 01 θ , r , t + w 03 θ , r , t ) φ r θ , r , t = 1 h 2 ( w 01 θ , r , t w 03 θ , r , t )
Kinematic relations of the ith layer of the beam structure for cylindrical coordinates can be expressed by the Von Karman nonlinear theory for large displacements for i = 1,3 , where superscript represents the layer number.
ε θ i = 1 R i u i θ + w i + 1 2 1 R i w i θ 2 γ θ r i = u i r + 1 R i w i θ u 0 i ε r 2 = w 2 r + 1 2 w 2 r 2
By using the generalized Hooke Law, the stress-strain relation of the isotropic materials can be written for beams as follows for i = 1,3 :
σ θ i τ θ r i = E 1 0 0 1 ν 2 ε θ i γ θ r i σ θ 2 σ r 2 τ θ r 2 = E 1 ν 2 1 ν 0 ν 1 0 0 0 1 ν 2 ε θ 2 ε r 2 γ θ r 2
Although shear strain is the most dominant for viscoelastic materials, strains in radial and circumferential directions are not neglected. The Lagrangian equation is energy conservation, which is the sum of kinetic energy K, potential energy P, internal strain energy U, and external force work V in this case. By using Hamilton’s principle, the equation of motion can be obtained, as shown below:
t 1 t 2 δ K δ U + δ V d t = 0
By inserting the constitutive relations (5) into Equation (6), the first variation of the energy terms can be obtained for i = 1,2 , 3 :
δ U = b 0 θ 0 h i 2 h i 2 ( σ θ i δ ε θ i + τ θ r i δ γ θ r i ) d r R i d θ + b 0 θ 0 h 2 2 h 2 2 ( σ r 2 δ ε r 2 ) d r R 2 d θ   δ K = b 0 θ 0 h i 2 h i 2 ρ i ( u ˙ i   δ u ˙ i + w ˙ i   δ w ˙ i ) d r R i d θ δ V = 0 θ 0 q R 1 d θ
Here b , ρ , q , R , and o v e r d o t represent width of the beam, density, distributed load, radius of curvature, and time derivative for the ith layer. Thus, governing equations and natural boundary conditions can be derived by inserting the kinematic relations (4) into Equation (7). From constitutive relations together with energy equations, the Euler equations are obtained for the dynamic analysis of sandwich beams in terms of cross-section forces and moments for u 01 , u 03 , w 01 , w 03 , φ 1 and φ 3 , respectively, as presented in Appendix A.

2.2. Generalized Differential Quadrature Method

The Differential Quadrature Method is a numerical method for solving partially differential equations that was founded in 1971 by Bellman and his associates. It basically converges the partial derivatives of a discrete point with respect to spatial derivatives to the functional values of the same discrete point, which is in the domain. The main motive to develop this method is that the necessity of the grid points to solve the problem is far less compared to other numerical methods, such as the central difference method. When GDQM is first developed, the results are partially bad for solving greater systems, but then Shu and his associates developed another method, the generalized differential quadrature method, to overcome this problem.
The derivation of a function at a discrete point with respect to a spatial variable is predicted by the weighted linear summation of the functional values of all discrete points in the domain as follows:
r θ r w θ j = m = 1 M d j m r w θ m
Here M represents grid numbers in the θ direction and d i k represents weighted coefficients. The first-degree partial derivative is obtained as weighted coefficients as follows:
d j m 1 = M ¯ θ 1 θ l θ k M ¯ θ 1 θ k j m ; j , m = 1 , 2 , , M s = 1 , s j M d j s 1 j = m ; j , m = 1 , 2 , , M
The derivation of M ¯ with respect to θ is described as follows:
M ¯ θ 1 θ j = s = 1 , s j M θ j θ s
To obtain weighted coefficients for higher-degree partial derivatives, the below formulation is used as follows:
d j m r = r d j j r 1 d j m 1 d j m r 1 θ f θ m j m   ; j , m = 1 , 2 , , M s = 1 , s i M d j s r j = m   ; j , m = 1 , 2 , , M
In the method, grid numbers and distribution have a significant effect on the solution. When the grids are determined to be equally spaced, Runge’s phenomena occur. To overcome this issue, the Chebyshev–Gauss–Lobatto grid distribution is applied as follows:
θ j = b 2 1 cos j 1 M 1 π           j = 1 , 2 , , M

2.3. Time Integration, Linearization of the Equation of Motion, and Viscoelasticity

2.3.1. Newmark-Beta Method and Newton–Raphson Method with GDQM

The Newmark-Beta Method is an integration scheme that is used for solving time-dependent differential equations implicitly. It is based on the two equations below.
u ˙ i + 1 = u ˙ i + 1 γ Δ t u ¨ i + γ Δ t u ¨ i + 1
u i + 1 = u i + Δ t u ˙ i + 1 2 β ( Δ t ) 2 u ¨ i + β ( Δ t ) 2 u ¨ i + 1
Here, β and γ describe the variation of the acceleration with respect to the time step and determine the accuracy and stability of the method. While β = 0.25 and γ = 0.5 , it is named the average acceleration method, stability is not affected by the size of the time step. In the solution of the initial value problem, the equation of motion for the i + 1 ’st time step can be written as follows:
m u ¨ i + 1 + c u ˙ i + 1 + k u i + 1 = F i + 1
Rearranging the equation above, acceleration and velocity at the i + 1 ’st time step can be described as follows:
u ¨ i + 1 = 1 β ( Δ t ) 2 u i + 1 u i 1 β Δ t u ˙ i 1 2 β 1 u ¨ i
u ˙ i + 1 = γ β Δ t u i + 1 u i + 1 γ β u ˙ i + Δ t 1 γ 2 β u ¨ i
Thus,
k ^ u i + 1 = F ^ i + 1
where
k ^ = k + γ β Δ t c + 1 β ( Δ t ) 2 m
F ^ i + 1 = F i + 1 + 1 β Δ t 2 m + γ β Δ t c u i + 1 β Δ t m + γ β 1 c + 1 2 β 1 m + Δ t γ 2 β 1 c u ¨ i
Here, k ^ and F ^ represent the effective stiffness matrix and effective load vector, respectively.
Unlike linear systems, the stiffness matrix k depends on the displacement value, which can be described as u i + 1 . While u i + 1 is calculated iteratively, the k matrix needs to be linearized. The Newton–Raphson method is widely used for these systems. For more detailed algebra, chapter 5 of the book “Dynamic of Structures” by Chopra [33] can be seen.
For this problem, the curved beam domain is discrete to N points. While sacrificing the boundary grids for the boundary condition equations, the inner grids for the governing equations would be N 2 for one degree of freedom. Thus, three degrees of freedom are used for this problem, namely. Thus, the beam domain would have 3   N points in total, and the solution domain would have 3 N 2 points. The governing equation for the dynamic analysis can be written in matrix form as follows:
K D ¯ + F D ¯ D ¯ + M D ¨ ¯ = P
where K is the linear stiffness matrix and F D ¯ is the nonlinear stiffness matrix with the dimensions of 3 N 2 × 3 N 2 . D ¯ T and P T are displacement and load vectors with the dimensions 3 N 2 × 1 .
D ¯ T = u 2 , u 3 , , u n 1 , φ 2 , φ 3 , , φ n 1 , w 2 , w 3 , , w n 1 P T = 0 , 0 , , 0 , 0 , 0 , , 0 , q 2 , q 3 , , q n 1
The nonlinear Equation (33) is solved by the Newton–Raphson method. The nonlinear F D ¯ matrix can be derived from the displacement vector and the tangential stiffness matrix is obtained as follows:
G D ¯ D ¯ = F D ¯ D ¯ D ¯
Applying the Newton–Raphson method to the governing equation results in the following:
K Δ D ¯ k + F D ¯ k 1 Δ D ¯ k + M D ¨ ¯ = Δ P k
where the residual { P } k matrix is as follows:
Δ P k = P K + F D ¯ k 1 D ¯ k 1
Iterations continue until the error ε is minimized to a certain value, i.e., 10 5 at k . The iteration is described as follows:
| | D ¯ k | | ε  
For more detailed algebra, chapter 9 of the book titled “Differential Quadrature and Differential Quadrature Based Element Methods: Theory and Applications” by Wang [34] can be seen.

2.3.2. Viscoelasticity

For the generalized Maxwell model, stress can be written in the following form:
σ t = C 0   ε t C 1 ε i t
where C 0 and C 1 are creep compliance coefficients and ε t and ε i t are elastic and inelastic strain terms of the generalized Maxwell model. The stress tensor of the material is described in terms of the creep compliance function, which means the specimen would have a stress value at time t when a constant strain is applied as a Heaviside step function at time t.
σ t = C   ε t σ i t = C 1   ε t ε i t = ξ   ε ˙ i t σ t = σ i t + σ t
where C 0 = C + C 1 are the creep compliance coefficients of the generalized Maxwell model and ξ is the viscous damping constant. Using these equations, the first-order differential equation is obtained as follows:
ε t r ε i t r ε i t = 0
where r = C 1 ξ , and by using the Laplace transformation, the first-order differential equation is solved for inelastic strain, and it gives the following:
ε i t = t 1 r   e t s / r   ε s   d s
Substituting this solution into Equation (A5), stress can be described in terms of elastic and inelastic strain as follows:
σ t = C 0 ε t C 1 r t   e t s / r   ε s   d s
where s represents the integrant time variable for the Laplace transformation rule.
To model the viscoelastic material in this problem, the multiaxial case of the time-dependent generalized Maxwell model is used. As a convolution equation over time that requires a high computational cost, the below equation is rewritten in incremental form, in which time integration is used numerically.
σ i j t + Δ t = C i j k l   [ ε i j t + Δ t C 1 C 0 ξ t + Δ t e t + Δ t s r ε k l s d s ]
where ξ and C i j k l represent the viscous damping constant and constitutive matrix, respectively. Dividing the integral into two parts, in which the first part is from - to t and the second part is from t to t + t . Then, by applying the trapezoidal rule to the second part and rearranging the equation, the stress term can be obtained as follows:
[ σ i j ] t + Δ t = [ σ i j ] t   . e Δ t / r C i j k l . ε k l t . e Δ t r . 1 + Δ t 2 C 1 C 0 ξ + C i j k l   . ε k l t + Δ t . 1 Δ t 2 C 1 C 0 ξ
The first and second terms of the right-hand side of Equation (33) represent internal forces for the viscoelastic layer at the previous time step, which are t and are calculated by using the kinematic and constitutive relations (4-5) by using calculated displacement values as shown in Equation (14) at the previous time step. Then, these vectors are extracted from the external force vector in Equation (22) as shown in Equation (A3) in Appendix A and are represented as F 1 and F 2 . The third term of the right-hand side of Equation (33) represents the multiplication of stress terms and the first variation of strain terms at the initial time step is t + t for the viscoelastic layer in energy (Equation (7)). For more detailed algebra, chapter 6 of “Nonlinear Dynamics of Structures” by Oller can be seen [35].

3. Validation Studies

In this section, three different problems are investigated to determine the accuracy and validity of modeling geometric nonlinearities, moving loads, and viscoelasticity, as follows:

3.1. Example 1

The problem setup with simple supported boundary conditions and moving loads can be seen in Figure 2.
In this section, analyses are made to determine the accuracy and validity of the present method. The results are compared with the literature in which the forced vibration of a beam due to a moving load is solved by Eftekhari [10]. In the problem, a concentrated F load moves from one end to the other with a constant velocity for each case in t time, and the response of the beam center is investigated. In order to model point load as a distributed load, an exponential function behaving similar to the Dirac-delta function is introduced:
F t = F α π exp [ ( x i v t ) 2 α 2 ] exp [ ( x n v t ) 2 α 2 ] T
where F, v, t and α represent the amplitude of the point load, velocity of the point load, time, and regularization parameter, respectively, for the beam with n grid points. For a more detailed expression of the function, an article by Eftekhari [10] can be seen.
In the graphs below, w c d and w c s represent the dynamic response of the beam center and static displacement of the beam center, respectively. Dividing w c d with w c s results in the normalization of the results. Also, the time required for the point load to travel through the beam is divided by the period for normalization as well. The results are compared to the exact solution values in the article. Data digitization is used to determine values. The geometric and mechanical properties of the sandwich beam are considered in the problem under a moving load with an amplitude of 4.45 N, as given in Table 1 below:
In Figure 3, a convergence study is also made for 31, 41, and 51 grid points. However, instabilities are observed for 21 grid points, which are not shown in the graphs. Different from the validation article, the Dirac-delta coefficient and time interval values are 0.005 and 200, respectively. It can be seen that convergence of the grid points occurred earlier for higher velocity values.

3.2. Example 2

Now we consider Example 1 by adding geometrical nonlinearities. The results are compared with the explicit FEM analyses performed by ANSYS. In GDQM analysis, 51 grid points are used, and a forced vibration response is obtained. Equivalent mass damping with ζ = 0.2 is used to simulate viscoelastic behavior. However, normalization of time and displacement is not made in this case. The time history of the beam center is observed again. As can be seen in Figure 4, a good agreement is observed.

3.3. Example 3

In this example, the nonlinear transient response of a semi-circular viscoelastic arch subjected to pressure load in the reverse normal direction, as shown in Figure 5, is analyzed. Again, 51 grid points are used to investigate the mid-point displacement of the arch. A load is applied to the model with the Heaviside step function. In addition, load is applied in two steps as half amplitude for each 12,501-time interval. The circular-arch modeling parameters are presented in Table 2.
The finite element model is built by solid elements, and the external load is applied as a pressure load. The model consists of 600 solid elements, 1212 nodes, and 7272 degrees of freedom, and it is constrained so that only in-plane motion is allowed. On the other hand, the pressure load is multiplied by the beam’s width and converted to a distributed load acting on the circular arch in a reverse normal direction.
Here, ST and LT represent short-term and long-term, respectively. Transverse deflection of the central point and damping are observed under the load. An acceptable agreement is observed from the results, as shown in Figure 6.

4. Case Studies

In this section, nonlinear dynamic analyses of a three-layered curved sandwich beam are presented. Aluminum material is used for the face layers, while viscoelastic material is used for the core layer. Fixed boundary conditions are applied on both ends. The moving point load is applied in two load steps on the top layer towards the arc center, as shown in Figure 7, and the response of the beam center on the top layer is calculated.

4.1. Sandwich Beam Subjected to Moving Load (Forced Vibration)

For this case, the thicknesses of the top, core, and bottom layers are fixed at 0.3 mm, 0.4 mm, and 0.3 mm, respectively, which results in a total 1 mm arch thickness. For all parametric analyses, the Dirac-delta coefficient is 0.005. The constant average acceleration method is used in the Newmark-Beta formulation. Viscoelastic material, which is introduced in Example 3. The geometry and material properties of the problem and the parameters of the solution procedure are listed in Table 3 below:
In total, 2146 parametric analyses are conducted in order to investigate the effect of the curvature of the beam under certain load amplitudes and velocity values. Maximum displacement and oscillation in the transverse direction of the beam center are mainly investigated. The length of the beam arc is held constant, while the radius of curvature and arc angle are parametrically varied. The results have converged for 55 grid points for the problem considered. The amplitude of the loads is chosen to satisfy the requirement that the maximum displacement-total thickness ratio of 1 is exceeded to observe the nonlinear effects.
Increasing displacement is expected with increasing radius of curvature, i.e., the straighter the beam, the larger the displacement in linear analyses. However, this is not observed for all cases considered, as seen in Figure 8, due to the nonlinear effects. Moreover, with decreasing curvature and increasing load amplitude, the frequency of oscillations increases in all cases. Beams recover to the undeformed position when the load exits the beams, only in some cases. It is predicted that the load travel frequency does not match the natural frequency of the beam.
From Figure 9a,b, it can be clearly seen that the displacement of the beam center first increases with the increasing value of the load velocity, then decreases. On the other hand, displacements increase with the load amplitude, as expected. It can be said that there is a nonlinear relationship between the displacement, velocity, and load amplitude.

4.2. Sandwich Beam Subjected to Moving Load (Free Vibration)

In this section, the free oscillation of the beam is investigated after the load is removed. Different from Section 4.1, the thicknesses of the top, core, and bottom layers are fixed at 1 mm, which results in a total of 3 mm of arch thickness. The radius of curvature is 1 m, and the long-term shear modulus of the viscoelastic material is 3.5 mPa. In total, 30,576 parametric analyses were done to observe cancellation and resonance phenomena for cantilever beams. In the literature, simply supported beams are widely studied for this phenomenon [14].
Figure 10 shows the variation of maximum mid-point displacement with load amplitude and velocity after the load leaves the beam. It is observed that for loads with relatively lower amplitude, both the cancellation and resonance velocity values decrease with the increasing load. In contrast, the resonance and cancellation velocity values increase with increasing load for the loads with relatively higher amplitude. It can be said that the ratio of load effect frequency to structural natural frequency is not the only parameter that determines the cancellation and resonance phenomena but also the amplitude of the load. This result can be supported by the effect of force amplitude on nonlinear frequency response functions due to geometrical nonlinear effects, which can be found in the literature. The displacements for minimum and maximum load cases, i.e., 1 N and 26 N, respectively, are presented in Figure 11 for a closer inspection.
Time histories of the beam mid-point under load with a 26 N amplitude where cancellation and responance are observed are also presented in Figure 12 and Figure 13, respectively. Time history variation with the first, second, third, fourth, fifth, and sixth cancellation velocity values in Figure 11 corresponds to Figure 12a–f, respectively. Furthermore, as seen in Figure 12, there are n number of oscillations during forced vibration in which the time history is symmetrical, which corresponds to the nth value of cancellation velocity. The travel time of the load from one end to the other is one-third of the total analysis time; hence, the free vibration duration is twice as long as the forced vibration duration. However, instead of a complete cancellation, residual responses are observed, as in [14]. Another interesting point is that in the free vibration zone, higher modes are observed with increasing cancellation velocity values, as Figure 12 shows. Time history variation with the first, second, third, fourth, and fifth resonance velocity values in Figure 11 corresponds to Figure 13a–e, respectively. It can be seen that higher modes are not observed for all resonance velocity values.

4.3. Sandwich Beam Subjected to Multiple Moving Loads (Forced Vibration)

However, when the beam is excitated by eight moving loads with a 26 N load amplitude and a distance of the beam length between all eight loads, the time history of the mid-point displacement is shown in Figure 14. Resonance phenomena are observed at 118.4 m/s velocity, which does not coincide with the resonance condition for a single moving load. The resonance velocity for the forced vibration response for eight consecutive moving loads is represented by a vertical line in Figure 11 that does not coincide with the resonance velocity of the free vibration response of a single moving load. This is possibly due to the nonlinear effects associated with the large strain analyses of the structure.

5. Conclusions

Vibration analyses due to the moving point load are carried out by both GDQM. Also, validation analyses are conducted by ANSYS and literature, and good agreement is observed even though beam geometry is modeled in polar coordinates and first-order shear deformation theory is used for the kinematic results.
  • Increasing the radius of curvature does not cause an increase in mid-point displacement for all cases within the forced vibration.
  • Resonance and cancellation velocity values change with the load amplitude due to the nonlinear effects.
  • The load travels with the nth cancellation velocity, causing n number of oscillations during the forced vibration, in which the time history of the mid-point is symmetrical.
  • The resonance velocity for the beam is not the same for the free vibration condition after leaving a single moving load and the forced vibration condition with multiple moving loads.
For a wider range of parametric analyses, i.e., deeper curved beams, curvilinear coordinates and higher-order shear deformation theory can be used for improvement.

Author Contributions

All authors contributed to the study’s conception and design. Material preparation, data collection, and analysis were performed by M.M.S., O.D. and A.A. The first draft of the manuscript was written by M.M.S., O.D. and A.A., and all authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The equations of motion mentioned in Section 2.1 are given as follows:
Q θ r 1 R 1 N θ θ 1 R 1 + 1 h 2 1 2 R 2 Q θ r 2 1 h 2 R 2 M θ θ 2 θ 1 2 R 2 N θ θ 2 θ = R 1 u ¨ 01 I ¯ 0 1 + R 2 4 u ¨ 01 + 2 u ¨ 03 2 h 1 φ ¨ 1 + h 3 φ ¨ 3 12 I ¯ 0 2
Q θ r 3 R 3 N θ θ 3 R 3 1 h 2 + 1 2 R 2 Q θ r 2 + 1 h 2 R 2 M θ θ 2 θ 1 2 R 2 N θ θ 2 θ = R 3 u ¨ 03 I ¯ 0 3 + R 2 4 u ¨ 03 + 2 u ¨ 01 h 1 φ ¨ 1 + 2 h 3 φ ¨ 3 12 I ¯ 0 2
1 R 1 Q θ r 1 θ + N θ θ 1 R 1 1 R 1 2 N θ θ 1 w 01 θ 2 N θ θ 1 R 1 2   2 w 01 θ 2 + M θ θ 2 h 2 R 2 + N θ θ 2 2 R 2 + N r r 2 h 2 1 2 R 2 Q θ r 2 θ 1 h 2 R 2 R θ r 2 θ + N r r 2 w 01 h 2 2 N r r 2 w 03 h 2 2 1 h 2 R 2 2 M θ θ 2 w 01 θ 2 1 4 R 2 2 N θ θ 2 w 01 θ 2 1 h 2 2 R 2 2 S θ θ 2 w 01   θ 2 1 4 R 2 2 N θ θ 2 w 03 θ 2 + 1 h 2 2 R 2 2 S θ θ 2 w 03 θ 2 1 h 2 R 2 2 M θ θ 2 2 w 01 θ 2 1 4 R 2 2 N θ θ 2 2 w 01 θ 2 1 4 R 2 2 N θ θ 2 2 w 03 θ 2 + S θ θ 2 h 2 2 R 2 2 2 w 03 θ 2 S θ θ 2 h 2 2 R 2 2   2 w 01 θ 2 = R 1 w ¨ 01 I ¯ 0 1 + R 2 2 w ¨ 01 + w ¨ 03 6 I ¯ 0 2 q + F 1 + F 2
1 R 3 Q θ r 3 θ + N θ θ 3 R 3 1 R 3 2   N θ θ 3 w 03   θ 2 N θ θ 3     R 3 2   2 w 03 θ 2 M θ θ 2 h 2 R 2 + N θ θ 2 2 R 2 N r r 2 h 2 1 2 R 2 Q θ r 2 θ + 1 h 2 R 2 R θ r 2 θ N r r 2 w 01 h 2 2 + N r r 2 w 03 h 2 2 + 1 h 2 R 2 2 M θ θ 2 w 03 θ 2 1 4 R 2 2 N θ θ 2 w 03 θ 2 1 h 2 2 R 2 2 S θ θ 2 w 03 θ 2 1 4 R 2 2 N θ θ 2 w 01 θ 2 + 1 h 2 2 R 2 2 S θ θ 2 w 01 θ 2 + M θ θ 2 h 2 R 2 2 2 w 03 θ 2 N θ θ 2   4 R 2 2 2 w 03 θ 2 + S θ θ 2 h 2 2 R 2 2 2 w 01 θ 2 N θ θ 2 4 R 2 2 2 w 01 θ 2 S θ θ 2 h 2 2 R 2 2 2 w 03 θ 2 = R 3 w ¨ 03 I ¯ 0 3 + R 2 2 w ¨ 03 + w ¨ 01 6 I ¯ 0 2
Q θ r 1 h 2 1 R 1 M θ θ 1 θ 1 2 h 2 1 4 R 2 h 1 Q θ r 2 + h 1 2 h 2 R 2 M θ θ 2 θ + h 1 4 R 2 N θ θ 2 θ = R 1 φ ¨ 1 I ¯ 2 1 R 2 h 1 + 4 u ¨ 01 + 2 u ¨ 03 2 h 1 φ ¨ 1 + h 3 φ ¨ 3 24 I ¯ 0 2
Q θ r 3 h 2 1 R 3 M θ θ 3 θ 1 2 h 2 + 1 4 R 2 h 3   Q θ r 2 + h 3 2 h 2 R 2 M θ θ 2 θ h 3 4 R 2 N θ θ 2 θ = R 3 φ ¨ 3 I ¯ 2 3 + R 2 h 3 + 2 u ¨ 01 + 4 u ¨ 03 h 1 φ ¨ 1 + 2 h 3 φ ¨ 3 24 I ¯ 0 2
In addition, the natural boundary conditions for u 01 ,   u 03 ,   w 01 ,   w 03 ,   φ 1 and φ 3 are derived respectively as below:
N θ θ 1 R 1 + N θ θ 2 2   R 2 + M θ θ 2 h 2   R 2 = 0
N θ θ 3 R 3 + N θ θ 2 2   R 2 M θ θ 2 h 2 R 2 = 0
M θ θ 1 R 1 h 1 M θ θ 2 2 h 2 R 2 h 1 N θ θ 2 4 R 2 + M θ θ 2 h 2 R 2 2 w 01 θ + N θ θ 2 4 R 2 2 w 01 θ + S θ θ 2 h 2 2 R 2 2   w 01 θ + N θ θ 2   4 R 2 2 w 03 θ S θ θ 2 h 2 2 R 2 2 w 03 θ = 0
M θ θ 3 R 3 h 3 M θ θ 2 2 h 2 R 2 h 1 N θ θ 2 4 R 2 M θ θ 2 h 2 R 2 2 w 03 θ + N θ θ 2 4 R 2 2 w 03 θ + S θ θ 2 h 2 2 R 2 2   w 03 θ + N θ θ 2   4 R 2 2 w 01 θ S θ θ 2 h 2 2 R 2 2 w 01 θ = 0
Q θ r 1 R 1 + N θ θ 1     R 1 2   w 01 θ + Q θ r 2 R 2 + R θ r 2 h 2   R 2 = 0
Q θ r 3 R 3 + N θ θ 3     R 3 2   w 03 θ + Q θ r 2 R 2 + R θ r 2 h 2   R 2 = 0
where m and n represent θ and r for i = 1,2 , 3 . Forces and moment resultants can be described as:
N m m i = h i 2 h i 2 σ m m i   d z i ,         M m n i = h i 2 h i 2 z i σ m m i   d z i S m m i = h i 2 h i 2 z 2 i σ m m i   d z i ,       Q m n i = h 2 2 h 2 2 σ m n i   d z i R m n i = h 2 2 h 2 2 σ m n i   d z i ,         I j i = h i 2 h i 2 ρ i z j d z i

References

  1. Singh, G.; Sharma, A.K.; Venkateswara Rao, G. Large-amplitude free vibrations of beams—A discussion on various formulations and assumptions. J. Sound Vib. 1990, 142, 77–85. [Google Scholar] [CrossRef]
  2. Gafsi, W.; Najar, F.; Choura, S.; El-Borgi, S. Confinement of Vibrations in Variable-Geometry Nonlinear Flexible Beam. Shock Vib. 2014, 2014, 687340. [Google Scholar] [CrossRef]
  3. Bert, C.W.; Malik, M. Differential Quadrature Method in Computational Mechanics: A Review. Appl. Mech. Rev. 1996, 49, 1–28. [Google Scholar] [CrossRef]
  4. Jang, S.K.; Bert, C.W.; Striz, A.G. Application of differential quadrature to static analysis of structural components. Int. J. Numer. Meth. Eng. 1989, 28, 561–577. [Google Scholar] [CrossRef]
  5. Feng, Y.; Bert, C.W. Application of the quadrature method to flexural vibration analysis of a geometrically nonlinear beam. Nonlinear Dyn. 1992, 3, 13–18. [Google Scholar] [CrossRef]
  6. Karami, G.; Malekzadeh, P. A new differential quadrature methodology for beam analysis and the associated differential quadrature element method. Comput. Methods Appl. Mech. Eng. 2002, 191, 3509–3526. [Google Scholar] [CrossRef]
  7. Zhong, H.; Guo, Q. Nonlinear Vibration Analysis of Timoshenko Beams Using the Differential Quadrature Method. Nonlinear Dyn. 2003, 32, 223–234. [Google Scholar] [CrossRef]
  8. Ghasemi, A.R.; Mohandes, M. The effect of finite strain on the nonlinear free vibration of a unidirectional composite Timoshenko beam using GDQM. Adv. Aircr. Spacecr. Sci. 2016, 3, 379–397. [Google Scholar] [CrossRef]
  9. Mahmoud, A.A.; Esmaeel, R.A.; Nassar, M.M. Application of the generalized differential quadrature method to the free vibrations of delaminated beam plates. Eng. Mech. 2007, 14, 431–441. [Google Scholar]
  10. Eftekhari, S.A. A Differential Quadrature Procedure with Regularization of the Dirac-delta Function for Numerical Solution of Moving Load Problem. Lat. Am. J. Solids Struct. 2015, 12, 1241–1265. [Google Scholar] [CrossRef]
  11. Eftekhari, S.A. A simple and accurate mixed Ritz-DQM formulation for free vibration of rectangular plates involving free corners. Ain Shams Eng. J. 2016, 7, 777–790. [Google Scholar] [CrossRef]
  12. Eftekhari, S.A. A simple and systematic approach for implementing boundary conditions in the differential quadrature free and forced vibration analysis of beams and rectangular plates. J. Solid Mech. 2015, 7.4, 7.4,374–399. [Google Scholar]
  13. Eftekhari, S.A. A differential quadrature procedure for linear and nonlinear steady state vibrations of infinite beams traversed by a moving point load. Meccanica 2016, 51, 2417–2434. [Google Scholar] [CrossRef]
  14. Yang, Y.-B.; Wu, C.-M.; Yau, J.-D. Dynamic response of a horizontally curved beam subjected to vertical and horizontal moving loads. J. Sound Vib. 2001, 242, 519–537. [Google Scholar] [CrossRef]
  15. Hajianmaleki, M.; Qatu, M.S. Static and vibration analyses of thick, generally laminated deep curved beams with different boundary conditions. Compos. Part B Eng. 2012, 43, 1767–1775. [Google Scholar] [CrossRef]
  16. Poojary, J.; Roy, S.K. In-plane vibration of curved beams subjected to moving loads using finite element method. J. Phys. Conf. Ser. 2019, 1240, 012048. [Google Scholar] [CrossRef]
  17. Kurtaran, H. Large displacement static and transient analysis of functionally graded deep curved beams with generalized differential quadrature method. Compos. Struct. 2015, 131, 821–831. [Google Scholar] [CrossRef]
  18. Kurtaran, H. Geometrically nonlinear transient analysis of thick deep composite curved beams with generalized differential quadrature method. Compos. Struct. 2015, 128, 241–250. [Google Scholar] [CrossRef]
  19. Nikkhoo, A. Application of differential quadrature method to investigate dynamics of a curved beam structure acted upon by a moving concentrated load. Indian J. Sci. Technol. 2012, 5, 8. [Google Scholar] [CrossRef]
  20. Kang, K.; Bert, C.W.; Striz, A.G. Vibration Analysis of Horizontally Curved Beams with Warping Using DQM. J. Struct. Eng. 1996, 122, 657–662. [Google Scholar] [CrossRef]
  21. Akgün, G.; Kurtaran, H. Geometrically nonlinear transient analysis of laminated composite super-elliptic shell structures with generalized differential quadrature method. Int. J. Non-Linear Mech. 2018, 105, 221–241. [Google Scholar] [CrossRef]
  22. Malik, M.; Bert, C.W. Implementing multiple boundary conditions in the dq solution of higher-order pdes: Application to free vibration of plates. Int. J. Numer. Meth. Eng. 1996, 39, 1237–1258. [Google Scholar] [CrossRef]
  23. Fung, T.C. Stability and accuracy of differential quadrature method in solving dynamic problems. Comput. Methods Appl. Mech. Eng. 2002, 191, 1311–1331. [Google Scholar] [CrossRef]
  24. Bacciocchi, M.; Eisenberger, M.; Fantuzzi, N.; Tornabene, F.; Viola, E. Vibration analysis of variable thickness plates and shells by the Generalized Differential Quadrature method. Compos. Struct. 2016, 156, 218–237. [Google Scholar] [CrossRef]
  25. Soares da Costa Azevêdo, A.; Soares, A.S.C. Dynamic analysis of a viscoelastic Timoshenko beam. In Proceedings of the 24th ABCM International Congress of Mechanicl Engineering, ABCM, Curitiba, Brazil, 15 February 2017. [Google Scholar]
  26. Chen, W.-H.; Lin, T.-C. Dynamic analysis of viscoelastic structures using incremental finite element method. Eng. Struct. 1982, 4, 271–276. [Google Scholar] [CrossRef]
  27. Shafei, E.; Faroughi, S.; Rabczuk, T. Nonlinear transient vibration of viscoelastic plates: A NURBS-based isogeometric HSDT approach. Comput. Math. Appl. 2021, 84, 1–15. [Google Scholar] [CrossRef]
  28. Tolpekina, T.; Pyckhout-Hintzen, W.; Persson, B.N.J. Linear and Nonlinear Viscoelastic Modulus of Rubber. Lubricants 2019, 7, 22. [Google Scholar] [CrossRef]
  29. Arikoglu, A. A new fractional derivative model for linearly viscoelastic materials and parameter identification via genetic algorithms. Rheol. Acta 2014, 53, 219–233. [Google Scholar] [CrossRef]
  30. Demir, O.; Balkan, D.; Peker, R.C.; Metin, M.; Arikoglu, A. Vibration analysis of curved composite sandwich beams with viscoelastic core by using differential quadrature method. J. Sandw. Struct. Mater. 2020, 22, 743–770. [Google Scholar] [CrossRef]
  31. Arikoglu, A.; Ozkol, I. Vibration Analysis of Composite Sandwich Plates by the Generalized Differential Quadrature Method. AIAA J. 2012, 50, 620–630. [Google Scholar] [CrossRef]
  32. Taskin, M.; Arikoglu, A.; Demir, O. Vibration and Damping Analysis of Sandwich Cylindrical Shells by the GDQM. AIAA J. 2019, 57, 3040–3051. [Google Scholar] [CrossRef]
  33. Chopra, A.K. Dynamics of Structures: Theory and Applications to Earthquake Engineering; Prentice Hall: Upper Saddle River, NJ, USA, 2015; ISBN 978-0-273-77424-2. [Google Scholar]
  34. Wang, X. Differential Quadrature and Differential Quadrature Based Element Methods: Theory and Applications; Elsevier Science: Amsterdam, The Netherlands, 2015; ISBN 978-0-12-803081-3. [Google Scholar]
  35. Oller, S. Nonlinear Dynamics of Structures; Lecture Notes on Numerical Methods in Engineering and Sciences; Springer: Berlin/Heidelberg, Germany, 2014; ISBN 978-3-319-05193-2. [Google Scholar]
Figure 1. Coordinate system, configuration, and representation of beam layers.
Figure 1. Coordinate system, configuration, and representation of beam layers.
Symmetry 16 00238 g001
Figure 2. Moving point load on the simply supported beam.
Figure 2. Moving point load on the simply supported beam.
Symmetry 16 00238 g002
Figure 3. Comparison of the results for two different velocity values: (a) 15.6 m/s and (b) 46.8 m/s.
Figure 3. Comparison of the results for two different velocity values: (a) 15.6 m/s and (b) 46.8 m/s.
Symmetry 16 00238 g003
Figure 4. Comparison of the results for 31.2 m/s.
Figure 4. Comparison of the results for 31.2 m/s.
Symmetry 16 00238 g004
Figure 5. Distributed load on the fixed circular arch.
Figure 5. Distributed load on the fixed circular arch.
Symmetry 16 00238 g005
Figure 6. Time history of mid-point displacements.
Figure 6. Time history of mid-point displacements.
Symmetry 16 00238 g006
Figure 7. Curved sandwich beam representation under moving load.
Figure 7. Curved sandwich beam representation under moving load.
Symmetry 16 00238 g007
Figure 8. Time history of mid-point displacements: (a) 8 m/s, (b) 16 m/s, (c) 32 m/s, and (d) 64 m/s.
Figure 8. Time history of mid-point displacements: (a) 8 m/s, (b) 16 m/s, (c) 32 m/s, and (d) 64 m/s.
Symmetry 16 00238 g008aSymmetry 16 00238 g008b
Figure 9. Present method results: (a) R c = 0.75   m and (b) R c = 2   m .
Figure 9. Present method results: (a) R c = 0.75   m and (b) R c = 2   m .
Symmetry 16 00238 g009aSymmetry 16 00238 g009b
Figure 10. Variation of maximum amplitude with force and velocity.
Figure 10. Variation of maximum amplitude with force and velocity.
Symmetry 16 00238 g010
Figure 11. Variation of maximum amplitude with velocity.
Figure 11. Variation of maximum amplitude with velocity.
Symmetry 16 00238 g011
Figure 12. Time history of mid-point displacement under 26 N: (a) 51 m/s, (b) 26.4 m/s, (c) 17.6 m/s, (d) 13.08 m/s, (e) 10.48 m/s, and (f) 8.72 m/s.
Figure 12. Time history of mid-point displacement under 26 N: (a) 51 m/s, (b) 26.4 m/s, (c) 17.6 m/s, (d) 13.08 m/s, (e) 10.48 m/s, and (f) 8.72 m/s.
Symmetry 16 00238 g012
Figure 13. Time history of mid-point displacement under 26 N: (a) 123.4 m/s, (b) 36.5 m/s, (c) 22 m/s, (d) 15.36 m/s, and (e) 11.76 m/s.
Figure 13. Time history of mid-point displacement under 26 N: (a) 123.4 m/s, (b) 36.5 m/s, (c) 22 m/s, (d) 15.36 m/s, and (e) 11.76 m/s.
Symmetry 16 00238 g013
Figure 14. Time history of mid-point displacement under eight moving loads under 26 N.
Figure 14. Time history of mid-point displacement under eight moving loads under 26 N.
Symmetry 16 00238 g014
Table 1. Description of the moving load problem.
Table 1. Description of the moving load problem.
Length 0.1016   m Density 10,686.9   k g / m 3
Width × Thickness 0.00635   ×   0.00635   m 2 Dirac-delta Coefficient 0.005
Young’s Modulus 2.068   ×   10 11   P a Time Intervals200
Table 2. Problem data for the circular arc.
Table 2. Problem data for the circular arc.
Length 0.1   m Creep Decay Exponent15.7
Width × Thickness 0.01   m   ×   0.0004   m Density 1500   k g / m 3
Radius of Curvature 1   m ST Poisson’s Ratio 0.25
ST Shear Modulus 9.5   ×   10 6   P a Distributed Load 0.01 N / m
LT Shear Modulus 0.55   ×   10 6   P a Timestep4 × 10−5
Table 3. Problem data for the present method.
Table 3. Problem data for the present method.
Beam Arc Length 0.1   m Width × Thickness 0.002   ×   0.001   m 2
Radius of Curvature 2   m 0.75   m Timestep4 × 10−6
Young’s Modulus (Faces) 70   ×   10 9   P a Density (Faces) 2700   k g / m 3
Poisson’s Ratio (Faces)0.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Serveren, M.M.; Demir, O.; Arikoglu, A. Nonlinear Dynamic Analysis of a Curved Sandwich Beam with a Time-Dependent Viscoelastic Core Using the Generalized Differential Quadrature Method (GDQM). Symmetry 2024, 16, 238. https://doi.org/10.3390/sym16020238

AMA Style

Serveren MM, Demir O, Arikoglu A. Nonlinear Dynamic Analysis of a Curved Sandwich Beam with a Time-Dependent Viscoelastic Core Using the Generalized Differential Quadrature Method (GDQM). Symmetry. 2024; 16(2):238. https://doi.org/10.3390/sym16020238

Chicago/Turabian Style

Serveren, Mehmet Mert, Ozgur Demir, and Aytac Arikoglu. 2024. "Nonlinear Dynamic Analysis of a Curved Sandwich Beam with a Time-Dependent Viscoelastic Core Using the Generalized Differential Quadrature Method (GDQM)" Symmetry 16, no. 2: 238. https://doi.org/10.3390/sym16020238

APA Style

Serveren, M. M., Demir, O., & Arikoglu, A. (2024). Nonlinear Dynamic Analysis of a Curved Sandwich Beam with a Time-Dependent Viscoelastic Core Using the Generalized Differential Quadrature Method (GDQM). Symmetry, 16(2), 238. https://doi.org/10.3390/sym16020238

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