• 제목/요약/키워드: Finite Element-Boundary Element Method

검색결과 1,342건 처리시간 0.029초

SU/PG 기법을 이용한 이송이 지배적인 흐름 수치모의 (Numerical Simulation of Convection-dominated Flow Using SU/PG Scheme)

  • 송창근;서일원
    • 대한토목학회논문집
    • /
    • 제32권3B호
    • /
    • pp.175-183
    • /
    • 2012
  • 본 연구에서는 천수방정식에 Galerkin법과 Petrov-Galerkin 기법의 일종인 SU/PG 기법을 적용하여 유한요소 천수흐름 해석 모형을 개발하고, 다양한 실험수로에서 이송이 지배적인 흐름을 수치 모의하였다. 수로 내부에 얇은 판 형태의 구조물이 존재하는 경우 Fr 수와 Re 수가 매우 낮은 흐름에서는 Galerkin 기법과 SU/PG 기법이 동일한 결과를 나타냈으나, Fr=1.58인 경우 Galerkin법은 발산하여 해를 얻을 수 없었다. 이 경우 SU/PG법은 Newton-Raphson법에 의한 5회의 반복에 의해 수렴된 유속결과를 구할 수 있었다. 사류와 상류가 혼재하여 천이류가 나타나는 단면확대 수로 모의에서 SU/PG 기법을 이용한 본 연구의 경우 상류단 수심조건이 잘 유지되며, 도수가 발생하는 지점 및 도수에 의한 수심 경사, 도수 후의 수심이 모두 Khalifa(1980)의 실험결과와 매우 근사하였다. 이송이 지배적인 사류(Fr=2.74)에 의한 사각도수 모의의 경우에도 Galerkin 기법은 최초 모의시간의 첫 번째 반복 이후 발산하였으나, SU/PG 기법은 도수 경계면을 수치진동 없이 잘 포착하였으며, 해석해와 비교한 수심 및 유속의 최대 오차는 0.2% 이내로 나타나 기존 연구(Levin 등, 2006; Ricchiuto 등, 2007)에 비해 더욱 정확한 결과를 도출하였다.

홍수터 식생을 고려한 불규칙한 단면에서의 수위-유량 곡선 및 단위유량 횡분포 예측 (Prediction of Stage Discharge Curve and Lateral Distribution of Unit Discharge in an Arbitrary Cross Section Channel with Floodplain Vegetation)

  • 김태범;장지연;신재국;최성욱
    • 한국수자원학회논문집
    • /
    • 제44권2호
    • /
    • pp.157-167
    • /
    • 2011
  • 본 연구에서는 복단면개수로및 불규칙한 하상을보이는 횡단면 상에서의 수위-유량곡선 및단위유량횡분포예측을 위한 유한요소모형을 개발하였다. 지배방정식은 정상류와 종방향 등류를 가정한 운동량방정식을 이용하며, 수면은 횡단면에 걸쳐 일정하다고 가정한다. 홍수터 식생의 영향을 반영하기 위해 식생항력이 지배방정식에 포함되었으며, 수치해를 구하기 위해 유한요소법을 적용하였다. 단면형상과 Manning의 조도계수, 식생정보, 종방향 하상경사를 입력자료로 수위-유량 곡선을 예측 가능하며, 주어진 수위에서의 흐름방향 단위유량의 횡방향 분포를 예측할 수 있다. 개발된 모형의 검증을 위해 실측자료 및 Darby and Thorne (1996)의 모형 결과, 그리고 비선형 k-$\epsilon$ 모형의 결과와도 비교하였다. 검증된 모형의 알고리즘을 2차원 모형의 상류단 경계조건 설정에 활용하여, 유입유량을 절점별 단위유량으로 분배시켰을 때와 그렇지 않았을 때의 결과를 비교하였다.

