• Geomechanics and Engineering
• /
• v.16 no.6
• /
• pp.619-626
• /
• 2018

A simple prediction procedure was investigated for calculating the stresses and displacements of a circular opening. Unlike existed approaches, the proposed approach starts each step with a radius increment. The stress for each annulus could be obtained analytically, while strain increments for each step can be determinate numerically from the compatility equation by finite difference approximation, flow rule and Hooke's law. In the successive manner, the distributions of stresses and displacements could be found. It should be noted that the finial radial stress and displacement were equal to the internal supporting pressure and deformation at the tunnel wall, respectively. By assuming different plastic radii, GRC and the evolution curve of plastic radii and internal supporting pressures could be obtained conveniently. Then the real plastic radius can be calculated by using linear interpolation in the evolution curve. Some numerical and engineering examples were performed to demonstrate the accuracy and validity for the proposed procedure. The comparisons results show that the proposed procedure was faster than that in Lee and Pietrucszczak (2008). The influence of annulus number and dilation on the accuracy of solutions was also investigated. Results show that the larger the annulus number was, the more accurate the solutions were. Solutions in Park et al. (2008) were significantly influenced by dilation.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Steel and Composite Structures
• /
• v.44 no.5
• /
• pp.633-649
• /
• 2022

This paper is first implemented with the bending behavior of three-dimensional functionally graded (3DFG) plates in the framework of level set-based topology optimization (LS-based TO). Besides, due to the suitable properties of the current design domain, the thin-plate spline (TPS) is recognized as a RBF to construct the LS function. The overall mechanical properties of the 3DFG plate are assessed using a power-law distribution scheme via Mori-Tanaka micromechanical material model. The bending response is obtained using the first-order shear deformation theory (FSDT). The mixed interpolation of four elements of tensorial components (MITC4) is also implemented to overcome a well-known shear locking problem when the thickness becomes thinner. The Hamilton-Jacobi method is utilized in each iteration to enforce the necessary boundary conditions. The mathematical formulas are expressed in great detail for the LS-based TO using 3DFG materials. Several numerical examples are exhibited to verify the efficiency and reliability of the current methodology with the previously reported literature. Finally, the influences of FG materials in the optimized design are explained in detail to illustrate the behaviors of optimized structures.

• Structural Engineering and Mechanics
• /
• v.88 no.5
• /
• pp.439-449
• /
• 2023

Non-linear vibration characteristics of functionally graded CNT-reinforced composite (FG-CNTRC) cylindrical shell panel on elastic foundation have not been sufficiently examined. In this situation, this study aims at the profound numerical investigation of the non-linear vibration response of FG-CNTRC cylindrical panels on Winkler-Pasternak foundation by introducing an accurate and effective 2-D meshfree-based non-linear numerical method. The large-amplitude free vibration problem is formulated according to the first-order shear deformation theory (FSDT) with the von Karman non-linearity, and it is approximated by Laplace interpolation functions in 2-D natural element method (NEM) and a non-linear partial derivative operator H_NL. The complex and painstaking numerical derivation on the curved surface and the crucial shear locking are overcome by adopting the geometry transformation and the MITC3+ shell elements. The derived nonlinear modal equations are iteratively solved by introducing a three-step iterative solving technique which is combined with Lanczos transformation and Jacobi iteration. The developed non-linear numerical method is estimated through the benchmark test, and the effects of foundation stiffness, CNT volume fraction and functionally graded pattern, panel dimensions and boundary condition on the non-linear vibration of FG-CNTRC cylindrical panels on elastic foundation are parametrically investigated.

• Smart Structures and Systems
• /
• v.29 no.1
• /
• pp.1-16
• /
• 2022

Despite recent breakthroughs in deep learning and computer vision fields, the pixel-wise identification of tiny objects in high-resolution images with complex disturbances remains challenging. This study proposes a modified U-net for tiny crack segmentation in real-world steel-box-girder bridges. The modified U-net adopts the common U-net framework and a novel Self-Attention-Self-Adaption (SASA) neuron as the fundamental computing element. The Self-Attention module applies softmax and gate operations to obtain the attention vector. It enables the neuron to focus on the most significant receptive fields when processing large-scale feature maps. The Self-Adaption module consists of a multiplayer perceptron subnet and achieves deeper feature extraction inside a single neuron. For data augmentation, a grid-based crack random elastic deformation (CRED) algorithm is designed to enrich the diversities and irregular shapes of distributed cracks. Grid-based uniform control nodes are first set on both input images and binary labels, random offsets are then employed on these control nodes, and bilinear interpolation is performed for the rest pixels. The proposed SASA neuron and CRED algorithm are simultaneously deployed to train the modified U-net. 200 raw images with a high resolution of 4928 × 3264 are collected, 160 for training and the rest 40 for the test. 512 × 512 patches are generated from the original images by a sliding window with an overlap of 256 as inputs. Results show that the average IoU between the recognized and ground-truth cracks reaches 0.409, which is 29.8% higher than the regular U-net. A five-fold cross-validation study is performed to verify that the proposed method is robust to different training and test images. Ablation experiments further demonstrate the effectiveness of the proposed SASA neuron and CRED algorithm. Promotions of the average IoU individually utilizing the SASA and CRED module add up to the final promotion of the full model, indicating that the SASA and CRED modules contribute to the different stages of model and data in the training process.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.223-233
• /
• 2009

