• 제목/요약/키워드: Lagrangian equation

검색결과 212건 처리시간 0.022초

유체-구조 연성해석 기반 해저케이블 위해인자의 수중낙하 특성 비교 (Comparison of Underwater Drop Characteristics for Hazard Apparatuses on Subsea Cable Using Fluid-Structure Interaction Analysis)

  • 장경호;김정훈;송창용
    • 한국해양공학회지
    • /
    • 제32권5호
    • /
    • pp.324-332
    • /
    • 2018
  • It is known that damages to the subsea cables used for electric power transmission between islands and countries, including renewable energy from offshore wind power, current, tides, etc., cost much to restore, which causes social and economic losses. Various types of fishing rigs and anchors have been reported to be the greatest hazards to subsea cables. It is possible to design and construct a suitable protection facility for a subsea cable by precisely estimating the underwater behavior of such hazardous apparatuses. In this study, numerical simulations of the underwater behaviors of various hazardous apparatuses were carried out using fluid-structure interaction (FSI) analysis as a basic study to simulate the actual behavior phenomena of hazardous apparatuses in relation to a subsea cable. In addition, the underwater drop characteristics according to the types of hazardous apparatuses were compared. In order to verify the accuracy of the FSI analysis method used in this study, we compared the test results for underwater drops of a steel ball bearing. Stock anchors, stockless anchors, and rocket piles, which were actually reported to be the cases of damage to subsea cables along the southwest coast of Korea, were considered as the hazardous apparatuses for the numerical simulations. Each hazardous apparatus was generated by a Lagrangian model and coupled with the fluid domain idealized by the Eulerian equation to construct the three-dimensional FSI analysis model. The accuracy of the numerical simulation results was verified by comparing them with the analytical solutions, and the underwater drop characteristics according to the types of hazard apparatuses were compared.

유한요소법을 이용한 방향성 시추의 굴진율 연구 (A study on the excavation rate of directional drilling using finite element method)

  • 정태준;신영기
    • 플랜트 저널
    • /
    • 제17권3호
    • /
    • pp.42-46
    • /
    • 2021
  • 드릴 스트링의 진동 문제는 수년간 드릴링 성능 저하의 주요 원인 중 하나로 인식되어 왔으며, 굴착 작업 시 발생하는 과도한 진동은 드릴링의 효율성, 파이프 피로, 비트의 수명 단축으로 인하여 고장을 초래할 수 있다. 이러한 진동의 원인은 드릴 스트링이 굴착 작업 중 궤적에 따라 파이프의 굽힘과 wellbore와의 접촉으로 인해 마찰이 발생하고 이러한 진동은 일반적으로 축 방향, 굴곡 및 비틀림 변형을 일으킨다. 본 연구에서는 Khulief와 Al-Naser가 제시한 모델을 바탕으로 드릴 스트링에 6자유도(DOF)의 구성요소를 갖는 모델을 적용하여 curved beam의 수치해석 값과 Analytical값을 비교하여 검증하고 드릴 스트링에 hookload와 WOB 경계조건을 주어 각 element 마다 동적 거동을 분석하였다. 실제 궤적을 적용하여 드릴 파이프의 굽힘이나 중력으로 인하여 Wellbore와 접촉되는 부분에 마찰을 적용하였고, 또한 마찰 작용 시 일정한 축 방향 속도를 유지하기 위한 PI제어 값을 설계하여 drillstring 전체의 각속도 변화와 실제 드릴 굴착 작업 중 발생하는 stick slip현상을 관찰하였다.

