• 대한기계학회논문집A
• /
• 제27권9호
• /
• pp.1579-1588
• /
• 2003

In this paper, the finite element equations for the transient linear viscoelastic stress analysis are presented in time domain, whose variational formulation is derived by using the Galerkin's method based on the equations of motion in time integral. Since the inertia terms are not included in the variational formulation, the time integration schemes such as the Newmark's method widely used in the classical dynamic analysis based on the equations of motion in time differential are not required in the development of that formulation, resulting in a computationally simple and stable numerical algorithm. The viscoelastic material is assumed to behave as a standard linear solid in shear and an elastic solid in dilatation. To show the validity of the presented method, two numerical examples are solved nuder plane strain and plane stress conditions and good results are obtained.

• 한국지반공학회논문집
• /
• 제22권9호
• /
• pp.45-59
• /
• 2006

탄소성 구성방정식은 주로 미분방적식(rate equation)으로 이루어져 있기 때문에 유한요소법 등을 이용한 지반구조물 해석시 미분방정식들에 대한 수치적분을 수행할 수 있는 방법이 필요하다. 구조물의 거동을 해석할시 미분방정식들을 위한 적분방법은 해석결과의 정확성과 유한요소법 모델링의 안전성에 큰 영향을 미치고 있다. 본 논문에서는 최근에 개발되어 사용되고 있는 흙에 관한 구성모델인 "Two-surface soil plasticity model (Manzari and Dafalias 1997)"을 Implicit return-mapping 수치적분방법을 이용하여 실행하는 과정을 제시한다. 본 연구에서 사용된 수치적분방법은 Closest-Point-Projection Method(CPPM) 방법으로 탄성 예측자-소성 교정자(elastic predictor-plastic corrector) 개념을 Implicit Backward Euler방법으로 체계화 시킨 알고리듬이다. 본 연구에서 수행한 "Two-surface soil plasticity model"은 조립토의 비선형거동을 해석하며, Bounding surface 개념 및 비선형 등방경화와 이동경화법칙을 사용하는 모델이다. 본 연구는 CPPM 방법이 정확하고 안정되며 유용한 수치적분을 수행할 수 있는 알고리듬이라는 것을 제시한다. 또한, CPPM 알고리듬은 구성방정식의 해를 반복적으로 해석하는 동안 "Consistent tangent operator $d{\sigma}/d{\varepsilon}$"를 제공하므로, 비선형 유한요소 해석이 2차(quadratic convergence rate)의 수렴 조건을 만족하는데 기여한다는 것을 보여준다.

• 대한수학회보
• /
• 제59권2호
• /
• pp.507-528
• /
• 2022

We propose a path integral method to construct a time stepwise local volatility for the stock index market under Dupire's model. Our method is focused on the pricing with the Monte Carlo Method (MCM). We solve the problem of randomness of MCM by applying numerical integration. We reconstruct this task as a matrix equation. Our method provides the analytic Jacobian and Hessian required by the nonlinear optimization solver, resulting in stable and fast calculations.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2001년도 춘계학술대회논문집
• /
• pp.457-464
• /
• 2001

The aeroacoustic phenomena in the simple rectangular open cavity are well published by many researchers. But the geometry shapes of aircraft landing gear wells, weapon bays, etc. are more complicate than that of the simple retangular cavity. They are more similar to the cavity having cover-plates at adges, or Helmholtz resonator. Therefore, the effects of cover-plates existing on edges of rectangular open cavity are numerically investigated in this paper. The compressible Navier-Stokes equations are solved for two-dimensional cavities with laminar boundary layers upstream. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Physically correct numerical boundary conditions and buffer zone techniques are implemented to produce time-accurate solutions in the whole computation domain. The computational domain is large enough to directly resolve a portion of the radiated acoustic field. Results show that the cover-plates existing on edges of cavity reduce the noise convected from cavity, make the frequency of noise become higher, and change the directivity pattern. So these results can be used in the design of a low noise cavity.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2000년도 춘계학술대회논문집
• /
• pp.104-109
• /
• 2000