낮은 강도를 갖는 산 안드레아 단층의 열탄성 특성 (Thermoelastic Aspects of the San Andreas Faults under Very Low Strength)

  • 박무춘;한욱
    • 한국지구과학회지
    • /
    • 제21권3호
    • /
    • pp.315-322
    • /
    • 2000
  • 본 연구는 장기간에 걸쳐 산 안드레아 단층계 내에서 56개 지점의 단층이동률에 대한 지질학적인 측정자료를 기준으로 모델을 설정하였다. 모델은 산안드레아 단층을 중심으로 한 수렴대에서 낮은 마찰(${\mu}$=0.3)을 갖는 단층군에 대해 최적의 결과를 보여주고 있다. 저강도를 갖는 단층에 대해 국지적인 이상값이나 대표값을 결정하는 것은 분명히 중요한 의미를 갖는다. 더욱이 이러한 연구는 지구조적인 체계에서의 단층의 강도를 결정하는 데 도움이 될 것으로 보인다. 예상치 못한 원인에 의한 공극압력이나 마찰법칙의 적절성에 대한 의문을 고려하지 않을 수 없을 것이다. 이러한 문제를 해결하기 위해서, 다른 가설하에서 단층의 유체학적인 모의 실험이 가능한 유한요소법을 적용하기 위해 세 가지의 단층분석 모델을 시도하였다. 계산된 모델은 추정된 유체역학적 특성과 판구조경계 조건을 만족하며, 현재의 지진파 표면속도, 변형률과 강도의 예측값을 나타내고 있다. 모델 연구의 결과는 평균 단층이동률, 강도의 방향과 측지학적인 자료의 예측값 범위 내에서 실제 측정치에 접근하고 있음을 보여준다. 본 연구는 저강도를 갖는 산 안드레아 단층계에서의 상호관련성을 해석하기 위한 열탄성 특성의 적용 결과를 잘 제시하고 있다.

  • PDF

도로 폭이 좁은 굴착공사에서 연결부가 적용되는 장지간 주형의 수치해석적 거동 평가 (Numerical Analysis on the Structure Behavior of the Connected Long-span Beam during Excavation in Narrow Streets)

  • 최광수;하상봉;이환우
    • 한국전산구조공학회논문집
    • /
    • 제33권4호
    • /
    • pp.263-270
    • /
    • 2020
  • 본 연구의 대상은 도로 폭이 좁은 시가지에서 굴착공사 시 적용되는 장지간 주형의 연결부이다. 일반적으로 적용되고 있는 연결부에서 상부 플랜지의 단차와 피로균열 등의 문제로 연결부의 신뢰도가 저하된다. 연결부의 결함을 보완하고 안전성을 향상시킨 개선형 연결부를 개발하였다. 유한요소 기반의 상용프로그램(ABAQUS)를 이용하여 개선형 연결부의 거동을 평가하였다. 먼저, 개선형 연결부에 적용되는 고장력 볼트 연결 및 강재와 콘크리트의 합성거동을 구현하기 위한 수치해석 방법을 제안하였다. 비교논문의 실험결과와 수치해석 결과의 비교를 통하여 개선형 연결부를 해석하는데 있어 수치해석 방법의 적합성을 검증하였다. 본 연구에서 제안하는 수치해석 방법을 적용하여 개선형 연결부와 일반형 연결부가 적용된 장지간 주형을 해석하였다. 장지간 주형의 탄소성 거동과 연결부의 응력분포를 수치 해석적으로 비교분석하였다. 개선형 연결부의 도입으로 25%의 압축응력이 감소되며 구조적 성능 개선효과 및 안전성을 확인하였다.

Avoidance of Internal Resonances in Hemispherical Resonator Assemblies from Fused Quartz Connected by Indium Solder

  • 세르게이 사라플로프;이희남;박상진
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2013년도 춘계학술대회 논문집
    • /
    • pp.835-841
    • /
    • 2013
  • Modern solid-state gyroscopes (HRG) with hemispherical resonators from high-purity quartz glass and special surface superfinishing and ultrathin gold coating become the best instruments for precise-grade inertial reference units (IRU) targeting long-term space missions. Designing of these sensors could be a notable contribution into development of Korea as a space nation. In participial, 40mm diameter thin-shell resonator from high-purity fused quartz, fabricated as a single-piece with its supporting stem has been designed, machined, etched, tuned, tested, and delivered by STM Co. (ATS of Ukraine) several years ago; an extremely-high Q-factor (upto 10~20 millions) has been shown. Understanding of the best way how to match such a unique sensor with inner glass assembly of the gyro means how to use the high potential in a maximal extent; and this has become the urgent task. Inner quartz glass assembly has a very thin indium (In) layer soldered the resonator and its silica base (case), but effects of internal resonances between operational modal pair of the shell-cup and its side (parasitic) modes can notable degrade the potential of the sensor as a whole, instead of so low level of resonator's intrinsic losses. Unfortunately, there are special combinations of dimensions of the parts (so-called, "resonant sizes"), when intensive losses of energy occurs. The authors proposed to use the length of stem's fixture as an additional design parameter to avoid such cases. So-called, a cyclic scheme of finite element method (FEM) and ANSYS software were employed to estimate different combinations of gyro assembly parameters. This variant has no mismatches of numerical origin due to FEM's discrete mesh. The optimum length and dangerous "resonant lengths" have been found. The special attention has been paid to analyses of 3D effects in a cup-stem transient zone, including determination of a difference between the positions of geometrical Pole of the resonant hemisphere and of its "dynamical Pole", i.e., its real zone of oscillation node. Boundary effects between the shell (cup) and 3D short "beams" (inner and outer stems) have been ranged. The results of the numerical experiments have been compared with the classic model of a quasi-hemispherical shell band with inextensional midsurface, and the solution using Rayleigh's functions of the $1^{st}$ and $2^{nd}$ kinds. To guarantee the truth of the recommended sizes to a designer of the real device, the analytical and FEM results have been compared with experimental data for a party of real resonators. The consistency of the results obtained by different means has been shown with errors less than 5%. The results notably differ from the data published earlier by different researchers.

  • PDF

