• 제목/요약/키워드: Non-linear Numerical model

검색결과 423건 처리시간 0.026초

2차원 멀티패치 모델의 아이소-지오메트릭 해석 (Isogeometric Analysis for Two-dimensional Multipatch Model)

  • 김민근;구본용
    • 한국전산구조공학회논문집
    • /
    • 제30권6호
    • /
    • pp.515-522
    • /
    • 2017
  • 본 논문에서는 2개 이상의 기하형상이 순응 또는 비순응 경계면에서 접합된 멀티패치 문제에 대한 아이소-지오메트릭 해석에 대해서 연구하였다. 패치 경계면에서 응력의 연속성을 표현하는 방법으로 Nitsche 방법론과 마스터-슬레이브 방법에 기반한 방법론에 대해서 지배방정식을 유도하고 아이소-지오메트릭 이산화를 수행하였다. 멀티패치 문제에 대해서 두 방법론의 차이점을 간단하게 비교하였으며, 후처리 과정에서 사용되는 NURBS 곡면 기반의 응력 복원법에 대해서 기술하였다. 수치예제에서 비순응 경계면을 가지는 멀티패치 빔 문제를 통해 Nitshce 방법론을 검증하였으며, 응력집중을 가지는 문제에서 소개된 두 방법론이 유사한 결과를 보이는 것을 확인하였다. 소개된 NURBS 곡면 기반의 응력 복원법을 후처리에서 도입할 경우 멀티패치 문제의 경계면에서 개선된 연속적인 응력을 보임을 알 수 있다.

PFC를 이용한 입자와 비파쇄 평면과의 접촉면에서의 전단 거동 (Shear behavior at the interface between particle and non-crushing surface by using PFC)

  • 김은경;이정학;이석원
    • 한국터널지하공간학회 논문집
    • /
    • 제14권4호
    • /
    • pp.293-308
    • /
    • 2012
  • 암반 절리면과 같이 입자와 연속체 평면의 접촉면에서의 전단거동은 전체 구조물의 거동을 지배할 수 있다. 암반설계의 효율을 높이기 위해서는 입자와 연속체 평면의 접촉면 전단거동 메커니즘에 대한 기초적인 이해와 접촉면 전단강도를 정확하게 산정하는 것이 필수적이다. 본 연구에서는 연속체 평면의 표면 거칠기의 영향을 알아보기 위하여 개별요소법 수치해석 프로그램인 $PFC^{2D}$를 사용하였다. 이 때 표면 거칠기는 매끄러운 평면, 중간 거칠기 평면, 거친 평면의 세 가지로 구분하였다. 접촉면 전단시험에서 입자 형상 및 입자 파쇄의 영향을 연구하기 위하여 one ball 모델과 clump 및 cluster 모델을 형성하여 상대적으로 비교 분석하였다. 이 때 입자의 형상은 원형, 삼각형, 직사각형, 정사각형으로 모델링하였다. 수치해석 결과, 표면 거칠기가 클수록 접촉면 전단강도 및 마찰각이 크게 나타났으며, 또한 간극률이 크게 증가하였다. 표면이 매끄러운 one ball 모델보다 작은 입자들의 결합으로 이루어져 표면이 굴곡진 모델, 즉 clump 모델의 접촉면 전단강도 및 마찰각이 크게 나타났다. 입자의 결합이 파괴되는 모델, 즉 cluster 모델의 접촉면 전단 강도 및 마찰각이 같은 형상의 clump 모델보다 작게 나타났으며, 파괴포락선은 비선형으로 나타났다. 이러한 결과로부터 연속체 평면의 거칠기 및 입자의 형상이 입자와 평면의 접촉면 전단거동 특성에 미치는 영향을 확인하였다.