The modes of oscillation and radiated acoustic fields of compressible flows over open cavities are investigated computationally. The compressible Navier-Stokes equations are solved for two-dimensional cavities with laminar boundary layers upstream. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Physically correct numerical boundary conditions are implemented to produce time-accurate solutions in the whole computation domain. The computational domain is large enough to directly resolve a portion of the radiated acoutic field. The results show a transition from a shear layer mode, for shorter cavities and lower Mach numbers, to a wake mode for longer cavities and higher Mach numbers. The shear layer mode is well characterized by Rossiter modes and these oscillations lead to intense upstream acoustic radiation dominated by a single frequency. The wake mode is characterized instead by a large-scale vortex shedding. Acoustic radiation is more intense, with multiple frequencies present.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2005년도 춘계 학술대회논문집
• /
• pp.123-127
• /
• 2005

The numerical analyses of the complicated flows are widely attempted in these days. Because of the enormous demanding memory and calculation time, parallel processing is used for these problems. In order to obtain calculation efficiency, it is important to choose proper domain decomposition technique and numerical algorithm. In this research we enhanced the efficiency of the CFDS code developed by ADD, using parallel computation and newly developed numerical algorithms. For the huge amount of data transfer between blocks non-blocking method is used, and newly developed data transfer algorithm is used for non-aligned block interface. Recently developed RoeM scheme is adpoted as a spatial difference method, and AF-ADI and LU-SGS methods are used as a time integration method to enhance the convergence of the code. Analyses of the flows around the ONERA M6 wing and the high angle of attack missile configuration are performed to show the efficiency improvement.

• 대한기계학회논문집
• /
• 제16권2호
• /
• pp.336-347
• /
• 1992

The three dimensional inviscid transonic cascade flow was investigated numerically, incorporation a four stage Runge-Kutta integration method proposed by Jameson. Time marching to the steady state was accelerated by using optimum time step and enthalpy damping. In describing the boundary conditions at inlet and outlet, Riemann invariants are considered. By adding a second and a fourth order artificial viscocities, the numerical instability due to the propagation of undamped disturbance or the rapid change of state near the shock has been prevented. The numerical results for are bump cascade, cambered two dimensional turbine cascade and three dimensional stator cascade agreed reasonably well with previous results. It has been known that the accuracy of the solution depended a lot on the modeling of the leading or trailing edge.

• 대한기계학회논문집
• /
• 제18권5호
• /
• pp.1134-1143
• /
• 1994

A study of adiabatic shear band formation and propagation of 4340 steel was done using the stepped speciment which was subjected to high velocity impact. The high velocity impact was performed on compression Hopkinson bar impact machine. After the controlled impact, the specimen was prepared for visual inspection. Numerical simulation was also performed with same geometrical dimension using explicit time integration finite element code. Experimental results were then compared with the numerical prediction. It was found that the numerical prediction is quite accurate, average thickness of adiabatic shear band is about $10{\mu}m$, the macro crack around shoulder is due to folding, and the deformation control ring is effective to freeze the propagation of adiabatic shear band.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2002년도 금형가공 심포지엄
• /
• pp.274-284
• /
• 2002

Deep drawing process for rectangular drawn section is different with that for axisymmetric circular one. Therefore deep drawing process for rectangular drawn section requires several intermediate steps to generate the final configuration without any significant defect. In this study, finite element analysis for multi-stage deep drawing process for high precision rectangular cases is carried out especially for an extreme aspect ratio. The analysis is performed using rigid-plastic finite element method with an explicit time integration scheme of the commercial program, LS-DYNA3D. The sheet blank is modeled using eight-node continuum brick elements. The results of analysis show that the irregular contact condition between blank and die affects the occurrence of failure, and the difference of aspect ratio in the drawn section leads to non-uniform metal flow, which may cause failure. A series of experiments for multi-stage deep drawing process for the rectangular cases are conducted, and the deformation configuration and the thickness distribution of the drawn rectangular cases are investigated by comparing with the results of the numerical analysis. The numerical analysis with an explicit time integration scheme shows good agreement with the experimental observation.