• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 한국소음진동공학회 2008년도 추계학술대회논문집
• /
• pp.71-81
• /
• 2008

This paper proposes a spectral element which can represent dynamic responses in high frequency domain such as Lamb waves on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by piezoelectric layer (PZT layer) bonded on a base plate. In the two layer beam model, a PZT layer is assumed to be rigidly bonded on a base beam. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with electro mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are formulated through equations of motions converted into frequency domain. A detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through comparison results with the conventional 2-D FEM and the previously developed spectral elements.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 한국해양공학회 2004년도 학술대회지
• /
• pp.204-209
• /
• 2004

This paper describes dynamic manipulability analysis of robotic arms moving in viscous fluid. The Manipulability is a functionality of manipulator system in a given configuration and under the limits of joint ability with respect to the tasks required to bt performed. To investigate the manipulability of underwater robotic arms, a modeling and analysis method are presented. The dynamic equation of motion of underwater manipulator is derived from the Lagrange - Euler equation considering with the hydraulic forces caused by added mass, buoyancy and hydraulic drag. The hydraulic drag term in the equation: is established as analytical form using Denavit - Hartenberg (D-H) link coordination of manipulator. Two analytical approaches based on Manipulability Ellipsoid are presented to visualize the manipulability of robotic arm moving in viscous fluid. The one is scaled ellipsoid which transforms the boundary of joint torque to acceleration boundary of end-effector by normalizing the torque in joint space while the other is shifted ellipsoid which depicts total acceleration boundary of end-effector by shifting the ellipsoid in work space. An analysis example of 2-link manipulator with proposed analysis scheme is presented to validate the method.

• Proceedings of the KSME Conference
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.588-594
• /
• 2001

본 연구를 통해 초음속 전투기 날개의 공력-구조를 동시에 고려한 다학제간 설계를 수행하였다. 공력해석을 위해 사용된 3 차원 Euler Code는 수렴 속도를 개선하기 위해 Multigrid를 적용하였으며, 3차의 transfinite interpolation을 사용하여 O-H type의 공력해석 격자계를 생성하였다. 구조 분야는 절점당 54개의 자유도를 가지는 9 절점 쉘 혼합 유한요소(9-node shell mixed finite element)를 사용하여 해석을 수행하였다. 설계변수는 공력쪽으로 날개의 평면형상에 관련된 변수 3개, 구조쪽은 날개 윗면과 아래면의 표피두께에 관련된 4개의 설계변수 사용하였으며, D-optimality 조건을 만족시키는 실험점들에 대해 공력해석과 구조해석이 연동된 정적 공탄성 해석을 수행한 후, 반응면 기법을 이용하여 목적함수와 제약조건에 대한 반응면을 구성하였다. 단일점 설계를 수행한 후 이를 바탕으로 3개의 설계점을 동시에 고려한 다점 설계를 수행하였으며, 공력만을 고려한 설계 결과와 공력-구조를 동시에 고려한 다학제간 설계결과의 비교를 통해 다학제간 설계의 타당성과 우수성을 입증하였다.

• Structural Monitoring and Maintenance
• /
• 제4권2호
• /
• pp.115-131
• /
• 2017

Classical flutter of wind turbine blades indicates a type of aeroelastic instability with fully attached boundary layer where a torsional blade mode couples to a flapwise bending mode, resulting in a mutual rapid growth of the amplitudes. In this paper the monitoring problem of onset of flutter is investigated from a detection point of view. The criterion is stated in terms of the exceeding of a defined envelope process of a specific maximum torsional vibration threshold. At a certain instant of time, a limited part of the previously measured torsional vibration signal at the tip of blade is decomposed through the Empirical Mode Decomposition (EMD) method, and the 1st Intrinsic Mode Function (IMF) is assumed to represent the response in the flutter mode. Next, an envelope time series of the indicated modal response is obtained in terms of a Hilbert transform. Finally, a flutter onset criterion is proposed, based on the indicated envelope process. The proposed online flutter monitoring method provided a practical and direct way to detect onset of flutter during operation. The algorithm has been illustrated by a 907-DOFs aeroelastic model for wind turbines, where the tower and the drive train is modelled by 7 DOFs, and each blade by means of 50 3-D Bernoulli-Euler beam elements.

• The Journal of the Acoustical Society of Korea
• /
• 제25권1E호
• /
• pp.32-35
• /
• 2006

The screech tone of underexpanded jet is numerically calculated without any specific modeling for the screech tone itself. Fourth-order optimized compact scheme and fourth-order Runge-Kutta method are used to solve the 2D axisymmetric Euler equation. Adaptive nonlinear artificial dissipation model and generalized characteristic boundary condition are also used. The screech tone, generated by a closed loop between instability waves and quasi-periodic shock cells at the near field, is reasonably analyzed with present numerical methods for the underexpanded jet having Mach number 1.13. First of all, the centerline mean pressure distribution is calculated and compared with experimental and other numerical results. The instantaneous density contour plot shows Mach waves due to mixing layer convecting supersonically, which propagate downstream. The pressure signal and its Fourier transform at upstream and downstream shows the directivity pattern of screech tone very clearly. Most of all, we can simulate the axisymmetric mode change of screech tone very precisely with present method. It can be concluded that the basic phenomenon of screech tone including the frequency can be calculated by using high-order and high-resolution schemes without any specific numerical modeling for screech tone feedback loop.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.197-202
• /
• 1995

