• Geomechanics and Engineering
• /
• v.18 no.3
• /
• pp.225-234
• /
• 2019

A numerical stepwise approach for cavity expansion problem in strain-softening rock or soil mass is investigated, which is compatible with Mohr-Coulomb and generalized Hoek-Brown failure criteria. Based on finite difference method, plastic region is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium and compatibility equations, the material parameters of the rock or soil mass are assumed to be the same in each ring. For the strain-softening behavior, the strength parameters are assumed to be a linear function of deviatoric plastic strain (${\gamma}p^*$) for each ring. Increments of stress and strain for each ring are calculated with the finite difference method. Assumptions of large-strain for soil mass and small-strain for rock mass are adopted, respectively. A new numerical stepwise approach for limited pressure and plastic radius are obtained. Comparisons are conducted to validate the correctness of the proposed approach with Vesic's solution (1972). The results show that the perfectly elasto-plastic model may underestimate the displacement and stresses in cavity expansion than strain-softening coefficient considered. The results of limit expansion pressure based on the generalised H-B failure criterion are less than those obtained based on the M-C failure criterion.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 1998.11a
• /
• pp.27-31
• /
• 1998

Most models for the simulation of multi-phase flow and multi-species transport employ the capillary approach which uses the Darcy's law for the representation of mass flux of each phase. The capillary approach based on the Darcy's law require many empirical coefficients with complex functional dependencies rather than rigrous mathematical and physical formulation. The shortcoming of the capillary approach cause the numerical errors in the simulations by the multi-phase flow and transport models. This study discuss some of the problems related with the use of models.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.18 no.1
• /
• pp.143-149
• /
• 1993

An approach is presented for efficient numerical integration of the Sormnerfeld type integrals pertinent to the microstrip surface Green's function arising in the problem of an electric current point source on an infinite planar grounded dielectric substrate. This approach, valid for both lossless and lossy dielectric substrates, is based on the deformation of the integration contour via a coordinate transformation and Cauchy's residue theory, and identifies clearly the effects of surface waves. I ts useful application is in a rigorous moment method analysis of micros trip antenna arrays and microstrip guided wave structures. The efficiency and the usefulness of the present approach are emphasized through some numerical calculations of the impedance matrix elements with associated CPU times.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.89-95
• /
• 2006

This study has numerically modelled the combustion processes of the turbulent swirling premixed lifted flames in the low-swirl burner (LSB). In these turbulent swirling premixed flames, the four tangentially-injected air jets induce the turbulent swirling flow which plays the crucial role to stabilize the turbulent lifted flame. In the present approach, the turbulence-chemistry interaction is represented by the level-set based flamelet model. Two-dimensional and three-dimensional computations are made for the various swirl numbers and nozzle length. In terms of the centerline velocity profiles and flame liftoff heights, numerical results are compared with experimental data The three-dimensional approach yields the much better conformity with agreements with measurements without any analytic assumptions on the inlet swirl profiles, compared to the two-dimensional approach. Numerical clearly results indicate that the present level-set based flamelet approach has realistically simulated the structure and stabilization mechanism of the turbulent swirling stoichiometric and lean-premixed lifted flames in the low-swirl burner.

• Smart Structures and Systems
• /
• v.15 no.5
• /
• pp.1311-1327
• /
• 2015

In this paper, we have applied a new class of approximate analytical methods called Variational Approach (VA) for high nonlinear vibration equations. Three examples have been introduced and discussed. The effects of important parameters on the response of the problems have been considered. Runge-Kutta's algorithm has been used to prepare numerical solutions. The results of variational approach are compared with energy balance method and numerical and exact solutions. It has been established that the method is an easy mathematical tool for solving conservative nonlinear problems. The method doesn't need small perturbation and with only one iteration achieve us to a high accurate solution.

• Steel and Composite Structures
• /
• v.48 no.6
• /
• pp.641-648
• /
• 2023

In this article, we explore the issue concerning semiconductors half-space comprised of materials with varying thermal conductivity. The problem is within the framework of the generalized thermoelastic model under one thermal relaxation time. The half-boundary space's plane is considered to be traction free and is subjected to a thermal shock. The material is supposed to have a temperature-dependent thermal conductivity. The numerical solutions to the problem are achieved using the finite element approach. To find the analytical solution to the linear problem, the eigenvalue approach is used with the Laplace transform. Neglecting the new parameter allows for comparisons between numerical findings and analytical solutions. This facilitates an examination of the physical quantities in the numerical solutions, ensuring the accuracy of the proposed approach.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.5-9
• /
• 2003

Stress and deflection of Active Fiber Composite(AFC) embedded and/or attached composite structures are numerically investigated at the constituent level by the Direct Numerical Simulation(DNS). The DNS approach which models and simulates the fiber and matrix directly using 3D finite elements need to be solved by efficient way. To handle this large scale problem, parallel program for solving piezoelectric behavior was developed and run on the parallel computing environment. Also, the stress result from DNS approach is compared with that from uniform field model.

• Structural Engineering and Mechanics
• /
• v.87 no.5
• /
• pp.461-473
• /
• 2023

This study aims to solve for nonlinear cylindrical shell systems with a semi-analytical and numerical approach implementing the P-T method. The procedures and conditions for such a study are presented in practically solving and analyzing the cylindrical shell systems. An analytical model for a nonlinear thick cylindrical shell (TCS) is established on the basis of the stress function and Reddy's higher-order shear deformation theory (HSDT). According to Reddy's HSDT, Hooke's law in three dimensions, and the von-Kármán equation, the stress-strain relations are developed for the thick cylindrical shell systems, and the three coupled nonlinear governing equations are thus established and discretized as per the Galerkin method, for implementing the P-T method. The solution generated with the approach is continuous everywhere in the entire time domain considered. The approach proposed can also be used to numerically solve and analyze the nonlinear shell systems. The procedures and recurrence relations for numerical solutions of shell systems are presented. To demonstrate the application of the approach in numerically solving for nonlinear cylindrical shell systems, a specific nonlinear cylindrical shell system subjected to an external excitation is solved numerically. In numerically solving for the system, the present approach shows higher efficiency, accuracy, and reliability in comparison with that of the Runge-Kutta method. The approach with the P-T method presented is practically sound especially when continuous and high-quality numerical solutions for the shell systems are considered.

• Computers and Concrete
• /
• v.23 no.1
• /
• pp.69-80
• /
• 2019

Modeling approach for mesoscopic model of concrete depicting mass transportation and physicochemical reaction is important since there is growing demand for accuracy and computational efficiency of numerical simulation. Mesoscopic numerical simulation considering binder, aggregate and Interfacial Transition Zone (ITZ) generally produces huge number of DOFs, which is inapplicable for full structure. In this paper, a three-dimensional multiscale approach describing three-phase structure of concrete was discussed numerically. An effective approach generating random aggregate in polygon based on checking centroid distance was introduced. Moreover, ITZ elements were built by parallel expanding the surface of aggregates on inner side. By combining mesoscopic model including full-graded aggregate and macroscopic model, cases related to diffusivity and thickness of ITZ, volume fraction and grade of aggregate were studied regarding the consideration of multiscale compensation. Results clearly showed that larger analysis model in multiscale model expanded the diffusion space of chloride ion and decreased chloride content in front of rebar. Finally, this paper addressed some worth-noting conclusions about the chloride distribution and rebar corrosion regarding the configuration of, rebar diameter, concrete cover and exposure period.