• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.3
• /
• pp.361-369
• /
• 2001

This paper deals with numerical efficiency of nonlinear solution technique for space trusses. It will propose the combined Arc-length method to trace structural behavior after reaching buckling load as opposed to the current Arch-length method. The combined Arc-length method uses the current stiffness parameter as a control variable. It uses Secant-Newton method in stable path and applies Arc-length method in unstable path. To evaluate efficiency of solution technique, the accuracy of solution, convergence, and computing time concerning illustrative numerical examples are compared with the current Arc-length method. It show that the combined Arc-length method, as proposed in this paper, is superior to the current Arc-length method in numerical nonlinear analysis.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.418-423
• /
• 2008

Among various sensitivity evaluation techniques, semi-analytical method is quite popular since this method is more advantageous than analytical method and global finite difference method. However, SAM reveals severe inaccuracy problem when relatively large rigid body motions are identified for individual elements. Such errors result from the numerical differentiation of the pseudo load vector calculated by the finite difference scheme. In the present study, the adjoint variable method combined with complex variable is proposed to obtain the shape and size sensitivity for structural optimization. The complex variable can present accurate results regardless of the perturbation size as well as easy to be implemented. Through a few numerical examples of the static problem for the structural sensitivity, the efficiency and reliability of the adjoint variable method combined with complex variable is demonstrated.

• Geomechanics and Engineering
• /
• v.29 no.4
• /
• pp.391-405
• /
• 2022

A newly proposed grouting simulation method, the sequential diffusion solidification method was introduced into the numerical simulation of combined bi-borehole grouting. The traditional, critical and difficult numerical problem for the temporal and spatial variation simulation of the slurry is solved. Thus, numerical simulation of grouting and blocking of flowing water in karst conduits is realized and the mechanism understanding of the combined bi-borehole technology is promoted. The sensitivity analysis of the influence factors of combined bi-borehole grouting was investigated. Through orthogonal experiment, the influences of proximal and distal slurry properties, the initial flow velocity of the conduit and the proximal and distal slurry injection rate on the blocking efficiency are compared. The velocity variation, pressure variation and slurry deposition phenomenon were monitored, and the flow field characteristics and slurry outflow behavior were analyzed. The interaction mechanism between the proximal and distal slurries in the combined bi-borehole grouting is revealed. The results show that, under the orthogonal experiment conditions, the slurry injection rate has the greatest impact on blocking. With a constant slurry injection rate, the blocking efficiency can be increased by more than 30% when using slurry with weak time-dependent viscosity behavior in the distal borehole and slurry with strong time-dependent viscosity behavior in the proximal borehole respectively. According to the results of numerical simulation, the grouting scheme of "intercept the flow from the proximal borehole by quick-setting slurry, and grout cement slurry from the distal borehole" is put forward and successfully applied to the water inflow treatment project of China Resources Cement (Pingnan) Limestone Mine.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.709-714
• /
• 2000

A splitting method for the direct numerical simulation of solid-liquid mixtures is presented, where a symmetric pressure equation is newly proposed. Through numerical experiment, it is found that the newly proposed splitting method works well with a matrix-free formulation fer some bench mark problems avoiding an erroneous pressure field which appears when using the conventional pressure equation of a splitting method. When deriving a typical pressure equation of a splitting method, the motion of a solid particle has to be approximated by the 'intermediate velocity' instead of treating it as unknowns since it is necessary as a boundary condition. Therefore, the motion of a solid particle is treated in such an explicit way that a particle moves by the known form drag (pressure drag) that is calculated from the pressure equation in the previous step. From the numerical experiment, it was shown that this method gives an erroneous pressure field even for the very small time step size as a particle velocity increases. In this paper, coupling the unknowns of particle velocities in the pressure equation is proposed, where the resulting matrix is reduced to the symmetric one by applying the projector of the combined formulation. It has been tested over some bench mark problems and gives reasonable pressure fields.

• Structural Engineering and Mechanics
• /
• v.42 no.6
• /
• pp.867-882
• /
• 2012

