본 논문에서는 역전파 방법 기반 자동미분법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용 하였다. 응력제한조건이 있는 위상최적화문제는 특이점(singularity)과 응력의 국부성(local nature of stress constraint)문제, 그리고 설 계 변수에 대한 비선형성의 문제를 포함하고 최적해를 얻기가 매우 힘들다. 특이점 문제를 해결하기 위해서 응력 완화(stress relaxation) 기법을 사용하였고, 응력의 국부성을 해결하기 위해 p-norm을 이용한 전역 응력치를 제한조건에 사용하였다. 설계 변수에 대한 비선 형성을 극복하기 위해 해석적인 방법으로 정확한 설계민감도를 구하는 것이 중요하다. 위상최적설계에서 기존에는 보조변수방법 (adjoint variable method)을 사용하여 빠르고 정확한 설계민감도를 구했지만, 설계민감도를 해석적으로 구해야 하고, 보조평형방정식 을 추가로 풀어야 하는 어려움이 있다. 이를 해결하기 위해서 인공신경망에서 최적 가중치(weights)와 편차(biases)를 구할 때 쓰이는 역전파 기법을 이용하여 설계민감도를 구하고 이를 응력제한조건을 고려한 위상최적설계에 적용하였다. 역전파 기법은 자동미분에 쓰이는 기법으로 목적함수나 제한조건에 대한 설계민감도를 별도의 수식유도 없이 간단하게 구할 수 있는 장점이 있다. 또한, 미분값 을 구하는 역전파의 과정이 보조평형방정식을 푸는 것보다 계산시간이 빠르고 해석적 방법으로 구한 설계민감도와 같은 정확도를 보 여준다

본 연구의 대상은 신축이음장치를 설치하지 않고 상부구조와 교대를 일체화하는 무조인트 교량이다. 무조인트 교량은 국내에 2009년 이후 본격적으로 도입되었다. 공용기간이 짧고 설계, 시공 및 유지관리 경험이 부족하여 장기거동에 대한 신뢰도가 아직은 부족하다. 수치해석을 통해 다수의 교량을 분석하는 경우 수치해석 모델은 안정적인 정확성 유지와 모델 구축의 편의성이 확보되어야 한다. 본 연구에서는 다양한 형식을 가진 무조인트 교량의 수치해석 모델을 선정하기 위해 민감도 해석을 수행하였다. 민감도 해석은 상용 유한요소 프로그램인 MIDAS Civil과 ABAQUS를 사용하여 수행하였다. solid 요소 기반인 모델을 기준으로 하여 구조모델간 평균 및 최대 상대오차를 분석하였다. 해석결과 beam 요소 기반인 모델은 상대오차가 크게 발생하였고 shell 요소 기반인 모델은 상대오차가 아주 미소하였다. 따라서 무조인트 교량의 최적 수치해석 모델은 상대오차에 의한 변위 형상의 유사성과 정밀도를 유지하면서 실용적인 모델인 shell 요소 기반 모델이 가장 적합한 것으로 판단하였다.

PURPOSES: In this study, a numerical parametric study was performed to evaluate the effect of angular velocity and weight of wheel, and density of road-bed particles on corrugation development. METHODS : Discrete element method coupled with rigid body dynamics was applied to simulate a wheel-running circular table with variations in independent parameters, such as wheel angular velocity, wheel weight, and particle density. The position profiles for travel distance from origin were compared and analyzed to confirm if the trend from numerical analysis agrees with the analytical solution. RESULTS: The angular velocity of the wheel exhibits a clear inverse relationship with the development of corrugation even though the weight of wheel does not demonstrate clear trends for both long-wave and short-wave corrugation. The density of road-bed particles is observed to have clear proportional effect on corrugation development. The movement of corrugation to the running direction, which was observed in previous research, is also observed for various conditions. CONCLUSIONS : The parametric study using discrete element method with rigid body dynamics clearly exhibits good agreement with analytical solution for initiation of corrugation. The coupled method is confirmed to supply additional information that cannot be delivered by analytical solution only.

