• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.335-341
• /
• 2017

In this paper, beam finite elements with higher-order derivatives' continuity are formulated and evaluated for various boundary conditions. All the beam elements are based on Euler-Bernoulli beam theory. These higher-order beam elements are often required to analyze structures by using newly developed higher-order beam theories and/or non-classical beam theories based on nonlocal elasticity. It is however rare to assess the performance of such elements in terms of boundary and loading conditions. To this end, two higher-order beam elements are formulated, in which $C^2$ and $C^3$ continuities of the deflection are enforced, respectively. Three different boundary conditions are then applied to solve beam structures, such as cantilever, simply-support and clamped-hinge conditions. In addition to conventional Euler-Bernoulli beam boundary conditions, the effect of higher-order boundary conditions is investigated. Depending on the boundary conditions, the oscillatory behavior of deflections is observed. Especially the geometric boundary conditions are problematic, which trigger unstable solutions when higher-order deflections are prescribed. It is expected that the results obtained herein serve as a guideline for higher-order derivatives' continuous finite elements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.9
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.36 no.5
• /
• pp.355-369
• /
• 2018

LiDAR (Light Detection And Ranging) strip adjustment is process to improve geo-referencing of the ALS (Airborne Laser Scanner) strips that leads to seamless LiDAR data. Multiple strips are required to collect data over the large areas, thus the strips are overlapped in order to ensure data continuity. The LSA (LiDAR Strip Adjustment) consists of identifying corresponding features and minimizing discrepancies in the overlapping strips. The corresponding features are utilized as control features to estimate transformation parameters. This paper applied SURF (Speeded Up Robust Feature) to identify corresponding features. To improve determination of the corresponding feature, false matching points were removed by applying three schemes: (1) minimizing distance of the SURF feature vectors, (2) selecting reliable matching feature with high cross-correlation, and (3) reflecting geometric characteristics of the matching pattern. In the strip adjustment procedure, corresponding points having large residuals were removed iteratively that could achieve improvement of accuracy of the LSA eventually. Only a few iterations were required to reach reasonably high accuracy. The experiments with simulated and real data show that the proposed method is practical and effective to airborne LSA. At least 80 % accuracy improvement was achieved in terms of RMSE (Root Mean Square Error) after applying the proposed schemes.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.21-25
• /
• 2007

A general purpose program NUFLEX has been extended for two-phase flows with topologically complex interface and cavitation flows with liquid-vapor phase change caused by large pressure drop. In analysis of two-phase flow, the phase interfaces are tracked by employing a LS(Level Set) method. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. Also, the cavitation process is computed by including the effects of evaporation and condensation for bubble formation and collapse as well as turbulence in flows. The volume-faction and continuity equations are adapted for cavitation models with phase change. The LS and cavitation formulation are implemented into a general purpose program for 3-D flows and verified through several test problems.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.659-663
• /
• 2010

In recent, flood inundation damages by hydraulic structure failures have increased drastically and thus a variety of countermeasures were needed to minimize such damages. A real-time flood inundation prediction technique is essential to protect and mitigate flood inundation damages. In the context of real time flood inundation modeling, this study aims to develop a grid based two-dimensional numerical method for flood inundation modeling using globally-available DEM data: SRTM with $90m{\times}90m$ spatial resolution. The newly-developed model guarantees computational efficiency in terms of geometric data processing by direct application of DEM for flood inundation modeling and also have good compatibility with various types of raster data when compared to a commercial model such as FLUMEN. The model, which employed the leap-frog algorithm to solve shallow water and continuity equations, can simulate inundating flow from channel to lowland and also returning flow from lowland to channel by comparing water levels between channel and lowland in real time. We applied the model to simulate the BaekSan levee break in the Nam river during a flood period from August 10 to 13, 2002. The simulation results had good agreements with the field-surveyed data in terms of inundated area and also showed physically-acceptable velocity vector maps with respect to inundating and returning flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.32-40
• /
• 1994

