• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2006년도 춘계 학술발표회 논문집
• /
• pp.286-300
• /
• 2006

Nonlinear cross-anisotropic properties of soil is critical for exact numerical simulation. Theoretically, initial cross-anisotropic properties can be evaluated from triaxial tests with bender element tests, and nonlinear cross-anisotropic properties over initial strain level cannot be evaluated from triaxial tests. In this study, a supposed condition among nonlinear cross-anisotropic properties is suggested to calculate nonlinear cross-anisotropic properties from triaxial tests. Maximum strain and incremental strain energy are applied to combine triaxial test results and theoretical normalized shear modulus curve, respectively Based on combined results, nonlinear cross-anisotropic properties are calculated. Numerical simulation for triaxial tests Is carried out to verify the applicability of the supposed condition with calculated cross-anisotropic properties and simplified nonlinear cross-anisotropic model.

• Structural Engineering and Mechanics
• /
• 제40권2호
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Nuclear Engineering and Technology
• /
• 제53권9호
• /
• pp.3100-3111
• /
• 2021

Numerical modeling for the safety-related equipment used in a nuclear power plant (i.e., cabinet facilities) plays an essential role in seismic risk assessment. A full finite element model is often time-consuming for nonlinear time history analysis due to its computational modeling complexity. Thus, this study aims to generate a simplified model that can capture the nonlinear behavior of the electrical cabinet. Accordingly, the distributed plasticity approach was utilized to examine the stiffness-degradation effect caused by the local buckling of the structure. The inherent dynamic characteristics of the numerical model were validated against the experimental test. The outcomes indicate that the proposed model can adequately represent the significant behavior of the structure, and it is preferred in practice to perform the nonlinear analysis of the cabinet. Further investigations were carried out to evaluate the seismic behavior of the cabinet under the influence of the constitutive law of material models. Three available models in OpenSees (i.e., linear, bilinear, and Giuffre-Menegotto-Pinto (GMP) model) were considered to provide an enhanced understating of the seismic responses of the cabinet. It was found that the material nonlinearity, which is the function of its smoothness, is the most effective parameter for the structural analysis of the cabinet. Also, it showed that implementing nonlinear models reduces the seismic response of the cabinet considerably in comparison with the linear model.

• 한국해안·해양공학회논문집
• /
• 제23권5호
• /
• pp.335-344
• /
• 2011

곡면의 경계를 가지는 수로에서 비선형 파랑의 상호작용을 모의하기 위한 비정수압 자유수면 모형이 개발되었다. 제안된 모형은 비선형의 3차원 Navier-Stokes 방정식을 곡선좌표 영역에서 계산단계 분리법의 일종인 압력수정법에 의하여 수치적으로 해석된다. 특히, 연직방향으로 변형된 형태의 엇갈린 격자를 이용하여 상대적으로 간단하게 압력방정식과 자유수면 경계조건을 구성하였다. 개발된 모형의 수치해석 정확도는 2차원의 수치 파수조에서 파랑의 비선형 정도에 대하여 5차의 스토스우크스 해석해와 비교하였다. 본 모형의 실제적 적용은 원형의 수로에서 회절과 반사에 의해 변형되는 비선형 파의 변형에 초점을 맞추어 수행하였다. 두개의 파를 중첩한 고비선형의 파에 대한 경우를 제외하고 수치해석 결과는 비선형적인 영향을 고려하지 않은 해석해의 선형적인 중첩과 일치하였다. 두개의 파를 중첩한 고비선형의 파에 대한 모의를 통하여 본 모형은 원형의 수로에서 비선형 군파의 변형에관한 수치적인 모의 가능성을 제시하였다.

• 한국지구과학회지
• /
• 제34권1호
• /
• pp.51-58
• /
• 2013

대기 모델링 연구에서 시간 간격을 적절하게 결정하는 것은 중요한 문제이다. 본 연구에서는 비선형 대기 모형에서 수치 해의 시간 간격에 대한 민감도를 조사하였다. 이를 위해 간단한 무차원화된 역학 모형을 사용하여 시간 간격과 비선형성 인자를 바꾸어가며 수치 실험을 수행하였다. 실험 결과, 비선형성 인자가 영향을 줄 만큼 크지 않고 절단 오차를 무시할 수 있는 경우에는 수치 해가 시간 간격에 민감하지 않았다. 그러나 비선형성 인자가 큰 경우에는 수치 해가 시간 간격에 민감한 것으로 밝혀졌다. 이 경우, 시간 간격이 감소할수록 공간 필터의 강도가 증가하여 작은 규모의 현상이 약하게 모의되었다. 이는 일반적으로 시간 간격이 감소하면 절단 오차가 감소하여 더 정확한 수치 해가 도출된다는 사실과 상충한다. 이러한 충돌은 비선형 모형의 수치 해를 안정하게 하기 위해 공간 필터가 반드시 필요하기 때문에 피할 수 없다.

