• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.5
• /
• pp.275-282
• /
• 2018

Using isogeometric analysis(IGA) gives more accurate results for higher order mode in eigenvalue problem than using the finite element method(FEM). This is because the FEM has $C^0$ continuity between elements, whereas IGA guarantee $C^{P-1}$ between elements for p-th order basis functions. In this paper, a mode based reduced model is constructed by using IGA and dynamic behavior analysis is performed using this advantage. Craig-Bampton(CB) method is applied to construct the reduced model. Several numerical examples were performed to compare the eigenvalue analysis results for various order of element basis function by applying the IGA and FEM to simple rod analysis. We have confirmed that numerical error increases in the higher order mode as the continuity between elements decreases in the IGA by allowing internal knots multiplicity. The accuracy of the solution can be improved by using the IGA with high inter-element continuity when high-frequency external force acts on the reduced model for dynamic behavior analysis.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.04a
• /
• pp.344-349
• /
• 1995

셀구조물은 선박, 항공기, 우주선 등의 산업분야에 자주 사용되는 중요한 구조요소로서, 이에 대한 연구도 공학적 필요성에 의해 많은 연구가 진행되어 왔다. Maddox등의 Rayleigh-Ritz법, Srinivasan and Bobby의 적분방정식 기법, Petyt의 유한요소법, Irie, Henderson등의 전달매트릭스법, 그 외 Webster, Petyt and Deb Nath, Blevins, Srinivasan and Krishnan등에 의해 여러가지 해석방법의 진동해석이 행해져 왔다. 현재 주로 사용되고 있는 매트릭스 구조해석 및 진동해석 방법으로는 유한요소법과 전달매트릭스법을 들 수 있다. 유한요소법은 범용적인 해석 프로그램의 개발은 수월하지만 대규모의 선형 연립방정식의 해법에 귀착되므로, 기억용량이 큰 대형컴퓨터가 필요하고, 반면 전달매트릭스법은 기억용량이 적은 퍼스널컴퓨터로도 계산수행이 가능하나 고차의 고유진동수를 구할 경우나, 중간에 단단한 경탄성지지가 존재할 경우 등에는 수치계산상의 문제점이 지적되고 있다. 이에 대한 대책으로, 상태변수를 바꾸어 넣은 Riccati방법, 각 절점의 기지의 상태변수를 제거하는 Frontal법, 전달매트릭스와 강성매트릭스를 결합시키는 방법 등이 보고 되고 있다. 이에 자자들은 퍼스널컴퓨터 이용에 적합한 고속, 고정도의 구조해석 및 진동해석 기법으로 전달영향계수법을 제안하여, 원판구조물, 사각판구조물, 원통형 셀구조물 등의 여러가지 형상구조물의 자유진동 해석에 적용해서, 종래의 전달매트릭스법에 비해 계산정도 및 계산속도의 양면에서 매우 우수함을 보고한바 있다. 본 연구에서는 동적영향계수의 축차전달에 그 개념을 두고 있는 전달영향계수법을 좀 더 일반화시키고 체계화시키기 위하여 탄성지지를 갖는 셀구조물에 적용하여 자유진동 해석 알고리즘을 정식화한 후, 간단한 모델에 대한 수치실험을 통해서 전달 영향계수법으로 구한 해를 전달매트릭스법의 결과와 비교.검토하여 본 방법의 유용성을 확인하였다.

• Journal of the Korea Society for Simulation
• /
• v.30 no.1
• /
• pp.11-20
• /
• 2021

As unmanned underwater systems have recently emerged, a high-speed underwater channel modeling technique, which is one of the most important techniques in the system, has received a lot of attention. In this paper, we proposed a high-speed sound propagation model and verified the applicability through quantitative performance analyses. We used a high-order finite difference method (FDM) for wave propagation modeling in the water, and a domain decomposition method was adopted using multiple general-purpose graphics processing units (GPUs) to increase the calculation efficiency. We compared the results of the model we proposed with the analytic solution in the half-infinite media and results of the Virtual Timeseries Experiment (VirTEX) model, which is based on the ray method. Finally, we analyzed the performance of the model quantitatively using numerical examples. Through quantitative analyses of the improvement in computational performance, we confirmed that the computational speed increases linearly as the number of GPUs increases. The computation times are increased by 2 times and 8 times, respectively, when the domain size of computation and the maximum frequency are doubled. We expect that the proposed high-speed underwater channel modeling technique is able to contribute to the enhancement of national defense as an underwater communication channel model and analysis tool to develop the underwater communication technique for the unmanned underwater system.

• Journal of Internet Computing and Services
• /
• v.24 no.2
• /
• pp.27-36
• /
• 2023

Waveguides are transmission lines that guide electromagnetic waves in the desired direction and are utilized in various fields such as medical devices, radar systems, and satellite communications. Computational electromagnetics (CEM) is essential for designing and optimizing waveguides. The finite element method (FEM), which is one of the numerical analysis techniques, is efficient in solving closed problems such as waveguides. In order to apply FEM for waveguide analysis, boundary conditions that truncate the computational domain are required. This paper performs electromagnetic simulations using absorbing boundary conditions (ABC) and waveguide port boundary conditions (WPBC) in 2/D and 3/D waveguides using the finite element method and compared their performances. The accuracy of the analysis was verified by comparing the results with HFSS, a representative commercial electromagnetic simulation software. Simulation results confirm that applying WPBC allows for smaller analysis domains than ABC.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.11a
• /
• pp.1494-1498
• /
• 2007