초기 처짐을 갖는 Spherical Shell의 동적 특성에 관한 연구(II) - 초기 처짐에 따른 동적 특성 - (Dynamic Characteristics Analysis of Spherical Shell with Initial Deflection(II) - Effects of Initial Deflection -)

  • 조진구
    • 한국농공학회지
    • /
    • 제40권5호
    • /
    • pp.91-99
    • /
    • 1998
  • The widespread use of thin shell structures has created a need for a systematic method of analysis which can adequately account for arbitrary geometric form and boundary conditions as well as arbitrary general type of loading. Therefore, the stress and analysis of thin shell has been one of the more challenging areas of structural mechanics. A wide variety of numerical methods have been applied to the governing differential equations for spherical and cylindrical structures with a few results applicable to practice. The analysis of axisymmetric spherical shell is almost an every day occurrence in many industrial applications. A reliable and accurate finite element analysis procedure for such structures was needed. Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometry changes on the response is also significant in many cases. Therefore both material and geometric nonlinear effects should be considered. In general, the shell structures designed according to quasi-static analysis may fail under conditions of dynamic loading. For a more realistic prediction on the load carrying capacity of these shell, in addition to the dynamic effect, consideration should also include other factors such as nonlinearities in both material and geometry since these factors, in different manner, may also affect the magnitude of this capacity. The objective of this paper is to demonstrate the dynamic characteristics of spherical shell. For these purposes, the spherical shell subjected to uniformly distributed step load was analyzed for its large displacements elasto-viscoplastic static and dynamic response. Geometrically nonlinear behaviour is taken into account using a Total Lagrangian formulation and the material behaviour is assumed to elasto-viscoplastic model highly corresponding to the real behaviour of the material. The results for the dynamic characteristics of spherical shell in the cases under various conditions of base-radius/central height(a/H) and thickness/shell radius(t/R) were summarized as follows : The dynamic characteristics with a/H. 1) AS the a/H increases, the amplitude of displacement in creased. 2) The values of displacement dynamic magnification factor (DMF) were ranges from 2.9 to 6.3 in the crown of shell and the values of factor in the mid-point of shell were ranged from 1.8 to 2.6. 3) As the a/H increases, the values of DMF in the crown of shell is decreased rapidly but the values of DMF in mid-point shell is increased gradually. 4) The values of DMF of hoop-stresses were range from 3.6 to 6.8 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.3 to 2.6, and the values of DMF of stress were larger than that of displacement. The dynamic characteristics with t/R. 5) With the thickness of shell decreases, the amplitude of the displacement and the period increased. 6) The values of DMF of the displacement were ranged from 2.8 to 3.6 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.1 to 2.2.

  • PDF