Inelastic behavior of systems with flexible base

  • Fernandez-Sola, Luciano R.;Huerta-E catl, Juan E.
    • Earthquakes and Structures
    • /
    • 제14권5호
    • /
    • pp.411-424
    • /
    • 2018
  • This study explores the inelastic behavior of systems with flexible base. The use of a single degree of freedom system (ESDOF) with equivalent ductility to represent the response of flexible base systems is discussed. Two different equations to compute equivalent ductility are proposed, one which includes the contribution of rigid body components, and other based on the overstrength of the structure. In order to asses the accuracy of ESDOF approach with the proposed equations, the behavior of a 10-story regular building with reinforced concrete (RC) moment resisting frames is studied. Local and global ductility capacity and demands are used to study the modifications introduced by base flexibility. Three soil types are considered with shear wave velocities of 70, 100 and 250 m/s. Soil-foundation stiffness is included with a set of springs on the base (impedance functions). Capacity curves of the building are computed with pushover analysis. In addition, non linear time history analysis are used to asses the ductility demands. Results show that ductility capacity of the soil-structure system including rigid body components is reduced. Base flexibility does not modify neither yield and maximum base shear. Equivalent ductility estimated with the proposed equations is fits better the results of the numerical model than the one considering elastoplastic behavior. Modification of beams ductility demand due to base flexibility are not constant within the structure. Some elements experience reduced ductility demands while other elements experience increments when flexible base is considered. Soil structure interaction produces changes in the relation between yield strength reduction factor and structure ductility demand. These changes are dependent on the spectral shape and the period of the system with fixed and flexible base.

RC교각을 이용한 실시간 하이브리드 실험 시스템의 적용성 연구 (Verification of Real-time Hybrid Test System using RC Pier Model)

  • 이진행;박민석;채윤병;김철영
    • 한국지진공학회논문집
    • /
    • 제22권4호
    • /
    • pp.253-259
    • /
    • 2018
  • Structure behaviors resulting from an earthquake are experimentally simulated mainly through a shaking table test. As for large-scale structures, however, size effects over a miniature may make it difficult to assess actual behaviors properly. To address this problem, research on the hybrid simulation is being conducted actively. This method is to implement numerical analysis on framework members that affect the general behavior of the structure dominantly through an actual scale experiment and on the rest parts by applying the substructuring technique. However, existing studies on hybrid simulation focus mainly on Slow experimental methods, which are disadvantageous in that it is unable to assess behaviors close to the actual level if material properties change depending on the speed or the influence of inertial force is significant. The present study aims to establish a Real-time hybrid simulation system capable of excitation based on the actual time history and to verify its performance and applicability. The hybrid simulation system built up in this study utilizes the ATS Compensator system, CR integrator, etc. in order to make the target displacement the same with the measured displacement on the basis of MATLAB/Simulink. The target structure was a 2-span bridge and an RC pier to support it was produced as an experimental model in order for the shaking table test and Slow and Real-time hybrid simulations. Behaviors that result from the earthquake of El Centro were examined, and the results were analyzed comparatively. In comparison with the results of the shaking table test, the Real-time hybrid simulation produced more similar maximum displacement and vibration behaviors than the Slow hybrid simulation. Hence, it is thought that the Real-time hybrid simulation proposed in this study can be utilized usefully in seismic capacity assessment of structural systems such as RC pier that are highly non-linear and time-dependent.

제방 성토시 지반거동에 대한 점토의 비선형 투수성의 영향 (Effect of the Non-linear Permeability of Clays on the Behavior of Soils in Embankment Construction)

  • 김태훈;한태곤;류기정;이송
    • 한국지반공학회논문집
    • /
    • 제20권6호
    • /
    • pp.61-73
    • /
    • 2004
  • 일반 자연 점토지반의 투수계수는 다양한 응력수준에 따라 비선형적인 특성을 보이며, 이러한 지반 투수계수의 비선형성은 상재하중 재하시 지반의 압밀시간 및 지반 변형 특성에 영향을 주게 된다. 따라서, 본 연구에서는 지반의 비선형 투수계수의 영향성을 고찰하기 위해 수정 Cam-clay 모델과 Biot의 압밀이론을 접목한 연립이론(Coupled theory)을 이용하여 FEM 수치해석을 실시하였다. 검토 결과 비선형 투수성은 비선형 특성을 나타내는 투수계수 변화비와 전체 간극비에 대한 평균 투수계수의 크기에 따라 지반 변형과 지중 과잉간극수압에 미치는 영향성이 상이하게 나타났다. 따라서, 보다 현장에 근접한 해석을 실시하기 위해서는 지반 투수계수의 변화비와 평균 투수계수의 크기에 따라 투수계수의 비선형 특성을 고려해야할 것인지를 판단해야 할 것이다.