This paper describes a smooth surface interpolating method of ship hull using a three-dimensional currie net that comes from the mesh curve fairing process. Geometric continuity(($GC^1$) is preserved across the boundary curve between patches. The three-dimensional curve net can have nonrectangular topologies, such as triangular and pentagonal topology. Among the boundary curve interpolation methods, Hermite blended Coons patch, Convex combination, and Gregory patch interpolation method are used to generate the ship hull surface. To check the fairness of the surface, the numerical method of surface/surface intersection problem is adopted. An application to an actual ship hull is given as an example.

• Advances in nano research
• /
• v.13 no.4
• /
• pp.351-368
• /
• 2022

Recently, promising structural technologies like multi-function, ultra-load bearing capacity and tailored structures have been put up for discussions. Finite Element (FE) modelling is probably the best-known option capable of treating these superior properties and multi-domain behavior structures. However, advanced materials such as Functionally Graded Material (FGM) and nanocomposites suffer from problems resulting from variable material properties, reinforcement aggregation and mesh generation. Motivated by these factors, this research proposes a unified shape function for FGM, nanocomposites, graded nanocomposites, in addition to traditional isotropic and orthotropic structural materials. It depends not only on element length but also on the beam's material properties and geometric characteristics. The systematic mathematical theory and FE formulations are based on the Timoshenko beam theory for beam structure. Furthermore, the introduced element achieves C1 degree of continuity. The model is proved to be convergent and free-off shear locking. Moreover, numerical results for static and free vibration analysis support the model accuracy and capabilities by validation with different references. The proposed technique overcomes the issue of continuous properties modelling of these promising materials without discarding older ones. Therefore, introduced benchmark improvements on the FE old concept could be extended to help the development of new software features to confront the rapid progress of structural materials.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.10
• /
• pp.55-65
• /
• 2014

Recent researches on debris flow is focused on understanding its movement mechanism and building a numerical simulator to predict its behavior. However, previous simulators emulating fluid-like debris flow have limitations in numerical stability, geometric modeling and application of various boundary conditions. In this study, depth integration is applied to continuity equation and force equilibrium for debris flow. Thickness of sediment, and average velocities in x and y flow direction are chosen for main variables in the analysis, which improve numerical stability in the area with zero thickness. Petrov-Galerkin formulation uses a discontinuous test function of the weighted matrix from DG scheme. Presented mechanical constitutive model combines fluid and granular behaviors for debris flow. Effects on slope angle, inducing debris height, and bottom friction resistance are investigated for a simple slope. Numerical results also show the effect of embankment at the bottom of the slope. Developed numerical simulator can assess various risk factors for the expected area of debris flow, and facilitate embankment design in order to minimize damage.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.291-309
• /
• 2023

This paper addresses the finite element modeling of functionally graded porous (FGP) beams for free vibration and buckling behaviour cases. The formulated finite element is based on simple and efficient higher order shear deformation theory. The key feature of this formulation is that it deals with Euler-Bernoulli beam theory with only three unknowns without requiring any shear correction factor. In fact, the presented two-noded beam element has three degrees of freedom per node, and the discrete model guarantees the interelement continuity by using both C⁰ and C¹ continuities for the displacement field and its first derivative shape functions, respectively. The weak form of the governing equations is obtained from the Hamilton principle of FGP beams to generate the elementary stiffness, geometric, and mass matrices. By deploying the isoparametric coordinate system, the derived elementary matrices are computed using the Gauss quadrature rule. To overcome the shear-locking phenomenon, the reduced integration technique is used for the shear strain energy. Furthermore, the effect of porosity distribution patterns on the free vibration and buckling behaviours of porous functionally graded beams in various parameters is investigated. The obtained results extend and improve those predicted previously by alternative existing theories, in which significant parameters such as material distribution, geometrical configuration, boundary conditions, and porosity distributions are considered and discussed in detailed numerical comparisons. Determining the impacts of these parameters on natural frequencies and critical buckling loads play an essential role in the manufacturing process of such materials and their related mechanical modeling in aerospace, nuclear, civil, and other structures.

• Bulletin of the Korean Mathematical Society
• /
• v.60 no.4
• /
• pp.1085-1100
• /
• 2023

For any given function f, we focus on the so-called prescribed mean curvature problem for the measure e^-f(|x|2)dx provided thate^-f(|x|2) ∈ L¹(ℝⁿ⁺¹). More precisely, we prove that there exists a smooth hypersurface M whose metric is ds² = dρ² + ρ²d𝜉² and whose mean curvature function is ${\frac{1}{n}}(\frac{u^p}{{\rho}^{\beta}})e^{f({\rho}^2)}{\psi}(\xi)$ for any given real constants p, β and functions f and ψ where u and ρ are the support function and radial function of M, respectively. Equivalently, we get the existence of a smooth solution to the following quasilinear equation on the unit sphere 𝕊ⁿ, $${\sum_{i,j}}({{\delta}_{ij}-{\frac{{\rho}_i{\rho}_j}{{\rho}^2+|{\nabla}{\rho}|^2}})(-{\rho}ji+{\frac{2}{{\rho}}}{\rho}j{\rho}i+{\rho}{\delta}_{ji})={\psi}{\frac{{\rho}^{2p+2-n-{\beta}}e^{f({\rho}^2)}}{({\rho}^2+|{\nabla}{\rho}|^2)^{\frac{p}{2}}}}$$ under some conditions. Our proof is based on the powerful method of continuity. In particular, if we take $f(t)={\frac{t}{2}}$, this may be prescribed mean curvature problem in Gauss measure space and it can be seen as an embedded result in Gauss measure space which will be needed in our forthcoming papers on the differential geometric analysis in Gauss measure space, such as Gauss-Bonnet-Chern theorem and its application on positive mass theorem and the Steiner-Weyl type formula, the Plateau problem and so on.