• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.1-6
• /
• 2005

Analogue physical modelling using granular materials (i.e., sandbox experiments) has been applied with great success to a number of geological problems at various scales. Such physical experiments can also be simulated numerically with the Discrete Element Method (DEM). In this study, we apply the DEM simulation to the collision between the Indian subcontinent and the Eurasian Plate, one of the most significant current tectonic processes in the Earth. DEM simulation has been applied to various kinds of dynamic modelling, not only in structural geology but also in soil mechanics, rock mechanics, and the like. As the target of the investigation is assumed to be an assembly of many tiny particles, DEM simulation makes it possible to treat an object with large and discontinuous deformations. However, in DEM simulations, we often encounter difficulties when we examine the validity of the input parameters, since little is known about the relationship between the input parameters for each particle and the properties of the whole assembly. Therefore, in our previous studies (Yamada et al.,2002a,2002b,2002c), we were obliged to tune the input parameters by trial and error. To overcome these difficulties, we introduce a numerical biaxial test with the DEM simulation. Using the results of this numerical test, we examine the validity of the input parameters used in the collision model. The resulting collision model is quite similar to the real deformation observed in eastern Asia, and compares well with GPS data and in-situ stress data in eastern Asia.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.22 no.1
• /
• pp.138-144
• /
• 2013

The automation technology for overlay welding is needed due to the occurrence of severe corrosion and abrasion on the surface of internal contact in different shape of fittings. In Korea, different shapes of fittings have been manufactured by using the imported equipment of overlay welding automation at some companies. Thus the research on the development of overlay welding automation system (in short, OWAS) for a large L-type tube is urgently needed. In this paper, the investigation is focused on the optimal design of a supporting base for the (currently developing) OWAS of large L-type tube. Specifically we assume that the base which supports the equipment during the process of overlay welding is loaded as self-weight in the direction of gravity through static analysis especially when it is rotated 180 degree on the OWAS. For optimal design of a supporting base for OWAS of large L-type tube, Solidworks(R) (for 3-dimensional modelling) and ANASYS Workbench(R) (for structural analysis) are incorporated so as to proceed an optimization routines based on Response Surface Method (RSM) and Design of Experiment (DOE). In more specific, DOE finds out major factors (or dimensions) of the supporting base by using MINITAB(R). Then the regression equations between design variables (the major factors of supporting base) and response variables (deformation, stress and safety factor for the supporting base), which will be resulted in by RSM, verify the major factors of DOE. In the next step, Central Composite Design (CCD) plans 20 simulations of ANASYS Workbench(R) and then figures out the optimal values of design variables which will be reflected on the manufacturing of supporting base. Finally welding experiment is conducted to figure out the influence of overlay welding quality in applying the optimized design values of supporting base to the actual OWAS.

• Journal of Korean Society of societal Security
• /
• v.2 no.2
• /
• pp.63-68
• /
• 2009

Seismic design of newly built bridges can be considered and carried out during construction process according to the revised road bridge design standard issued recently. While for the existing reinforced concrete bridge priers under service before new standard implements, their resistance capacity against lateral seismic loading is inferior. In this research, seismic reinforcing for existing bridge piers by nickel-chrome alloy bar has been analyzed. Based on the established model by MIDAS program, the behaviors of bridge piers including deformation and stress with and without nickel-chrome alloy reinforcing bars have been compared and discussed under lateral seismic loading. And the advantages of using nickel-chrome alloy bar as seismic reinforcement over other materials, such as good performance, good economy etc. have been demonstrated by comparison with other researches. Also the anti-seismic efficiency of nickel-chrome alloy reinforcing bars has been confirmed by MIDAS modeling analysis.

• Journal of the Korea Convergence Society
• /
• v.6 no.6
• /
• pp.241-246
• /
• 2015

In cold winter season, the apparatuses of heating and heater which warm up the interior of a room are necessary element and the used amount of these apparatuses from the year 2000 has been increased abruptly. But, the fire accident and the danger of fire are also increased. Therefore, 3D modelling is done by referring three kinds of the electric heaters as the heaters of ceramic, carbon and near infrared ray sold in the city for the design of more safe heating apparatuses in this study. The thermal analyses with these models are carried out and the durabilities due to the thermal deformation and stress are studied. By the background of the study results derived in this study ultimately, the durabilities of electric heater models due to each shape can be anticipated and contributed to the development of new heating apparatus with more safe resistance to fire. And it is possible to be grafted onto the convergence technique at design and show the esthetic sense.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.3
• /
• pp.187-195
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory), considering the effect of compression softening in cracked concrete, and macro-scopic and rotating crack models etc., was presented for the nonlinear behaviour of structural concrete. Considering the computational efficency and the ability of modelling the post-ultimate behaviour as major concerns, the Incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Non1inear Analysis of Structural Concrete by FEM : Monotonic Loading) developed in this way enables the predictions of strength and deformation capacities in a full range, crack patterns and their corresponding widths, and yield extents of reinforcement. As the verification purpose of NASCOM, the predictions were made for Bhide's Panel(PB21) and Leonhardt's deep beam tests. The predicted results shows somewhat stiff behaviour for the panel test, and vice versa for deep beam tests. More refining process would be necessary hereafter in terms of more accurately simulating the effects of tension-stiffening and compression softening in concrete.