본 논문은 아이소-지오메트릭 해석에서 h-세분화를 이용한 국부 세분화법과 이에 따른 설계 민감도 해석의 방법론을 연 구하였다. 다중 조밀도 방식을 이용하여 경계면에서 변위 적합조건을 만족하였고, 기존의 아이소-지오메트릭 해석의 텐서곱 으로 인해 발생하는 원치 않는 자유도 증가의 문제를 극복하였다. 해석에서의 변위 적합조건과 마찬가지로, 설계 민감도 해석에서도 변위 결과와 마찬가지로 똑같은 적합조건을 만족하도록 하는 방법론을 제시하였다. 수치 예제를 통하여 본 방법론의 효율성을 입증하였고, 특별히 응력 집중 문제에서의 결과와 민감도 값을 비교하며 경계면에서의 적합조건을 확인하였다.

국내 고속도로 콘크리트 중앙분리대는 SB5-B(270kJ)의 충돌등급에 저항하도록 설계되어 있다. 그러나 최근 대형 화물차 량의 충돌사고가 지속적으로 증가하는 경향을 보이고 있어 SB6(420kJ) 등급으로의 상향이 필요하다. 충돌등급 상향을 위한 새로운 중앙분리대 단면을 제시하기 위해서는 실제 충돌시험을 수행하여 기준 통과여부를 결정하며, 충돌시험 수행을 위한 적정 단면을 제시하기 위해서는 충돌해석을 통해 선정한다. 이러한 충돌해석의 정확도 향상을 위해서는 차량 모델, 콘크리트 단면 열화상태, 콘크리트 피복 두께 등 다양한 변수에 대한 정확한 변수 선정이 필요하다. 따라서 본 연구에서는 공차 중량, 단면 열화, 콘크리트 피복 두께에 대한 변수연구를 수행하여 충돌저항성능 저감을 확인하였다. 전체 중량뿐만 아닌 공차 중 량에 따라 중앙분리대의 충돌저항성능에 차이가 있는 것으로 확인되었으며, 10cm 이하의 콘크리트 피복 두께에서는 충돌저 항성능이 민감하게 증가 또는 감소한다. 단면 열화가 발생할 경우 중앙분리대의 충돌저항성능의 감소가 발생하여 열화정도 에 따른 보수 및 보강이 이루어져야 하는 것으로 판단된다. 따라서 콘크리트 구조물과 차량의 충돌해석을 수행할 경우 트럭 의 공차중량 비율, 콘크리트의 피복두께 및 열화에 대한 영향이 상세하게 고려될 필요가 있음을 확인하였다.

페리다이나믹스 이론과 이진분해 기법의 병렬연산을 이용하여 동적 균열진전 문제에 대한 애조인 형상 설계민감도 해석법을 개발하였다. 페리다이나믹스에서는 균열의 연속적인 분기를 다룰 수 있으며, Explicit 시간적분법을 채택한다. 설계민감도 해석은 애조인 변수법은 경로의존성 문제에는 적합하지 않으나 여기서는 응답해석의 경로를 이미 알고 있으므로 채택하여 사용할 수 있었다. 얻어진 해석적 설계민감도는 유한차분과 비교하여 그 정확성을 검증하였다. 유한차분법은 설계섭동량에 민감하여 비선형성이 강한 페리다이나믹스 문제에서 부정확한 설계민감도를 제시할 수 있다. 정확한 설계민감도 해석을 위해서는 이산화과정에서 C1 연속성을 가지는 체적율이 필요함을 알 수 있었다.

본 논문에서는 브리징 스케일 분해를 기반으로 멀티스케일 문제에 대한 설계민감도 해석법을 개발하였다. 나노 기술의 급속한 발전으로 인해 나노 수준의 해석의 필요성이 지속적으로 증가하고 있다. 최근 분자동역학과 연속체역학의 연성문제에서 많은 해석 방법들이 개발되었다. 본 논문에서는 연성시스템 해석을 위해 브리징 스케일 기법을 사용한다. 전체 영역의 분자동역학 시스템의 해석은 많은 양의 계산 비용이 들기때문에 분자동역학과 연속체 시뮬레이션의 연성시스템을 선호한다. 분자동역학과 연속체 수준 사이의 정보 교환은 분자동역학과 연속체의 경계에서 일어난다. 브리징 스케일 법에서 일반화된 랑지벵 방정식은 축소된 영역의 분자동역학 시스템 해석을 위하여 요구되고, 시간이력 커널을 사용하여 구한 GLE 힘은 분자동역학 시스템에서 경계에 있는 원자들에 작용한다. 그러므로 분자동역학과 연속체 수준의 시뮬레이션을 분리해서 해석할 수 있으며 계산 과정을 가속시킬 수 있다. 연성문제의 시뮬레이션 이후에는 설계의 최적화를 위해 설계민감도 해석의 필요성이 자연스럽게 나타나며 전체 시스템의 성능은 나노 스케일의 효과를 고려해서 최적화된다. 설계구배 기반 최적화에서 설계민감도가 요구되지만 유한차분법으로 구한 민감도는 문제가 대형화될 때 계산비용의 제한때문에 비실용적이나 해석적 설계민감도는 효율적인 강점을 갖는다. 본 연구에서는 연성된 분자동역학-연속체 멀티스케일 문제에서 해석적 설계민감도를 유도하여 정확성과 향후 최적설계로의 활용 가능성을 확인하였다.

