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

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

Three-Dimensional Vibration Analysis of Cantilevered Laminated Composite Plates

한국전산구조공학회 한국전산구조공학회 2001 한국전산구조공학회논문집 Vol.14 No.3

본 논문은 캔틸레버 복합 적층판의 고유진동에 대하여 3차원으로 해석한 연구를 제시하고 있다. 본 연구에서는 고정단에서의 경계조건을 엄밀히 만족하고 수학적으로 완전한 다항식으로 표현되는 근사 변위와 Ritz 방법을 이용하여 Lagrangian범함수의 정상값을 구하였다. 3차원 모델의 정확도는 무차원 진동수의 수렴연구를 통하여 이루어졌으며, 또한 기존 문헌상의 해석 및 실험 결과와의 비교를 통하여 본 연구 결과의 정확성을 검토하였다. 본 논문에서 제시된 3차 원 진동수의 결과를 이용하여 캔틸레버 복합 적층판의 기하학 및 재료 매개변수 즉, 형상비(a/b), 폭두께비(a/h), 재료의 직교 이방성, 플라이 수(NP), 섬유 배향각(), 및 적층 순서가 미치는 효과를 설명하였다. This paper presents the three-dimensional (3-D) study of the natural vibration of cantilevered laminated composite plates. The Ritz method is used to obtain stationary values of the associated Lagrangian functional with displace- ments approximated by mathematically complete polynomials satisfying the boundary conditions at the clamped edge exactly. The accuracy of the 3-D model is established through a convergence study of non-dimensional frequencies followed by a comparison of the converged 3-D solutions with analytical and experimental findings in the existing literature. A wide scope of 3-D frequency results explain the influence of a number of geometrical and material parameters for cantilevered laminated plates, namely aspect ratio (a/b), width-to-thickness ratio (a/h), orthotropy of material, number of plies (NP), fiber orientation angle (), and stacking sequence.

경계 요소법을 이용한 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

본 연구에서는 직선 강상자형 거더의 단면변형에 의한 변형 및 응력계산을 위한 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.

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

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

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

최근, 유한요소해석견과의 신뢰도를 향상시키기 위하여 활발하게 연구되고 있는 적응유한요소해석은 반복계산을 통해서 해석결과의 오차가 사용자에 의해 지정된 허용오차와 같아지도록 하는 해석방법이다. 이와 간은 적응유한요소해석은 해석결과의 오차평가와 이에 따른 유한요소의 재구성과정으로 나누어진다. 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.

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

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

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

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

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

동적 불안정 좌굴현상에 관한 연구는 다소 발표되고 있으나, 주기성을 가진 하중하에서의 동적 좌굴을 다룬 연구는 그리 많지 않은 편이다. 주기성을 가진 하중하에서의 거동은 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).