• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.2
• /
• pp.45-56
• /
• 1994

A computationally efficient numerical method based on potential flow is developed for time-domain simulation of the nonlinear free-surface flows around a 2-dimensional hydrofoil. This numerical method, namely, spectral/boundary-element method, is a mixed one of the high-order spectral method and the boundary-element method in time-domain. The high-order spectral method is used to calculate the nonlinear evolution of free-surface, and the boundary-element method is used to calculate the effects of the hydrofoil and the shed vortex. As application examples, nonlinear free-surface flows around a 2-dimensional hydrofoil which starts from the rest and translates near the free-surface with or without harmonic oscillations are calculated. Nonlinear/unsteady results of free-surface waves and hydrodynamic farces are shown and discussed. Particularly, the results of steady-state which are obtained as a special case of the present unsteady solution are compared with others' calculated and experimental results, and good agreements are observed.

• Smart Structures and Systems
• /
• v.23 no.6
• /
• pp.589-613
• /
• 2019

Cables are prone to vibration due to their low inherent damping characteristics. Recently, negative stiffness dampers have gained attentions, because of their promising energy dissipation ability. The viscous inertial mass damper (termed as VIMD hereinafter) can be viewed as one realization of the inerter. It is formed by paralleling an inertial mass part with a common energy dissipation element (e.g., viscous element) and able to provide pseudo-negative stiffness properties to flexible systems such as cables. A previous study examined the potential of IMD to enhance the damping of stay cables. Because there are already models for common energy dissipation elements, the key to establish a general model for IMD is to propose an analytical model of the rotary mass component. In this paper, the characteristics of the rotary mass and the proposed analytical model have been evaluated by the numerical and experimental tests. First, a series of harmonic tests are conducted to show the performance and properties of the IMD only having the rotary mass. Then, the mechanism of nonlinearities is analyzed, and an analytical model is introduced and validated by comparing with the experimental data. Finally, a real-time hybrid simulation test is conducted with a physical IMD specimen and cable numerical substructure under distributed sinusoidal excitation. The results show that the chosen model of the rotary mass part can provide better estimation on the damper's performance, and it is better to use it to form a general analytical model of IMD. On the other hand, the simplified damper model is accurate for the preliminary simulation of the cable responses.

• Wind and Structures
• /
• v.27 no.6
• /
• pp.381-397
• /
• 2018

Nonlinear behavior in fluid-structure interaction (FSI) of bridge decks becomes increasingly significant for modern bridges with increasing spans, larger flexibility and new aerodynamic deck configurations. Better understanding of the nonlinear aeroelasticity of bridge decks and further development of reduced-order nonlinear models for the aeroelastic forces become necessary. In this paper, the amplitude-dependent and neutral angle dependent nonlinearities of the motion-induced loads are further highlighted by series of computational fluid dynamics (CFD) simulations. An effort has been made to investigate a semi-analytical time-domain model of the nonlinear motion induced loads on the deck, which enables nonlinear time domain simulations of the aeroelastic responses of the bridge deck. First, the computational schemes used here are validated through theoretically well-known cases. Then, static aerodynamic coefficients of the Great Belt East Bridge (GBEB) cross section are evaluated at various angles of attack, leading to the so-called nonlinear backbone curves. Flutter derivatives of the bridge are identified by CFD simulations using forced harmonic motion of the cross-section with various frequencies. By varying the amplitude of the forced motion, it is observed that the identified flutter derivatives are amplitude-dependent, especially for $A^*_2$ and $H^*_2$ parameters. Another nonlinear feature is observed from the change of hysteresis loop (between angle of attack and lift/moment) when the neutral angles of the cross-section are changed. Based on the CFD results, a semi-analytical time-domain model for describing the nonlinear motion-induced loads is proposed and calibrated. This model is based on accounting for the delay effect with respect to the nonlinear backbone curve and is established in the state-space form. Reasonable agreement between the results from the semi-analytical model and CFD demonstrates the potential application of the proposed model for nonlinear aeroelastic analysis of bridge decks.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.785-802
• /
• 2018