목포해역 낙조류 우세현상의 수치모의 (Numerical Modeling of Ebb-Dominant Tidal Flow in the Mokpo Coastal Zone)

  • 정태성;최종화
    • 한국해안·해양공학회논문집
    • /
    • 제22권5호
    • /
    • pp.333-343
    • /
    • 2010
  • 목포해역에서 조위와 조류 관측자료를 분석하여, 낙조류 우세 조류특성을 보이는 것을 확인하였으며 낙조류 우세를 발생하는 물리적 요인에 대해 검토하였다. 노출한계수심, 바닥마찰응력 산정방법, 해안 매립, 조석의 진폭, 비선형 조석, 와점성계수와 목포해역의 낙조류 우세현상의 관련성에 대해서 2차원 해수유동 모형을 적용하여 자세히 조사하였다. 다양한 조건에 대한 모의결과로부터 노출한계수심과 와점성계수는 모의결과에 거의 영향을 주지 않으며, 바닥마찰응력이 증가함에 따라 조류의 비선형성이 증가하여 낙조류 우세현상이 나타나는 것을 확인할 수 있었다. 간사지가 미치는 영향은 간사지가 육지로 바뀌면 낙조류 우세가 심해지고, 간사지가 바다로 바뀌면 낙조류 우세가 사라지는 것을 확인할 수 있었다. 그동안 목포해역에서 진행되어온 해안매립사업은 낙조우세현상을 심화시켜 왔음을 확인할 수 있었으며, 비선형 조석은 비대칭 조류를 생성하는 데 결정적인 역할을 하여 조류의 비대칭성인 강한 목포해역에서 조류를 모의할 때는 개방경계 조위에 비선형조석을 포함해야 함을 확인할 수 있었다.

기계학습 기반의 메타모델을 활용한 ZnO 바리스터 소결 공정 최적화 연구 (Sintering process optimization of ZnO varistor materials by machine learning based metamodel)

  • 김보열;서가원;하만진;홍연우;정찬엽
    • 한국결정성장학회지
    • /
    • 제31권6호
    • /
    • pp.258-263
    • /
    • 2021
  • ZnO 바리스터는 다결정구조를 가지는 반도체 소자로 결정립과 입계의 미세구조 제어를 통해 비선형적인 전류/전압 특성을 가지기 때문에 서지(surge)전압으로부터 회로를 보호하는 역할을 한다. 이러한 ZnO 바리스터에서 원하는 전기적 물성을 얻기 위해서는 소결 공정에서 미세구조의 제어가 중요하다. 따라서 소결 공정에서 중요한 변수들과 소결체의 전기적 물성인 유전율로 구성된 데이터셋을 정의한 후 실험계획법 기반으로 데이터를 수집했다. 수집된 실험데이터셋을 기계학습 알고리즘에 학습하여 메타모델을 개발했고, 개발된 메타모델에 수치기반 최적화 알고리즘인 HMA(Hybrid Metaheuristic Algorithm)를 적용하여 최대 유전율을 가질 수 있는 공정조건을 도출했다. 이러한 메타모델 기반의 최적화를 다변수 시스템인 세라믹공정에 적용한다면 최소한의 실험만으로 최적 공정조건 탐색이 가능할 것으로 판단된다.

Direct numerical simulations of viscoelastic turbulent channel flows at high drag reduction

  • Housiadas Kostas D.;Beris Antony N.
    • Korea-Australia Rheology Journal
    • /
    • 제17권3호
    • /
    • pp.131-140
    • /
    • 2005
  • In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

