• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• Journal of the Korean earth science society
• /
• v.38 no.2
• /
• pp.105-114
• /
• 2017

Two dimensional finite element method with quadrilateral basis functions was applied to the spherical high order filter on the spherical surface limited area domain. The basis function consists of four shape functions which are defined on separate four grid boxes sharing the same gridpoint. With the basis functions, the first order derivative was expressed as an algebraic equation associated with nine point stencil. As the theory depicts, the convergence rate of the error for the spherical Laplacian operator was found to be fourth order, while it was the second order for the spherical Laplacian operator. The accuracy of the new high order filter was shown to be almost the same as those of Fourier finite element high order filter. The two-dimension finite element high order filter was incorporated in the weather research and forecasting (WRF) model as a hyper viscosity. The effect of the high order filter was compared with the built-in viscosity scheme of the WRF model. It was revealed that the high order filter performed better than the built in viscosity scheme did in providing a sharper cutoff of small scale disturbances without affecting the large scale field. Simulation of the tropical cyclone track and intensity with the high order filter showed a forecast performance comparable to the built in viscosity scheme. However, the predicted amount and spatial distribution of the rainfall for the simulation with the high order filter was closer to the observed values than the case of built in viscosity scheme.

• Structural Engineering and Mechanics
• /
• v.12 no.5
• /
• pp.507-526
• /
• 2001

In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

• Nuclear Engineering and Technology
• /
• v.26 no.2
• /
• pp.212-224
• /
• 1994

A numerical method has been developed for studying the dynamics of a flexible cylinder in a coaxial cylindrical duct, immersed in inviscid flow. The unsteady inviscid fluid-dynamic force acting on the oscillating cylinder has been estimated more rigorously by means of a spectral collocation method without simplification of governing equations. This numerical approach is applicable to the system haying wider annular gap and/or shorter length of cylinder as compared to existing potential theory. The governing equation of the unsteady flow was obtained from Laplace equation. The equation of cylinder motion coupled with the fluid motion was discretized by Galerkin's method, from which the dynamic behaviour of the system has been evaluated. The effect of the length of the cylinder and the annular gap on the critical flour velocity, where the system loses stability by buckling, was investigated. To validate the numerical method, the potential flow theory developed by Hobson based on thin film approximation has been improved. Typical results of the present numerical theory on the dynamics and stability of the system are compared with those of available existing theory and the present approximate results. Good agreement was found between the results. It was also found that a nondimensional critical flow velocity becomes larger as increasing the annular gap and decreasing the length of cylinder.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.6
• /
• pp.692-710
• /
• 2018

The Three-dimensional (3-D) dynamical behaviors of a fluid-conveying pipe subjected to vortex-induced vibration are investigated with different internal flow velocity ${\nu}$. The values of the internal flow velocity are considered in both subcritical and supercritical regimes. During the study, the 3-D nonlinear equations are discretized by the Galerkin method and solved by a fourth-order Runge-Kutta method. The results indicate that for a constant internal flow velocity ${\nu}$ in the subcritical regime, the peak Cross-flow (CF) amplitude increases firstly and then decrease accompanied by amplitude jumps with the increase of the external reduced velocity. While two response bands are observed in the In-line (IL) direction. For the dynamics in the lock-in condition, 3-D periodic, quasi-periodic and chaotic vibrations are observed. A variety of CF and IL responses can be detected for different modes with the increase of ${\nu}$. For the cases studied in the supercritical regime, the dynamics shows a great diversity with that in the subcritical regime. Various dynamical responses, which include 3-D periodic, quasi-periodic as well as chaotic motions, are found while both CF and IL responses are coupled while ${\nu}$ is beyond the critical value. Besides, the responses corresponding to different couples of ${\mu}_1$ and ${\mu}_2$ are obviously distinct from each other.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.5
• /
• pp.829-837
• /
• 1992

This paper presents a new approach for the analysis of hot strip rolling processes. The approach is based on the finite element method and capable of predicting velocity field in the strip, temperature field in the strip, temperature field in the roll, and roll pressure. Basic finite element formulations for heat transfer analysis are described with emphasis on the treatment of numerical instability resulting from a standard Galerkin formulation. Comparison with the theoretical solutions found in the literature is made for the evaluation of the accuracy of the temperature solutions. An iterative scheme is developed for dealing with strong correlations between the metal flow characteristics and the thermal behavior of the roll-strip system. A series of process simulations are carried out to investigate the effect of various process parameters including interface friction, interface heat transfer coefficient, roll speed, reduction in thickness, and spray zone. The results are shown and discussed.

• Smart Structures and Systems
• /
• v.16 no.1
• /
• pp.141-162
• /
• 2015

The magneto-rheological visco-elastomer (MRVE) is used as a smart core to control the stochastic micro-vibration of a sandwich plate with supported mass. The micro-vibration response of the sandwich plate with MRVE core and supported mass under stochastic support motion excitations is studied and compared to evaluate the vibration suppression capability. The effects of the supported mass and localized magnetic field on the stochastic micro-vibration response of the MRVE sandwich plate are taken into account. The dynamic characteristics of the MRVE core in micro-vibration are described by a non-homogeneous complex modulus dependent on vibration frequency and controllable by applied magnetic fields. The partial differential equations for the coupled transverse and longitudinal motions of the MRVE sandwich plate with supported mass are derived from the dynamic equilibrium, constitutive and geometric relations. The simplified ordinary differential equations are obtained for the transverse vibration of the MRVE sandwich plate under localized magnetic fields. A frequency-domain solution method for the stochastic micro-vibration response of sandwich plates with supported mass is developed based on the Galerkin method and random vibration theory. The expressions of frequency-response functions, response power spectral densities and root-mean-square velocity responses of the plate in terms of the one-third octave frequency band are obtained for micro-vibration evaluation. Finally, numerical results are given to illustrate the large response reduction capacity of the MRVE sandwich plate with supported mass under stochastic support motion excitations, and the influences of MRVE parameters, supported mass and localized magnetic field placement on the micro-vibration response.