In this paper, two different computational methods, called Rayleigh-Ritz and collocation are developed to estimate the ultimate strength of composite plates. Progressive damage behavior of moderately thick composite laminated plates is studied under in-plane compressive load and uniform lateral pressure. The formulations of both methods are based on the concept of the principle of minimum potential energy. First order shear deformation theory and the assumption of large deflections are used to develop the equilibrium equations of laminated plates. Therefore, Newton-Raphson technique will be used to solve the obtained system of nonlinear algebraic equations. In Rayleigh-Ritz method, two degradation models called complete and region degradation models are used to estimate the degradation zone around the failure location. In the second method, a new energy based collocation technique is introduced in which the domain of the plate is discretized into the Legendre-Gauss-Lobatto points. In this new method, in addition to the two previous models, the new model named node degradation model will also be used in which the material properties of the area just around the failed node are reduced. To predict the failure location, Hashin failure criteria have been used and the corresponding material properties of the failed zone are reduced instantaneously. Approximation of the displacement fields is performed by suitable harmonic functions in the Rayleigh-Ritz method and by Legendre basis functions (LBFs) in the second method. Finally, the results will be calculated and discussions will be conducted on the methods.

• Steel and Composite Structures
• /
• v.30 no.1
• /
• pp.31-46
• /
• 2019

In the present study, the various dynamic properties of MWCNT embedded fiber reinforced polymer uniform and tapered composite (MWCNT-FRP) plates are investigated. Various configurations of a tapered composite plate with ply-drop off and uniform composite plate have been considered for the development of the finite element formulation and experimental investigations. First order shear deformation theory (FSDT) has been used to derive the kinetic and potential energy equations of the hybrid composite plates by including the effect of rotary inertia, shear deformation and non-uniformity in thickness of the plate. The governing equations of motion of FRP composite plates without and with MWCNT reinforcement are derived by considering a nine- node rectangular element with five degrees of freedom (DOF) at each node. The effectiveness of the developed finite element formulation has been demonstrated by comparing the natural frequencies and damping ratio of FRP composite plates without and with MWCNT reinforcement obtained experimentally. Various parametric studies are also performed to study the effect of CNT volume fraction and CNT aspect ratio of the composite plate on the natural frequencies of different configurations of CNT reinforced hybrid composite plates. Further the forced vibration analysis is performed to compare the dynamic response of the various configurations of MWCNT-GFRP composite plate with GFRP composite plate under harmonic excitations. It was observed that the fundamental natural frequency and damping ratio of the GFRP composite plate increase approximately 8% and 37% respectively with 0.5wt% reinforcement of MWCNT under CFCF boundary condition. The natural frequencies of MWCNT-GFRP hybrid composite plates tend to decrease with the increase of MWCNT volume fraction beyond 2% due to agglomeration of CNT's. It is also observed that the aspect ratio of the CNT has negligible effect on the improvement of dynamics properties due to randomly orientation of CNT's.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.399-413
• /
• 2020

The numerical simulation of liquid sloshing in the three-dimensional tanks under horizontal excitation and roll excitation was carried out, and the inhibition effect of different baffles on the sloshing phenomenon was investigated. The numerical calculations were carried out by the nonlinear Boundary Element Method (BEM) with Green's theorem based on the potential flow, which was conducted with the governing equation corresponding to the boundaries of each region. The validity of the method was verified by comparing with experimental values and published literatures. The horizontal baffle, the vertical baffle and the T-shaped baffle in the sloshing tanks were investigated respectively, and the baffles' position, dimension and the liquid depth were provided and discussed in detail. It is drawn that the baffle shape plays a non-negligible role in the tank sloshing. The vertical baffle is a more effective way to reduce the sloshing amplitude when the tank is under a horizontal harmonic excitation while the horizontal baffle is a more effective way when the tank is under a roll excitation. The amplitude of free surface elevation at right tank wall decreases with the increasing of the horizontal baffle length and the vertical baffle height. Although the T-shaped baffle has the best suppression effect on tank sloshing under horizontal excitation, it has limited suppression effect under roll excitation and will complicate the sloshing phenomenon when changing baffle height.

• Steel and Composite Structures
• /
• v.45 no.2
• /
• pp.175-191
• /
• 2022