항공 라이다 자료를 이용한 토석류 발생지역의 지형복원기법 개발 (Development of the Topography Restoration Method for Debris Flow Area Using Airborne LiDAR Data)

  • 우충식;윤호중;이창우;이규성
    • 한국지리정보학회지
    • /
    • 제14권3호
    • /
    • pp.174-187
    • /
    • 2011
  • 항공 LiDAR 측량으로 토석류 발생 전 후의 지형자료를 취득하는 경우 토석류로 인하여 유출된 토사량을 알 수 있다. 그러나 토석류 발생지를 미리 예측하여 촬영하기가 힘들고, 토석류 발생 지역의 과거 항공 LiDAR 자료는 존재가능성이 낮아 토석류 발생이전 지형자료를 이용하는 것은 어렵다. 따라서 본 연구에서는 토석류 발생지역의 토사량 추정을 위해 발생전 지형을 복원하고, 토사유출의 공간적 범위를 파악할 수 있는 지형복원기법을 개발하였다. 지형복원기법은 토석류 발생지역에서 추출한 선형 및 비선형 횡단면을 가우시안혼합모델로 수식화하고 중심점 추정방법과 근사정확도로 근사결과를 평가하여 토석류 발생이전의 지형을 복원한다. 지형복원기법은 토석류 발생 전 후의 항공 LiDAR 자료를 이용하여 두 가지 방법으로 검증하였다. 먼저 토석류 발생구간에서 추출한 각 횡단면을 지형복원하여 발생전 항공 LiDAR 자료와 비교하였다. 또한 토석류 발생지역에 지형복원기법을 적용한 뒤 지형자료를 제작하여 토석류 발생전 항공 LiDAR DEM과 비교하여 검증하였다. 지형복원기법의 검증한 결과 전반적으로 근사정확도가 0.5m에 가까운 높은 정확도를 나타냈다.

Behaviour of steel-fibre-reinforced concrete beams under high-rate loading

  • Behinaein, Pegah;Cotsovos, Demetrios M.;Abbas, Ali A.
    • Computers and Concrete
    • /
    • 제22권3호
    • /
    • pp.337-353
    • /
    • 2018
  • The present study focuses on examining the structural behaviour of steel-fibre-reinforced concrete (SFRC) beams under high rates of loading largely associated with impact problems. Fibres are added to the concrete mix to enhance ductility and energy absorption, which is important for impact-resistant design. A simple, yet practical non-linear finite-element analysis (NLFEA) model was used in the present study. Experimental static and impact tests were also carried out on beams spanning 1.3 meter with weights dropped from heights of 1.5 m and 2.5 m, respectively. The numerical model realistically describes the fully-brittle tensile behaviour of plain concrete as well as the contribution of steel fibres to the post-cracking response (the latter was allowed for by conveniently adjusting the constitutive relations for plain concrete, mainly in uniaxial tension). Suitable material relations (describing compression, tension and shear) were selected for SFRC and incorporated into ABAQUS software Brittle Cracking concrete model. A more complex model (i.e., the Damaged Plasticity concrete model in ABAQUS) was also considered and it was found that the seemingly simple (but fundamental) Brittle Cracking model yielded reliable results. Published data obtained from drop-weight experimental tests on RC and SFRC beams indicates that there is an increase in the maximum load recorded (compared to the corresponding static one) and a reduction in the portion of the beam span reacting to the impact load. However, there is considerable scatter and the specimens were often tested to complete destruction and thus yielding post-failure characteristics of little design value and making it difficult to pinpoint the actual load-carrying capacity and identify the associated true ultimate limit state (ULS). To address this, dynamic NLFEA was employed and the impact load applied was reduced gradually and applied in pulses to pinpoint the actual failure point. Different case studies were considered covering impact loading responses at both the material and structural levels as well as comparisons between RC and SFRC specimens. Steel fibres were found to increase the load-carrying capacity and deformability by offering better control over the cracking process concrete undergoes and allowing the impact energy to be absorbed more effectively compared to conventional RC members. This is useful for impact-resistant design of SFRC beams.