• Computers and Concrete
• /
• v.31 no.2
• /
• pp.163-171
• /
• 2023

Face gear transmission is widely used in aerospace shunt-confluence transmission system. Tooth wear is one of the main factors affecting its bearing transmission performance. Furthermore, the installation errors of face gear are inevitable. In order to study the wear mechanism of face gear tooth surface with installation errors, based on tooth contact analysis numerical method and Archard wear theory, the UMESHMOTION subroutine in ABAQUS is developed.Combining with Arbitrary Lagrangian-Eulerian adaptive mesh technology, the finite element mesh wear model of abraded face gear pair is established.The preprocessing conditions are set to generate the inp files.Then,the inp files for each corner are imported and batch processed in ABAQUS.The loading tooth contact problem at each rotation angle is solved and the load distribution coefficient among gear tooth, tooth root bending stress, tooth surface contact stress and loaded transmission error are obtained. Results show that the tooth root wear is the most serious and the wear at the pitch cone is close to 0.The wear law of tooth surface along tooth width direction is convex parabola and the wear law along tooth height direction is concave parabola.

• Nuclear Engineering and Technology
• /
• v.56 no.7
• /
• pp.2498-2515
• /
• 2024

In this study, we propose a fully parallelized adaptive particle refinement (APR) algorithm for smoothed particle hydrodynamics (SPH) to construct a stable and efficient multi-resolution computing system for nuclear safety analysis. The APR technique, widely employed by SPH research groups to adjust local particle resolutions, currently operates on a serialized algorithm. However, this serialized approach diminishes the computational efficiency of the system, negating the advantages of acceleration achieved through high-performance computing devices. To address this drawback, we propose a fully parallelized APR algorithm designed to enhance both efficiency and computational accuracy, facilitated by a new adaptive smoothing length model. For model validation, we simulated both hydrostatic and hydrodynamic benchmark cases in 2D and 3D environments. The results demonstrate improved computational efficiency compared to the conventional SPH method and APR with a serialized algorithm, and the model's accuracy was confirmed, revealing favorable outcomes near the resolution interface. Through the analysis of jet breakup, we verified the performance and accuracy of the model, emphasizing its applicability in practical nuclear safety analysis.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.26 no.1
• /
• pp.1-17
• /
• 2024

Analyzing the tunnel excavation behavior and its effect on the surrounding ground involves large deformation behavior. Therefore, in order to properly simulate the tunnel excavation process and rigorously investigate the actual effect of excavation on surrounding ground and tunnel structure large deformation analysis method is required. In this study, two major numerical approaches capable of considering large deformations behavior were applied to investigate the effect of tunnel boring machine excavation on the surrounding ground: coupled Eulerian-Lagrangian (CEL) and the automatic remeshing (AR) method. Relative performance of both approaches was evaluated through the ground response due to TBM excavation. The ground response will be quantified by estimating the range of the excavation damaged zone (EDZ). By comparing the results, the range of the EDZ will be suggested on the vertical and horizontal direction along the TBM excavation surface. Based on the computed results, it was found that the size of EDZ around the excavation surface and the tendencies was in good agreement among the two approaches. Numerical results clearly show that the size of the EDZ around the tunnel tends to be larger for rock with higher RMR rating. The size of the EDZ is found to be direct proportional to the tunnel diameter, whereas the depth of the tunnel is inversely proportional due to higher confinement stress around the excavation surface.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.6
• /
• pp.556-563
• /
• 2019

The ice keel gouge and seabed interaction is one of the major considerations in the design of an Arctic pipeline system. Ice keel and seabed interaction engineering models based on experimental data, which give an explicit equation for estimating the ice gouging depth, have been suggested. The suggested equations usually overestimate the ice keel gouging depth. In addition, various types of numerical analyses have been carried out to verify the suggested engineering model equations in comparison to the experimental data. However, most of numerical analysis results were also overestimated compared with the laboratory experimental data. In this study, a numerical analysis considering the contact condition and geostatic stress was carried out to predict the ice keel gouging depth and compared with the previous studies. Considering the previously mentioned conditions, more accurate results were produced compared with the laboratory experiment results and the error rate was reduced compared to previous numerical analysis studies.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2014.10a
• /
• pp.411-414
• /
• 2014