Impulsive shooting noise is basically complex phenomenon which contains the linear and non-linear characteristics. For those reasons, numerical analysis of impulsive shooting noise has the difficulties in control of the numerical stability and accuracy on the simulation. In this research, Wave-number Extended Finite Volume Scheme (WEFVS) is applied to the numerical analysis of impulsive shooting noise. In the muzzle blast flow simulation, the generation of the precursor wave and the induced vortex ring are observed. Consequently, blast wave. vortex ring interaction and vortex ring. bow shock wave interaction are evaluated on the shooting process using the accurate and stable scheme. The sound generation in the interactions can be explained by the vorticity transport theorem. The shear layer is evolved behind the projectiles due to the jet flow. In these computations, the impulsive shooting noise is generated by the complex interaction with shooting process and is propagated to the far-field boundary. The impulsive shooting noise generation can be observed by the applications of WEFVS and analyzed by the physical phenomena.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.76-86
• /
• 1998

Structural intensity of plates experiencing bending vibration is analytically evaluated using the modal analysis based on assumed mode method. To evaluate the convergence of structural intensity according to the number of superposition modes, the power obtained by structural intensity integration over the closed curve containing the excitation source is compared with the power injected into plates. The erect of power reduction due to the material internal loss is evaluated using the intensity around a localized damping point, In addition, the dominant component among internal forces in the power transfer by the bending vibration of plates and the change of power flow due to stiffener are also investigated.

• Korean Journal of Optics and Photonics
• /
• v.16 no.5
• /
• pp.423-427
• /
• 2005

In this paper, we propose easily tooling method for Seidel third order aberration, which are not well utilized in actual design process due to the complication of mathematical operation and the difficulty of understanding Seidel third order aberration theory, even though most insightful and systematic means in pre-designing for the initial data of optimization. First, using paraxial ray tracing and Seidel third order aberration theory, spherical aberration coefficient is derived for a two-mirror system with a finite object distance. The coefficient, which is expressed as a higher-order nonlinear equation, consists of design parameters(object distance, two curvatures, and inter-mirror distance) and effective focal length(EFL). Then, the expressed analytical equation is solved by using a computer with numerical analysis method. From the obtained numerical solutions satisfying the nearly zero coefficient condition($<10^{-6}$), linear fitting process offers a linear relationship called the curvature linear equation between two mirrors. Consequently, this linear equation has two worthy meanings: the equation gives a possibility to obtain initial design data for optimization easily. And the equation shows linear relationship to a two-mirror system with a finite object distance under the condition of corrected third order spherical aberration.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.7
• /
• pp.18-25
• /
• 2005

In this study, the main focus is the simulation of acoustic wave scattering in trailing edge and the analysis of the generation mechanism of instability wave by the interaction of trailing edge, shear flow and initial disturbance. The numerical algorithm is based on CFD/CAA hybrid method with high-order computational aeroacoustic method. It is found that steady mean flow gradient terms play a crucial role on the generation of instability wave through the comparison of simulations of Simple Linearized Euler Equation and Full Linearized Euler Equation. Through the comparison with the results of Full Navier-Stokes Equation, it is reasonable and efficient to use the Full Linearized Euler Equation in the initial generation mechanism of the instability wave near the trailing edge.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.37-43
• /
• 1996

Breaking waves were generated in a 2-D flume. A piston-type wavemaker was operated in accordance with signals which consist of elementary harmonics with appropriate phase differences. These phase differences were estimated by using a linear wave theory so that wave crests were to be concentrated at the same position. The stroke of wavemaker was controlled to create plunging-type breaking waves. The signal with small amplitude could not generate breaking waves. In the case of moderate amplitudes, various breaking waves could be obtained. Most of breaking waves were spilling type. Only when the wavemaker was operated with appropriate amplitude, plunging-type breaking waves were generated. The parameters of breaking waves are the wave steepness and the frequency bandwidth. If the central frequency was low, breaking waves were not generated. Based on experimental data, we found that the wave height of breaking inception was H = 0.0113 gT$^2$. We also made computations by using a mixed Euler-Lagrangian scheme under the assumption of potential flow. The numerical results show good agreements with tank measurements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.3
• /
• pp.187-195
• /
• 2024

This paper presents two-dimensional boundary value problems of the stress-based gradient elasticity within the smoothed finite element method (S-FEM) framework. Gradient elasticity is introduced to address the limitations of classical elasticity, particularly its struggle to capture size-dependent mechanical behavior at the micro/nano scale. The Ru-Aifantis theorem is employed to overcome the challenges of high-order differential equations in gradient elasticity. This theorem effectively splits the original equation into two solvable second-order differential equations, enabling its incorporation into the S-FEM framework. The present method utilizes a staggered scheme to solve the boundary value problems. This approach efficiently separates the calculation of the local displacement field (obtained over each smoothing domain) from the non-local stress field (computed element-wise). A series of numerical tests are conducted to investigate the influence of the internal length scale, a key parameter in gradient elasticity. The results demonstrate the effectiveness of the proposed approach in smoothing stress concentrations typically observed at crack tips and dislocation lines.