• Structural Engineering and Mechanics
• /
• v.88 no.1
• /
• pp.83-94
• /
• 2023

An accurate and highly efficient inverse element labelled iPCB is developed based on the inverse finite element method (iFEM) for real-time shape estimation of plane-curved structures (such as arch bridges) utilizing onboard strain data. This inverse problem, named shape sensing, is vital for the design of smart structures and structural health monitoring (SHM) procedures. The iPCB formulation is defined based on a least-squares variational principle that employs curved Timoshenko beam theory as its baseline. The accurate strain-displacement relationship considering tension-bending coupling is used to establish theoretical and measured section strains. The displacement fields of the isoparametric element iPCB are interpolated utilizing nonuniform rational B-spline (NURBS) basis functions, enabling exact geometric modelling even with a very coarse mesh density. The present formulation is completely free from membrane and shear locking. Numerical validation examples for different curved structures subjected to different loading conditions have been performed and have demonstrated the excellent prediction capability of iPCBs. The present formulation has also been shown to be practical and robust since relatively accurate predictions can be obtained even omitting the shear deformation contributions and considering polluted strain measures. The current element offers a promising tool for real-time shape estimation of plane-curved structures.

• Computers and Concrete
• /
• v.33 no.1
• /
• pp.91-102
• /
• 2024

Beam-shaped components commonly rotate along a fixed axis when massive mechanical structures like rotors, jet engine blades, motor turbines, and rotating railway crossings perform their functions. For these structures to be useful in real life, their mechanical behavior is essential. Therefore, this is the first article to use the modified shear deformation theory type hyperbolic sine functions theory and the FEM to study the static bending response of rotating functionally graded GPL-reinforced composite (FG-GPLRC) beams with initial geometrical deficiencies in thermal media. Graphene platelets (GPLs) in three different configurations are woven into the beam's composition to increase its strength. By comparing the numerical results with those of previously published studies, we can assess the robustness of the theory and mechanical model employed in this study. Parameter studies are performed to determine the effect of various geometric and physical variables, such as rotation speed and temperature, on the bending reactions of structures.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• The Journal of Engineering Geology
• /
• v.3 no.2
• /
• pp.125-148
• /
• 1993

Groundwater flow in fractured rock masses is controlled by combined effects of fracture networks, state of geostafic stresses and crossflow between fractures and rock matrix. Furthermore the scaie dependent, anisotropic properties of hydraulic parameters results mainly from irregular paftems of fracture system, which can not be evaluated properly with the methods available at present. The basic assumpfion of discrete fracture network model is that groundwater flows only along discrete fractures and the flow paths in rock mass are determined by geometric paftems of interconnected fractures. The characteristics of fracture distribution in space and fracture hydraulic parameters are represented as the probability density functions by stochastic simulation. The discrete fracture network modelling was aftempted to characterize the groundwater flow in the vicinity of existing large cavems located in Wonjeong-ri, Poseung-myon, Pyeungtaek-kun. The fracture data of $1\textrm{km}^2$ area were analysed. The result indicates that the fracture sets evaluated from an equal area projection can be grouped into 6 sets and the fracture sizes are distributed in longnormal. The conductive fracture density of set 1 shows the highest density of 0.37. The groundwater inflow into a carvem was calculated as 29ton/day with the fracture transmissivity of $10^{-8}\textrm{m}^2/s$. When the fracture transmissivity increases in an order, the inflow amount estimated increases dramatically as much as fold, i.e 651 ton/day. One of the great advantages of this model is a forward modelling which can provide a thinking tool for site characterization and allow to handle the quantitative data as well as qualitative data.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Journal of Korean Society for Geospatial Information Science
• /
• v.16 no.3
• /
• pp.117-126
• /
• 2008