Analytical (Rayleigh-Ritz method) and numerical studies are carried out and buckling interaction curves are developed for simply supported plates of varying aspect ratios ranging from 1 to 5, under the combined action of in-plane shear and tension. A multi-step buckling procedure is employed in the Finite Element (FE) model instead of a regular single step analysis in view of obtaining the buckling load under the combined forces. Both the analytical (classical) and FE studies confirm the delayed shear buckling characteristics of thin plate under the combined action of shear and tension. The interaction curves are found to be linear and are found to vary with plate aspect ratio. The interaction curve developed using Rayleigh-Ritz method is found to deviate in an increasing trend from that of validated FE model as plate aspect ratio is increased beyond value of 1. It is found that the observed deviation is due to the insufficient number of terms that is been considered in the assumed deflection function of Rayleigh-Ritz method and a convergence study is suggested as a solution.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.21 no.6
• /
• pp.1453-1462
• /
• 1996

A numerical simulation is carried out on the directiveity synthesis of ultrasonic transducers by point source array. Directive method with combined LMS(Least-Mean-Square) method is practiced by means of a iterative method to realize the desired directivity. The directiviey of quasi-ideal beam with a beam width and a directive arbitrary specified was chosen. On the numerical resut, Proposed algorithm shows higher speed of clculating simulation than that of LMS method, and make adaptive control, which enables the desired directivity. Numerical simulations are carried out by PC(CPU:80486 DX2, RAM 16Mbyte).

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.152-157
• /
• 1999

For the simulation of the blood flow in coronary artery, the system modeling of coronary hemodynamics is combined with CFD technique. The blood flow in coronary artery interacts with the global coronary circulation. Especially in case of the coronary artery with stenosis, the interaction plays an important role in the hemodynamics of the circulation. In this study we present a combined numerical approach using both the CFD technique for flow simulation and the global system model of coronary circulation. We use a lumped parameter model for the global simulation of coronary circulation whereas the finite element method is employed to compute the viscous flow field in stenosed coronary artery, The time variation of the pressure drop due to stenosis is obtained from the proposed numerical method. Numerical results shows that the flow resistance and pressure drop due to stenosis has a relatively large value in systole.

• Journal of electromagnetic engineering and science
• /
• v.9 no.4
• /
• pp.211-217
• /
• 2009

A three-dimensional locally one-dimensional multiresolution time-domain(LOD-MRTD) method is introduced and unconditional stability is proved analytically. The updating formulations have fewer terms on the right-hand side than those of an alternating direction implicit MRTD(ADI-MRTD). The validation of the method is presented using the resonance frequency problem of an empty cavity. The reduction of the numerical dispersion technique is also combined with the proposed method. The numerical examples show that the combined method can improve the accuracy significantly.

• Structural Engineering and Mechanics
• /
• v.16 no.3
• /
• pp.341-360
• /
• 2003

The goal of this paper is to determine the eigenvalues of a uniform rectangular plate carrying any number of spring-damper-mass systems using an analytical-and-numerical-combined method (ANCM). To this end, a technique was presented to replace each "spring-damper-mass" system by a massless equivalent "spring-damper" system with the specified effective spring constant and effective damping coefficient. Then, the mode superposition approach was used to transform the partial differential equation of motion into the matrix equation, and the eigenvalues of the complete system were determined from the associated characteristic equation. To verify the reliability of the presented theory, all numerical results obtained from the ANCM were compared with those obtained from the conventional finite element method (FEM) and good agreement was achieved. Since the order of the property matrices for the equation of motion obtained from the ANCM is much lower than that obtained from the FEM, the CPU time required by the ANCM is much less than that by the FEM.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.2
• /
• pp.125-136
• /
• 2018

In this article, a homotopy perturbation transform method (HPTM) and the Laplace transform combined with Taylor expansion method are presented for solving Volterra integral equations with a convolution kernel. The (HPTM) is innovative in Laplace transform algorithm and makes the calculation much simpler while in the Laplace transform and Taylor expansion method we first convert the integral equation to an algebraic equation using Laplace transform then we find its numerical inversion by power series. The numerical solution obtained by the proposed methods indicate that the approaches are easy computationally and its implementation very attractive. The methods are described and numerical examples are given to illustrate its accuracy and stability.