• Steel and Composite Structures
• /
• v.53 no.1
• /
• pp.1-27
• /
• 2024

This work presents the effectiveness of differential quadrature shape functions (i.e., Lagrange interpolation polynomial, Cardinal sine function, Delta Lagrange kernel and Regularized Shannon kernel) in the solution of nonlinear vibration of multilayers piezoelectric plates with nonlinear elastic support. A piezoelectric composite laminated plate is rested on nonlinear Winkler and Visco-Pasternak elastic foundations problems. Based on 3D elasticity theory and piezoelectricity, the governing equations of motion are derived. Differential quadrature methods based on four shape functions are presented as numerical techniques for solving this problem. The perturbation method is implemented to solve the obtained nonlinear eigenvalue problem. A MATLAB code is written for each technique for solving this problem and extract the numerical results. To validate these methods, the computed results are we compare with the previous exact results. In addition, parametric analyses are offered to investigate the influence of length to thickness ratio, elastic foundation parameters, various boundary conditions, and piezoelectric layers thickness on the natural frequencies and mode shapes. Consequently, it is discovered that the obtained results via the proposed schemes can be applied in structural health monitoring.

• Journal of the Korean Mathematical Society
• /
• v.58 no.6
• /
• pp.1433-1448
• /
• 2021

In the present work, we investigate a viscoelastic equation involving a logarithmic nonlinear source term. After proving the existence of solutions, we establish a general decay estimate of the solution using energy estimates and theory of convex functions. This result extends and complements some previous results of [9, 21].

• Structural Engineering and Mechanics
• /
• v.35 no.4
• /
• pp.387-407
• /
• 2010

In this paper, using perturbation and Galerkin method, the response of a resonant viscoelastic microbeam to an electric actuation is obtained. The microbeam is under axial load and electrical load. It is assumed that midplane is stretched, when the beam is deflected. The equation of motion is derived using the Newton's second law. The viscoelastic model is taken to be the Kelvin-Voigt model. In the first section, the static deflection is obtained using the Galerkin method. Exact linear symmetric mode shape of a straight beam and its deflection function under constant transverse load are used as admissible functions. So, an analytical expression that describes the static deflection at all points is obtained. Comparing the result with previous research show that using deflection function as admissible function decreases the computation errors and previous calculations volume. In the second section, the response of a microbeam resonator system under primary and secondary resonance excitation has been obtained by analytical multiple scale perturbation method combined with the Galerkin method. It is shown, that a small amount of viscoelastic damping has an important effect and causes to decrease the maximum amplitude of response, and to shift the resonance frequency. Also, it shown, that an increase of the DC voltage, ratio of the air gap to the microbeam thickness, tensile axial load, would increase the effect of viscoelastic damping, and an increase of the compressive axial load would decrease the effect of viscoelastic damping.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.15 no.6
• /
• pp.3435-3441
• /
• 2014

A dynamic modeling and prediction of polyurethane foam material, which is used as the seat in vehicles is very important for improving the ride quality of vehicle occupants. In this study, parameters to define the nonlinear stiffness and time-variant characteristics of the viscoelasticity of polyurethane foam were obtained using a static compression test. Polynomial functions and convolution integral were used to model the nonlinear and viscoelastic characteristics of polyurethane foam mathematically. The dynamic behaviors excited by the seat floor displacement were analyzed using a numerical integration method for the nonlinear vibration model. As a result, the viscoelastic characteristics of polyurethane foam was found to be an important parameter for improving the ride quality.

• Journal of the Korean Society for Railway
• /
• v.18 no.3
• /
• pp.194-202
• /
• 2015

In this paper, dynamic modeling with viscoelastic properties of a human body resting on a seat is presented to quantitatively analyze ride quality of passengers exposed to vertical vibrations. In describing the motions of a seated body, a 5 degree-of-freedom multibody model from the literature is investigated. The viscoelastic characteristics of seats used in railway vehicles are mathematically formulated with nonlinear stiffness characteristics and convolution integrals representing time delay terms. Transfer functions for the floor input are investigated and it is found that these are different in accordance with the input magnitude due to nonlinear characteristics of the seat. Measured floor input at the railway vehicle is used to analyze realistic human vibration characteristics. Frequency weighted RMS acceleration values are calculated and the effects of the seat design parameters on the frequency weighted RMS acceleration values are presented.

