• Journal of the Korea Concrete Institute
• /
• v.16 no.6 s.84
• /
• pp.841-850
• /
• 2004

The modeling of uplift pressure within dam, on the foundation on which it was constructed, and on the interface between the dam and foundation is a critical aspect in the analysis of concrete gravity dams, i.e. crack stability in concrete dam can correctly be predicted when uplift pressures are accurately modelled. Current models consider a uniform uplift distribution, but recent experimental results show that it varies along the crack faces and the procedures for modeling uplift pressures are well established for the traditional hand-calculation methods, but this is not the case for finite element (FE) analysis. In large structures, such as dams, because of smaller size of the fracture process zone with respect to the structure size, limited errors should occur under the assumptions of linear elastic fracture mechanics (LEFM). In this paper, the fracture behaviour of concrete gravity dams mainly subjected to uplift Pressure at the crack face was studied. Triangular type, trapezoidal type and parabolic type distribution of the uplift pressure including uniform type were considered in case of evaluating stress intensity factor by surface integral method. The effects of body forces, overtopping pressures are also considered and a parametric study of gravity dams under the assumption of LEFM is performed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.588-593
• /
• 2006

자연하천은 일반적으로 만곡수로나 사행수로 형태를 보이고 있으며, 직선수로에서와 달리 원심력에 기인한 이차류 영향을 받게 된다. 이차류에 의해서 수면에서는 만곡부 바깥쪽으로, 하상에서는 만곡부 안쪽으로의 흐름특성을 보이게 된다. 만곡부 안쪽으로 가해지는 하상 전단응력에 기인하여 하상에서의 입자가 만곡부 안쪽으로 이송되며, 결과적으로 만곡부 안쪽에는 점사주가, 바깥쪽에는 소(pool)가 생성된다. 또한 지형경사의 생성으로 입자에 가해지는 중력효과도 변화된다. 따라서 이와 같은 자연하천의 흐름과 하상변동을 수치모의 하기 위해서는 만곡부 이차류 특성을 고려한 모형이 필요하다. 본 연구에서는 수심 적분된 흐름방정식과 하상토 보존방정식 (Exner equation)을 이용한 하상변동을 위한 비연계 수치모형을 위해서 하상토 보존방정식의 유한요소 알고리즘을 개발하였다. 하상토 보존방정식은 흐름 특성에 따른 평형 유사량의 공간변화율을 이용하여 일정 기간 동안의 하상 변화량을 계산한다. 이 때 이차류에 의한 하상 전단응력의 편각 및 지형경사 변화에 따른 실제 입자의 이송방향을 보정하여 평형 유사량이 계산된다. 이러한 보정식을 적용시키기 위해서는 유속성분의 공간변화량 및 지형경사의 공간성분이 필요하다. 유한요소법은 연속성 변수를 이산화시켜 근사해를 구하는 수치기법의 일종이기 때문에, 요소망이 불규칙적으로 구성되었을 경우 임의의 절점에서 연속성을 지닌 변수의 공간변화율을 계산하는데 어려움이 있다. 따라서 본 연구에서는 평형 유사량 계산 시에 절점이 아닌 요소 내부에서 평형 유사량을 계산하는, 하상토 보존방정식의 새로운 유한요소 알고리즘을 개발하고, 새로운 알고리즘을 적용시킨 수치모형의 검증을 행하였다. 경계조건 알고리즘의 검증으로 위해서 Soni 등 (1980)이 행한 상류 유입 유사량에 따른 하상변동을 수치 모의하고 실험치와 비교하였으며, Sutmuller와 Glerum (1980)이 수행한 만곡수로에서의 하상변동을 모의하고 실험과 비교하였다. 새로운 알고리즘을 적용시킨 하상토 보존방정식의 유한요소 수치모형의 결과는 매우 안정적이며, 실험과 매우 유사한 결과를 얻을 수 있었다. 본 수치모델은 현재 균일한 입자의 하상토만을 고려하므로, 입자분급이나 하상 장갑화 현상 등은 무시한다.

• Journal of Astronomy and Space Sciences
• /
• v.14 no.1
• /
• pp.142-149
• /
• 1997

This paper applied one of the well-known optimal control theory, namely, linear quadratic regulator(LQR), to the station-keeping maneuvers(SKM) for a geostationary satellite. The boundary conditions to transfer the system with a good accuracy at a terminal time were based upon the predicted orbital data which are created due to the Earth's non-uniform mass distribution's effect during 14 days and due to luni-solar effect during 28 days. Through the linearization of the nonlinear system equation with respect to reference orbit and the numerical integration of Riccati equation, the optimal trajectories and the corresponding control law have been obtained by using LQR. From the comparison of ${\Delta}V$ obtained by LQR with the ${\Delta}V$ obtained anatically by geometric method, Station Keeping Maneuvers(SKM) via LQR may provide comparable results to a real system. Furthermore it will demonstrate the possibility in fuel optimization and life extension of geostationary satellite.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.14 no.3
• /
• pp.11-22
• /
• 2010

