경계요소법에 의한 대규모 3차원 지하구조물의 동적해석

한국전산구조공학회 한국전산구조공학회 1995 전산구조공학 Vol.8 No.3

3차원 대규모 지하구조물의 동적응답을 결정하기 위한 일반적인 수치해석이 제안되었다. 지반과 구조물을 해석하기 위하여 Laplace 변환을 적용한 경계요소법을 설명하였고, 지반-구조물계에 작용하는 외부 동적하중과 지진파를 고려할 수 있도록 공식화하였다. 동적교란이 전파되는 경우에 시간영역의 응답을 얻기 위하여는 구해진 변화된 해를 수치적인 Laplce 역변환을 수행하여야 하지만 동적교란이 조화적인 경우에는 응답이 주파수 영역으로부터 직접 얻어지며, 역변환이 필요하지 않다. 이 방법의 특징은 높은 정확도와 효율성이며, 지반-구조물계에 대하여 초기조건 및 점탄성 재료의 거동을 쉽게 고려할 수 있다는 것이다. 그러므로 이 방법은 다양한 지하구조물의 동적거동과 지진에 대한 취약함을 연구하기 위한 적절한 도구로 사용되어 질 수 있다.

순차설계영역을 이용한 근사 형상최적에 관한 연구

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

기계부품의 설계는 초기 설계, 해석, 성능 평가의 반복 과정을 통하여 수행된다. 설계자는 각 과정에서 특성에 맞는 프로그램을 사용하고 있다. 본 연구에서는 순차 설계 영역을 이용한 형상최적화를 수행하였다. 순차 설계영역의 근사함수를 구하기 위하여 Pro/Engineer 와 ANSYS 실행의 자동화를 수행하였다. 전체 설계영역을 근사식으로 표현하기에는 어려움이 있다. 정확도가 높은 근사식을 만들기 위하여 순차설계영역을 설정하여 각 단계에서 수렴한 해로 이동량을 결정하고, 두 번 연속하여 순차설계 영역에 존재하면 수렴조건을 만족하는 것으로 하였다. 각 단계의 해는 순차이차 계획법인 PLBA(Pshenichny-Lim-Belegundu Arora)알고리즘을 이용하여 구하였다. Mechanical design process is generally accomplished by design, analysis, and test. Designers use programs fitting purpose, and obtain repeatedly a response of a simulation program, a sub-program for optimization. In this paper, shape optimization using approximate optimization technique is carried out with sequential design domain(SDD). In addition, algorithm executing Pro/Engineer and ANSYS automatically are adopted in the approximate optimization program by SDD. It is difficult for design problem to be approximated accurately for the whole range of design space. However, more or less accurate approximation is constructed if SDD is applied to that case. SDD starts with a certain range which is off-seted from midpoint of an initial design domain and then SDD of the next step is determined by a move limited. Convergence criterion is defined such that optimal point must be located within SDD during the two steps. Also, the PLBA(Pshenichny-Lim-Belegundu-Arora) algorithm is used to solve approximate optimization problems. This algorithm uses the second-order information and the active set strategy, in order to seek the direction of design variables.

경계 요소법을 이용한 2차원 비등방성 복합재료 탄성체의 비파괴 결함 추정

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