• Smart Structures and Systems
• /
• v.19 no.1
• /
• pp.21-31
• /
• 2017

To control the stochastic vibration of a vibration-sensitive instrument supported on a beam, the beam is designed as a sandwich structure with magneto-rheological visco-elastomer (MRVE) core. The MRVE has dynamic properties such as stiffness and damping adjustable by applied magnetic fields. To achieve better vibration control effectiveness, the optimal bounded parametric control for the MRVE sandwich beam with supported mass under stochastic and deterministic support motion excitations is proposed, and the stochastic and shock vibration suppression capability of the optimally controlled beam with multi-mode coupling is studied. The dynamic behavior of MRVE core is described by the visco-elastic Kelvin-Voigt model with a controllable parameter dependent on applied magnetic fields, and the parameter is considered as an active bounded control. The partial differential equations for horizontal and vertical coupling motions of the sandwich beam are obtained and converted into the multi-mode coupling vibration equations with the bounded nonlinear parametric control according to the Galerkin method. The vibration equations and corresponding performance index construct the optimal bounded parametric control problem. Then the dynamical programming equation for the control problem is derived based on the dynamical programming principle. The optimal bounded parametric control law is obtained by solving the programming equation with the bounded control constraint. The controlled vibration responses of the MRVE sandwich beam under stochastic and shock excitations are obtained by substituting the optimal bounded control into the vibration equations and solving them. The further remarkable vibration suppression capability of the optimal bounded control compared with the passive control and the influence of the control parameters on the stochastic vibration suppression effectiveness are illustrated with numerical results. The proposed optimal bounded parametric control strategy is applicable to smart visco-elastic composite structures under deterministic and stochastic excitations for improving vibration control effectiveness.

• Steel and Composite Structures
• /
• v.44 no.2
• /
• pp.155-169
• /
• 2022

In the present paper an analytical model was developed to study the non-linear vibrations of Functionally Graded Carbon Nanotube (FG-CNT) reinforced doubly-curved shallow shells using the Multiple Scales Method (MSM). The nonlinear partial differential equations of motion are based on the FGM shallow shell hypothesis, the non-linear geometric Von-Karman relationships, and the Galerkin method to reduce the partial differential equations associated with simply supported boundary conditions. The novelty of the present model is the simultaneous prediction of the natural frequencies and their mode shapes versus different curvatures (cylindrical, spherical, conical, and plate) and the different types of FG-CNTs. In addition to combining the vibration analysis with optimization algorithms based on the genetic algorithm, a design optimization methode was developed to maximize the natural frequencies. By considering the expression of the non-dimensional frequency as an objective optimization function, a genetic algorithm program was developed by valuing the mechanical properties, the geometric properties and the FG-CNT configuration of shallow double curvature shells. The results obtained show that the curvature, the volume fraction and the types of NTC distribution have considerable effects on the variation of the Dimensionless Fundamental Linear Frequency (DFLF). The frequency response of the shallow shells of the FG-CNTRC showed two types of nonlinear hardening and softening which are strongly influenced by the change in the fundamental vibration mode. In GA optimization, the mechanical properties and geometric properties in the transverse direction, the volume fraction, and types of distribution of CNTs have a considerable effect on the fundamental frequencies of shallow double-curvature shells. Where the difference between optimized and not optimized DFLF can reach 13.26%.

• Wind and Structures
• /
• v.34 no.2
• /
• pp.173-183
• /
• 2022

The galloping of iced conductors has long been a severe threat to the safety of overhead transmission lines. Compared with normal transmission lines, the ultra-high-voltage (UHV) transmission lines are more prone to galloping, and the damage caused is more severe. To control the galloping of UHV lines, it is necessary to conduct a comprehensive analysis of galloping characteristics. In this paper, a large-span 1000-kV UHV transmission line in China is taken as a practical example where an 8-bundled conductor with D-shaped icing is adopted. Galerkin method is employed for the time history calculation. For the wind attack angle range of 0°~180°, the galloping amplitudes in vertical, horizontal, and torsional directions are calculated. Furthermore, the vibration frequencies and galloping shapes are analyzed for the most severe conditions. The results show that the wind at 0°~10° attack angles can induce large torsional displacement, and this range of attack angles is also most likely to occur in reality. The galloping with largest amplitudes in all three directions occurs at the attack angle of 170° where the incoming flow is at the non-iced side, due to the strong aerodynamic instability. In addition, with wind speed increasing, galloping modes with higher frequencies appear and make the galloping shape more complex, indicating strong nonlinear behavior. Based on the galloping amplitudes of three directions, the full range of wind attack angles are divided into five galloping regions of different severity levels. The results obtained can promote the understanding of galloping and provide a reference for the anti-galloping design of UHV transmission lines.