본 논문에서는 아이소-지오메트릭 기법을 기반으로 민들린 후판에 대한 형상 설계민감도 해석법을 제시하였다. 아이소-지오메트릭 기법은 정확한 기하학적 형상의 표현, 요소 사이의 높은 연속성 등 바람직한 강점들을 가지고 있으며 궁극적으로는 해석해로의 빠른 수렴성과 정확한 설계민감도를 제공한다. 선형 형상함수를 사용하는 유한요소법과는 달리 아이소-지오메트릭 기법에서는 높은 차수의 NURBS 기저함수를 활용하여 CAD 형상의 법선벡터와 곡률을 정확하게 고려한다. 전단 잠김(Shear locking) 현상을 극복하기 위해서 선택적 감소적분(Selective reduced integration) 기법을 사용하였다. 이 간단한 방법은 복잡한 정식화 과정 없이 정확한 아이소-지오메트릭 형상 설계민감도 해석을 수행한다. 굽힘 문제에 대한 수치예제를 통하여 제안된 아이소-지오메트릭 해석과 유한요소 해석을 비교하였으며, 유한차분 설계민감도와 비교하여 아이소-지오메트릭 형상 설계민감도는 매우 정확함을 확인하였다.

Finite element analysis is to approximate a geometry model developed in computer-aided design(CAD) to a finite element model, thus the conventional shape design sensitivity analysis and optimization using the finite element method have some difficulties in the parameterization of geometry. However, isogeometric analysis is to build a geometry model and directly use the functions describing the geometry in analysis. Therefore, the geometric properties can be embedded in the NURBS basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. In this study, the isogeometric structural analysis and shape design sensitivity analysis in the generalized curvilinear coordinate(GCC) systems are discussed for the curved geometry. Representing the higher order geometric information, such as normal, tangent and curvature, yields the isogeometric approach to be the best way for generating exact GCC systems from a given CAD geometry. The developed GCC isogeometric structural analysis and shape design sensitivity analysis are verified to show better accuracy and faster convergency by comparing with the results obtained from the conventional isogeometric method.

최근 구조물의 국부적인 손상이 전체적인 붕괴로 이어지는 연쇄붕괴 현상에 관한 연구가 활발히 진행되고 있다. 연쇄붕괴에 관한 기존의 연구는 대부분 해석변수의 불확실성을 포함하지 않는 확정론적인 방법이므로, 해석결과에 대한 신뢰도를 알 수 없다. 본 논문에서는 재료의 항복강도, 활하중의 크기, 감쇠비, 탄성계수 등치 설계변수들이 기둥이 제거됨에 따라 발생하는 수직변위에 영향을 미치는 민감도를 분석하였다. 이를 위하여 몬테카를로 시뮬레이션, 일계이차법, 토네이도 다이어그램의 세 가지 해석기법을 적용하였다. 비선형정적 해석결과에 의하면 난수로 설정한 해석변수들 중에서 보의 항복강도가 수직변위의 변동폭이 가장 컸으며, 비선형동적해석의 경우 보의 항복강도와 감쇠비가 서로 유사한 변동폭을 가지는 것으로 나타났다.