본 연구에서는 경계 요소법을 이용하여 2차원 비등방성 탄성체의 손상 규명을 수행한다. 경계에서의 적분항만을 포함하는 본 수치모델은 1차원으로 줄게된다. 이러한 장점은 특히 균열 역학과 같은 문제에 있어서 중요한 의미를 갖는다. 또한 일정 영역을 분할하는 기법을 피함으로서 수치해석 과정을 간편하고 효율적으로 수행할수 있기 때문에 역문제 해결에 있어서 장점을 갖는다. 본 연구에서는 기존의 등방성 재료에 대한 비파괴 추정기법을 복합신소재와 같은 비등방성 재료로 이루어진 탄성체의 해석에 대하여 확장한다. 먼저 경계요소법에 의한 수치모델의 타당성을 기존의 문헌과 비교 검증하며, 서로 다른 특성을 보이는 비등방성 형식의 변화에 따라 실제 측정시 발생하는 노이즈 영향을 분석한다. 수치예제는 적층 형태 및 하중조건에 대하여 수행하며, 결함 추정에 미치는 적층 형태의 영향을 시험한다. In this paper, the defects of two-dimensional anisotropic elastic bodies are identified by using the boundary element method. The use of numerical models that contain only boundary integral terns reduces the dimensionality of the problem by one. This advantage is particularly important in problems such as crack mechanics. Avoiding domain meshing is also particularly advantageous in the solution of inverse problems since it overcomes mesh perturbations and simplifies the procedure. In this paper, nondestructive approaches for the existing isotropic materials are extended to analyze the elastic bodies made of anisotropic materials such as composites. After verifying that the proposing boundary element model is in good agreement with numerical results reported by other investigators, the effect of noise in the measurements on the identifiability is studied with respect to different design parameters of layered composites. Sample studies are carried out for various layup configurations and loading conditions. The effects of the layup sequences in detecting flaw of composites is explored in this paper.

고차전단변형 판이론을 이용한 채널단면을 갖는 복합적층 절판 구조물의 유한요소 진동 해석

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

본 연구에서는 유한요소법을 이용한 채널단면을 갖는 복합재료 적층 구조물의 자유진동을 다룬다. 복합적층 절판구조물에 고차항 판이론을 적용하기 위하여 개발된 유한요소 프로그램은 Lagrangian 및 Hermite 보간함수를 병용하여 면내회전각 자유도를 포함한 절점 당 8개의 자유도를 갖는다. 전단보정계수의 가정을 필요로 하지 않고 전단변형의 3차항 비선형 특성이 고려된 본 논문의 절판 요소는 국부좌표계와 전체좌표계에 대한 좌표변환행렬에 의하여 요소 당 32의 국부요소행렬로 구성된다. 본 해석 프로그램의 결과는 기존의 고전적 이론 및 일차항 이론에 의한 문헌 결과와 비교ㆍ분석하였으며, 화이버 보강각도, 길이-두께비, 기하학적 형상 변화 등의 다양한 매개변수 연구를 수행하였다. 본 연구에서는 특히 경계조건 및 길이-두께비 변화에 따라 예측하기 힘든 복잡한 거동을 보이는 복합적층 채널단면 구조물의 자유진동에 대하여 정밀한 고차항 이론 적용에 의한 엄밀 해석의 필요성을 제기하였다. This study deals with free vibrations of laminated composite structures with a channel section using finite element method. In this paper, the mixed finite element method using Lagrangian and Hermite interpolation functions is adopted and a high-order plate theory is used to analyze laminated composite non-prismatic folded plates with a channel section more accurately for free vibration. The theory accounts for parabolic distribution of the transverse shear stress and requires no shear correction factors supposed in the first-order plate theory. An 3232 matrix is assembled to transform the system element matrices from the local to global coordinates using a coordinate transformation matrix, in which an eighth drilling degree of freedom (DOF) per node is appended to the existing 7-DOF system. The results in this study are compared with those of available literatures for the conventional and first-order plate theory. Sample studies are carried out for various layup configurations and length-thickness ratio, and geometric shapes of plates. The significance of the high-order plate theory in analyzing complex composite structures with a channel section is enunciated in this paper.

횡하중과 연직하중을 받는 고층건물의 변위설계를 위한 재분배기법 개발

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

