• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.2
• /
• pp.114-120
• /
• 2004

Propeller blades have complex airfoil section type geometry and the thickness is continuously varied to both its length and cord-wise direction. in the present research, the finite element analysis program PROSTEC (Propeller Stress Evaluation Code) is developed to calculate the structural responses of propeller blades in irregular ship wake field. To represent the curved and skewed geometry of propeller blades accurately, 20-node curved solid element using the quadratic shape function is adopted. Input data for the analysis including the geometry and pressure distribution of propeller blades can be generated automatically from the propeller design program. And to visualize the results of analysis on windows system conveniently, the post processor PROSTEC-POST is developed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.3
• /
• pp.403-414
• /
• 1994

Nonlinear analysis of RC containment structure under thermal load and pressure is presented to trace the behaviour after an assumed LOCA. The temperature distribution varying with time through the wall thickness is determined by transient finite element analysis with the two time level scheme in time domain. The layered shell finite elements are used to represent the containment structures in nuclear power plants. Both geometric and material nonlinearities are taken into account in the finite element formulation. The constitutive relation of concrete is modeled according to Drucker-Prager yield criteria in compression. Tension stiffening model is used to represent the tensile behaviour of concrete including bond effect. The reinforcing bars are modeled by smeared layer at the location of reinforcements accounting elasto-plastic axial behaviors. The steel liner model under Von Mises yield criteria is adopted to represent elastic-perfect plastic behaviour. Geometric nonlinearity is formulated to consider the large displacement effect. Thermal stress components are determined by the initial strain concept during each time step. The temperature differential between any two consecutive time steps is considered as a load incremental. The numerical results from this study reveal that nonlinear temperature gradient based on transient thermal analysis will produces excessive large displacement. Nonlinear behavior of containment structures up to ultimate stage can be traced reallistically. The present study allows more realistic analysis of concrete containment structures in nuclear power plants.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Advances in nano research
• /
• v.10 no.2
• /
• pp.151-163
• /
• 2021

The main target of this study is to investigate nonlinear transient responses of moving polymer nano-size plates fortified by means of Graphene Platelets (GPLs) and resting on a Winkler-Pasternak foundation under a transverse pressure force and a temperature variation. Two graphene spreading forms dispersed through the plate thickness are studied, and the Halpin-Tsai micro-mechanics model is used to obtain the effective Young's modulus. Furthermore, the rule of mixture is employed to calculate the effective mass density and Poisson's ratio. In accordance with the first order shear deformation and von Karman theory for nonlinear systems, the kinematic equations are derived, and then nonlocal strain gradient scheme is used to reflect the effects of nonlocal and strain gradient parameters on small-size objects. Afterwards, a combined approach, kinetic dynamic relaxation method accompanied by Newmark technique, is hired for solving the time-varying equation sets, and Fortran program is developed to generate the numerical results. The accuracy of the current model is verified by comparative studies with available results in the literature. Finally, a parametric study is carried out to explore the effects of GPL's weight fractions and dispersion patterns, edge conditions, softening and hardening factors, the temperature change, the velocity of moving nanoplate and elastic foundation stiffness on the dynamic response of the structure. The result illustrates that the effects of nonlocality and strain gradient parameters are more remarkable in the higher magnitudes of the nanoplate speed.

• Smart Structures and Systems
• /
• v.30 no.3
• /
• pp.221-237
• /
• 2022