본 논문에서는 등기하 해석법을 이용하여 평면 탄성문제의 변분식을 유도하였다. 등기하 해석법은 새로이 부각되고 있는 해석법으로서 기저 함수가 NURBS(Non-Uniform Rational B-Splines) 로부터 직접 생성되므로 해 공간은 CAD 모델을 구성하는 함수로써 표현된다. 또한 CAD 모델의 B-Spline 기저 함수를 직접 사용하므로 기하학적으로 엄밀한 형상을 표현할 수 있고 요소망의 재구성 없이 해석모델을 정밀화(Refinement)할 수 있는 강점이 있다. 본 논문에서는 이를 확장하여 연속체 기반의 애드조인트 설계 민감도 해석법을 사용하는 등기하 설계민감도 해석법을 유도하였다. 기존의 유한요소 기반형상 최적설계는 형상의 매개화에 어려움을 겪었으나 등기하 기반 최적설계에서는 기하학적 정보가 이미 B-spline 기저함수와 조정점에 포함되어 있으므로 이러한 어려움을 피할 수 있는 잠재력을 가지고 있다. 몇몇 수치 예제를 통해서 등기하 해석법을 사용한 설계 민감도 해석을 수행하였으며 유한차분 민감도와 비교하여 정확성을 확인하였다.

본 연구에서는 방사성폐기물의 화학처리공정에 자주 사용되는 유동관식 장치 중 튜브형 반응기, 다단식 용매추출 장치, 흡착탑 등 물질전달이 수반되는 장치에 있어 각종 매개변수들이 반응수율이나 물질전달수율에 미치는 영향과 민감도를 살펴보았다. 먼저 각 장치에 대한 거동을 묘사하기 위하여 수학적 모델링을 수행하였고 전산모사를 통하여 해당 장치의 거동을 예측하였다. 그리고 그 결과로부터 해당 공정의 고유한 매개변수들이 반응수율 또는 물질전달수율에 미치는 영향과 민감도를 분석하였다. 튜브형 반응기에서는 확산계수, 반응속도상수 등이 반응수율에 미치는 영향을, 다단식 용매추출 장치에서는 연속상과 분산상의 분배계수, 연속상 흐름의 역혼합 등이 추출수율 및 장치 내 농도 분포에 미치는 영향을 고찰하였다. 또 흡착탑에 있어서는 흡착평형상수 및 유체-흡착재간 물질전달계수 등이 흡착 속도에 미치는 영향을 조사하였다.

본 연구에서는 지진하중을 받는 고층 RC 골조구조물의 횡변위를 정량적으로 제어할 수 있는 방안을 제시한다. 이를 위해 수학적인 일반성을 가지면서 큰 규모의 문제도 효율적으로 다룰 수 있는 근사화 개념을 도입하여 횡변위 구속조건식을 설정한다. 아울러 구조부재의 단면특성 관계식을 설정함으로써 설계변수의 수를 줄여주고, 초기에 주어진 단면형상이 최적설계 과정동안 계속 유지된다는 가정을 이용하여 최적설계결과에서 구해진 단면특성에 따라 부재단면크기를 산출하는 방안을 강구한다. 특히 근사화된 횡변위구속조건식을 정식화 하기 위해 동적 변위민감도해석 방안이 고려된다. 이와 같이 제시된 동적 강성최적설계 기법의 효용성을 검토하기 위해 10층과 50층 규모의 삼차원 RC 골조구조물 모델이 고려된다.

본 연구에서는 민감도 해석을 이용하여 전단벽-골조 구조시스템의 횡변위를 정량적으로 제어할 수 있는 강성최적설계방안을 제시하고자 한다. 이를 위해 먼저 골조와 전단벽요소 사이의 변위자유도 적합성 문제를 해결하기 위한 요소강성행렬을 구성하며, 또한 수학적계획법의 일반성을 유지하면서도 큰 규모의 문제도 효율적으로 다를 수 있는 근사화 재념을 도입하여 횡변위 구속조건식을 설정한다. 아울러 전단벽 및 골조부재의 단면특성 관계식을 설정함으로써 설계변수의 수를 줄여주고, 이를 이용하여 강성행렬도함수의 산정을 용이하게 한다. 특히 골조의 경우 초기에 주어진 단면형상이 최적설계 과정동안 계속 유지된다는 가정을 이용하여 최적설계결과에서 구해진 단면특성에 따라 부재단면크기를 산출하고, 전단벽은 사용자의 의도에 따라 두께 또는 부재길이를 재산정하는 방안을 강구한다. 이와 같이 제시된 강성최적설계기법의 효용성을 검토하기 위해 두 가지 형태의 20층 전단벽-골조 구조물의 예제가 고려된다.

