• Steel and Composite Structures
• /
• 제29권6호
• /
• pp.785-802
• /
• 2018

In this paper, two different computational methods, called Rayleigh-Ritz and collocation are developed to estimate the ultimate strength of composite plates. Progressive damage behavior of moderately thick composite laminated plates is studied under in-plane compressive load and uniform lateral pressure. The formulations of both methods are based on the concept of the principle of minimum potential energy. First order shear deformation theory and the assumption of large deflections are used to develop the equilibrium equations of laminated plates. Therefore, Newton-Raphson technique will be used to solve the obtained system of nonlinear algebraic equations. In Rayleigh-Ritz method, two degradation models called complete and region degradation models are used to estimate the degradation zone around the failure location. In the second method, a new energy based collocation technique is introduced in which the domain of the plate is discretized into the Legendre-Gauss-Lobatto points. In this new method, in addition to the two previous models, the new model named node degradation model will also be used in which the material properties of the area just around the failed node are reduced. To predict the failure location, Hashin failure criteria have been used and the corresponding material properties of the failed zone are reduced instantaneously. Approximation of the displacement fields is performed by suitable harmonic functions in the Rayleigh-Ritz method and by Legendre basis functions (LBFs) in the second method. Finally, the results will be calculated and discussions will be conducted on the methods.

• Structural Engineering and Mechanics
• /
• 제68권5호
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• 지구물리와물리탐사
• /
• 제17권4호
• /
• pp.231-241
• /
• 2014

최근 엔지니어링 목적의 단일채널 탄성파 탐사가 많이 수행되고 있다. 단일채널 탄성파 탐사는 일반적으로 특별한 자료처리 없이 효율적으로 지하 지질구조를 파악할 수 있는 장점이 있지만, 복잡한 지질구조에 대한 정확한 영상화에는 한계를 가진다. 자원개발 목적의 다중채널 탄성파 탐사에서는 최근 파형역산에 대한 연구가 활발히 진행되고 있으며, 복잡한 지하구조에 대해서도 정확한 지하영상화 결과를 제시하고 있다. 이에 본 논문에서는 단일채널 탄성파 탐사 자료를 이용하여 지하 속도모델을 구하기 위한 탄성파 파형역산 알고리즘을 제안하였다. 제안한 알고리즘은 단일채널 탄성파 탐사를 고려하여 지하 매질을 1차원으로 가정하였으며, 벌림에 의한 지연시간을 제거하여 벌림에 의한 효과를 보정하였다. 파형역산은 안정적인 해의 계산이 가능한 가우스-뉴턴법을 이용하였다. 알고리즘은 수정된 Marmousi2 모델에 적용하여 검증하였으며, 부산항에서 얻은 현장자료에 적용해 보았다.

• 한국지반공학회지:지반
• /
• 제8권2호
• /
• pp.5-20
• /
• 1992

The physical-based and lumped-parameter hydrologic groundwater flow model for predicting the rainfall-triggered rise of groundwater levels in hillside slopes is developed in this paper to assess the risk of landslides. The developed model consists of a vertical infiltration model for unsaturated zone linked to a linear storage reservoir model(LSRM) for saturated zone. The groundwater flow model has uncertain constants like soil depttL slope angle, saturated permeability, and potential evapotranspiration and four free model parameters like a, b, c, and K. The free model parameters could be estimated from known input-output records. The BARD algorithm is uses as the parameter estimation technique which is based on a linearization of the proposed model by Gauss -Newton method and Taylor series expansion. The application to examine the capacity of prediction shows that the developed model has a potential of use in forecast systems of predicting landslides and that the optimal estimate of potential 'a' in infiltration model is the most important in the global optimum analysis because small variation of it results in the large change of the objective function, the sum of squares of deviations of the observed and computed groundwater levels. 본 논문에서는 가파른 산사면에서 산사태의 발생을 예측하기 위한 수문학적 인 지하수 흐름 모델을 개발하였다. 이 모델은 물리적인 개념에 기본하였으며, Lumped-parameter를 이용하였다. 개발된 지하수 흐름 모델은 두 모델을 조합하여 구성되어 있으며, 비포화대 흐름을 위해서는 수정된 abcd 모델을, 포화대 흐름에 대해서는 시간 지체 효과를 고려할 수 있는 선형 저수지 모델을 이용하였다. 지하수 흐름 모델은 토층의 두께, 산사면의 경사각, 포화투수계수, 잠재 증발산 량과 같은 불확실한 상수들과 a, b, c, 그리고 K와 같은 자유모델변수들을 가진다. 자유모델변수들은 유입-유출 자료들로부터 평가할 수 있으며, 이를 위해서 본 논문에서는 Gauss-Newton 방법을 이용한 Bard 알고리즘을 사용하였다. 서울 구로구 시흥동 산사태 발생 지역의 산사면에 대하여 개발된 모델을 적용하여 예제 해석을 수행함으로써, 지하수 흐름 모델이 산사태 발생 예측을 위하여 이용할 수 있음을 입증하였다. 또한, 매개변수분석 연구를 통하여, 변수 a값은 작은 변화에 대하여 목적함수값에 큰 변화를 일으키므로 a의 값에 대한 최적값을 구하는 것이 가장 중요한 요소라는 결론을 얻었다.