Three-dimensional (3-D) optimally-interpolated sea surface temperature (SST) field was produced by using AQUA/AMSR-E satellite data, and its limitations were described by comparing the temporal average of sea surface temperatures. The 3-D OI (Optimum Interpolation) SST showed a small error of less than $0.05^{\circ}C$ in the central North Pacific, but yielded large errors of greater than $0.4^{\circ}C$ at the coastal area where the satellite microwave data were not available. OI SST composite around pixels with no observation due to heavy rainfall or cloudy pixels had estimation errors of $0.1-0.15^{\circ}C$. Comparison with temporal means showed a tendency that overall OI SSTs were underestimated around heavy cloudy pixels and smoothed out by reducing the magnitude of SST fronts. In the low-latitude areas near the equator, OI SST field produced discontinuity, originated from the window size for the OI procedure. This was mainly caused by differences in the spatial scale of oceanic features. Infernal Rossby deformation radius, as a measure of spatial stale, showed dominant latitudinal variations with O(1) difference in the North Pacific. This study suggests that OI SST methodology should consider latitudinally-varying size of window and the characteristics of spatial scales of oceanic phenomena with substantial dependency on latitude and vertical structure of density.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.21-30
• /
• 2004

This study deals with free vibrations of laminated composite structures with a channel section using finite element method. In this paper, the mixed finite element method using Lagrangian and Hermite interpolation functions is adopted and a high-order plate theory is used to analyze laminated composite non-prismatic folded plates with a channel section more accurately for free vibration. The theory accounts for parabolic distribution of the transverse shear stress and requires no shear correction factors supposed in the first-order plate theory. An 32×32 matrix is assembled to transform the system element matrices from the local to global coordinates using a coordinate transformation matrix, in which an eighth drilling degree of freedom (DOF) per node is appended to the existing 7-DOF system. The results in this study are compared with those of available literatures for the conventional and first-order plate theory. Sample studies are carried out for various layup configurations and length-thickness ratio, and geometric shapes of plates. The significance of the high-order plate theory in analyzing complex composite structures with a channel section is enunciated in this paper.

• Journal of Biomedical Engineering Research
• /
• v.24 no.4
• /
• pp.259-266
• /
• 2003

Non-rigid registration between different modality images with shape deformation can be used to diagnosis and study for inter-patient image registration, longitudinal intra-patient registration, and registration between a patient image and an atlas image. This paper proposes a hierarchical registration method using bodily tissue based segmentation for registration between color images and CT images of the Visible Human leg areas. The cross-sectional color images and the axial CT images are segmented into three distinctive bodily tissue regions, respectively: fat, muscle, and bone. Each region is separately registered hierarchically. Bounding boxes containing bodily tissue regions in different modalities are initially registered. Then, boundaries of the regions are globally registered within range of searching space. Local boundary segments of the regions are further registered for non-rigid registration of the sampled boundary points. Non-rigid registration parameters for the un-sampled points are interpolated linearly. Such hierarchical approach enables the method to register images efficiently. Moreover, registration of visibly distinct bodily tissue regions provides accurate and robust result in region boundaries and inside the regions.

• The Journal of Korean Institute of Next Generation Computing
• /
• v.15 no.5
• /
• pp.64-74
• /
• 2019

Although a non-rigid registration has high demands in clinical practice, it has a high computational complexity and it is very difficult for ensuring the accuracy and robustness of registration. This study proposes a method of applying a non-rigid registration to 3D magnetic resonance images of brain in an unsupervised learning environment by using a deep-learning network. A feature vector between two images is produced through the network by receiving both images from two different patients as inputs and it transforms the target image to match the source image by creating a displacement vector field. The network is designed based on a U-Net shape so that feature vectors that consider all global and local differences between two images can be constructed when performing the registration. As a regularization term is added to a loss function, a transformation result similar to that of a real brain movement can be obtained after the application of trilinear interpolation. This method enables a non-rigid registration with a single-pass deformation by only receiving two arbitrary images as inputs through an unsupervised learning. Therefore, it can perform faster than other non-learning-based registration methods that require iterative optimization processes. Our experiment was performed with 3D magnetic resonance images of 50 human brains, and the measurement result of the dice similarity coefficient confirmed an approximately 16% similarity improvement by using our method after the registration. It also showed a similar performance compared with the non-learning-based method, with about 10,000 times speed increase. The proposed method can be used for non-rigid registration of various kinds of medical image data.