• Coupled systems mechanics
• /
• 제8권5호
• /
• pp.393-414
• /
• 2019

This paper presents seismic analysis of concrete gravity dams considering soil-structure-fluid interaction. Displacement based plane strain finite element formulation is considered for the dam and foundation domain whereas pressure based finite element formulation is considered for the reservoir domain. A direct coupling method has been adopted to obtain the interaction effects among the dam, foundation and reservoir domain to obtain the dynamic responses of the dam. An efficient absorbing boundary condition has been implemented at the truncation surfaces of the foundation and reservoir domains. A parametric study has been carried out considering each domain separately and collectively based on natural frequencies, crest displacement and stress at the neck level of the dam body. The combined frequency of the entire coupled system is very less than that of the each individual sub-system. The crest displacement and neck level stresses of the dam shows prominent enhancement when coupling effect is taken into consideration. These outcomes suggest that a complete coupled analysis is necessary to obtain the actual responses of the concrete gravity dam. The developed methodology can easily be implemented in finite element code for analyzing the coupled problem to obtain the desired responses of the individual subdomains.

• Structural Engineering and Mechanics
• /
• 제85권1호
• /
• pp.89-104
• /
• 2023

This work presents a proposal for employing reduced numerical integration in the formulation of the 4-node quadrilateral solid finite element. The use of these low-order integration rules leads to numerical instabilities such as those producing the hourglass effect. The proposed procedure allows evaluating a given constitutive model only in one integration point, achieving an attractive computational cost reduction and, also, successfully controls the hourglass effect. A validation of the proposal is included and discussed throughout the paper. To show the efficiency of the proposal, several application examples of masonry structures are studied and discussed. To represent the non-linear mechanical behaviour of masonry a plastic-damage model is implemented within the application of this sub-integration scheme. Also, in order to have a full and computationally efficient strategy to determine the behaviour of masonry structures, involving its evolution to collapse, a homogenization technique with a macro-modeling approach is used. The methodology discussed throughout this paper demonstrates a substantial computational cost reduction and an improved approximation of the non-linear problem evidenced by a reduction of up to 85% of the computational time for some cases.

• Earthquakes and Structures
• /
• 제26권1호
• /
• pp.59-75
• /
• 2024

This paper presents a dam-reservoir interaction model that includes, water compressibility, sloshing of surface water, and radiation damping at the far-end reservoir, to investigate the influence of concrete deterioration on seismic behavior along with seismic performance of gravity dams. Investigations on seismic performance of the dam body have been conducted using the linear time-history responses obtained under six real and 0.3 g normalized earthquake records with time durations from 10 sec to 80 sec. The deterioration of concrete is assumed to develop due to mechanical and chemical actions over the dam lifespan. Several computer programs have been developed in FORTRAN 90 and MATLAB programming languages to analyze the coupled problem considering two-dimensional (2D) plane-strain condition. According to the results obtained from this study, the dam structure shows critical responses at the later ages (75 years) that could cause disastrous consequences; the critical effects of some earthquake loads such as Chi-Chi with 36.5% damage and Loma with 56.2% damage at the later ages of the selected dam body cannot be negligible; and therefore, the deterioration of concrete along with its effects on the dam response should be considered in analysis and design.

• 한국전기전자재료학회:학술대회논문집
• /
• 한국전기전자재료학회 2010년도 하계학술대회 논문집
• /
• pp.31-31
• /
• 2010

