• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.1
• /
• pp.77-85
• /
• 1988

A stochchastic analysis procedure of generating floor response spectra for proportionally damped linear systems subject to earthquake loading is presented. Theories of random vibration and mode acceleration method are used in the formulation of governing equations. The structure-oscillator interaction is not considered. It is assumed that the input motions and oscillator responses are stationary Gaussian processes with mean zero. The nonstationary characteristics of earthquake motion are incorporated in the peak factor which is based on Vanmarcke's theory. Floor response spectra for both resonance and non-resonance cases are calculated under the assumption that the peak factors for structure and oscillator are equal to that for ground response spectrum. The validity of this method is demonstrated by comparing the results obtained by proposed method with those by time history analyses. The results obtained by this method are conservative and accurate with tolerable precision. This method saves much computing time compared with time history analysis method.

• International Journal of Safety
• /
• v.1 no.1
• /
• pp.36-42
• /
• 2002

This paper is aimed at predicting the life of rubber insulating gloves under normal operating stresses from relatively rapid test performed at higher stresses. Specimens of rubber insulating gloves are subject to multiple stress conditions, i.e. combined electrical and thermal stresses. Two modes of electrical stress, step voltage stress and constant voltage stress are used in specimen aging. There are two types of test for electrical stress in this experiment: the one is Breakdown Voltage (BDV) test under step voltage stress and thermal stress and the other is lifetime test under constant voltage stress and temperature stress. The ac breakdown voltage defined as the break-down point of insulation that leakage current excesses a limit value, l0mA in this experiment, is determined. Because the very high variability of aging data requires the application of statistical model, Weibull distribution is used to represent the failure times as the straight line on Weibull probability paper. Weibull parameters are deter-mined by three statistical methods i.e. maximum likelihood method, graphical method and least squares method, which employ SAS package, Weibull probability paper and FORTRAN, respectively. Two chosen models for predicting the life under simultaneous electrical and thermal stresses are inverse power model and exponential model. And the constants of life equation for multistress aging are calculated using numerical method, such as Gauss Jordan method etc.. The completion of life equation enables to estimate the life at normal stress based on the data collected from accelerated aging test. Also the comparison of the calculated lifetimes between the inverse power model and the exponential model is carried out. And the lifetimes calculated by three statistical methods with lower voltage than test voltage are compared. The results obtained from the suggested experimental method are presented and discussed.

• Structural Engineering and Mechanics
• /
• v.34 no.6
• /
• pp.779-799
• /
• 2010

In order to consider high-order effects on the actual limit state function, a new response surface method is proposed for structural reliability analysis by the use of high-order approximation concept in this study. Hermite polynomials are used to determine the highest orders of input random variables, and the sampling points for the determination of highest orders are located on Gaussian points of Gauss-Hermite integration. The cross terms between two random variables, only in case that their corresponding percent contributions to the total variation of limit state function are significant, will be added to the response surface function to improve the approximation accuracy. As a result, significant reduction in computational cost is achieved with this strategy. Due to the addition of cross terms, the additional sampling points, laid on two-dimensional Gaussian points off axis on the plane of two significant variables, are required to determine the coefficients of the approximated limit state function. All available sampling points are employed to construct the final response surface function. Then, Monte Carlo Simulation is carried out on the final approximation response surface function to estimate the failure probability. Due to the use of high order polynomial, the proposed method is more accurate than the traditional second-order or linear response surface method. It also provides much more efficient solutions than the available high-order response surface method with less loss in accuracy. The efficiency and the accuracy of the proposed method compared with those of various response surface methods available are illustrated by five numerical examples.

• Journal of the Korean earth science society
• /
• v.37 no.5
• /
• pp.277-285
• /
• 2016

We introduce a depth scaling strategy to improve the accuracy of frequency-domain elastic full waveform inversion (FWI) using the new pseudo-Hessian matrix for seismic data without low-frequency components. The depth scaling strategy is based on the fact that the damping factor in the Levenberg-Marquardt method controls the energy concentration in the gradient. In other words, a large damping factor makes the Levenberg-Marquardt method similar to the steepest-descent method, by which shallow structures are mainly recovered. With a small damping factor, the Levenberg-Marquardt method becomes similar to the Gauss-Newton methods by which we can resolve deep structures as well as shallow structures. In our depth scaling strategy, a large damping factor is used in the early stage and then decreases automatically with the trend of error as the iteration goes on. With the depth scaling strategy, we can gradually move the parameter-searching region from shallow to deep parts. This flexible damping factor plays a role in retarding the model parameter update for shallow parts and mainly inverting deeper parts in the later stage of inversion. By doing so, we can improve deep parts in inversion results. The depth scaling strategy is applied to synthetic data without lowfrequency components for a modified version of the SEG/EAGE overthrust model. Numerical examples show that the flexible damping factor yields better results than the constant damping factor when reliable low-frequency components are missing.

• Journal of Astronomy and Space Sciences
• /
• v.5 no.1
• /
• pp.31-43
• /
• 1988

