• Structural Engineering and Mechanics
• /
• 제50권4호
• /
• pp.487-500
• /
• 2014

An analytical solution for the nonlinear in-plane free oscillations of the suspended cable which contains the quadratic and cubic nonlinearities is investigated via the homotopy analysis method (HAM). Different from the existing analytical technique, the HAM is indeed independent of the small parameter assumption in the nonlinear vibration equation. The nonlinear equation is established by using the extended Hamilton's principle, which takes into account the effects of the geometric nonlinearity and quasi-static stretching. A non-zero equilibrium position term is introduced due to the quadratic nonlinearity in order to guarantee the rule of the solution expression. Therefore, the mth-order analytic solutions of the corresponding equation are explicitly obtained via the HAM. Numerical results show that the approximate solutions obtained by using the HAM are in good agreement with the numerical integrations (i.e., Runge-Kutta method). Moreover, the HAM provides a simple way to adjust and control the convergent regions of the series solutions by means of an auxiliary parameter. Finally, the effects of initial conditions on the linear and nonlinear frequency ratio are investigated.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2005년도 춘계학술대회 논문집
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.

• Nuclear Engineering and Technology
• /
• 제49권5호
• /
• pp.953-967
• /
• 2017

In this paper, a complete station blackout (SBO) or complete loss of electrical power supplies is simulated and analyzed in a downward cooling 5-MW pool-type Material Testing Reactor (MTR). The scenario is traced in the absence of active cooling systems and operators. The code nodalization is successfully benchmarked against experimental data of the reactor's operating parameters. The passive heat removal system includes downward water cooling after pump breakdown by the force of gravity (where the coolant streams down to the unfilled portion of the holdup tank), safety flapper opening, flow reversal from a downward to an upward cooling direction, and then the upward free convection heat removal throughout the flapper safety valve, lower plenum, and fuel assemblies. Both short-term and long-term natural core cooling conditions are simulated and investigated using the RELAP5 code. Short-term analyses focus on the safety flapper valve operation and flow reversal mode. Long-term analyses include simulation of both complete SBO and long-term operation of the free convection mode. Results are promising for pool-type MTRs because this allows operators to investigate RELAP code abilities for MTR thermal-hydraulic simulations without any oscillation; moreover, the Tehran Research Reactor is conservatively safe against the complete SBO and long-term free convection operation.

• 한국산학기술학회논문지
• /
• 제18권2호
• /
• pp.71-76
• /
• 2017

선박 및 해양구조물 설계에 있어 롤댐핑에 관한 문제는 유체점성과 현상자체의 비선형성으로 인해 공학자들에게 있어 난제로 남아있다. 본 연구에서는 강제동요방법을 이용하여 원형실린더의 점성 롤감쇠에 관한 연구를 수행하였다. 토크 센서를 이용하여 강제동요 시 발생하는 롤모멘트(roll moment)를 강제동요 주기별로 계측하였고 이를 실험식(empirical formula)과 비교 검토 하였다. 점성에 의한 전단력으로부터 계측된 토크의 크기가 상대적으로 작은 값임에도 불구하고 실험식으로부터 얻어진 계산값들과 정성적으로 유사한 결과를 보였고, 일부 주기에서는 정량적으로도 잘 일치하는 결과를 보였다. 또한, PIV 계측기법을 통해 원형실린더 벽면주위의 유동을 면밀히 관찰하였으며, 유체점성으로 인해 원형실린더 벽면근처에서 경계층이 형성되고 자유수면에서 원형실린더의 주기적인 강제회전동요로 인해 미소한 크기의 파가 생성(wave making)됨을 PIV 계측결과의 분석을 통해 확인하였다. 본 연구에서는 점성 롤감쇠로 인한 실험식의 적합성을 모형시험을 통해 확인하였고, PIV 계측기법을 통해 벽면주위의 유체점성으로 인한 조파현상을 입증하였다.

• 대한조선학회지
• /
• 제14권3호
• /
• pp.13-22
• /
• 1977

The hydrodynamic forces acting on a forced oscillating 2-dimensional cylinder on a free surface of a fluid of a finite depth are calculated by distributing singularities on the immersed body surface. And the Haskind-Newman relation in a fluid of a finite depth is derived. The wave exciting force of the cylinder to an oscillation is also calculated by using the above relation. The method is applied to a circular cylinder swaying in a water of finite depth, and then, to a rectangular cylinder heaving, swaying, and rolling. The results of above cases give a good agreement with those by earlier investigators such as Bai, Keil, and Yeung. Also, this method is applied to a Lewis form cylinder with a half beam-to-draft ratio of 1.0 and a sectional area coefficient of 0.941, and to a bulbous section cylinder which is hard to represent by a mapping function. The results reveal that the hydrodynamic forces in heave increase as the depth of a water decrease, but in sway or roll, the tendency of the hydrodynamic forces is difficult to say in a few words. The exciting force to heave for a bulbous section cylinder becomes zero at two frequencies. The added mass moment of inertia for roll is seemed to mainly depend on the sectional shape than the water depth.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• 한국연소학회지
• /
• 제16권3호
• /
• pp.41-51
• /
• 2011

In this paper, we describe the behavior of two self-excitations in laminar attached free-jet flames under the influence of DC electric fields, one of buoyancy-driven and the other of diffusion-thermal instability, established from the horizontal and vertical injection. In the horizontal injection with removed buoyancy effect, oscillating flames with the frequency of 1.3 - 7.4 Hz were observed in a certain condition with Lewis number more than unity. On the other hand, it was appeared Lewis number induced self-excitation as well as buoyancy-driven self-excitation in the vertical upward injection with DC electric fields. This behavior had frequency range of 1.6 - 9.4 Hz and was exhibited to attribute the buoyancy effect. Finally, a well-defined division about two self-excitations having similar frequency range is briefly discussed.

• International Journal of Fluid Machinery and Systems
• /
• 제2권3호
• /
• pp.197-205
• /
• 2009

Since large-scale commercial wind turbine generator systems such as MW-class wind turbines are becoming widely operated, the vibration and distortion of the blade are becoming larger and larger. Therefore the soft structure design instead of the solid-design is one of the important concepts to reduce the structural load and the cost of the wind turbine rotors. The objectives of the study are development of the fluid-structure coupled analysis code and evaluation of soft rotor-blade design to reduce the unsteady structural blade load. In this paper, fluid-structure coupled analysis for the HAWT rotor blade is performed by free wake panel method coupled with hinge-spring blade model for the flapwise blade motion. In the model, the continuous deflection of the rotor blade is represented by flapping angle of the hinge with one degree of freedom. The calculation results are evaluated by comparison with the database of the NREL unsteady aerodynamic experiment. In the analysis the unsteady flapwise moments in yawed inflow conditions are compared for the blades with different flapwise eigen frequencies.

• 한국해양공학회지
• /
• 제20권6호
• /
• pp.67-74
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with large amplitude under the free~surface are solved in time-domain. Through the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.