비정상(unsteady) 압축성(compressible) 유동에 의한 공력음향(aeroacoustics)을 모사하여 공력소음원을 해석하기 위해서는 고차(high order)의 정확도와 높은 해상도(resolution)를 가지며, 상대적으로 계산시간을 많이 필요로 하지 않는 외재적(explicit) 유한차분법이 필수적으로 요구된다. 이것은 주어진 차분방식과 격자계로써 공간과 시간상에 존재하는 미소크기의 파동성분들을 충분히 구현하여야 만족할 만한 수치해를 얻을 수 있기 때문이다. 본 연구에서는, 그러한 유한차분법 중 최근에 관심의 대상이 되고있는 삼각(tridiagonal)또는 오각(pentadiagonal) 집적유한차분법(compact finite difference scheme)이 최대의 해상도를 갖도록 하는 수학적인 방법을 개발하고, 이 방법으로써 새롭게 집적유한차분법을 최적화하였다. 개발된 최적화 방법은, 푸리에 해석법(Fourier analysis)을 통하여 파동수(wavenumber) 영역에서 수학적으로 계산된 위상오차(phase error)를 최소화하는 것이며, 이러한 개념과 방법은 본 연구에서 처음으로 집적유한차분법에 적용되었다. 여러가지 절단정확도(truncation order)에 대해서 최적화 된 집적유한차분법들이 실제 공간과 시간상에서 보여주는 정확도와 오차특성을 알아보기 위하여, 이 방법들을 1차원 선형파동방정식에 적용하였고, 이 결과를 통하여 가장 정확하고 효과적인 절단정확도의 집적유한차분법을 선별하였다. 특히, 오각(pentadiagonal)법에 비해 더욱 효율적인 6차 삼각(tridiagonal)법을 1차원 Euler방정식에 적용하여, 비선형 파동에 대한 모사를 수행할 수 있었다.

• Proceedings of the KSME Conference
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• pp.387-394
• /
• 2000

Theoretical and experimental study is conducted on the contact pressure distribution of a radial tire in motion under various camber angles. Tire construction is modelled by a spring bedded elastic ring, consisted of sidewall springs and a composite belt ring. The contact area is assumed to be a trapezoidal shape varying with camber angles and weighted load. The basic equation in a quasi-static form is derived for the deformation of a running belt with a constant velocity by the aid of Lagrange-Euler transformation. Galerkin's method and stepwise calculation are applied for solving the basic equation and the mechanical boundary condition along both sides of the contact belt part subjected to shearing forces transmitted from the sidewall spring. Experimental results on the contact pressure, measured by pressure sensors embedded in the surface of the drum tester, correspond well with the calculated ones for the test tire under various camber angles, running velocities and weighted loads. These results indicate that a buckling phenomenon of the contact belt in the widthwise direction occurs due to the effect of camber angle.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1998년도 강연회 및 연구개발 발표회 논문집
• /
• pp.128-135
• /
• 1998

Procedures and results of aerodynamic design of a centrifugal compressor are presented for development of a 40kW class turbogenerator gas turbine. Specification of higher level of total pressure ratio of 4 and total efficiency of $80\%$ requires advanced methods of design and analysis. In the meanline design/analysis, a method with conventional loss modeling and a method with the two-zone model are alternately used for more reliable prediction. In the impeller blade generation, a series of Bezier curve are combined to produce meridional contours and distributions of blade camber angle and blade thickness. Intermediate profiles of blades are repeatedly produced and changed to be finally fixed through quasi-three dimensional Euler flow analysis. Three dimensional compressible turbulent flow analysis is then performed for the impeller to be confirmed in the final step of design. Satisfactory results in the aerodynamic performance are obtained, which assures that there is no need of aerodynamic re-design.

• Annual Conference of KIPS
• /
• 한국정보처리학회 2016년도 춘계학술발표대회
• /
• pp.189-191
• /
• 2016

실시간 물리 기반 3D 시뮬레이션에서 연산속도는 매우 중요한 요소이다. 객체의 움직임이나 변형과 같은 현상들은 복잡한 연산을 통해서 계산되기 때문에 일반적으로 시뮬레이션의 정확도와 연산속도는 반비례 관계에 있다. 현재 출시되고 있는 대부분의 게임에서는 물체의 움직임을 정확하게 표현하기보다 연산량을 줄이기 위해 물체의 움직임이나 변형을 비슷하게 표현하는데 중점을 두고 있다. 본 논문에서는 이러한 문제를 해결하기 위하여 OpenGL 4.3의 Compute shader를 사용하여 다이내믹 시뮬레이션의 연산 작업을 GPU 병렬처리로 처리하였다. Compute shader에서 파티클의 움직임을 계산하고 Shader storage buffer object에 저장하고 파티클들의 작업량을 적절한 Workgroup의 크기로 나누어 할당하여 최적의 처리속도를 제공하도록 구현하였다. Compute shader에서 파티클의 움직임을 표현하기 위해서 수치해법 중의 하나인 Euler method를 사용하였으며 실험 결과 파티클의 수가 4,194,304개일 때 CPU 방법에 비해 약 182배 빠른 연산속도 결과를 보였다. 추후 Compute shader를 활용하여 연산량이 많은 분야에 적용 가능할 수 있을 것으로 기대한다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• 제14권5호
• /
• pp.1244-1253
• /
• 1990

Time-Marching methods solving Euler equations are used for calculation of two-dimensional, steady, inviscid flow through a stationary compressor cascade. Calculation method is based on the Denton`s opposed difference scheme. A smoothing in the axial direction is used to increase the stability of solution. The computational grid consists of quadrilateral elements, one of which has four nodes at each corner and the grid points on the upper periodic boundaries are located one pitch away from those on the lower boundaries to satisfy the periodicity condition. Results of calculation show good agreement with other computational and experimental results, proving that the present method of calculation should be applied with confidence for the cascade flow with shock wave.