The growth of the high-quality GaN epilayers is of significant technological importance because of their commercializedoptoelectronic applications as high-brightness light-emitting diodes (LEDs) and laser diodes (LDs) in the visible and ultraviolet spectral range. The GaN-based heterostructural epilayers have the polar c-axis of the hexagonal structure perpendicular to the interfaces of the active layers. The Ga and N atoms in the c-GaN are alternatively stacked along the polar [0001] crystallographic direction, which leads to spontaneous polarization. In addition, in the InGaN/GaN MQWs, the stress applied along the same axis contributes topiezoelectric polarization, and thus the total polarization is determined as the sum of spontaneous and piezoelectric polarizations. The total polarization in the c-GaN heterolayers, which can generate internal fields and spatial separation of the electron and hole wave functions and consequently a decrease of efficiency and peak shift. One of the possible solutions to eliminate these undesirable effects is to grow GaN-based epilayers in nonpolar orientations. The polarization effects in the GaN are eliminated by growing the films along the nonpolar [$11\bar{2}0$] ($\alpha$-GaN) or [$1\bar{1}00$] (m-GaN) orientation. Although the use of the nonpolar epilayers in wurtzite structure clearly removes the polarization matters, however, it induces another problem related to the formation of a high density of planar defects. The large lattice mismatch between sapphiresubstrates and GaN layers leads to a high density of defects (dislocations and stacking faults). The dominant defects observed in the GaN epilayers with wurtzite structure are one-dimensional (1D) dislocations and two-dimensional (2D) stacking faults. In particular, the 1D threading dislocations in the c-GaN are generated from the film/substrate interface due to their large lattice and thermal coefficient mismatch. However, because the c-GaN epilayers were grown along the normal direction to the basal slip planes, the generation of basal stacking faults (BSFs) is localized on the c-plane and the generated BSFs did not propagate into the surface during the growth. Thus, the primary defects in the c-GaN epilayers are 1D threading dislocations. Occasionally, the particular planar defects such as prismatic stacking faults (PSFs) and inversion domain boundaries are observed. However, since the basal slip planes in the $\alpha$-GaN are parallel to the growth direction unlike c-GaN, the BSFs with lower formation energy can be easily formed along the growth direction, where the BSFs propagate straightly into the surface. Consequently, the lattice mismatch between film and substrate in $\alpha$-GaN epilayers is mainly relaxed through the formation of BSFs. These 2D planar defects are placed along only one direction in the cross-sectional view. Thus, the nonpolar $\alpha$-GaN films have different atomic arrangements along the two orthogonal directions ($[0001]_{GaN}$ and $[\bar{1}100]_{GaN}$ axes) on the $\alpha$-plane, which are expected to induce anisotropic biaxial strain. In this study, the anisotropic strain relaxation behaviors in the nonpolar $\alpha$-GaN epilayers grown on ($1\bar{1}02$) r-plane sapphire substrates by metalorganic chemical vapor deposition (MOCVO) were investigated, and the formation mechanism of the abnormal zigzag shape PSFs was discussed using high-resolution transmission electron microscope (HRTEM).

• 한국전산구조공학회논문집
• /
• 제28권5호
• /
• pp.441-449
• /
• 2015

본 논문에서는 커코프 판이론과 폰-칼만 비선형 변형율-변위 관계를 이용하여 서형화된 좌굴해석을 수행하였다. 평면응력과 좌굴문제에서 영률과 두께에 관한 설계민감도식을 유도하였고, 고유치를 최대화하면서 컴플라이언스를 최소화하는 위상최적설계 기법을 정식화하였다. 좌굴해석에서의 프리스트레스를 이용하여 판 좌굴문제에 적용할 수 있는 위상최적설계 기법을 개발하였다. 폰-칼만 비선형 변형률을 사용하여 좌굴문제의 응력행렬을 구성하는데 프리스트레스가 필요하므로 면외로의 운동을 도입하였다. 위상최적설계를 위하여 정규재료밀도를 설계변수로 하고, 목적함수는 최소 컴플라이언스와 최대 고유진동수로 하였으며 제한조건은 허용되는 재료량이다. 여러 수치예제를 통하여 개발된 설계민감도 해석법은 유한차분 민감도와 비교하여 매우 정확한 값을 가지고, 위상최적설계는 물리적으로 의미있는 결과를 제공함을 확인하였다.

• 한국전산구조공학회논문집
• /
• 제13권4호
• /
• pp.449-459
• /
• 2000

기존의 Mindlin이론을 바탕으로 한 감절점 쉘 요소는 두께가 두꺼운 쉘 구조를 해석할 경우 매우 좋은 결과가 얻어지는 것으로 알려져 있다. 그러나 두께가 얇은 쉘의 경우 전단구속 및 막구속 현상으로 인해 요소의 강성이 과대하게 평가되는 것으로 나타났다. 이러한 문제를 해결하기 위해 본 연구에서는 서로 다른 두 가지 막 및 전단 구속 해결방법을 조합하여 쉘 요소에 적용하였다. 첫 번째 요소는 전단 변형률에 가정된 전단변형률을 적용하였고, 막 변형률에는 감차적분법을 적용하였다. 두 번째 요소의 경우, 막 변형률과 전단 변형률 모두에 가정된 전단 변형률을 적용하여 쉘 구조해석을 수행하였다. 여러 가지 수치 해석을 수행한 결과, 두 가지 쉘 요소가 빠른 수렴성을 가지고 있는 것으로 나타났으며, 막 구속과 전단 구속 현상이 제거된 것을 나타났다.

• Structural Engineering and Mechanics
• /
• 제74권2호
• /
• pp.227-241
• /
• 2020

