• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.87-94
• /
• 1997

The circular plate resting on Boussinesq's half-space model under axisymmetric loading is studied by a finite element procedure to evaluate the distribution of contact pressure between plate and elastic half-space. The displacement of half-space due to axisymmetric surface loading can be evaluated by double integration of Boussinesq's solution. On that case the analytical integration can be executed for the radial direction but the analytical integration for the circumferential direction is impossible and the numerical integration should be considered. With the radial integration we can get non-dimensional function. Then the numerical integration for the formula is executed for the circumferential direction and the results are approximated 5th order Polynomials by using the least square method. With these 5th order approximate formula, the flexibility matrix of half-space is constructed as the coefficient matrix of nodal contact pressure by the finite element procedures. Iteration procedures are attempted by using this method to determine the separated region.

• Earthquakes and Structures
• /
• v.8 no.6
• /
• pp.1325-1348
• /
• 2015

Real-time hybrid testing (RTHT) involves virtual splitting of the structure into two parts: physical substructure that contains the key region of interest which is tested in a laboratory and numerical substructure that contains the remaining part of the structure in the form of a numerical model. This paper numerically assesses four step-by-step integration methods (Central difference method (CDM), Operator splitting method (OSM), Rosenbrock based method (RBM) and CR-integration method (CR)) which are widely used in RTHT. The methods have been assessed in terms of stability and accuracy for various realistic damping ratios of the physical substructure. The stability is assessed in terms of the spectral radii of the amplification matrix while the accuracy in terms of numerical damping and period distortion. In order to evaluate the performance of the methods, five carefully chosen examples have been studied - undamped SDOF, damped SDOF, instantaneous softening, instantaneous hardening and hysteretic system. The performance of the methods is measured in terms of a non-dimensional error index for displacement and velocity. Based on the error indices, it is observed that OSM and RBM are robust and performs fairly well in all the cases. CDM performed well for undamped SDOF system. CR method can be used for the system showing softening behaviour. The error indices indicate that accuracy of OSM is more than other method in case of hysteretic system. The accuracy of the results obtained through time integration methods for different damping ratios of the physical substructure is addressed in the present study. In the presence of a number of integration methods, it is preferable to have criteria for the selection of the time integration scheme. As such criteria are not available presently, this paper attempts to fill this gap by numerically assessing the four commonly used step-by-step methods.

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.

• Journal of computational fluids engineering
• /
• v.9 no.2
• /
• pp.43-51
• /
• 2004

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

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

Interaction between fluid flow and thermal radiation has received considerable attention due to its numerous applications in engineering field. In this case the thermofluid properties of radiating fluid vary with the variation of temperature field caused by absorption and emission of radiant heat. To analyze the radiation heat transfer in radiating fluid, the simultaneous solution of the radiative transfer equation (RTE) and the fluid dynamics equations is required. This means that the numerical procedure used for the RTE must be computationally efficient to permit its inclusion in the other submodels, and must be compatible with the other transport equations. The finite volume method (FVM) and the discrete ordinates method (DOM) are usually employed to simulate radiation problems in generalized coordinates. These two representative methods are examined and compared, especially in view of the numerical integration of the radiation intensity over solid angle. The FVM shows better accuracy than the DOM owing to less constraints of the selection of control angle.

• Structural Engineering and Mechanics
• /
• v.9 no.6
• /
• pp.557-568
• /
• 2000

As time-variant reliability approaches become increasingly used for service life prediction of the aging infrastructure, the demand for computer solution methods continues to increase. Effcient computer techniques have become well established for the reliability analysis of structural systems. Thus far, however, this is largely limited to time-invariant reliability problems. Therefore, the requirements for time-variant reliability prediction of deteriorating structural systems under time-variant loads have remained incomplete. This study presents a computer program for $\underline{REL}$iability of $\underline{T}$ime-Variant $\underline{SYS}$tems, RELTSYS. This program uses a combined technique of adaptive importance sampling, numerical integration, and fault tree analysis to compute time-variant reliabilities of individual components and systems. Time-invariant quantities are generated using Monte Carlo simulation, whereas time-variant quantities are evaluated using numerical integration. Load distribution and post-failure redistribution are considered using fault tree analysis. The strengths and limitations of RELTSYS are presented via a numerical example.

• Structural Engineering and Mechanics
• /
• v.38 no.3
• /
• pp.283-300
• /
• 2011

A new conception of fundamental tasks in dynamics of the bridge-track-train systems (BTT), with the aim to evaluate moving load's models adequacy, has been developed. The 2D physical models of BTT systems, corresponding to the fundamental tasks, have been worked out taking into account one-way constraints between the moving unsprung masses and the track. A method for deriving the implicit equations of motion, governing vibrations of BTT systems' models, as well as algorithms for numerical integration of these equations, leading to the solutions of high accuracy and relatively short times of simulations, have been also developed. The derived equations and formulated algorithms constitute the basis for numerical simulation of vibrations of the considered systems.

• Journal of applied mathematics & informatics
• /
• v.24 no.1_2
• /
• pp.65-79
• /
• 2007

An optimal 3-point quadrature formula of closed type is derived. It is shown that the optimal quadrature formula has a better error bound than the well-known Simpson's rule. A corrected formula is also considered. Various error inequalities for these formulas are established. Applications in numerical integration are given.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.2
• /
• pp.117-136
• /
• 2001

A generalization of Ostrowski integral inequality for mappings whose derivatives are bounded and applications for general quadrature formulae are given.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.52-53
• /
• 2003