• 한국지진공학회논문집
• /
• 제26권6호
• /
• pp.227-235
• /
• 2022

For a member model in nonlinear structural analysis, a lumped plastic model that idealizes its flexural bending, shear, and axial behaviors by springs with the nonlinear hysteretic model is widely adopted because of its simplicity and transparency compared to the other rigorous finite element methods. On the other hand, a challenging task in its numerical solution is to satisfy the equilibrium condition between nonlinear flexural bending and shear springs connected in series. Since the local forces between flexural and shear springs are not balanced when one or both springs experience stiffness changes (e.g., cracking, yielding, and unloading), the additional unbalanced force due to overshooting or undershooting each spring force is also generated. This paper introduces an iterative scheme for numerical solutions satisfying the equilibrium conditions between flexural bending and shear springs. The effect of equilibrium iteration on analysis results is shown by comparing the results obtained from the proposed method to those from the conventional scheme, where the equilibrium condition is not perfectly satisfied.

• 대한수학회논문집
• /
• 제22권3호
• /
• pp.453-464
• /
• 2007

A computationally efficient numerical scheme is presented for the phase-field model of two-phase systems for anisotropic interfacial energy. The scheme is solved by using a nonlinear multigrid method. When the coefficient for the anisotropic interfacial energy is sufficiently high, the interface of the system shows corners or missing crystallographic orientations. Numerical simulations with high and low anisotropic coefficients show excellent agreement with exact equilibrium shapes. We also present spinodal decomposition, which shows the robustness of the pro-posed scheme.

• Structural Engineering and Mechanics
• /
• 제56권1호
• /
• pp.27-47
• /
• 2015

In this study, the accuracy and reliability of fully nonlinear method against equivalent linear method for dynamic analysis of soil-structure interaction is investigated comparing the predicted results of both numerical procedures with the results of experimental shaking table tests. An enhanced numerical soil-structure model has been developed which treats the behaviour of the soil and the structure with equal rigour. The soil-structural model comprises a 15 storey structural model resting on a soft soil inside a laminar soil container. The structural model was analysed under three different conditions: (i) fixed base model performing conventional time history dynamic analysis, (ii) flexible base model (considering full soil-structure interaction) conducting equivalent linear dynamic analysis, and (iii) flexible base model performing fully nonlinear dynamic analysis. The results of the above mentioned three cases in terms of lateral storey deflections and inter-storey drifts are determined and compared with the experimental results of shaking table tests. Comparing the experimental results with the numerical analysis predictions, it is noted that equivalent linear method of dynamic analysis underestimates the inelastic seismic response of mid-rise moment resisting building frames resting on soft soils in comparison to the fully nonlinear dynamic analysis method. Thus, inelastic design procedure, using equivalent linear method, cannot adequately guarantee the structural safety for mid-rise building frames resting on soft soils. However, results obtained from the fully nonlinear method of analysis fit the experimental results reasonably well. Therefore, this method is recommended to be used by practicing engineers.

• Geomechanics and Engineering
• /
• 제21권2호
• /
• pp.103-109
• /
• 2020

The inflow rate is of interest in the design of underground structures such as tunnels and buried pipes below the groundwater table. Soil permeability governing the inflow rate significantly affects the hydro-geological behavior of soils but is difficult to estimate due to its wide range of distribution, nonlinearity and anisotropy. Volume changes induced by stress can cause nonlinear stress-strain behavior, resulting in corresponding permeability changes. In this paper, the nonlinearity and anisotropy of permeability are investigated by conducting Rowe cell tests, and a nonlinear permeability model considering anisotropy was proposed. Model modification and parameter evaluation for field application were also addressed. Significance of nonlinear permeability was illustrated by carrying out numerical analysis of a tunnel. It is highlighted that the effect of nonlinear permeability is significant in soils of which volume change is considerable, and particularly appears in the short-term flow behavior.

• Nuclear Engineering and Technology
• /
• 제51권2호
• /
• pp.600-606
• /
• 2019

Studies on the application of the lead rubber bearing (LRB) isolation system to nuclear power plants are being carried out as one of the measures to improve seismic performance. Nuclear power plants with isolation systems require seismic probabilistic safety assessments, for which the seismic fragility of the structures, systems, and components needs be calculated, including for beyond design basis earthquakes. To this end, seismic response analyses are required, where it can be seen that the behaviors of the isolation system components govern the overall seismic response of an isolated plant. The numerical model of the LRB used in these seismic response analyses plays an important role, but in most cases, the extreme performance of the LRB has not been well studied. The current work therefore develops an extreme nonlinear numerical model that can express the seismic response of the LRB for beyond design basis earthquakes. A full-scale LRB was fabricated and dynamically tested with various input conditions, and test results confirmed that the developed numerical model better represents the behavior of the LRB over previous models. Subsequent seismic response analyses of isolated nuclear power plants using the model developed here are expected to provide more accurate results for seismic probabilistic safety assessments.