경계형 모델을 사용한 초고강도 섬유보강 콘크리트거더의 파괴역학적 해석 (Fracture Simulation of UHPFRC Girder with the Interface Type Model)

  • 궈이홍;한상묵
    • 한국전산구조공학회논문집
    • /
    • 제23권1호
    • /
    • pp.81-94
    • /
    • 2010
  • 본 연구는 선형 상보법으로 초고강도 섬유보강 콘크리트 I형보의 파괴역학적 해석을 수치해석으로 수행하였다. 기존의 보통강도 콘크리트에 대한 유사 취성 파괴역학적 수치해석을 기반으로 초고강도 섬유보강 콘크리트 재료역학적 구성모델파괴 면에 인장경화 관계를 도입함으로써 초고강도 섬유보강 콘크리트 I형 거더 해석을 개선시켰다. 상수변형률 삼각형 요소에 꼭지점 또는 요소의 중앙점 절점을 배제하고 요소의 변에 절점을 배치한 결합된 삼각형 요소를 사용하였다. 인장영역에서는 경화/연화 파괴역학적 구성모델을, 전단영역에서는 연화 파괴역학적 구성모델을, 경계절점의 압축에 대해서는 연화파괴역학적 구성모델을 사용하여 파괴역학적 해석을 수행하였다. Non-holonomic rate 형태로 경로에 의존적인 경화연화거동을 LCP로 방정식을 구성하였으며, 그 해는 PATH를 사용해서 구하였다. Piece-wise 비탄성 항복-파괴면은 두 개의 압축 caps, 두 개의 Mohr-Coulomb 파괴면, 인장항복면과 인장파괴면 등으로 구성하였다. 초고강도 섬유보강 콘크리트 거더의 변형거동과 파괴 상태와 비교하여 이 수치해석 방법에 대한 유효성을 검증하였다.

Step-down Piezoelectric Transformer Using PZT PMNS Ceramics

  • Lim Kee-Joe;Park Seong-Hee;Kwon Oh-Deok;Kang Seong-Hwa
    • KIEE International Transactions on Electrophysics and Applications
    • /
    • 제5C권3호
    • /
    • pp.102-110
    • /
    • 2005
  • Piezoelectric transformers(PT) are expected to be small, thin and highly efficient, and which are attractive as a transformer with high power density for step down voltage. For these reasons, we have attempted to develop a step-down PT for the miniaturized adaptor. We propose a PT, operating in thickness extensional vibration mode for step-down voltage. This PT consists of a multi-layered construction in the thickness direction. In order to develop the step-down PT of 10 W class and turn ratio of 0.1 with high efficiency and miniaturization, the piezoelectric ceramics and PT designs are estimated with a variety of characteristics. The basic composition of piezoelectric ceramics consists of ternary yPb(Zr$_{x}$Ti$_{1-x}$)O$_{3}$-(1-y)Pb(Mn$_{1/3}$Nb1$_{1/3}$Sb$_{1/3}$)O$_{3}$. In the piezoelectric characteristics evaluations, at y=0.95 and x=0.505, the electromechanical coupling factor(K$_{p}$) is 58$\%$, piezoelectric strain constant(d$_{33}$) is 270 pC/N, mechanical quality factor(Qr$_{m}$) is 1520, permittivity($\varepsilon$/ 0) is 1500, and Curie temperature is 350 $^{\circ}C$. At y = 0.90 and x = 0.500, kp is 56$\%$, d33 is 250 pC/N, Q$_{m}$ is 1820, $\varepsilon$$_{33}$$^{T}$/$\varepsilon$$_{0}$ is 1120, and Curie temperature is 290 $^{\circ}C$. It shows the excellent properties at morphotropic phase boundary regions. PZT-PMNS ceramic may be available for high power piezoelectric devices such as PTs. The design of step-down PTs for adaptor proposes a multi-layer structure to overcome some structural defects of conventional PTs. In order to design PTs and analyze their performances, the finite element analysis and equivalent circuit analysis method are applied. The maximum peak of gain G as a first mode for thickness extensional vibration occurs near 0.85 MHz at load resistance of 10 .The peak of second mode at 1.7 MHz is 0.12 and the efficiency is 92$\%$.