Earth geopotential models were used to determine the reference surface in geoid modelling and many global geopotential models were developed after 1980's. Nowadays, EGM96 and EIGEN-GL04C geopotential models have been most widely used in the world, but what so called EGM2008 earth geopotential model were developed in 2008 by NGA. In this paper, we intended to compare the results of spherical harmonic analyses using the three geopotential model, EGM96, EIGEN-GL04C and EGM2008. So, the spherical harmonic analyses were performed up to degree and order 360(in case of EGM2008, up to degree and order 720, 1440, 2190 in addition), on each $1'{\times}1'$grid point in and around Korean peninsula. Geometric geoid were calculated at 464 GPS/leveling points for accuracy evaluation and then the results of three geopotential models were compared to geometric geoid. The results show that the accuracy of EGM2008 is improved considerablely compared to EGM96 and EIGEN-GL04C and it is possible to calculate geoidal heights within 14cm standard deviation and 5.5cm standard deviation after LSC fitting in and around Korean peninsula using EGM2008 geopotential model.

• Journal of the Institute of Convergence Signal Processing
• /
• v.7 no.3
• /
• pp.108-115
• /
• 2006

In this paper, a recognition algorithm of the straight line components of lane markings and an obstacle in the travelling lane region is proposed. This algorithm also involve the pitching error correction algorithm due to traveling vehicle's fluctuation. In order to reduce their error a practical road image modelling algorithm using V.F. model and camera calibration procedure are suggested to adapt the geometric variations. It is obtained the 3D world coordinate data by the 2D road images. In experimental test, we showed that this algorithm is available to recognize lane markings and an obstacle in the traveling lane.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.5
• /
• pp.1479-1484
• /
• 2014

Based on the full optimized molecular geometric structures at B3LYP/cc-pvtz method, a new designed compound, 4,6,8-trinitro-4,5,7,8-tetrahydro-6H-furazano[3,4-f ]-1,3,5-triazepine was investigated in order to look for high energy density compounds (HEDCs). The analysis of the molecular structure indicates that the seven-membered ring adopts chair conformation and there exist intramolecular hydrogen bond interactions. IR spectrum and heat of formation (HOF) were predicted. The detonation velocity and pressure were evaluated by using Kamlet-Jacobs equations based on the theoretical density and condensed HOF. The bond dissociation energies and bond orders for the weakest bonds were analyzed to investigate the thermal stability of the title compound. The results show that $N_1-N_6$ bond is the trigger bond. The crystal structure obtained by molecular mechanics belongs to $Pna2_1$ space group, with lattice parameters Z = 4, a = 15.3023 ${\AA}$, b = 5.7882 ${\AA}$, c = 11.0471 ${\AA}$, ${\rho}=2.06gcm^{-3}$. In addition, the analysis of frontier molecular orbital shows the title compound has good stability and high chemical hardness.

• Computers and Concrete
• /
• v.22 no.4
• /
• pp.395-405
• /
• 2018

In this paper, the performance of post-installed adhesive bonded anchor embedded in concrete is assessed using numerical simulations. This study aims at studying the influence of parameters on the performance of a chemically bonded anchorage system. Non-linear finite element modelling and simulations are carried out by properly using the material properties and phenomenon. Materials parameters such as characteristic length, fracture energy, damage criteria, tension retention and crack width of concrete and interface characteristics are carefully assigned so as to obtain a most realistic behaviour of the chemical anchor system. The peak strength of two different anchor systems obtained from present numerical studies is validated against experimental results. Furthermore, validated numerical models are used to study the load transferring mechanism and damage progression characteristics of various anchors systems where strength of concrete, strength of epoxy, and geometry and disposition of anchors are the parameters. The process of development of strain in concrete adjacent to the anchor and energy dissipated during the course of damage progression are analysed. Results show that the performance of the considered anchorage system is, though a combined effect of material and geometric parameters, but a clear distinction could be made on the parameters to achieve a desired performance based on strength, slip, strain development or dissipated energy. Inspite the increase in anchor capacity with increase in concrete strength, it brings some undesirable performance as well. Furthermore, the pullout capacity of the chemical anchor system increases with a decrease in disparity among the strength of concrete and epoxy.