• Structural Engineering and Mechanics
• /
• v.45 no.4
• /
• pp.521-542
• /
• 2013

In this study a truss model is used for the geometrically nonlinear static and dynamic analysis of a thin shallow arch subject to snap-through. Thanks to the very simple geometry of a truss, the equilibrium conditions can be easily written and the global stiffness matrix can be easily updated with respect to the deformed structure, within each step of the analysis. A very coarse discretization is applied; so, in a very simple way, the high frequency modes are suppressed from the beginning and there is no need to develop a complicated reduced-order technique. Two short computer programs have been developed for the geometrically nonlinear static analysis by displacement control of a plane truss model of a structure as well as for its dynamic analysis by the step-by-step time integration algorithm of trapezoidal rule, combined with a predictor-corrector technique. These two short, fully documented computer programs are applied on the geometrically nonlinear static and dynamic analysis of a specific thin shallow arch subject to snap-through.

• Structural Engineering and Mechanics
• /
• v.5 no.5
• /
• pp.491-505
• /
• 1997

This paper deals with the formulation and solution of a wide class of structures, in the presence of both geometric and material nonlinearities, as a particular mathematical programming problem. We first present key ideas for the nonholonomic (path dependent) rate formulation for a suitably discretized structural model before we develop its computationally advantageous stepwise holonomic (path independent) counterpart. A feature of the final mathematical programming problem, known as a nonlinear complementarity problem, is that the governing relations exhibit symmetry as a result of the introduction of so-called nonlinear "residuals". One advantage of this form is that it facilitates application of a particular iterative algorithm, in essence a predictor-corrector method, for the solution process. As an illustrative example, we specifically consider the simplest case of plane trusses and detail in particular the general methodology for establishing the static-kinematic relations in a dual format. Extension to other skeletal structures is conceptually transparent. Some numerical examples are presented to illustrate applicability of the procedure.

• Journal of Institute of Control, Robotics and Systems
• /
• v.12 no.7
• /
• pp.671-678
• /
• 2006

This article proposes a smooth path generation and time scheduling method for general tasks defined by non-smooth path segments in industrial robotic applications. This method utilizes a simple 3rd order polynomial function for smooth interpolation between non-smooth path segments, so that entire task can effectively maintain constant line speed of operation. A predictor-corrector type numerical mapping technique, which correlates time based speed profile to the smoothed path in Cartesian space, is also provided. Finally simulation results show the feasibility of the proposed algorithm.

• Computers and Concrete
• /
• v.4 no.1
• /
• pp.67-82
• /
• 2007

The discrete crack-concept is applied to study the 3D propagation of tensile-dominated failure in plain concrete. To this end the Partition of Unity Finite Element Method (PUFEM) is utilized and the strong discontinuity approach is followed. A consistent linearized implementation of the PUFEM is combined with a predictor-corrector algorithm to track the crack path, which leads to a robust numerical description of concrete cracking. The proposed concept is applied to study concrete failure during the PCT3D test and the predicted numerical results are compared to experimental data. The proposed numerical concept provides a clear interface for constitutive models and allows an investigation of their impact on concrete cracking under 3D conditions, which is of significant scientific interests to interpret results from 3D experiments.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.2 no.4 s.7
• /
• pp.97-116
• /
• 2002

In this papers, a finite element formulation is proposed for dynamic analysis of vehicle-bridge interaction problems under realistic loading conditions. Although the formulation presented in this paper is based on the consideration of only a single traversing vehicle, it can be extended to include several different bridge configurations. The traversing vehicle and the vibrating bridge superstructure are considered as an integrated system. Hence, although material and geometric nonlinearities are excluded, this introduces nonlinearity into the problem. Various vehicle models, including those with suspension systems, are considered. Traveling speed of the vehicle can be varied. The finite element discretization of the bridge structure permits the inclusion of arbitrary geometrical configurations, and surface and boundary conditions. To obtain accurate solutions, time integration of the equation of vehicle-bridge motion is carried out by using the Newmark method in connection with a predictor-corrector algorithm.