An analysis method is proposed for the transient linear or nonlinear analysis of dynamic interactions between a flexible dam body and reservoir impounding compressible water under earthquake loadings. The coupled dam-reservoir system consists of three substructures: (1) a dam body with linear or nonlinear behavior; (2) a semi-infinite fluid region with constant depth; and (3) an irregular fluid region between the dam body and far field. The dam body is modeled with linear and/or nonlinear finite elements. The far field is formulated as a displacement-based transmitting boundary in the frequency domain that can radiate energy into infinity. Then the transmitting boundary is transformed for the direct coupling in the time domain. The near field region is modeled as a compressible fluid contained between two substructures. The developed method is verified and applied to various earthquake response analyses of dam-reservoir systems. Also, the method is applied to a nonlinear analysis of a concrete gravity dam. The results show the location and severity of damage demonstrating the applicability to the seismic evaluation of existing and new dams.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.26 no.2
• /
• pp.25-31
• /
• 1989

Two-dimensional large-amplitude motions in regular harmonic wave are treated in time domain, by satisfying the exact body boundary condition and the linear free surface condition. For the present numerical calculation, the method of free-surface spectral representation with simple source distribution on the instantaneous body surface has been extended to include the effect of the incident wave. Calculations of the wave exciting force are performed for a submerged circular cylinder fixed or oscillating with large amplitude. Especially, nonlinear effects on the time-mean forces are studied in detail. It is shown that relative motion between the body and the fluid particle gives a significant effect on the lift and drift forces. Also, large-amplitude motion of a submerged circular cylinder and that of a floating Lewis-form cylinder are directly simulated in time domain. In the calculation results, some nonlinear effects are shown.

• Journal of Radiation Protection and Research
• /
• v.26 no.2
• /
• pp.59-66
• /
• 2001

It is a hot issue to determine the spatial location and shape of tumor boundary in fractionated stereotactic radiotherapy (FSRT). We could get consecutive transaxial plane images from the phantom (paraffin) and 4 patients with brain tumor using helical computed tomography(HCT). K-means classification algorithm was adjusted to change raw data pixel value in CT images into classified average pixel value. The classified images consists of 5 regions that ate tumor region (TR), normal region (NR), combination region (CR), uncommitted region (UR) and artifact region (AR). The major concern was how to separate the normal region from tumor region in the combination area. Relative average deviation analysis was adjusted to alter average pixel values of 5 regions into 2 regions of normal and tumor region to define maximum point among average deviation pixel values. And then we drawn gross tumor volume (GTV) boundary by connecting maximum points in images using semi-automatic contour method by IDL(Interactive Data Language) program. The error limit of the ROI boundary in homogeneous phantom is estimated within ${\pm}1%$. In case of 4 patients, we could confirm that the tumor lesions described by physician and the lesions described automatically by the K-mean classification algorithm and relative average deviation analyses were similar. These methods can make uncertain boundary between normal and tumor region into clear boundary. Therefore it will be useful in the CT images-based treatment planning especially to use above procedure apply prescribed method when CT images intermittently fail to visualize tumor volume comparing to MRI images.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.28-37
• /
• 2002

In this study nonlinear dynamic behaviors of towed tow-tension cables are numerically analysed. In the case of a taut cable analysis, a bending stiffness term is usually neglected due to its minor effect but it plays an important role in a low-tension cable analysis. A low-tension cable may experience large displacements due to relatively small restoring forces and thus the effects of fluid and geometric non-linearities become predominant. The bending stiffness and non-linearity effects are considered in this work. In order to obtain dynamic behaviors of a towed low-tension cable, three-dimensional nonlinear dynamic equation is described and discretized by employing a finite difference method. An implicit method and Newton-Raphson iteration are adopted for the time integration and nonlinear solutions. 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 those of a in-house program of WHOI Cable with good agreements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.463-474
• /
• 2009

Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

• Journal of the Korean Institute of Telematics and Electronics T
• /
• v.36T no.2
• /
• pp.21-31
• /
• 1999

This paper presents a new possibility of calculating optimal control for large scale which consist of non-linear dynamic sub-systems using two level hierarchical structures method. And the proposed method is based on the idea of block pulse transformation to simplify the algorithm and its calculation. This algorithm used an expansion around the equilibrium point of the system to fix the second and higher order terms. These terms are compensated for iteratively at the second level by providing a prediction for the states and controls which form of a part of the higher order terms. In this new approach the quadratic penalty terms are not used in the cost function. This allows convergence over a longer time horizon and also provides faster convergence. And the method is applied to the problem of optimization of the synchronous machine. Results show that the new approach is superior to conventional numerical method or other previous algorithm.