지반조사를 위한 loop-loop 전자탐사 기법의 적용 (An Application of loop-loop EM Method for Geotechnical Survey)

  • 유진상;송윤호;설순지;송영수
    • 지구물리와물리탐사
    • /
    • 제4권2호
    • /
    • pp.25-33
    • /
    • 2001
  • 토목이나 건설 그리고 환경 분야에서 요구하는 지반 조사에 대한 공학적인 적용을 위하여 주파수 영역 루프-루프 전자탐사 연구를 수행하였다. 이를 위하여 송신 루프와 수신 루프가 동일면을 이루는 수평동일면 배열과 수직동일면 배열을 이용하였고, 각 배열에 대해 세 개의 이상성분을 측정 자료로 하여 영상을 구성하였다. 본 연구에서의 모형 반응 계산은 확장된 Born 근사를 이용한 2.5차원 적분방정식을 이용하였고, 역산 과정에서 분해능의 향상을 위해 ACB(Active Constraint Balancing)를 채택하였다. 모형 반응을 통한 1차원 및 2차원 역산 알고리듬을 적용시켜 본 결과, 층서구조에서는 비교적 1층의 전기전도도가 높을수록 모델 변수의 추정이 용이하였으며, 2차원 고립 이상체에 대한 역산 결과는 전도성 이상체와 비전도성 이상체의 위치를 잘 확인할 수 있었다. 또한 VCP배열 자료의 역산 결과보다 HCP배열 자료의 역산 결과에서 보다 나은 해상도를 보였으며, HCP배열과 VCP배열의 자료를 동시에 역산하였을 경우에 분해능의 향상을 확인하였다. 실제 루프-루프 전자탐사의 현장 자료 획득을 위해서 캐나다 Geonics사의 전자탐사 장비 EM34-3XL을 이용하였으며, 전기비저항 탐사 자료와 동일 측선상에서 비교 분석하였다. 역산 결과, 1차원 역산 보다는 2차원 역산 알고리듬을 이용하여 구성한 영상에서 전기비저항 탐사 결과와 매우 유사한 만족할만한 전기비저항 분포를 확인할 수 있었다. 따라서 현장에서의 개략적인 지반 조사를 위한 루프-루프 전자탐사의 응용이 기대된다.

  • PDF

위상민감도를 이용한 선형구조물의 레벨셋 기반 형상 최적설계 (Level Set Based Shape Optimization of Linear Structures using Topological Derivatives)

  • 윤민호;하승현;김민근;조선호
    • 한국전산구조공학회논문집
    • /
    • 제27권1호
    • /
    • pp.9-16
    • /
    • 2014
  • 레벨셋 기법과 위상민감도를 이용하여 선형 탄성 구조물에 대하여, 초기 설계형상에 의존성이 없는 위상 및 형상 최적설계 기법을 개발하였다. 레벨셋 기법에서는 복잡한 위상 형상변화를 쉽게 다루기 위해 초기 영역은 고정한 채 레벨셋 함수로 표현되는 암시적 이동경계로 경계를 표현한다. 해밀턴-자코비(H-J) 방정식과 수치적으로 강건한 기법인 'up-wind scheme'은 컴플라이언스 목적함수를 최소화시키고 허용체적 제약조건을 만족시키면서, 초기 암시적 경계를 법선 속도장에 따라 최적의 형상으로 이끌어 낸다. 점근적인 정규화 개념에 근거하여, 구멍의 반지름을 0으로 접근시켜 형상 미분의 극한을 취한 위상민감도를 고려하였다. 최적조건으로부터 유도된 라그란지안의 감소 방향을 이용하여 H-J 방정식을 갱신하기 위한 속도장을 결정하였다. 개발한 방법에서는 위상민감도로부터 얻어지는 지표를 이용하여 구멍을 언제든지 어디에서나 생성가능하기 때문에 초기 구멍이 최적 형상을 얻기 위해 요구되지 않는다는 사실을 확인하였다. 또한 효율적인 최적화 과정을 위해서는 구멍 생성을 위한 조정변수의 적절한 선택이 중요함을 확인하였다.