• Proceedings of the Korean Society of Rheology Conference
• /
• 2006.06a
• /
• pp.143-147
• /
• 2006

• Korea-Australia Rheology Journal
• /
• v.19 no.1
• /
• pp.35-42
• /
• 2007

The effects of particle size distribution on the magnetorheological response of inverse ferrofluids was investigated using controlled mixtures of two monodisperse non-magnetisable powders of sizes $4.6\;{\mu}m\;and\;80{\mu}m$ at constant volume fraction of 30%, subjected to large amplitude oscillatory shear flow. In the linear viscoelastic regime (pre-yield region), it was found that the storage and loss moduli were dependent on the particle size as well as the proportion of small particles, with the highest storage modulus occurring for the monodisperse small particles. In the nonlinear regime (post yield region), Fourier analysis was used to compare the behaviour of the $1^{st}\;and\;3^{rd}$ harmonics ($I_{1}\;and\;I_{3}\;respectively$) as well as the fundamental phase angle as functions of the applied strain amplitude. The ratio of $I_{3}/I_{1}$ was found to become more pronounced with decreasing particle size as well as with increasing proportion of small particles in the bidisperse mixtures. Furthermore, the phase angle was able to clearly show the transition from solid-like to viscous behaviour. The results suggested that the nonlinear response of a bidisperse IFF is dependent on particle size as well as the proportion of small particles in the system.

• Communications of the Korean Mathematical Society
• /
• v.38 no.3
• /
• pp.943-966
• /
• 2023

In this paper, we study the initial-boundary value problem for viscoelastic wave equations of Kirchhoff type with Balakrishnan-Taylor damping terms in the presence of the infinite memory and external time-varying delay. For a certain class of relaxation functions and certain initial data, we prove that the decay rate of the solution energy is similar to that of relaxation function which is not necessarily of exponential or polynomial type. Also, we show another stability with g satisfying some general growth at infinity.

• Korea-Australia Rheology Journal
• /
• v.15 no.3
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

• Structural Monitoring and Maintenance
• /
• v.6 no.4
• /
• pp.317-346
• /
• 2019

Performance assessment of pavements proves useful, in terms of handling the ride quality, controlling the travel time of vehicles and adequate maintenance of pavements. Roughness profiles provide a good measure of the deteriorating condition of the pavement. For the accurate estimates of pavement roughness from dynamic vehicle responses, vehicle parameters should be known accurately. Information on vehicle parameters is uncertain, due to the wear and tear over time. Hence, condition monitoring of pavement requires the identification of pavement roughness along with vehicle parameters. The present study proposes a scheme which estimates the roughness profile of the pavement with the use of accurate estimates of vehicle parameters computed in parallel. Pavement model used in this study is a two-layer Euler-Bernoulli beam resting on a nonlinear Pasternak foundation. The asphalt topping of the pavement in the top layer is modeled as viscoelastic, and the base course bottom layer is modeled as elastic. The viscoelastic response of the top layer is modeled with the help of the Burgers model. The vehicle model considered in this study is a half car model, fitted with accelerometers at specified points. The identification of the coupled system of vehicle-pavement interaction employs a coupled scheme of an unbiased minimum variance estimator and an optimization scheme. The partitioning of observed noisy quantities to be used in the two schemes is investigated in detail before the analysis. The unbiased minimum variance estimator (MVE) make use of a linear state-space formulation including roughness, to overcome the linearization difficulties as in conventional nonlinear filters. MVE gives estimates for the unknown input and fed into the optimization scheme to yield estimates of vehicle parameters. The issue of ill-posedness of the problem is dealt with by introducing a regularization equivalent term in the objective function, specifically where a large number of parameters are to be estimated. Effect of different objective functions is also studied. The outcome of this research is an overall measure of pavement condition.