고층 건물의 구조시스템 및 물량은 최대변위 및 층간 변위 제어를 포함하는 강성설계에 의해서 결정되는 경향이 있으나 강도설계와는 달리 특정 소수 부재의 단면 변화로 강성제한 조건을 만족시키기는 어렵다. 이러한 이유로 이제까지의 고층건물 구조설계는 설계자의 경험과 직관에 의존한 반복적인 설계과정을 피할 수 없었다. 그러므로 본 연구에서는 민감도 해석이나 반복적인 구조해석 없이 변위를 제어할 수 있는 실용적인 고층건물 변위제어기법을 제안하였다. 기존 연구에서 개발된 재분배기법은 구조물이 횡하중을 받는 경우만을 고려하였으나 실제 구조물은 횡하중과 연직하중을 동시에 받고 있으므로 이를 고려한 재분배기법을 개발하였다. 개발된 횡하중과 연직하중을 고려한 재분배기법을 강도설계모듈과 결합하여 고층건물 변위설계 모델을 제안하였으며 제안된 변위설계모델을 20층 강접골조 전단벽 구조물과 60층 아웃리거 구조물의 변위설계에 적용하여 기존의 횡하중 만을 고려한 변위설계모델의 결과와 비교 분석하였다. Drift design of a high-rise building is a governing factor in the determination of structural weights and lateral resisting systems. However, high-rise buildings are composed of tens of thousands of structural member, designer can not know which members are active to lateral drift control and how much they contribute to lateral drifts. Resizing technique was proved to be a practical method for drift design of high-rise buildings. However, no resizing algorithm has been considered the effect of vertical loads in drift designs. Thus, in this paper, a resizing algorithm has been developed for drift designs of high-rise buildings subjected to both lateral and vertical loads. The drift design model has been applied to drift designs of two high-rise building examples.

정현형 조화하중에 의한 아치 구조물의 동적 좌굴 특성에 관한 연구

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

동적 불안정 좌굴현상에 관한 연구는 다소 발표되고 있으나, 주기성을 가진 하중하에서의 동적 좌굴을 다룬 연구는 그리 많지 않은 편이다. 주기성을 가진 하중하에서의 거동은 STEP 하중하에서의 거동과는 다르리라 예상된다. 본 논문에서는 동적 불안정의 기본 메커니즘을 파악하기 위하여 양단 핀으로 고정된 정현형 아치가 정현형 조화하중을 받았을 때의 얕은 아치를 대상으로 한다. Newmark- 법에 의한 수치적분을 이용하여 비선형 운동 방정식의 변위 응답을 구하고, 얻어진 비선형 변위 응답으로 FFT(Fast Fourier Transform)에 의한 연속 응답 스펙트럼을 구해 동적 불안정 특성에 관해서 분석한다. The dynamic instability for snapping phenomena has been studied by many researchers. Few paper deal with the dynamic buckling under the load with periodic characteristics, and the behavior under periodic excitation is expected the different behavior against STEP excitation. We investigate the fundamental mechanisms of the dynamic instability when the sinusoidal shaped arch structures are subjected to sinusoidal harmonic excitation with pin-ends. By using Newmark- method, we can get the nonlinear displacement response, and using this analyze characteristics of the dynamic instability through the running response spectrum by FFT(Fast Fourier Transform).

응력파 전파 수치모의를 위한 일차원 유한요소모형의 분산 특성 및 제어

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

본 논문에서는 응력파 전파를 수치모의할 때 발생하는 수치적인 분산효과를 제거하기 위해 파동방정식에 기초한 일차원 유한요소모형을 이용하여 수치분산오차의 특성을 분석하고 분산오차를 제어할 수 있는 방법을 제안하였다. 질량행렬을 그대로 사용하는 경우와 집중질량행렬을 사용하는 경우에 대한 수치분산오차를 분석하였다. 개발된 분산제어기법은 공간미분항의 시간단계 가중치 및 질량집중도를 조정하는 음해법과 인위적인 분산항을 추가하는 양해법의 두가지 방법이다. 제안된 분산보정기법을 이용하여 계산한 수치해와 파동방정식의 해석해를 비교한 결과 본 연구에서 제안한 분산보정기법의 타당성을 확인하였다. With an aim at eliminating the numerical dispersion error arising from the numerical simulation of stress wave propagation, numerical dispersion characteristics of the wave equation based one-dimensional finite element model are analyzed and some dispersion control scheme are proposed in this paper The dispersion analyses are carried out for two types of mass matrix, namely the consistent and the lumped mass matrices. Based on the finding of the analyses, dispersion correction techniques are developed for both the implicit and explicit schemes. For the implicit scheme, either the weighting factor for the spatial derivatives of each time level or the lumping coefficient for mass matrix is adjusted to minimize the numerical dispersion. In the case of the explicit scheme an artificial dispersion term is introduced in the governing equation. The validity of the dispersion correction techniques proposed in this study is demonstrated by comparing the numerical solutions obtained using the Present techniques with the analytical ones.

