• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2003년도 추계학술대회논문집
• /
• pp.486-491
• /
• 2003

The paper presents the relationship between the eigenvalue branches and the corresponding flutter modes of cantilevered pipes conveying fluid. The pipes are located on elastic foundations which can be regarded as a soil model. In this paper, elastic foundations are assumed as linear distributed translational springs. Governing equations of motion are derived by extended Hamilton's principle, and the numerical scheme using finite element method is applied to obtain the discretized equations. The critical How velocity and stability maps of the pipe are investigated according to the variation of elastic foundation parameters, mass ratios of the pipe and internal damping Parameter. Also, the vibrational modes associated with flutter are shown.

• Structural Engineering and Mechanics
• /
• 제54권3호
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• 대한기계학회논문집A
• /
• 제25권7호
• /
• pp.1139-1146
• /
• 2001

The reactor pressure vessel(RPV) is usually cladded with stainless steel to prevent corrosion and radiation embrittlement, and a number of subclad cracks have been found during an in-service-inspection. These subclad cracks should be assured for a safe operation under normal conditions and faulted conditions such as pressurized thermal shock(PTS). Currently available integrity assessment procedure for an RPV, ASME Code Sec. XI, are built on the basis of linear fracture mechanics (LEFM). In PTS condition, however, thermal stress and mechanical stress give rise to high tensile stress at the cladding and elastic-plastic behavior is expected in this area. Therfore, ASME Code Sec. XI is overly conservative in assessing the structural integrity under PTS condition. In this paper, the fracture parameter (stress intensity factor, K, and RT(sub)NDT) from elastic analysis using ASME Sec. XI and finite element method were validated against 3-D elastic-plastic finite element analyses. The difference between elastic and elastic-plastic analysis became significant with increasing crack depth. Therfore, it is recommended to perform elastic-plastic analysis for the accurate assessment of subclad cracks under TPS which causes plastic deformation at the cladding.

• Korea-Australia Rheology Journal
• /
• 제15권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.

• 한국도로학회논문집
• /
• 제17권3호
• /
• pp.77-83
• /
• 2015

PURPOSES : In this study, a fracture-based finite element (FE) model is proposed to evaluate the fracture behavior of fiber-reinforced asphalt (FRA) concrete under various interface conditions. METHODS : A fracture-based FE model was developed to simulate a double-edge notched tension (DENT) test. A cohesive zone model (CZM) and linear viscoelastic model were implemented to model the fracture behavior and viscous behavior of the FRA concrete, respectively. Three models were developed to characterize the behavior of interfacial bonding between the fiber reinforcement and surrounding materials. In the first model, the fracture property of the asphalt concrete was modified to study the effect of fiber reinforcement. In the second model, spring elements were used to simulated the fiber reinforcement. In the third method, bar and spring elements, based on a nonlinear bond-slip model, were used to simulate the fiber reinforcement and interfacial bonding conditions. The performance of the FRA in resisting crack development under various interfacial conditions was evaluated. RESULTS : The elastic modulus of the fibers was not sensitive to the behavior of the FRA in the DENT test before crack initiation. After crack development, the fracture resistance of the FRA was found to have enhanced considerably as the elastic modulus of the fibers increased from 450 MPa to 900 MPa. When the adhesion between the fibers and asphalt concrete was sufficiently high, the fiber reinforcement was effective. It means that the interfacial bonding conditions affect the fracture resistance of the FRA significantly. CONCLUSIONS : The bar/spring element models were more effective in representing the local behavior of the fibers and interfacial bonding than the fracture energy approach. The reinforcement effect is more significant after crack initiation, as the fibers can be pulled out sufficiently. Both the elastic modulus of the fiber reinforcement and the interfacial bonding were significant in controlling crack development in the FRA.

• Structural Engineering and Mechanics
• /
• 제7권6호
• /
• pp.533-550
• /
• 1999