저층수 배사관 내 유입된 사석 배출능력에 대한 연구 (A Experimental Study on Exclusion Ability of Riprap into Bypass Pipe)

  • 정석일;이승오
    • 대한토목학회논문집
    • /
    • 제37권1호
    • /
    • pp.239-246
    • /
    • 2017
  • 국내 중소하천의 횡단구조물인 보 또는 낙차공은 대부분 고정식 콘크리트 구조물이며, 저층수의 배제가 쉽지 않다. 횡단구조물로 인해 유사가 퇴적되며, 유사에 흡착한 오염물들이 그대로 하천의 바닥을 오염시키고 있다. 이에 저층수 및 퇴적유사에 대한 관심이 증가하고 있는 실정이며, 이러한 대안의 하나로 횡단구조물 상류와 하류를 하상 아래로 연결시키는 구조물로써, 보 상류 저층의 물 및 유사 배제를 목적으로 저층수배출관을 설치하는 방안이 있다. 그러나 사석이 유입되고, 배제 되지 않을 경우 효율성이 크게 저하될 가능성이 있다. 이에 본 연구에서는 저층수 배출관 내 사석을 배제할 수 있는 능력에 대한 연구를 수행하였다. 사석과 거동이 유사한 유사(sediment)의 한계조건(critical condition) 중 한계전단력(critical shear stress) 유도과정과 달랑베르의 원리(d'Alembert principle)를 응용하여 이동 중인 사석이 배제될 수 있는 조건(${\tau}_c{^*}$)을 유도하였다. 그러나 저층수 배출관 내 유입된 사석은 정지상태가 아닌 이동 중이므로, Lagrangian 기법을 활용하여 수리실험에서 도출된 유속으로 상대속도(relative velocity)를 제시하였다. 수리실험은 축척효과(scale effect)를 최소화하기 위해 폭이 5.0 m이고, 높이가 1.0 m인 광폭 개수로를 제작하였으며, 사용된 사석은 가공된 완전 구형을 사용하였다. 실험 결과 유속과 구형 입자 속도와의 비가 0.5~0.7 사이로 나타났으며, 이러한 결과를 유도된 식에 적용하여, 최종적으로 사석이 배제되는 조건을 도출하게 되었다. 구간은 입자레이놀즈수($Re_p$)와 무차원 한계 전단력(${\tau}_c{^*}$)에 따라 크게 3가지로 구분되었다. 배제 구간(exclusion section), 확률적 배제 구간(probabilistic exclusion section), 비배제 구간(no exclusion section)이다. 본 연구결과는 횡단구조물의 저층수 배출관 설계시 유용한 기초 정보를 제공할 수 있을 것이다.

유수단면 흐름장에서 Shear 효과를 갖는 Stokes Drift의 계산법 (New Method of Computing the Stokes Drift Including Shear Effect in the Cross-Sectional Flow Field)

  • 김종화;박병수
    • 수산해양기술연구
    • /
    • 제33권1호
    • /
    • pp.9-26
    • /
    • 1997
  • Stokes drift(SD) and Lagrangian discharge(LD) are important factors for analysis of flushing time, tidal exchange, solute transport and pollutant dispersion. The factors should be calculated using the approached method to flow phenomena. The aim of this paper re-examines the previous procedures for computing the SD and LD, and is to propose the new method approached to stratified flow field in the cross-section of coastal region, e.g. Masan Bay. The intensity of velocity near the bottom boundary layer(BBL) depends on the sea-bed irregularity in the coastal estuaries. So we calculated the depth mean velocity(DMV) considering that of BBL omitted in Kjerfve's calculation method. It revealed that BBL effect resulting in application of the bay acts largely on DMV in half more among 1l stations. The new expression of SD and LD per unit width in the cross-section using the developed DMV and proposed decomposition procedure of current were derived as follow : $$Q=u_0+\frac{1}{2}H_1{U_1cos(\varphi_h-\varphi_u)+U_3cos(\varphi_h-\varphi{ud})} LD ED SD$(Q_{skim}+Q_{sk2}) The third term, $Q_{sk2}$, on the right-hand of the equation is showed newly and arise from vertical oscillatory shear. According to the results applied in 3 cross-sections including 11 stations of the bay, the volume difference between proposed and previous SD was founded to be almost 2 times more at some stations. But their mean transport volumes over all stations are 18% less than the previous SD. Among two terms of SD, the flux of second term, $Q_{skim}$, is larger than third term, $Q_{sk2}$, in the main channel of cross-section, so that $Q_{skim}$ has a strong dependence on the tidal pumping, whereas third term is larger than second in the marginal channel. It means that $Q_{sk2}$ has trapping or shear effect more than tidal pumping phenomena. Maximum range of the fluctuation in LD is 40% as compared with the previous equations, but mean range of it is showed 11% at all stations, namely, small change. It mean that two components of SD interact as compensating flow. Therefore, the computation of SD and LD depend on decomposition procedure of velocity component in obtaining the volume transport of temporal and spacial flow through channels. The calculation of SD and LD proposed here can separate the shear effect from the previous SD component, so can be applied to non-uniform flow condition of cross-section, namely, baroclinic flow field.

  • PDF

탄성지반상에 놓인 철근 콘크리트 축대칭 쉘의 정적 및 동적 해석(IV) -축대칭 쉘의 동적 응답에 대한 철근의 영향을 중심으로- (Static and Dynamic Analysis of Reinforced Concrete Axisymmetric Shell on the Elastic Foundation -Effect of Steel on the Dynamic Response-)

  • 조진구
    • 한국농공학회지
    • /
    • 제39권4호
    • /
    • pp.106-113
    • /
    • 1997
  • Dynamic loading of structures often causes excursions of stresses well into the inelastic range, and the influence of the geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. A structure in a nuclear power plant is a structure of importance which puts emphasis on safety. A nuclear container is a pressure vessel subject to internal pressure and this structure is constructed by a reinforced concrete or a pre-stressed concrete. In this study, the material nonlinearity effect on the dynamic response is formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion is numerically solved by a central difference scheme. The constitutive relation of concrete is modeled according to a Drucker-Prager yield criterion in compression. The reinforcing bars are modeled by a smeared layer at the location of reinforcements, and the steel layer model under Von Mises yield criteria is adopted to represent an elastic-plastic behavior. To investigate the dynamic response of a nuclear reinforced concrete containment structure, the steel-ratios of 0, 3, 5 and 10 percent, are considered. The results obtained from the analysis of an example were summarized as follows 1. As the steel-ratio increases, the amplitude and the period of the vertical displacements in apex of dome decreased. The Dynamic Magnification Factor(DMF) was some larger than that of the structure without steel. However, the regular trend was not found in the values of DMF. 2. The dynamic response of the vertical displacement and the radial displacement in the dome-wall junction were shown that the period of displacement in initial step decreased with the steel-ratio increases. Especially, the effect of the steel on the dynamic response of radial displacement disapeared almost. The values of DMF were 1.94, 2.5, 2.62 and 2.66, and the values increased with the steel-ratio. 3. The characteristics of the dynamic response of radial displacement in the mid-wall were similar to that of dome-wall junction. The values of DMF were 1.91, 2.11, 2.13 and 2.18, and the values increased with the steel-ratio. 4. The amplitude and the period of the hoop-stresses in the dome, the dome-wall junction, and the mid-wall were shown the decreased trend with the steel-ratio. The values of DMF were some larger than those of the structure without steel. However, the regular trend was not found in the values of DMF.

  • PDF

루프-루프 전자탐사자료의 효과적인 2.5차원 역산 (An efficient 2.5D inversion of loop-loop electromagnetic data)

  • 송윤호;김정호
    • 지구물리와물리탐사
    • /
    • 제11권1호
    • /
    • pp.68-77
    • /
    • 2008
  • 2.5차원 전자탐사 적분방정식의 확장된 Born 근사해 또는 국소 비선형 근사에 기초하여 루프-루프 전자탐사 역산 알고리듬이 개발되었다 송수신 배열은 수평 동일면(HCP) 또는 수직 동일면(VCP) 방식이고, 다중 주파수 및 다중 송수신 간격을 포함할 수 있으며 PC에서 작동된다. 안정적이고 고해상도를 유지하는 역산이 가능하도록 변수분해 행렬과 Backus-Gilbert 분산 함수 분석을 통해 감도 분포의 함수로서의 공간적으로 변화하는 최적 Lagrange 곱수 결정 알고리듬을 포함하였다. HCP와 VCP 배열 자료가 지하 전기비저항 구조에 따라 서로 다른 감도를 가짐에 따라 동시 역산에서 안정성과 해상도에 영향을 미치게 되므로, 계산값과 측정값 차의 분산에 따라 가중치를 적용하는 방식을 도입하였다. 모델링 코드의 정확성은 통상적으로 루프-루프 전자탐사에서 사용하는 주파수 및 송수신 간격 범위에서 유한차분법에 의해 계산된 결과와의 비교를 통하여 증명되었다. 개발된 역산 알고리듬은 먼저 반무한 공간내 전도체 및 저항체가 포함된 모델에 대한 계산자료에 적용되어 성능이 입증되었다. 현장자료에 적용하고 그 결과 영상을 전기비저항 탐사자료에 대한 역산 결과와 비교하여, 의미있는 지하구조의 영상을 얻을 수 있음을 확인하였다.

다중 슬라이딩 표면 제어 기법에 기반한 쿼드로터의 능동 결함 허용 제어 (Active Fault Tolerant Control of Quadrotor Based on Multiple Sliding Surface Control Method)

  • 황남웅;김병수
    • 한국전자통신학회논문지
    • /
    • 제17권1호
    • /
    • pp.59-70
    • /
    • 2022
  • 본 논문에서는 쿼드로터의 모터 하나가 완전히 고장이 발생한 경우 쿼드로터의 위치 제어를 위한 능동 결함 허용 제어 방법을 제안한다. 소각의 가정 없이 라그랑지 방정식을 사용하여 쿼드로터의 동적 방정식을 구한다. 제안한 방법에서는 모터의 결함 검출을 위해 고장 검출 및 진단(FDD) 모듈과 고장 검출 및 분리(FDI) 모듈로 구성되는 고장 검출모듈을 설계한다. FDD 모듈에서는 구해진 동력학에 기반하여 쿼드로터의 상태를 관측하는 비선형 관측기를 설계한다. 관측된 쿼드로터의 상태들를 이용하여, 유수 신호를 설계하고 결함을 검출하기 위한 유수 신호의 적절한 문턱 값을 설정한다. 또한 설계된 추가 조건을 사용하여 결함 위치를 알아내기 위한 FDI 모듈을 설계한다. 모터의 결함을 검출한 후 쿼드로터가 원하는 경로로 비행하기 위해 다중 슬라이딩 표면 제어 기법에 기반한 결함 허용 제어기를 설계한다. 마지막으로, 모의실험을 통해 제안한 능동 결함 허용 제어 방법이 효용성을 검증한다.

복소 전기비저항 3차원 역산 알고리듬 개발 (Development of Three-dimensional Inversion Algorithm of Complex Resistivity Method)

  • 손정술;신승욱;박삼규
    • 지구물리와물리탐사
    • /
    • 제24권4호
    • /
    • pp.180-193
    • /
    • 2021
  • 복소 전기비저항 탐사기법은 진동수 영역에서 전기비저항과 위상을 측정하여 지하 매질의 다양한 특성정보를 획득할 수 있는 탐사기법으로 최근 그 활용성이 증가하고 있다. 이 논문에서는 복소 전기비저항 탐사기법의 활용성을 높이기 위하여 획득한 자료에 대한 3차원 역산 알고리듬을 개발하였다. 이를 위한 모델링에는 전자기 커플링 효과를 무시하는 경우에 적용할 수 있는 포아송 방정식을 적용하였으며, 역산에는 기존의 평활화된 역산법을 복소수로 확장하는 방법으로 알고리듬을 개발하였다. 역산의 안정성 및 현장자료의 적용성을 높이기 위하여 라그랑지 곱수를 역산 과정에서 오차 벡터와 모델 증분 벡터의 크기에 따라 자동적으로 조정되도록 하는 기법을 도입하였다. 또한, 잡음이 많이 포함된 위상자료로 인한 자료의 손실을 보완하기 위하여 역산반복 단계에서 초반부는 전기비저항 자료만을, 후반부는 전기비저항 자료와 위상 자료를 모두 역산하는 두 단계로 구성된 역산기법을 제시하였다. 수치 모형실험에 대한 역산 시험결과 안정적인 역산 결과를 얻을 수 있었으며, 개발된 3차원 역산 알고리듬을 국내 천열수 광산 인근에서 수행한 복소 전기비저항 탐사자료 해석에 적용하여 그 타당성을 확인하였다.