• Steel and Composite Structures
• /
• 제19권1호
• /
• pp.191-208
• /
• 2015

The sensitivity-based finite element model updating method has received increasing attention in damage detection of structures based on measured modal parameters. Finding an optimization technique with high efficiency and fast convergence is one of the key issues for model updating-based damage detection. A new simple and computationally efficient optimization algorithm is proposed and applied to damage detection by using finite element model updating. The proposed method combines the Gauss-Newton method with region truncation of each iterative step, in which not only the constraints are introduced instead of penalty functions, but also the searching steps are restricted in a controlled region. The developed algorithm is illustrated by a numerically simulated 25-bar truss structure, and the results have been compared and verified with those obtained from the trust region method. In order to investigate the reliability of the proposed method in damage detection of structures, the influence of the uncertainties coming from measured modal parameters on the statistical characteristics of detection result is investigated by Monte-Carlo simulation, and the probability of damage detection is estimated using the probabilistic method.

• 대한기계학회논문집
• /
• 제18권10호
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• 한국공간구조학회논문집
• /
• 제21권4호
• /
• pp.49-56
• /
• 2021

This study aims to develop a form-finding algorithm for a single-layered pneumatic membrane. The initial shape of this pneumatic membrane, which is an air-supported type pneumatic membrane, is to find a state in which a given initial tension and internal pneumatic pressure are in equilibrium. The algorithm developed to satisfy these conditions is that a nonlinear optimization problem based on the force method considering the deformed shape is formulated, and, it's able to find the shape by iteratively repeating the process of obtaining a solution of the governing equations. An computational technique based on the Gauss-Newton method was used as a method for obtaining solutions of nonlinear equations. In order to verify the validity of the proposed form-finding algorithm, a single-curvature pneumatic membrane example and a double-curvature air pneumatic membrane example were adopted, respectively. In the results of these examples, it was possible to well observe the step-by-step convergence process of the shape of the pneumatic membrane, and it was also possible to confirm the change in shape according to the air pressure. In addition, the calculation results of the shape and internal force after deformation due to initial tension, air pressure, and self-weight were obtained.

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2003년도 Proceedings of the international symposium on the fusion technology
• /
• pp.390-397
• /
• 2003

A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.

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

• 한국지구물리탐사학회:학술대회논문집
• /
• 한국지구물리탐사학회 2008년도 공동학술대회
• /
• pp.51-56
• /
• 2008

등방성 매질에서의 파형역산에 대한 연구는 1980년대부터 꾸준히 이루어져 왔으나 이방성 매질에 대한 연구는 그렇지 못하다. 본 연구에서는 이방성 매질에 대한 시간영역 셀기반 유한 차분 모델링 기법을 이용해 2차원 TI 구조에서의 파형역산 알고리듬을 개발하였다. 반복적인 비선형 역산에서 최대 급경사 방향은 역시간 구조보정의 역전파 방법을 이용하여 간접적으로 계산하였고, 이를 정규화 시키기 위해 슈도-헤시안 행렬을 이용하였다. 본 연구에서 제시된 시간영역 파형역산 기법을 이방성 매질을 포함한 2층 구조와 이방성 Marmousi 모형 자료에 적용하고 이를 등방성 매질만을 고려한 기존의 파형역산 결과와 비교하였다. 본 연구의 결과를 통해 이방성 매질을 등방성 매질로 가정하고 파형역산을 수행할 경우 정확한 영상을 얻을 수 없기 때문에, 실제 탐사 자료의 파형역산을 수행할 경우 이방성 매질을 고려해야 좀 더 정확한 지하 구조를 파악할 수 있음을 확인하였다.