The bi-axial and shear buckling behavior of laminated hypar shells having rhombic planforms are studied for various boundary conditions using the present mathematical model. In the present mathematical model, the variation of transverse shear stresses is represented by a second-order function across the thickness and the cross curvature effect in hypar shells is also included via strain relations. The transverse shear stresses free condition at the shell top and bottom surfaces are also satisfied. In this mathematical model having a realistic second-order distribution of transverse shear strains across the thickness of the shell requires unknown parameters only at the reference plane. For generality in the present analysis, nine nodes curved isoparametric element is used. So far, there exists no solution for the bi-axial and shear buckling problem of laminated composite rhombic (skew) hypar shells. As no result is available for the present problem, the present model is compared with suitable published results (experimental, FEM, analytical and 3D elasticity) and then it is extended to analyze bi-axial and shear buckling of laminated composite rhombic hypar shells. A C⁰ finite element (FE) coding in FORTRAN is developed to generate many new results for different boundary conditions, skew angles, lamination schemes, etc. It is seen that the dimensionless buckling load of rhombic hypar increases with an increase in c/a ratio (curvature). Between symmetric and anti-symmetric laminations, the symmetric laminates have a relatively higher value of dimensionless buckling load. The dimensionless buckling load of the hypar shell increases with an increase in skew angle.

• 대한기계학회논문집
• /
• 제13권5호
• /
• pp.820-830
• /
• 1989

본 논문에서는 타이어의 각 부분의 물성치 계산을 위한 식을 유한요소법에 적용할 수 있도록 제안하였다. 이 식은 강철 코드의 굽힙효과를 고려 하였으며, 특히 각 요소에서 전단변형이 일어나는 동안의 굽힘효과를 고려하였다. 유한요소 공식화는 가상일의 원리에 의하여 평형 방정식으로부터 유도하였고, Updated refer- ence coordinate에 대해 증분해석을 적용하여 Updated Lagrangian공식화를 하였다. 그리고 차량하중에 의하여 타이어가 노면에 접지될때의 응력상태를 게산할 수 있도록 접촉문제 공식화를 유한요소 공식화에 첨가 하였다.

• 한국전산구조공학회논문집
• /
• 제24권5호
• /
• pp.513-520
• /
• 2011

본 연구에서는 스트라치 시스템의 긴장설치과정 해석을 수행하기 위한 개선된 명시적 해석법을 제안하였다. 스트라치 시스템은 Stressed-Arch에서 유래한 용어로 슬리브에 의한 갭이 도입된 유동 하현재가 내부의 긴장재에 의한 초기장력의 도입으로 서서히 닫히게 되고, 이에 따라 전체 구조물이 상승하여 최종적인 아치형태의 구조물을 형성하는 독창적인 구조시스템이다. 스트라치 시스템의 초기장력 도입과정을 긴장설치(stress-erection process) 과정이라 하며, 초기곡률의 도입에 따라 유동 상현재에는 과도한 초기변형이 발생하여 소성거동에 의한 강체회전이 발생하는 불안정 구조물이 된다. 본 연구에서는 이러한 스트라치 시스템의 불안정 거동특성을 해석하기 위해서 소성힌지가 적용된 보-기둥요소를 사용하여 유동상현재를 모델링하였고, 불안정 구조물의 해석법에 효과적으로 사용되는 동적이완법의 개선된 알고리즘을 개선하여, 실제 스트라치 구조물의 긴장설치과정 해석을 수행하였으며, 실제 프로젝트에 대한 해석결과의 분석을 통하여 제안된 해석법의 적용성을 검증하였다.

• 비파괴검사학회지
• /
• 제19권1호
• /
• pp.16-24
• /
• 1999

잔류응력은 공학적 구조물에 결함을 야기하는 유해한 응력이므로 산업현장에서는 이를 측정하기 위한 많은 관심을 가지고 있다. 현재 잔류응력 측정은 일반적인 방법을 사용하고 있지만 이러한 방법들은 어느정도 문제점을 내포하고 있다. 그러므로 레이저 스페클 간섭법에 유한요소법과 점 가열법을 혼합하여 산업 현장에 유해한 잔류응력을 측정할 수 있는 새로운 기법을 개발하였다. 이 기법에서 사용된 레이저 스페클 간섭계는 측정 시편의 응력이 풀릴 때의 면내 변형을 측정하기 위해서 사용되어 졌으며, 유한요소법은 응력 풀림을 위한 가열부의 온도와 그 밖의 변수들을 결정하기 위해서 사용되어 졌다. 잔류응력은 측정하고자 하는 부위의 가열 냉각에 의한 변형률에 의해 결정되어지며, 이에 대한 간단한 모델링을 제시하였다. 또한 본 실험은 변형의 3차원 정량화를 수행하기 위해서 레이저 스페클 간섭계의 전자 광학적 기법 중의 하나인 위상이동법을 이용하여 수행되어졌다.