본 논문에서는 다양한 문제들의 형상 설계 민감도 해석에 대한 효율적인 경계기반 기법을 제시하였다 우선 문제에서 정의된 일반적인 함수들에 대한 연속체 형태의 식에 근거하여, 경계 적분 형태의 해석적 민감도 식을 유도하였다. 이 식은 다양한 형상 설계 문제들의 경사를 계산하는데 편리하게 사용할 수 있다. 그리고 경계법은 형상 변분 벡터가 전체 도메인이 아닌 경계에서만 요구된다는 장점이 있는데, 여기서 경계 형상 변분은 형상 함수의 복잡한 해석적 미분 대신 형상을 미소 증분시킴으로써 편리하게 계산할 수 있다. 제시한 방법의 효율성을 보이기 위해 포텐셜 유동 문제와 필렛(fillet)에서의 응력 집중 문제에 이를 적용하였다.

평면 아치 구조물에 대해 선형 탄성 변분방정식에 기반을 둔 설계민감도 해석을 위한 일반적 이론을 개발하였다. 아치 구조물내의 임의 마디에 정의된 응력범함수를 고려하였고 이에 대한 설계민감도 공식을 유도하기 위해 전미분(material derivative) 개념과 보조(adjoint) 변수 방법을 도입하였다. 얻어진 민감도 공식은 구조해석 결과를 얻고 나면 이들로부터 단순 대수연산을 통해 계산이 되므로 적용이 간편할 뿐 아니라 해의 정확도가 높은 잇점이 있다. 본 방법은 아치의 형상을 매개변수를 통해 표현하므로 얕은 아치에 국한하지 않고 어떠한 형상도 고려가 가능하며, 나아가서 아치의 형상변화를 형상에 대해 수직뿐 아니라 접선방향도 포함하여 일반적으로 고려하므로 다양한 형상설계가 가능하다. 몇 가지 예제에서 민감도 계산을 수행함으로써 본 방법의 정확도와 효율성을 입증하였으며, 두 가지의 설계최적화 문제를 대상으로 실제로 두께 및 형상최적설계를 수행하였다.

It is main objective of this approach to present a method to analyse stochastic design sensitivity for problems of structural dynamics with randomness in design parameters. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element approach. An alternative form of the constant functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The terminal problem of the adjoint system is solved using equivalent homogeneous equations excited by initial velocities. The numerical procedures are shown to be much more efficient when based on the fold superposition method: the generalized co-ordinates are normalized and the correlated random variables are transformed to uncorrelated variables, whereas the secularities are eliminated by the fast Fourier transform of complex valued sequences. Numerical algorithms have been worked out and proved to be accurate and efficient : they can be readily adapted to fit into the existing finite element codes whose element derivative matrices can be explicitly generated. The numerical results of two cases -2 dimensional portal frame for the comparison with reference and 3-dimensional frame structure - for the deterministic sensitivity analysis are presented.

본 연구에서는 중복 고유치를 갖는 비비례 감쇠 진동계의 고유치와 고유벡터의 민감도를 계산하는 새로운 방법을 제시하였다. 제안 방법은 매우 간단하면서도 수치적 안정성이 보장되고 정확한 해를 주는 방법이다. 제안 방법에서는 (n+m)차의 대칭 행렬로 이루어진 대수방정식을 해석함으로써 n개의 자유도를 갖는 감쇠계에 있어서 m차의 중복도를 갖는 고유치와 고유벡터의 설계변수에 대한 미분을 구한다. 제안 방법의 검증을 위해 5자유도를 갖는 단순구조물의 민감도해석을 예제에서 다루고 있다. 예제에서의 설계변수는 모델의 부분강성으로 하였다.

본 연구에서는 중복되지 않는 고유치를 갖는 비비례 감쇠계의 고유치와 고유벡터의 민감도를 계산하는 새로운 방법을 제시하였다. 제안 방법에서는 (n+1)차의 대칭 행렬로 이루어진 대수방정식을 해석함으로써 n개의 자유도를 갖는 감쇠계의 고유치와 고유벡터의 설계변수에 대한 미분을 구한다. 제안 방법은 매우 간단하면서도 수치적 안정성이 보장되고 정확한 해를 주는 방법이다. 제안 방법의 검증을 위해 7자유도를 갖는 차량모델의 민감도해석을 예제에서 다루고 있다. 예제에서의 설계변수는 콘테이너의 질량으로 하였다.