A review of previous studies shows that although there is a considerable difference between buckling loads of structures under follower and non-follower lateral loads, only the buckling load of FGM elliptical cylindrical shell under non-follower lateral load was investigated in the literature. This study is the first to obtain the buckling load of elliptical FGM cylindrical shells under follower lateral load and also make a comparison between buckling loads of elliptical FGM cylindrical shells under follower and non-follower lateral loads. Moreover, this research is the first one to derive the load potential function of elliptical cylindrical shell. In this regard, the FGM cylindrical elliptical shell was modeled using the semi-analytical finite strip method and based on the First Shear Deformation Theory (FSDT). The shell is discretized by strip elements aligned in the longitudinal direction. The Lagrangian and harmonic shape functions were considered in the circumference and longitudinal directions, respectively. The buckling pressure of the shell under follower and non-follower lateral loads was obtained from eigenvalue problem. The results obtained from the model were compared with those presented in the literature to evaluate the validity of the model. A comparison index was defined to compare the buckling loads of the shell under follower and non-follower lateral load. A parametric study was carried out to investigate the effects of material properties and shell geometry characteristics on the comparison index. For the elliptical cylindrical shells with length-to-radius ratio greater than 16 and major-to-minor axis ratio greater than 0.6, the comparison index reaches to more than 20 percent which is significant. Moreover, the maximum difference is about 30 percent in some cases. The results obtained from the parametric study indicate that the buckling load of long elliptical cylindrical shell under non-follower load is not reliable.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.5
• /
• pp.1054-1065
• /
• 1993

A computer simulation is performed on the effectiveness of the active minimization of harmonically excited enclosed sound fields for producing global reduction in the amplitude of the pressure fluctuations. In this study for the appreciable reductions in total time averaged acoustic potential energy, $E_{pp}$, the transducer location strategies for three dimensional active noise control is presented based on a state space modal which approximates the closed acoustic field.In this study, the above theoretical basis is used to investigate the application of active control to sound fields of low modal density. By the used of room-like 3-dimensional rectangular enclosure it is demonstrated that the reductions in $E_{pp}$ can be achieved by using a single secondary source, provided that the source is placed within the half a wavelength from the primary source and placed away from nodal line of the sound field. Concerning the reductions in $E_{pp}$ by minimzing the pressure in sound fields by the use of 3-dimensional rectangular enclosure, the effects of the number of sensors and the locations of these sensors are investigated. When a few modes dominate the response it is found that if only a limited number of sensors are located away from nodal line and located at the pressure maxima of the sound field such as at each corner of a rectangular enclosure.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.12 no.3
• /
• pp.165-172
• /
• 2009

A set of equations for description of transformation of harmonic waves is proposed here. Velocity potential function and separation of variables are introduced for the derivation. The continuity equation is in a vertical plane is integrated through the water so that a horizontal one-dimensional wave equation is produced. The new equation composed of the complex velocity potential function, further be modified into. A set up of equations composed of the wave amplitude and wave phase gradient. The horizontally one-dimensional equations on the wave amplitude and wave phase gradient are the first and second-order ordinary differential equations. They are solved in a one-way marching manner starting from a side where boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient. Simple spatially-centered finite difference schemes are adopted for the present set of equations. The equations set is applied to three test cases, Booij's inclined plane slope profile, Massel's smooth bed profile, and Bragg's wavy bed profile. The present equations set is satisfactorily verified against existing theories including Massel's modified mild-slope equation, Berkhoff's mild-slope equation, and the full linear equation.

• Journal of the Korean Magnetics Society
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2017

Flywheel Energy Storage System (FESS) is composed by flywheel generating rotating potential energy and motor/generator set charging and discharging electric potential energy. The flywheel and motor/generator is connected by rotating shaft. And torque characteristics of motor/generator part can influence charging and mechanical traits of FESS. This paper analyze about motor/generator design method of 5 [kWh] FESS and torque ripple, harmonic effects by change of slots. At First, this paper proposes a method to estimate the flywheel size and the rotor size of the motor from the the rotational kinetic energy by inertia of FESS. The number of induction motor rotor slots for torque ripple reduction in the high speed operation region is selected. This paper performs to reduce the noise and vibration of the flywheel composed of coaxial with motor/generator and flywheel and realize the high efficiency.