• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.05a
• /
• pp.216-221
• /
• 2004

In two dimensional incompressible flaws, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint and hence achieve an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flaw analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flaw showed the present procedure to be stable, very efficient and useful in flaw analyses in conjunction with hierarchical elements.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.5 s.54
• /
• pp.11-18
• /
• 2003

In two-dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure (HIP) has been implemented on a SUPG finite element formulation. By using the SUPG formulation, one can escape from the LBB constraint hence, achieving an equal order formulation. In this paper, we increased the order of interpolation up to cubic. The conjugate gradient squared (CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper and efficient iterative procedure for higher order finite element formulation has not been available, thus far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient, and useful in flow analyses, in conjunction with hierarchical elements.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far. The numerical results by the present HIP for the lid driven cavity flow showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Structural Engineering and Mechanics
• /
• v.10 no.5
• /
• pp.427-440
• /
• 2000

The formulas for the perturbation analysis with modes of close eigenvalues are derived in this paper. Emphasis is made on the consistency of the straightforward perturbation process, given the complete terms of perturbations in the zeroth-order, which is a form of Rayleigh quotient, and in the higher-orders. By dividing the perturbation of eigenvector into two parts, the first-order perturbation with respect to the modes of close eigenvalues is moved into the zeroth-order perturbation. The normality condition is employed to compute the higher-order perturbations of eigenvector. The algorithm can be condensed to a single mode with a distinct eigenvalue, and this can accelerate the convergence of the perturbation analysis. The example confirms that the perturbation approximation obtained from the suggested procedure is in a good accuracy on the eigenvalues, eigenvectors, and normality.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.659-683
• /
• 2014

The static analysis of structures with arbitrary cross-section geometry and material lamination via a refined one-dimensional (1D) approach is presented in this paper. Higher-order 1D models with a variable order of expansion for the displacement field are developed on the basis of Carrera Unified Formulation (CUF). Classical Euler-Bernoulli and Timoshenko beam theories are obtained as particular cases of the first-order model. Numerical results of displacement, strain and stress are provided by using the finite element method (FEM) along the longitudinal direction for different configurations in excellent agreement with three-dimensional (3D) finite element solutions. In particular, a layered thin-walled cylinder is considered as first assessment with a laminated conventional cross-section. An atherosclerotic plaque is introduced as a typical structure with arbitrary cross-section geometry and studied for both the homogeneous and nonhomogeneous material cases through the 1D variable kinematic models. The analyses highlight limitations of classical beam theories and the importance of higher-order terms in accurately detecting in-plane cross-section deformation without introducing additional numerical problems. Comparisons with 3D finite element solutions prove that 1D CUF provides remarkable three-dimensional accuracy in the analysis of even short and nonhomogeneous structures with arbitrary geometry through a significant reduction in computational cost.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.199-211
• /
• 2024

The loading capacity of engineering structures/components reduces after the initiation and propagation of crack eventually leads to the final failure. Hence, it becomes essential to deal with the crack and its effects at the design and simulation stages itself, by detecting the prone area of the fracture. The phase-field (PF) method has been accepted widely in simulating fracture problems in complex geometries. However, most of the PF methods are formulated with second order continuity theoryinvolving C⁰ continuity. In the present study, PF method based on fourth-order (i.e., higher order) theory, maintaining C¹ continuity has been proposed for ductile fracture simulation. The formulation includes fourth-order derivative terms of phase field variable, varying between 0 and 1. Applications of fourth-order PF theory to ductile fracture simulation resulted in novelty in this area. The proposed formulation is numerically solved using a two-dimensional finite element (FE) framework in 3-layered manner system. The solutions thus obtained from the proposed fourth order theory for different benchmark problems portray the improvement in the accuracy of the numerical results and are well matched with experimental results available in the literature. These results are also compared with second-order PF theory and a comparison study demonstrated the robustness of the proposed model in capturing ductile behaviour close to experimental observations.

• Geomechanics and Engineering
• /
• v.13 no.3
• /
• pp.385-410
• /
• 2017

This work presents a simple and refined nth-order shear deformation theory for mechanical and thermal buckling behaviors of functionally graded (FG) plates resting on elastic foundation. The proposed refined nth-order shear deformation theory has a new displacement field which includes undetermined integral terms and contains only four unknowns. Governing equations are obtained from the principle of minimum total potential energy. A Navier type analytical solution methodology is also presented for simply supported FG plates resting on elastic foundation which predicts accurate solution. The accuracy of the present model is checked by comparing the computed results with those obtained by classical plate theory (CPT), first-order shear deformation theory (FSDT) and higher-order shear deformation theory (HSDT). Moreover, results demonstrate that the proposed theory can achieve the same accuracy of the existing HSDTs which have more number of variables.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.495-504
• /
• 2017

This paper proposes a new direct numerical integration algorithm for solving equation of motion in structural dynamics problems with nonlinear stiffness. The new implicit method's degree of accuracy is higher than that of existing methods due to the higher order of the acceleration. Two parameters are defined, leading to a new family of unconditionally stable methods, which helps to take greater time steps in integration and eliminate concerns about the duration of solving. The method developed can be utilized for a number of solid plane finite elements, examples of which are given to compare the proposed method with existing ones. The results indicate the superiority of the proposed method.

• JSTS:Journal of Semiconductor Technology and Science
• /
• v.9 no.4
• /
• pp.216-224
• /
• 2009

The endless scaling down of the semiconductor process made the impact of the power network noise on the performance of the state-of-the-art chip a serious design problem. This paper compares the performances of two popular circuit reduction approaches used to improve the efficiency of power network noise analysis: moment matching-based model order reduction (MOR) and node elimination-based MOR. As the benchmarks, we chose PRIMA and R2Power as the matching-based MOR and the node elimination-based MOR. Experimental results indicate that the accuracy, efficiency, and memory requirement of both methods very strongly depend on the structure of the given circuit, i.e., numbers of the nodes and sources, and the number of moments to preserve for PRIMA. PRIMA has higher accuracy in general, while the error of R2Power is also in the acceptable range. On the other hand, PRIMA has the higher efficiency than R2Power, only when the numbers of nodes and sources are small enough. Otherwise, R2Power clearly outperforms PRIMA in efficiency. In the memory requirement, the memory size of PRIMA increases very quickly as the numbers of nodes, sources, and preserved moments increase.