The differential correction process determining osculating orbital elements as correct as possible at a given instant of time from tracking data of artificial satellite was accomplished. Preliminary orbital elements were used as an initial value of the differential correction procedure and iterated until the residual of real observation (O) and computed observation(C) was minimized. Tracking satellite was NOAA-9 or TIROS-N series. Two types of tracking data were prediction data precomputed from mean orbital elements of TBUS and real data obtained from tracking 1.70 GHz HRPT signal of NOAA-9 using 5 meter auto-track antenna in Radio Research Laboratory. Accrding to thacking data either Gause method or Herrick-Gibbs method was applied to preliminary orbit determination. In the differential correction stage we used both of the Escobal(1975)'s analytical method and numerical method using f, g series for the comparision. The results between analytical and numerical ones are nearly consistent. And the differentially corrected orbit converged to the same value in spite of the differences between preliminary orbits of each time span.

• Journal of the Institute of Electronics and Information Engineers
• /
• v.54 no.8
• /
• pp.99-106
• /
• 2017

Electrical impedance tomography is a nondestructive imaging modality in which the internal resistivity distribution is reconstructed based on the injected currents and measured voltages inside a domain of interest. In this paper, an adaptive threshold value based region of interest (ROI) method is proposed to improve the spatial resolution of reconstructed images as well as to reduce the computational time of the inverse problem. Adaptive threshold value is calculated by INTERMODES method and ROI is determined from the domain based on this value. Moreover, the computational domain of image reconstruction is restricted within a ROI and iterative Gauss-Newton method is employed to estimate the resistivity distribution. To evaluate the performance of the proposed method, numerical experiments have been performed and the results are analyzed.

• Economic and Environmental Geology
• /
• v.26 no.2
• /
• pp.207-215
• /
• 1993

In modem petroleum reservoir engineering, the characterization of reservoir heterogeneities is very important to accurately understand and predict reservoir production performance. Formation evaluation for the description of reservoir is generally conducted by performing the analysis of well logging, core testing, and well testing. However, the measured data points by well logging or core testing are in general very sparse and hence reservoir properties should be interpolated and extrapolated from measured points to uncharacterized areas. In assigning the data for the unknown points, simple averaging technique is not feasible as optimum estimation method since this method does not account the spatial relationship between the data points. The main goal of this work is to develop PC-version of multi-purpose geostatistical model in which several stages are systematically proceeded. In the development of model, the simulator employs a automatic selection of semivariogram function such as exponential or spherical model with the best values of $R^2$. The simulator also implements a special algorithm for the fitting of semivariogram function to experimental sernivariogram. The special algorithm such as trial and error scheme is devised since this method is much more reliable and stable than Gauss-Newton method. The simulator has been tested under stringent conditions and found to be stable. Finally, the validity and the applicability of the developed model have been studied against some existing actual field data.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.69-76
• /
• 2003

Low-tension cables have been increasingly used in recent years due to deep-sea developments and the advent of synthetic cables. In the case of low-tension cables, large displacements may happen due to relatively small restoring forces of tension and thus the effects of fluid and geometric non-linearities and bending stiffness. A Fortran program is developed by employing a finite difference method. In the algorithm, an implicit time integration and Newton-Raphson iteration are adopted. For the calculation of huge size of matrices, block tri-diagonal matrix method is applied, which is much faster than the well-known Gauss-Jordan method in two point boundary value problems. Some case studies are carried out and the results of numerical simulations are compared with a in-house program of WHOI Cable with good agreements.

• Steel and Composite Structures
• /
• v.9 no.2
• /
• pp.131-158
• /
• 2009

This paper presents a modelling technique for the nonlinear analysis of composite steel-concrete beams with partial shear interaction. It extends the applicability of two stiffness elements previously derived by the authors using the direct stiffness method, i.e. the 6DOF and the 8DOF elements, to account for material nonlinearities. The freedoms are the vertical displacement, the rotation and the slip at both ends for the 6DOF stiffness element, as well as the axial displacement at the level of the reference axis for the 8DOF stiffness element. The solution iterative scheme is based on the secant method, with the convergence criteria relying on the ratios of the Euclidean norms of both forces and displacements. The advantage of the approach is that the displacement and force fields of the stiffness elements are extremely rich as they correspond to those required by the analytical solution of the elastic partial interaction problem, thereby producing a robust numerical technique. Experimental results available in the literature are used to validate the finite element proposed in the paper. For this purpose, those reported by Chapman and Balakrishnan (1964), Fabbrocino et al. (1998, 1999) and Ansourian (1981) are utilised; these consist of six simply supported beams with a point load applied at mid-span inducing positive bending moment in the beams, three simply supported beams with a point load applied at mid-span inducing negative bending moment in the beams, and six two-span continuous composite beams respectively. Based on these comparisons, a preferred degree of discretisation suitable for the proposed modelling technique expressed as a function of the ratio between the element length and depth is proposed, as is the number of Gauss stations needed. This allows for accurate prediction of the nonlinear response of composite beams.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.