• Journal of Korean Society of Water and Wastewater
• /
• v.13 no.3
• /
• pp.42-55
• /
• 1999

A mathematical model for biological nutrient removal in a sequencing batch reactor process, which is based on the IAWQ Activated Sludge Model No. 2 with a few modifications, has been developed. Twenty water quality components and twenty three kinetic equations are incorporated in the model. The model is structured in the matrix form based on the law of mass conservation using stoichiometry and kinetic equations. Stoichiometric coefficients and kinetic parameters included in the model equations are chosen from the literature. A multistep predictor-corrector algorithm of variable step-size is adopted for solving the vector nonlinear ordinary differential equations. The simulation for experimental results is conducted to evaluate the validity of the model and to calibrate coefficients and parameters. The simulation using the model well represents the experimental results from laboratory. The mathematical model developed in this study may be utilized for the design and operation of a sequencing batch reactor process under the steady and unsteady-state at various environmental conditions.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.3
• /
• pp.116-121
• /
• 2005

In the present paper, a new calculation algorithm far solving large scale environmental or geophysical flows with free surface is proposed where the non-hydrostatic pressure component is taken into consideration. Predictor-corrector fractional step approach with explicit, forward time marching scheme in the sigma coordinate system is employed. In order to validate the present calculation algorithm and to estimate the effects of non-hydrostatic pressure on resultant flow and free surface movements, example calculations are carried out for typical steady and unsteady flow problems. Present method can be applied to the meso-scale free surface flows with complex bottom topography where MAC-like 3-d hydrodynamic calculations are quite ineffective and uneconomic.

• Structural Engineering and Mechanics
• /
• v.60 no.5
• /
• pp.891-902
• /
• 2016

The beam structure with breathing cracks subjected to harmonic excitations was modeled by FEM based on Euler-Bernoulli theory, and a piecewise dynamical system was deduced. The precise integration method (PIM) was employed to propose an algorithm for analyzing the dynamic responses of the deduced system. This system was first divided into linear sub-systems, between which there are switching points resulted from the breathing cracks. The inhomogeneous terms due to the external excitations were tackled by introducing auxiliary variables to express the harmonic functions, hence the sub-systems are homogeneous. The PIM was then applied to solve the homogeneous sub-systems one by one. During the procedures, a predictor-corrector algorithm was presented to determine the switching points accurately. The presented method can provide solutions with an accuracy to a magnitude of $10^{-12}$ compared with exact solutions obtained by the theories of ordinary differential equations. The PIM results are much more accurate than Newmark ones with the same time step. Moreover, it is found that the PIM can maintain a high level of accuracy even when the time step increases within a relatively wide range.

• Journal of Korean Society of Steel Construction
• /
• v.23 no.5
• /
• pp.535-546
• /
• 2011

In this paper, semi-rigid elasto-plastic post-buckling analysis of a space frame was performed using various explicit arc-length methods. Various explicit arc-length methodsand a large-deformation and small-strain elasto-plastic 3D space frame element with semi-rigid connections and plastic hinges were developed. This element can be appliedto both explicit and implicit numerical algorithms. In this study, the Dynamic Relaxation method was adopted in the predictor and corrector processesto formulate an explicit arc-length algorithm. The developed "explicit-predictor" or "explicit-corrector" were used in the elasto-plastic post-buckling analysis. The Eulerian equations for a beam-column with finite rotation, which considers the bowing effects, were adopted for the elastic system and extended to theinelastic system with a plastic hinge concept. The derived tangent stiffness matrix was asymmetrical due to the finite rotation. The joint connection elements were introduced for semi-rigidity using a static condensation technique. Semi-rigid elasto-plastic post-buckling analyses were carried out to demonstrate the potential of the developed explicit arc-length method and advanced space frame element in terms of accuracy and efficiency.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.