• Journal of the Microelectronics and Packaging Society
• /
• v.20 no.4
• /
• pp.7-11
• /
• 2013

As thermosonic ball bonding is developed for more and more advanced applications in the electronic packaging industry, the control of process stresses induced on the integrated circuits becomes more important. If Cu bonding wire is used instead of Au wire, larger ultrasonic levels are common during bonding. For advanced microchips the use of Cu based wire is risky because the ultrasonic stresses can cause chip damage. This risk needs to be managed by e.g. the use of ultrasound during the impact stage of the ball on the pad ("pre-bleed") as it can reduce the strain hardening effect, which leads to a softer deformed ball that can be bonded with less ultrasound. To find the best profiles of ultrasound during impact, a numerical model is reported for ultrasonic bonding with capillary dynamics combined with a geometrical model describing ball deformation based on volume conservation and stress balance. This leads to an efficient procedure of ball bond modelling bypassing plasticity and contact pairs. The ultrasonic force and average stress at the bond zone are extracted from the numerical experiments for a $50{\mu}m$ diameter free air ball deformed by a capillary with a hole diameter of $35{\mu}m$ at the tip, a chamfer diameter of $51{\mu}m$, a chamfer angle of $90^{\circ}$, and a face angle of $1^{\circ}$. An upper limit of the ultrasonic amplitude during impact is derived below which the ultrasonic shear stress at the interface is not higher than 120 MPa, which can be recommended for low stress bonding.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.1C
• /
• pp.27-35
• /
• 2012

The relieving platform has the advantage of decreasing the total lateral earth pressure on the retaining wall and increasing the overall stability of the structure. Several modeling tests were performed to determine the earth pressure distribution on the retaining wall with a relieving platform and to compare it with that of the cantilever retaining wall. Different types of soil and angle of cutting surface were used to determine the effect of the soil characteristics and the backfill conditions on these earth pressure distributions. From the modeling tests, comparisons between the retaining wall with a relieving platform and the cantilever retaining wall show that the reduction of the lateral earth pressure and deformation of wall was indicated clearly on the retaining wall with a relieving platform. And the overall stability was increased by the relieving platform.

• Geomechanics and Engineering
• /
• v.12 no.2
• /
• pp.295-312
• /
• 2017

Cast in situ and grouted concrete helical piles with 150-200 mm diameter half cylindrical ribs have become an economical and effective choice in Shanghai, China for uplift piles in deep soft soils. Though this type of pile has been successful used in practice, the reinforcing mechanism and the contribution of the ribs to the total resistance is not clear, and there is no clear guideline for the design of such piles. To study the inclusion of ribs to the contribution of shear resistance, the shear behaviour between silty sand and concrete slabs with parallel ribs at different spacing and angles were tested in a large direct shear box ($600mm{\times}400mm{\times}200mm$). The front panels of the shear box are detachable to observe the soil deformation after the test. The tests were modelled with three-dimensional finite element method in ABAQUS. It was found that, passive zones can be developed ahead of the ribs to form undulated failure surfaces. The shear resistance and failure mode are affected by the ratio of rib spacing to rib diameter. Based on the shape and continuity of the failure zones at the interface, the failure modes at the interface can be classified as "punching", "local" or "general" shear failure respectively. With the inclusion of the ribs, the pull out resistance can increase up to 17%. The optimum rib spacing to rib diameter ratio was found to be around 7 based on the observed experimental results and the numerical modelling.

• Advances in aircraft and spacecraft science
• /
• v.3 no.1
• /
• pp.15-28
• /
• 2016

Although the aircraft industry was the first to use fibre composites, now they are increasingly used in a range of structural applications such as flooring, decking, platforms and roofs. Interlayer delamination is a major failure mode which threatens the reliability of composite structures. Delamination can grow in size under increasing loads with time and hence leads to severe loss of structural integrity and stiffness reduction. Delamination reduces the natural frequency and as a consequence may result in resonance. Hence, the study of the effects of delamination on the free vibration behaviour of multilayer composite structures is imperative. The focus of this paper is to develop a 3D FE model and investigate the free vibration behaviour of fibre composite multilayer sandwich panels with interlayer delaminations. A series of parametric studies are conducted to assess the influence of various parameters of concern, using a commercially available finite element package. Additionally, selected points in the delaminated region are connected appropriately to simulate bolting as a remedial measure to fasten the delamination region in the aim of reducing the effects of delamination. First order shear deformation theory based plate elements have been used to model each sandwich layer. The findings suggest that the delamination size and the end fixity of the plate are the most important factors responsible for stiffness reduction due to delamination damage in composite laminates. It is also revealed that bolting the delaminated region can significantly reduce the natural frequency variation due to delamination thereby improving the dynamic performance.