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

한국전산구조공학회 한국전산구조공학회 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.

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.

구조해석 및 설계용 프로그램의 현황조사

한국전산구조공학회 한국전산구조공학회 1990 전산구조공학 Vol.3 No.2

구조해석 및 설계분야에서의 컴퓨터 이용은 이미 일반화 되어 수많은 프로그램들이 개발되어 사용되고 있다. 70년대 초에 SAP, NASTRAN 같은 범용 해석프로그램이 개발되어 보급되면서, 구조해석 및 설계업무에서의 컴퓨터의 이용은 확산되었다. 이제는 국내의 소규모의 용역회사나 설계사무소에서도 컴퓨터를 이용하여 구조해석 및 설계를 하는 것이 상당히 보편화되었다고 할 수 있다. 구조해석 및 설계에 사용되는 컴퓨터프로그램은 SAP, STRUDUL 같은 대형 범용 프로그램에서부터 특수한 구조물을 해석 또는 설계하기 위하여 개발된 전용 프로그램까지 매우 다양하다. 또한 수많은 프로그램들이 개발되었기 때문에, 자기에게 당면해 있는 문제를 해결하기 위하여 기존의 프로그램 중 어떤 것을 선택하는 것이 가장 효율적인 것인가, 또는 새로운 프로그램을 개발해양 하는 것인가를 결정하기 어려운 경우가 많다. 이러한 관점에서 국내에서 현재 사용 또는 판매되고 있는 구조해석 및 설계용 프로그램의 현황을 조사하여, 그 결과를 정리하였다.

The Finite Element Analysis for Calculations of Equivalent Elastic Constants Using the Homogenization Method

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

본 논문은 구조물의 미시적 측면에서 유효평균탄성계수를 결정하기 위한 균질화기법인 점근적 방법을 적용하였고, 탄성값을 조사하기 위하여 유한요소법으로 정식화하였다. 수치 예로서 물성치가 각기 다른 등방성 재료를 적층한 부재의 임의 단면에서 단위요소를 해석영역으로 설정하고 산출된 탄성계수를 기존의 해석방법으로부터 산출된 값과 비교하였다. 균질화기법으로 산출된 탄성계수는 과소평가되어 나타나며, 이는 해석영역을 유한요소정식화하는 과정에서 수정항만큼 차이가 난다는 것을 증명하였다. 기존 해석방법으로는 복합재료의 탄성계수가 단순히 재료의 산술적 평균값으로 계산되는 것과는 달리, 미시적으로 복합재 단위요소의 반복성을 고려함으로써 제안된 해석방법이 보다 유용하다는 것을 보여 주었다. This paper discusses the homogenization method to determine effective average elastic constants of a linear structure by considering its microstructure. A detailed description on the homogenization method is given for the linear elastic material and then the finite element approximation is performed for an investigation of elastic properties. An asymptotic expansion is carried out in the cross- section area, or in the unit cell. Two and three lay - up structures made up of individual isotropic constituents are chosen for numerical exarnples to check discrepancies between results generated by this theoretical development and the conventional approach. Asymptotic characteristics of the process in extracting the stiffness of structure locally formed by spatial repetitions yield underestimated values of stiffness These discrepancies are detected by the asymptotic corrective term which is ascribed to considerations of microscopic perturbations and proved in the finite element formulation. The asymptotic analysis is the more reasonable in analysing the composite material, rather than the conventional approach to calculate the macroscopic average for elastic properties.

절점이동과 단항증가법에 의한 이차원 평면문제의 적응 유한요소 해석

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

최근, 유한요소해석견과의 신뢰도를 향상시키기 위하여 활발하게 연구되고 있는 적응유한요소해석은 반복계산을 통해서 해석결과의 오차가 사용자에 의해 지정된 허용오차와 같아지도록 하는 해석방법이다. 이와 간은 적응유한요소해석은 해석결과의 오차평가와 이에 따른 유한요소의 재구성과정으로 나누어진다. rp방법에서는 절점의 위치를 이동시켜 요소의 크기를 조절하는 r방법과 형상함수찻수를 증가시키는 p방법을 동시에 적용함으로써 적응해석의 유효성을 향상시키고자 하였다. 제안한 rp방법의 특성을 규명하고 적응해석의 유효성을 보이기 위하여 전형적인 이차원 평면문제들을 해석하고 그 결과를 검토하였다. Adaptive finite element analysis, in which its solution error meets with the user defined allowable error, is recently used to improve the reliability of finite element analysis results. This adaptive analysis is composed of two procedures; one is the error estimation of an analysis result and the other is the reconstruction of finite elements. In the (p-method, an element size is controlled by relocating of nodal positions (r-method) and the order of an element shape function is determined by the hierarchical polynomial (p-method) corresponding to the clement solution error by the enhanced SPR. In order to show the effectiveness and the accuracy of the suggested rp-method, various numerical examples were analyzed and these analysis results were examined by comparing with those obtained by the existed methods.

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

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

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

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

한국전산구조공학회 한국전산구조공학회 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

기계부품의 설계는 초기 설계, 해석, 성능 평가의 반복 과정을 통하여 수행된다. 설계자는 각 과정에서 특성에 맞는 프로그램을 사용하고 있다. 본 연구에서는 순차 설계 영역을 이용한 형상최적화를 수행하였다. 순차 설계영역의 근사함수를 구하기 위하여 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.

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

한국전산구조공학회 한국전산구조공학회 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).