• 대한토목학회논문집
• /
• 제8권1호
• /
• pp.77-85
• /
• 1988

본 연구(硏究)에서는 지진(地震)하중을 받는 선형비례감쇠(線形比例減衰) 시스템의 층응답(層應答) 스펙트럼을 random vibration 이론(理論)을 적용하여 계산(計算)하는 방법을 제시(提示)하였다. 해석(解析)방법으로는 모드가속도법(加速度法)을 사용하였으며 구조물(構造物)-기기(機器)의 상호작용(相互作用)은 고려하지 않았다. 입력지진운동(入力地震運動)과 기기(機器)의 응답(應答)을 평균(平均)값이 영(零)인 정상(定常) Gauss 과정(過程)으로 가정하였다. 입력지진(入力地震)의 천이특성을 Vanmarcke 방법(方法)에 따라 첨두계수(尖頭係數) 계산시 고려하였다. 층응답(層應答) 스펙트럼을 공진(共振)과 비공진(非共振)으로 구분(區分)하여 계산하였으며 응답(應答)계산시 구조물(構造物)과 진동체(振動體)의 첨두계수(尖頭係數)는 지반응답(地盤應答) 스펙트럼의 첨두계수(尖頭係數)와 동일(同一)하다고 가정하였다. 적용예(適用例)에서는 시간이력해석(時間履歷解析)의 결과와 비교(比較)함으로써 본 연구(硏究)의 타당성(妥當性)을 입증(立證)하였다. 본(本) 논문(論文)의 해석방법(解析方法)을 사용(使用)하면 비교적(比較的) 정확(正確)한 안전측(安全側)의 결과(結果)를 얻을 수 있으며 시간이력해석법(時間履歷解析法)에 비해 계산시간(計算時間)을 상당(相當)히 절약할 수 있다.

• International Journal of Safety
• /
• 제1권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
• /
• 제34권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.

• 한국지구과학회지
• /
• 제37권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
• /
• 제5권1호
• /
• pp.31-43
• /
• 1988

인공위성의 추적자료로부터 예비궤도를 결정하고 이것을 초기값으로 해서 계산한 관측치 (C) 와 실제의 관측치(O)의 차이가 최소가 될때까지 반복계산함으로써, 어떤 시각에 있어서 실제의 궤도에 가장 가까운 접융궤도 요소를 얻어내는 미분보정의 과정을 수행했다. 추적한 인공 위성은 TIROS-N계열 위성인 NOAA - 9이고 추척자료는 컴퓨터에 의한 프로그랭 추척을 위해서 TBUS의 평균궤도요소로부터 미리 계산한 예측자료와 위성으후부터 송신되는 1.707GHz의 HRPT신호를 전파연구소의 직경 5m인 자동추적안테나로 추적해서 얻은 실제 자료를 이용했다. 예비궤도의 결정은 추적 자료에 따라 Gauss 방법과 Herrick -Gibbs 방법을 사용했고 미분보정의 과정에 있어서 미분보정행렬의 각 원소들은 Escobal(1975)의 해석적인 방법과 f, g급수를 이용한 수치적인 방법으로 각각 구해서 그 결과를 비교해 보았다. 해석적인 방법파 수치적인 방법의 결과는 거의 일치하였고 각 시간간격에 따른 예비궤도가 다름에도 불구하고 미분보정된 궤도는 같은 값으로 수령했다.

• 전자공학회논문지
• /
• 제54권8호
• /
• pp.99-106
• /
• 2017

전기 임피던스 단층촬영법은 주입 전류와 측정 전압을 기반으로 관심 도메인 내부의 도전율/저항률 분포를 복원하는 비파괴 영상 복원 기법이다. 본 논문에서는 역문제 계산시간을 줄이고 더불어 공간 해상도도 향상시키기 위해, 적응 문턱치 기반의 ROI(region of interest) 방법을 제안하였다. INTERMODES 방법에 의해 적응 문턱치가 계산이 되고 이 값을 기반으로 전체 도메인으로부터 ROI가 결정된다. 그리고 영상 복원의 계산 도메인을 ROI 내로 국한시켜 반복적 가우스-뉴턴 방법을 적용하여 저항률 분포를 추정하였다. 제안한 방법의 성능을 평가하기 위해 수치실험을 수행하고 그 결과를 비교분석하였다.

• 자원환경지질
• /
• 제26권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.

• 한국전산구조공학회논문집
• /
• 제16권1호
• /
• pp.69-76
• /
• 2003

최근 들어 심해역에 대한 개발과 합섬섬유 재질 케이블의 발달로 인하여 저장력 케이블의 사용이 증가되었다. 저장력 케이블은 장력에 의한 복원력이 작기 때문에 대변위가 발생하게 되며, 따라서 기하학적 비선형이 강하게 나타나게 된다. 또한 해양환경에서는 유체 비선형도 작용하게 된다. 본 연구에서는 수치해석적 방법을 통하여 불균일하게 구성된 예인되는 저장력 케이블의 3차원 동적거동 해석을 수행한다. 수치해석에서는 유체 및 기하학적 비선형과 굽힘강성이 고려되며, 유한차분법(음해법)을 적용하여 풀이된다. 비선형 해를 구하기 위해서 뉴톤-랍슨 방법을 사용한다. 대형 행렬을 풀이하기 위하여 불록삼중대각행렬 풀이법이 적용되는데, 이 방법은 일반적인 행렬 풀이법인 가우스-조르단 방법에 비하여 계산시간을 상당히 줄일 수 있었다. 선배열 음탐 케이블에 대한 다양한 예제해석을 수행하였으며, 해석결과는 미국 우즈홀 해양연구소에서 개발된 프로그램 결과와 잘 일치하였다.

• Steel and Composite Structures
• /
• 제9권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
• /
• 제78권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.