Numerical modeling was conducted to estimate the amount of dislocation that may occur across a frictionless fracture during an earthquake using commercial code FLAC3D (Fast Lagrangian Analysis of Continua in 3 Dimensions). The applied motion was calculated to represent a Richter 6.0 magnitude earthquake at distances of 2 km from the fracture. The velocity-time history was generated from Svensk $K{\ddot{a}}arnbr{\ddot{a}}anslehantering$ AB report. In the report, The velocity field resulting from an earthquake on a fault located in the near-field (2 km distance) was modelled using a finite difference program, WAVE. The stress-time history was substituted for velocity-time history to perform dynamic analysis using FLAC3D. During the earthquake, the maximum dislocation and change of shear stress were about 1 cm and 2MPa, respectively. Because the fracture is frictionless in this study, all dislocations relax to zero after the earthquake motions have ceased.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.137-156
• /
• 2013

In this paper, a meshfree shell adaptive procedure is developed for the applications in the sheet metal forming simulation. The meshfree shell formulation is based on the first-order shear deformable shell theory and utilizes the degenerated continuum and updated Lagrangian approach for the nonlinear analysis. For the sheet metal forming simulation, an h-type adaptivity based on the meshfree background cells is considered and a geometric error indicator is adopted. The enriched nodes in adaptivity are added to the centroids of the adaptive cells and their shape functions are computed using a first-order generalized meshfree (GMF) convex approximation. The GMF convex approximation provides a smooth and non-negative shape function that vanishes at the boundary, thus the enriched nodes have no influence outside the adapted cells and only the shape functions within the adaptive cells need to be re-computed. Based on this concept, a multi-level refinement procedure is developed which does not require the constraint equations to enforce the compatibility. With this approach the adaptive solution maintains the order of meshfree approximation with least computational cost. Two numerical examples are presented to demonstrate the performance of the proposed method in the adaptive shell analysis.

• Journal of Ocean Engineering and Technology
• /
• v.19 no.1 s.62
• /
• pp.20-25
• /
• 2005

The authors performed the underwater explosion analysis for the liquified oxygen tank - a kind of fuel tank of a mid-size submarine, and tried to verify the structural safety for this structure. First, the authors reviewed the theory and application of underwater explosion analysis, using a Structure-Fluid Interaction technique and its finite element modeling scheme. Next, the authors modeled the explosive and sea water as fluid elements, the LOX tank as structural elements, and the interface between the two regions as the ALE scheme. The effect on shock pressure and impulse of fluid mesh size and shape are also investigated. Upon analysis, it was found that the shock pressure due to explosion propagated into the water region, and hit the structure region. The plastic deformation and the equivalent stress were apparent at the web frame and the shock mount of LOX structure, but these values were acceptable for the design criteria.

• Structural Engineering and Mechanics
• /
• v.22 no.3
• /
• pp.351-370
• /
• 2006

In this study, the earthquake damage response of the concrete arch dams was investigated including dam-reservoir interaction. A continuum damage model which is a second-order tensor and includes the strain softening behavior was selected for the concrete material. Fluid-structure interaction problem was modeled by Lagrangian approach. Sommerfeld radiation condition was applied to the truncated boundary of reservoir. The improved form of the HHT-${\alpha}$ time integration algorithm was used in the solution of the equations of motion. The arch dam Type 5 was selected for numerical application. For the dynamic input, acceleration records of the 10 December 1967 Koyna earthquake were chosen. These records were scaled with earthquake acceleration scale factor (EASF) and then used in the analyses. Solutions were obtained for empty and full reservoir cases. The effects of EASF and damping ratio on the response of the dam were studied.

• Journal of Environmental Impact Assessment
• /
• v.6 no.2
• /
• pp.81-92
• /
• 1997

This paper presents the results from the comparison analysis and evaluation between the air pollutant dispersion modeling results and the observation data in the area within a 10 km radius from the Boryong thermal power plants. The observation data used in this study were the air pollutant concentrations which had been continuously measured from 8 locations around the Boryong power plants by TMS(tele-monitoring system) for 3 months from September to November, 1996. The short-term and long-term predictions were carried out using ISC3 model and LPDM(Lagrangian Panicle Dispersion Model). The results of ISC3 modeling in a short-term showed highly as 0.7 in a correlation coefficient, but in a long-term showed just 0.54. On the other hand, LPDM showed 0.78 in a correlation coefficient for a long-term, but in a short-term showed highly value than the observation concentrations.