This paper is concerned with nonlinear modeling and efficient numerical simulation of electrostrictive materials and structures. Two types of such materials are considered: relaxor ferroelectric ceramics and electrostrictive polymers. For ceramics, a geometrically linear formulation is developed, whereas polymers are studied in a geometrically nonlinear regime. In the paper, we focus on constitutive modeling first. For the reversible constitutive response under consideration, we introduce the augmented Helmholtz free energy, which is composed of a purely elastic part, a dielectric part and an augmentation term. For the elastic part, we involve an additive decomposition of the strain tensor into an elastic strain and an electrostrictive eigenstrain, which depends on the polarization of the material. In the geometrically nonlinear case, a corresponding multiplicative decomposition of the deformation gradient tensor replaces the additive strain decomposition used in the geometrically linear formulation. For the dielectric part, we first introduce the internal energy, to which a Legendre transformation is applied to compute the free energy. The augmentation term accounts for the contribution from vacuum to the energy. In our formulation, the augmented free energy depends not only on the strain and the electric field, but also on the polarization and an internal polarization; the latter two are internal variables. With the constitutive framework established, a Finite Element implementation is briefly discussed. We use high-order elements for the discretization of the independent variables, which include also the internal variables and, in case the material is assumed incompressible, the hydrostatic pressure, which is introduced as a Lagrange multiplier. The elements are implemented in the open source code Netgen/NGSolve. Finally, example problems are solved for both, relaxor ferroelectric ceramics and electrostrictive polymers. We focus on thin plate-type structures to show the efficiency of the numerical scheme and its applicability to thin electrostrictive structures.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.2 no.3
• /
• pp.99-106
• /
• 1998

The performance of the turbofan engine, a medium scale civil aircraft which has been developing in Rep. of Korea, was analyzed and the control scheme for optimization the performance was studied. The dynamic and real-time linear simulation was performed in the previous study The result was that the fuel scedule of the step increase overshoot the limit temperature(3105 $^{\cire}R$) of the high pressure turbine and got small surge margine of the high pressure compressor. Therefore a control scheme such as the LQR(Linear Quadratic Regulator) was applied to optimizing the performance in this studies. The linear model was expected for designing controller and the real time linear model was developed to be closed to nonlinear simulation results. The system matrices were derived from sampling operating points in the scheduled range and then the least square method was applied to the interpolation between these sampling points, where each element of matrices was a function of the rotor speed. The control variables were the fuel flow and the low pressure compressor bleed air. The controlled linear model eliminated the inlet temperature overshoot of the high pressure turbine and obtained maximum surge margins within 0.55. The SFC was stabilized in the range of 0.355 to 0.43.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.169-172
• /
• 2002

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.1
• /
• pp.77-99
• /
• 2017

Offshore oscillating water columns (OWC) represent one of the most promising forms of wave energy converters. The hydrodynamic performance of such converters heavily depends on their interactions with ocean waves; therefore, understanding these interactions is essential. In this paper, a fully nonlinear 2D computational fluid dynamics (CFD) model based on RANS equations and VOF surface capturing scheme is implemented to carry out wave energy balance analyses for an offshore OWC. The numerical model is well validated against published physical measurements including; chamber differential air pressure, chamber water level oscillation and vertical velocity, overall wave energy extraction efficiency, reflected and transmitted waves, velocity and vorticity fields (PIV measurements). Following the successful validation work, an extensive campaign of numerical tests is performed to quantify the relevance of three design parameters, namely incoming wavelength, wave height and turbine damping to the device hydrodynamic performance and wave energy conversion process. All of the three investigated parameters show important effects on the wave-pneumatic energy conversion chain. In addition, the flow field around the chamber's front wall indicates areas of energy losses by stronger vortices generation than the rear wall.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2004.11a
• /
• pp.603-608
• /
• 2004

The paper analyzes a scheme of decontamination of radionuclides from concrete structures, in which rapid microwave heating is used to spall off a thin contaminated surface layer. The analysis is split in two parts: (1) The hygrothermal part of the problem, which consists in calculating the evolution of the temperature and pore pressure fields, and (2) the fracturing part, which consists in predicting the stresses, deformations and fracturing. The rate of the distributed source of heat due to microwaves in concrete is calculated on the basis of the standing wave normally incident to the concrete wall with averaging over both the time period and the wavelength because of the very short time period of microwaves compared to the period of temperature waves and the heterogeneity of concrete. The reinforcing bars parallel to the surface arc treated as a smeared steel layer. The microplane model M4 is used as the constitutive model for nonlinear deformation and distributed fracturing of concrete. The aim of this study is to determine the required microwave power and predict whether and when the contaminated surface layer of concrete spalls off. The effects of wall thickness, reinforcing bars, microwave frequencies and power are studied numerically. As a byproduct of this analysis, the mechanism of spalling of rapidly heated concrete is clarified.