• 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.

• 제어로봇시스템학회:학술대회논문집
• /
• 제어로봇시스템학회 1993년도 한국자동제어학술회의논문집(국제학술편); Seoul National University, Seoul; 20-22 Oct. 1993
• /
• pp.507-512
• /
• 1993

In this paper we discuss the nonlinear motion of a conservative two-link arm using first integrals, which includes one integral constant. In the analysis of the motion, the constant plays important role. First, we give some discussion on the free motion by focusing on the integral constant. As the result, the free motion can be classified into two types-the one is oscillation and the other is rotation. Second, we discuss the forced motion of the arm actuated only at the second joint. We take the first integral in a more general form, and show that the forced motion of the second link can be expressed as a variation of the integral constant. Also, the characteristic of the forced motion actuated by arbitrary constant torques is discussed.

• 한국해양공학회지
• /
• 제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.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2002년도 춘계 학술발표회 논문집
• /
• pp.177-184
• /
• 2002

Large amplitude sloshing can occur in contained fluid region due to the seismic ground motion. Also, The pressure by large amplitude sloshing damages the connections between the wall and roof of a fluid container and causes outflow of contained fluid. Therefore, to predict the dynamic behavior accurately, three dimensional analysis with the nonlinear boundary condition must be performed. In this study, the numerical solution procedure is developed using the boundary element method with the Lagrangian particle approach. In order to demonstrate the accuracy and validity of the developed method, the fluid motion for a free oscillation with small amplitude and a forced vibration are analyzed. And the numerical results are compared with the linear theory results and the previous studies with the nonlinear boundary condition.

• 대한설비공학회:학술대회논문집
• /
• 대한설비공학회 2005년도 동계학술발표대회 논문집
• /
• pp.335-340
• /
• 2005

Nonlinear and forced oscillations of supported viscous droplet were focused in the present study. The droplet has a free contact line with solid plate and inviscid fluid. Natural frequencies of a pendant droplet have been investigated experimentally by imposing the acoustic wave while the frequency is being increased at a fixed amplitude. The evaporation was observed at atmosphere pressure. The droplet was recorded throughout the entire evaporation process and transient variations of the volume was measured. The evaporation process of oscillating droplet with thermofoil has been also observed to investigate analyzing the resonance effect on the thermal characteristics of droplet. It is found that a pendant droplet shows the resonant behaviors at each mode similar to the theoretical analysis. During imposing the acoustic wave, the pendant droplet makes a rotating motion in its longitudinal axis which is a new shape oscillation mode. The evaporation rate of a pendant droplet at resonant frequency is significantly enhanced.

• 대한조선학회논문집
• /
• 제29권4호
• /
• pp.114-131
• /
• 1992

본 논문은 2차원 자유표면파문제에서 시간영역해법을 이용하여 2차원 운동문제에 적용할 수 있는 수치해석을 하였다. 경계조건으로는 엄밀한 물체표면 경계조건과 비선형 자유표면경계 조건을 부과했다. 수치해를 구하는데는 코시이론을 이용하여 제2종 프레드흘름 경계적분방정식을 도출하고 이를 이산화시켜 처리하였다. 수치계산을 위해 전영역을 유한한 영역으로 제한하여야 한다. 제한된 영역에서 방사해의 부과를 위해 전영역을 수치해석영역과 외부영역으로 나누고, 외부영역의 해는 그린 제2정리를 이용하고, 선형자유표면조건을 만족하는 과도그린함수를 사용한다. 위의 그린 제2정리를 이용한 식으로 부터 초기조건, 선형 자유표면조건, 무한원방조건을 이용하여 단순화시킨 다음 포텐셜과 유동함수의 관계식으로 치환하면 비선형해와 정합할 수 있는 정합행렬을 구할 수 있다. 본 논문에서 개발한 정합방법을 이용하여 적용할 문제로서 첫째는 주상체가 상하동요, 수평동요하는 경우 계산이고 두번째로는 수면하에서 타원형실린더가 일정속도로 항진할 때 계산을 수행한 결과를 고차스펙트럴방법과 비교하였다.

• 대한조선학회지
• /
• 제24권2호
• /
• pp.47-54
• /
• 1987

Linear and nonlinear hydrodynamic force, which acts on submerged circular and eilliptic cylinders in oscillations as well as in advancing motion, are investigated as an initial-boundary value problem using a numerical method, which makes use of the source distribution on the body surface and the spectral method for treating the free surface waves. In the numerical code developed here, the boundary condition at the body surface is linearized. Using the numerical code so attained, nonlinear effects for different forward speeds and of the large-amplitude motion are computed. One of the major findings is that, when the forward speed is large, the added mass has its minimum and the damping force change rapidly around the frequency corresponding to the speed-frequency parameter, $\tau$=0.25, Compared to the result of Grue's [10], who used linear theory to get abrupt changes in values of the added mass and the damping force at the frequency corresponding to $\tau$=0.25, the present study, which takes nonlinear effects into account, shows much smoother variations near the frequency.

• 대한기계학회논문집B
• /
• 제35권6호
• /
• pp.617-624
• /
• 2011

본 연구에서는 사각용기 내에서의 자유표면을 갖는 2차원 슬로싱 문제에 대하여 고찰하였다. 용기 내부의 유동은 수평방향의 조화운동을 통해 형성되며, 외부 가진 속도는 u=Asin($2{\pi}ft$)으로 정의 된다. 여기서 u, A�� 그리고 f는 외부로부터 작용하는 가진 속도, 변위 그리고 주파수를 각각 나타낸다. 큰 변위(A~O(1)) 슬로싱 문제의 해석을 위한 실험설비를 구축하였으며, 광범위한 영역에서의 PIV실험을 수행하였다. 실험을 통해 자유표면의 움직임(motion)을 각각 서로 다른 물리적 특성을 갖는 세 가지; 선형 슬로싱의 자유표면의 움직임과 유사한 standing wave motion, 사각용기의 속도가 0을 지나는 순간(the moment of turn-over) 측벽에서 발생하는 run-up motion 그리고 측벽으로부터 내부유체로 점차적으로 전파되는 bore motion으로 분류하였다.

• 대한조선학회논문집
• /
• 제30권2호
• /
• pp.86-97
• /
• 1993

수치방법은 포텐셜 유동의 가정하에서 Semi-Lagrangian 기법을 사용하여 2차원 쇄기의 비선형운동과 축대칭 물체의 강제 상하동요 운동에 대해서 개발되었다. 2차원에서 Cauchy 이론은 경계를 따라서 복소포텐셜과 그것의 미분치를 계산하기 위해 적용되었고, 3차원에서 Rankinering 쏘오스가 사용되고 대수방정식을 풀기위해서 그린 제2정리를 이용하였다. 해는 완전한 사유표면 조건을 수치적분함으로서 시간전진시킨다. 수치계산 예는 정속도로 입수하는 쇄기형 주상체와 정지 상태로 부터 강제상하동요하는 문제를 택하였다. 쇄기입수 문제는 Chapman [4], Kim[11]의 계산결과와 비교된다. 위에서 적용된 기법을 이용하여 구한 시간영역에서 힘을 Fourier 변환함으로서 부가질량계수, 감쇄계수, 2차조화력등이 얻어지고 Yamashita[5]의 실험치와 비교된다.