Optimal design of raft-pile foundations is examined by combining finite element technique and the optimization approach. The piles and soil medium are modeled by three dimensional solid elements while the raft is modelled by shell elements. Drucker-Prager criterion is adopted for the soil medium while the raft and the piles are assumed to be linear elastic. For the optimization process, the approximate semi-analytical method is used for calculating constraint sensitivities and a constraint approximation method which is a combination of the extended Bi-point approximation and Lagrangian polynomial approximation is used for predicting the behaviour of the constraints. The objective function of the problem is the volume of materials of the foundation while the design variables are raft thickness, pile length and pile spacing. The generalized reduced gradient algorithm is chosen for solving the optimization process. It is demonstrated that the method proposed in this study is promising for obtaining optimal design of raft-pile foundations without carrying out a large number of analyses. The results are also compared with those obtained from the previous study in which linear analysis was carried out.

• 한국시뮬레이션학회:학술대회논문집
• /
• 한국시뮬레이션학회 2001년도 The Seoul International Simulation Conference
• /
• pp.373-379
• /
• 2001

There is a big issue to generate 3D hexahedral finite element (FE) model, since a process to divide the whole domain into several simple-shaped sub-domains is required before generating a continuous mesh with mapped mesh generators. In general, it is nearly impossible to set up proper division numbers interactively to keep mesh connectivity between sub-domains on a complicated arbitrary-shaped domain. If mesh continuity between sub-domains is not required in an analysis, this complicated process can be omitted. Element-free Galerkin method (EFGM) can accept discontinuous meshes, which only requires nodal information. However it is difficult to choose a reasonable influenced domain in moving least squares scheme with non-uniformly distributed nodes in discontinuous FE models. A new FE scheme fur discontinuous mesh is proposed in this paper by applying improved EFGM with some modification to derive FE approximated function in discontinuous parts. Its validity is evaluated on linear elastic problems.

• 동력기계공학회지
• /
• 제15권3호
• /
• pp.18-24
• /
• 2011

In this study, the reflection of the L(0,2), axially symmetric guided elastic wave from defects in pipes above ground is examined using finite element method. Phase and group velocity dispersion curves for the pipe were presented for the selection of the excitation mode. Some simple signal processing was applied to determine the amplitude of each of the reflected waves and to calculate the reflection coefficient. The results show the reflection coefficient of this mode is very close to a linear function of the circumferential extent of the defect. The motivation for the work was the development of a technique for inspecting chemical plant pipelines, but the study addresses the nature of the reflection function and its general applicability.

• Journal of Advanced Research in Ocean Engineering
• /
• 제4권1호
• /
• pp.25-34
• /
• 2018

In this study, the behavior of the coupled mooring system and floating body is analyzed. The related works are introduced for the mooring analysis of the floating body. Equations motion are introduced for calculating mooring force connected with the floating body. For formulating the equations of motion, the concept of the constrained force is applied for compact expression of it. The input and output data of the module for calculating mooring force is defined. The static analysis and quasi-static analysis are performed. For the analysis, equilibrium equation for elastic catenary mooring line is used by employing finite element method, and the C# solver is developed in this research. The analysis results are validated by comparing with other research results.

• 대한용접접합학회:학술대회논문집
• /
• 대한용접접합학회 2002년도 Proceedings of the International Welding/Joining Conference-Korea
• /
• pp.378-383
• /
• 2002

In previous study, the decrease and recovery of total stiffness in welded structure was discussed on the basis of experimental examination through tensile loading and unloading test of welded specimen. The recovery of structure stiffness was caused by the release of welding residual stress through mechanical loading. In this study, analysis model that is indispensable for the effective application of MSR method was established on the basis of test and measurement result. Thermal elasto-plastic analysis for welding process was performed by non-coupled analysis. Analysis results of welding process were transfer to elasto-plastic model for tensile loading & unloading by restart technique. In elasto-plastic analysis model for mechanical loading & unloading, hardening appearance of weld metal was considered by rezoning technique and tying technique was used for JIG condition of test machine.