유한요소법에 의한 직방성재료 내 균열의 응력확대계수 결정

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

응력확대계수는 균열진전경로의 수치해석적 연구에서 널리 사용되고 있다. 그러나 많은 경우에서 균열선단주위 응력의 급수전개식 중 이어지는 항은 정량적으로 중요하다. 따라서 본 연구에서는 이의 항을 계산하기 위하여 등매개 2차특이요소를 이용한 유한요소해석을 수행하였다. 일례로 단축하중을 받는 무한 직방성체 내 경사균열에 대하여 균열요소크기와 균열경사각을 달리 주어 가며 해석을 수행하였으며, 수치해석결과는 이론해와 비교하여 잘 일치하고 있다. The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient

Path-Integral Solution을 이용한 랜덤동요된 조화가진력을 받는 임팩트시스템의 거동분석

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

랜덤동요된 조화가진력을 받는 임팩트시스템의 비선형거동을 개발된 반해석적절차에 의해 확률영역에서 분석하였다. 반해석적절차는 path-integral solution을 이용하여 임팩트시스템의 추계론적 미분방정식으로부터 구함으로 얻어진다. 결합확률밀도함수의 전개를 구하고 시스템의 비선형거동 특성인 혼돈거동에 대하여 분석하고 노이즈의 영향을 시간영역과 확률영역에서 알아보았다. 결과로부터 반해석적절차는 결합확률밀도함수를 통하여 임팩트시스템의 거동에 대한 정보를 제공하는 것을 알 수 있었다. 노이즈의 영향은 혼돈거동의 특성을 약화시키며 궁극적으로 사라지게 함을 알 수 있었으며 또한 혼돈거동의 특성이 상대적으로 높은 노이즈아래에서도 남아있는 것을 밝혔다. 결합확률밀도함수는 응답앙상블이 약정상과정임을 확인시켜 주었다. Nonlinear system responses of an impact system under randomly perturbed harmonic excitations are predicted in the probability domain by adopting the semi-analytical procedure previously developed. The semi-analytical procedure is obtained by solving the Fokker-Planck equation corresponding to the stochastic differential equation of the given impact system by utilizing the path-integral solution. The evolutionary joint probability density functions are generated by using the method, and the characteristics of nonlinear dynamic response behaviors of the system are examined. Noise effects on the responses are also examined. It Is found that the semi-analytical method can provides the accurate information of the responses via the joint probability functions for the impact system. It is found that the noises weaken and eventually terminate the chaos in the responses, but it is also found that the chaotic signatures reside in the presence of the external noise with relatively high intensity. The joint probability density function shows that the ensemble of the system responses are weakly stationary.

단면변형의 효과를 포함한 강상자형 거더의 엄밀한 해석

한국전산구조공학회 한국전산구조공학회 2004 한국전산구조공학회논문집 Vol.17 No.1

본 연구에서는 직선 강상자형 거더의 단면변형에 의한 변형 및 응력계산을 위한 Matlab 해석프로그램을 개발하고자 한다. 이를 위하여 단면변형이론을 요약하고 빔유사이론을 제시한다. 이후 탄성지반위의 보-기둥부재의 지배방정식을 제시하고, 일반화된 고유치해석을 통하여 집중 및 분포하중을 받는 보요소의 엄밀한 강성행렬을 계산한다. 본 연구의 효율성과 정확성을 입증하기 위하여 격벽을 갖는 상자형 거더의 뒤틀림응력을 계산하고 유한요소해와 비교한다. The main goal of this study is to develop MATLAB programming for an analysis of distortional deformations and stresses of the straight box girder. For this purpose, a distortional deformation theory is firstly summarized and then a BEF (Beam on Elastic Foundation) theory is presented using analogy of the corresponding variables. Finally, with governing equations of the beam-column element on elastic foundation, an exact element stiffness matrix of the beam element and nodal forces equivalent to concentrated and distributed loads are evaluated via a generalized linear eigenvalue problem. In order to verify the efficiency and accuracy of this method, distortional stresses of box girders with multiple diaphragms are presented and compared with results by FEA.