생물개스 발생시스템을 위한 지하매설콘크리트 다이제스터의 열전달에 관한 연구 (Study on the Heat Transfer Phenomenon around Underground Concrete Digesters for Bigas Production Systems)

  • 김윤기;고재균
    • 한국농공학회지
    • /
    • 제22권1호
    • /
    • pp.53-66
    • /
    • 1980
  • The research work is concerned with the analytical and experimental studies on the heat transfer phenomenon around the underground concrete digester used for biogas production Systems. A mathematical and computational method was developed to estimate heat losses from underground cylindrical concrete digester used for biogas production systems. To test its feasibility and to evaluate thermal parameters of materials related, the method was applied to six physical model digesters. The cylindrical concrete digester was taken as a physical model, to which the model,atical model of heat balance can be applied. The mathematical model was transformed by means of finite element method and used to analyze temperature distribution with respect to several boundary conditions and design parameters. The design parameters of experimental digesters were selected as; three different sizes 40cm by 80cm, 80cm by 160cm and l00cm by 200cm in diameter and height; two different levels of insulation materials-plain concrete and vermiculite mixing in concrete; and two different types of installation-underground and half-exposed. In order to carry out a particular aim of this study, the liquid within the digester was substituted by water, and its temperature was controlled in five levels-35。 C, 30。 C, 25。 C, 20。C and 15。C; and the ambient air temperature and ground temperature were checked out of the system under natural winter climate conditions. The following results were drawn from the study. 1.The analytical method, by which the estimated values of temperature distribution around a cylindrical digester were obtained, was able to be generally accepted from the comparison of the estimated values with the measured. However, the difference between the estimated and measured temperature had a trend to be considerably increased when the ambient temperature was relatively low. This was mainly related variations of input parameters including the thermal conductivity of soil, applied to the numerical analysis. Consequently, the improvement of these input data for the simulated operation of the numerical analysis is expected as an approach to obtain better refined estimation. 2.The difference between estimated and measured heat losses was shown to have the similar trend to that of temperature distribution discussed above. 3.It was found that a map of isothermal lines drawn from the estimated temperature distribution was very useful for a general observation of the direction and rate of heat transfer within the boundary. From this analysis, it was interpreted that most of heat losses is passed through the triangular section bounded within 45 degrees toward the wall at the bottom edge of the digesten Therefore, any effective insulation should be considered within this region. 4.It was verified by experiment that heat loss per unit volume of liquid was reduced as the size of the digester became larger For instance, at the liquid temperature of 35˚ C, the heat loss per unit volume from the 0. 1m$^3$ digester was 1, 050 Kcal/hr m$^3$, while at for 1. 57m$^3$ digester was 150 Kcal/hr m$^3$. 5.In the light of insulation, the vermiculite concrete was consistently shown to be superior to the plain concrete. At the liquid temperature ranging from 15。 C to 350 C, the reduction of heat loss was ranged from 5% to 25% for the half-exposed digester, while from 10% to 28% for the fully underground digester. 6.In the comparison of heat loss between the half-exposed and underground digesters, the heat loss from the former was fr6m 1,6 to 2, 6 times as much as that from the latter. This leads to the evidence that the underground digester takes advantage of heat conservation during winter.

  • PDF

해상풍력터빈의 고유진동수 예측을 위한 지반에 인입된 파일의 탄성지지보 모델 기반 수평 거동 해석 (An analysis of horizontal deformation of a pile in soil using a beam-on-spring model for the prediction of the eigenfrequency of the offshore wind turbine)

  • 유정수;백경민;김태룡
    • 한국음향학회지
    • /
    • 제35권4호
    • /
    • pp.261-271
    • /
    • 2016
  • 지반에 인입된 파일에 대한 거동 해석은 지반의 비선형 거동 특성으로 인해 일반적으로 유한요소에 기반한 수치적인 방법을 주로 이용한다. 그러나 수치 해석은 파일-지반 모델링 및 연산에 많은 노력과 시간을 요구하므로 파일의 제원과 지반 물성치가 확정되지 않은 초기 설계 단계에서는 활용에 많은 한계를 갖는다. 반면, 지반을 선형화한 이론해석의 경우 수치 해석에 비해 모델이 단순하고 연산 시간이 매우 짧으므로, 해석의 신뢰성이 확보된다면 지반-지지구조의 거동 특성을 초기에 예측하는데 유용할 것이다. 본 연구에서 대상으로 하는 풍력발전기의 경우, 초기 설계 단계에서 풍력터빈의 고유진동수 예측을 위한 동적 거동 해석이 요구되며, 이 때 지반에 인입된 풍력터빈의 지지구조는 탄성경계조건으로 단순화하여 동적 거동 해석에 반영할 수 있다. 이를 위해, 본 논문에서는 풍력터빈 지지구조의 선단부 탄성 계수를 도출하고자 지반에 인입된 파일에 대한 파일-지반 연성해석을 수행하였다. 해석 시 지반의 변형은 탄성범위 이내에 있다고 단순화하여 지반에 인입된 파일을 탄성지지된 보로 모델링하였다. 탄성지지보 모델을 이용해 파일 선단에 수평 하중 또는 모멘트가 작용할 때 발생하는 파일의 횡변형을 구하고, 이로부터 영향계수를 도출하였다. 풍력터빈의 지지구조에 대한 해석 예로써, 모노파일과 석션파일에 대해 파일 선단의 영향계수를 구하고, 이를 문헌의 결과와 비교함으로써 해석 결과의 신뢰도를 검증하였다. 또한 이 두 파일의 깊이에 따른 변형 및 선단부의 스프링 상수를 비교하여 지지구조의 강성 측면에서 모노파일과 석션파일의 특성을 살펴보았다.