• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.391-401
• /
• 2000

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and the results are compared to the existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results show that as cross-sectional area decreases with the reduction in downstream pressure, flow rate increases and reaches the maximum when the speed index (mean velocity divided by wave speed) is near the unity at the point of minimum cross-section area, indicating the flow limitation or choking (flow speed equals wave speed) in one-dimensional studies. for further reductions in downstream pressure, flow rate decreases. The flow limitation or choking consist of the main reasons of waterfall effect which occurs in the airways, capillaries of lung, and other veins. Cross-sectional narrowing is significant but localized. When the ratio of downstream-to-upstream wall thickness is 2, the area throat is located near the downstream end. As this ratio is increased to 3, the constriction moves to the upstream end of the tube.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.815-826
• /
• 1994

Skew bridges are found frequently in new bridge construction due to geographical conditions when new constructing bridges are put across the existing highways, railroads or rivers. This study is to investigate the static behaviors of the steel deck plates of skew bridges which are increasingly used in bridges due to outstanding quality of structural steels, development of welding techniques, in order to reduce dead loads and period of constructions. The static behaviours of steel deck plates are analyzed using general purpose FE code SAP90 by modeling the skewed deck plates with rigorous finite elements, as the skew angles vary. The results of finite element analysis for the behaviors of steel deck plates and concrete slabs in acute, obtuse corners and center of decks are compared and discussed as the skew angles vary from $90^{\circ}$ to $30^{\circ}$. Two types of decks are treated, as isotropic plates and orthotropic plates, respectively. From the results of finite element analysis, it is found that more moments, reactions, and deflections occur at the obtuse corners than at the center of skewed decks regardless of isotropy or orthotropy. Especially, in case of the skewed deck plates with skew angles less than 45 degrees, significantly large discrepancies for the values of those internal forces are shown between the skewed and right deck plates. This study estimates the characteristics of deck behaviors according to skew angles, and proposes limitations of skew angles and the ciritical regions of decks.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.18 no.3
• /
• pp.10-19
• /
• 2014

Corrugated steel plates are made by fabricating thin steel plates to have trapezoidal or sinusoidal corrugation, and the corrugated plates are able to maintain high out-of-plane rigidity even when they are used instead of thick flat plates. Also, corrugated steel plates have almost no axial rigidity due to the accordion effect. Thus, if they are applied to the webs of plate girders, designing can be easily conducted so that the webs bear only shear stresses. However, unlike flat plates, the shear buckling of corrugated steel plates has very complex characteristics where buckling occurs due to the interaction of local and global buckling, besides local buckling and global buckling. For the investigation of the cause and characteristics of this interactive buckling, studies on sinusoidal corrugated steel plates are fewer than studies on trapezoidal corrugated steel plates. Therefore, in this study, the shear buckling characteristics of sinusoidal corrugated steel plates and the occurrence pattern of interactive buckling were investigated. For the calculation of shear buckling strength, a finite element program was used, and the analysis results were compared with the exact solution. In addition, the characteristics of buckling stress change and the change of buckling mode shape depending on corrugation thickness and shape parameter were analyzed, and by comparing these results with the results of a theoretical equation, the timing of buckling mode change was analyzed.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.34 no.2
• /
• pp.148-154
• /
• 2014

Photoelasticity is a useful technique for obtaining the differences and directions of principal stresses in a model. In conventional photoelasticity, the photoelastic parameters are measured manually point by point. Identifying and measuring photoelastic data is time-consuming and requires skill. The fringe phase shifting method was recently developed and has been found to be convenient for measuring and analyzing fringe data in photo-mechanics. This paper presents an experimental study on the stress distribution along a horizontal line that passes the central point of a rhombus plate made of Photoflex (i.e., type of urethane rubber). The isoclinic fringe and/or principal stress direction is constant on this horizontal line, so a four-bucket phase shifting method can be applied. The method requires four photoelastic fringes that are obtained from a circular polariscope by rotating the analyzer at $0^{\circ}C$, $45^{\circ}C$, $90^{\circ}C$ and $135^{\circ}C$. Experimental measurements using the method were quantitatively compared with the results from FEM analysis; the results from the two methods showed comparable agreement.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5A
• /
• pp.649-656
• /
• 2008

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchic void element based on the integrals of Legendre polynomials is used to characterize the fracture behaviour of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of zero stiffness element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem. It is noted that the proposed hierarchical void element can be one of alternatives to analyze the patched crack problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.287-292
• /
• 2007

In general, the response of bulk material is independent of its size when it comes to considering classical elasticity theory. Because the surface to bulk ratio of the large solids is very small, the influence of surface can be negligible. But the surface effect plays important role as the surface to bulk ratio becomes larger, that is, the contribution of the surface effect must be considered in nano-size elements such as thin film or beam structure. Molecular dynamics computation has been a conventional way to analyze these ultra-thin structures but this method is limited to simulate on the order of $10^6{\sim}10^9$ atoms for a few nanoseconds, and besides, very time consuming. Analysis of structures in submicro to micro range(thin-film, wire etc.) is difficult with classical molecular dynamics due to the restriction of computing resources and time. Therefore, in this paper, the continuum-based method is considered to simulate the overall physical and mechanical properties of the structures in nano-scale, especially, for the thin-film.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1995.10a
• /
• pp.289-295
• /
• 1995

구조물의 동적부하에 대한 동적변형 응답을 정확히 예측하고, Over Design이나 Under Design이 아닌 합리적인 설계방안의 개발은 중요한 과제이다. 동적강도해석이나 소음 승차감과 같은 진동 및 충격에 기인하는 제반 문제를 복잡한 구조물을 대상으로 합리적으로 처리하기 위한 Dynamic Design Analysis는 높은 신뢰성의 추구와 더불어 필요불가결한 기술이 되고 있다. 동적해석 방법으로는 현재 유한요소법이 널리 사용되고 있으며 여러 종류의 범용 프로그램들이 보급되어 있는 실정이다. 그러나 특히 동적문제에 있어서는 형상이나 거동이 복잡한 구조물의 경우, 또는 차량의 차체와 같이 많은 장착물이 부착된 경우에는 유한요소법의 적용이 곤란하여, 지금까지 대처할 수 있는 유용한 방법이 없었다. 따라서 비교적 용이하고 간단하게 적용가능한 진동실험을 기초로 한 구조물의 동적 응답해석 및 설계 방안의 개발이 필요하다. 본 연구에서는 진동시험으로 얻어진 부분구조물의 응답특성과 결합특성으로부터 결합 후의 응답특성을 예측할 수 있는 방법을 전달함수합성이론을 기초로하여 프로그래밍 package화 한다. 그리고 평판구조물에 대하여 진동시험과 컴퓨터 시뮬레이션을 통하여 개발된 방법의 타당성을 검증한다. 또한 실제 차량에서 차체만의 진동시험과 엔진의 자유진동시험에 의한 시험데이터로부터 차체와 엔진이 마운트 결합된 후의 진동특성을 예측한다. 진동시험시에 입력과 출력에 노이즈가 필연적으로 혼입되어 주파수응답함수의 크기(magnitude)와 위상(phase)을 왜곡시킨다. 특히 위상의 왜곡은 복소수연산을 하는 전달함수합성법의 결과에 중요한 영향을 미치게 된다. 본 연구에서는 데이타 획득시 입력과 출력의 시간지연으로 생기는 위상왜곡을 보정하는 방법을 제시하고, 그 개선 정도를 조사한다.는 소견의 확실도로서 가능성을 표현한 것이다. 예를 들면, 진동진폭 스펙트럼상에 2X 성분이 상당히 크게 나타나 정렬불량의 가능성이 0.7 정도라고 판정하는 것 등은 이러한 수치적진리치를 이용하는 방법이다. 그러나 상기의 수치적 표현만으로는 확실도를 한개의 수치로서 대표하게 하는 것은 진단의 정밀도에 문제가 있을 것으로 생각된다. 따라서 언어적진리치가 도입되어 [상당히 확실], [확실], [약간 확실] 등의 언어적인 표현을 이용하여 애매성을 표현하게 되었다. 본 논문에서는 간이진단 결과로부터 추출된 애매한 진단결과중에서 가장 가능성이 높은 이상원인을 복수로 선정하고, 여러 종류의 수치화할 수 없는 언어적(linguistic)인 정보ㄷㄹ을 if-then 형식의 퍼지추론으로 종합하는 회전기계의 이상진단을 위한 정밀진단 알고리즘을 제안하고 그 유용성을 검토한다. 존재하여도 모우드 변수들을 항상 정확하게 구할 수 있으며, 또한 알고리즘의 안정성이 보장된 것이다.. 여기서는 실험실 수준의 평 판모델을 제작하고 실제 현장에서 이루어질 수 있는 진동제어 구조물에 대 한 동적실험 및 FRS를 수행하는 과정과 동일하게 따름으로써 실제 발생할 수 있는 오차나 error를 실험실내의 차원에서 파악하여 진동원을 있는 구조 물에 대한 진동제어기술을 보유하고자 한다. 이용한 해마의 부피측정은 해마경화증 환자의 진단에 있어 육안적인 MR 진단이 어려운 제한된 경우에만 실제적 도움을 줄 수 있는 보조적인 방법으로 생각된다.ofile whereas relaxivity at high field is not affected by τS. On the other hand, the change in τV does not affect low field profile but strongly in fluences on bot

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.5
• /
• pp.459-465
• /
• 2015

There exist different kinds of aircrafts, such as conventional airplane, rotorcraft, fighter, and unmanned aerial vehicle. Their shape and feature are dependent upon their own assigned mission. One of the fundamental analyses performed during the aircraft design is the structural analysis. It becomes more complicated and requires severe computations because of the recent complex trends in aircraft structure. In order for efficiency in the structural analysis, a simplified approach, such as equivalent beam or plate model, is preferred. However, it is not clear which analysis will be appropriate to analyze the realistic configuration, such as an aircraft wing, i.e., between an equivalent beam and plate analysis. It is necessary to assess the limitation for both the one-dimensional beam analysis and the two-dimensional plate theory. Thus, in this paper, the static structural analysis results obtained by EDISON solvers were compared with the three-dimensional results obtained from MSC NASTRAN. Before that, EDISON program was verified by comparing the results with those from MSC NASTRAN program and other analytic solutions.

• Journal of the Korean GEO-environmental Society
• /
• v.17 no.4
• /
• pp.5-15
• /
• 2016

For achieving stability and economic construction at a retaining wall construction site, quantitative parameters of soil properties with excavation steps coincides with the actual field site. The main parameters of retaining wall design such as deformation modulus and modulus of horizontal subgrade reaction are common with N value of standard penetration test. Therefore, this study is compared and analyzed about the mutual relationship which is SPT, PBT and PMT for overcoming inconsistency of the existing retaining wall design generalized. In addition, modulus of horizontal subgrade reaction and reduction factor with excavation steps are proposed through back analysis of elasto-plasticity and finite element method with actual field monitoring data. Finally, it is purpose that parameter errors are reduced for applying effective retaining wall design at a construction small and medium-sized.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.11
• /
• pp.1437-1446
• /
• 2000

A numerical analysis on the freeze coating process of a non-isothermal finite dimensional plate with a binary alloy is performed to investigate the growth and decay behavior of the solid and the mushy layer of the freeze coat and a complete procedure to calculate the process is obtained in this study. The continuously varying solid and mushy layers are immobilized by a coordinate transform and the resulting governing differential equations are solved by a finite difference technique. To account for the latent heat release and property change during solidification, proper phase change models are adopted. And the convection in the liquid melt is modeled as an appropriate heat transfer boundary condition at the liquid/mushy interface. The present results are compared with analytic solutions derived for the freeze coating of infinite dimensional plates and the discrepancy is found to be less than 0.5 percent in relative magnitude for all simulation cases. In addition the conservation of thermal energy is checked. The results show that the freeze coat grows proportional to the 1.2 square of axial position as predicted by analytic solutions ar first. But after the short period of initial growth, the growth rate of the freeze coat gradually decreases and finally the freeze coat starts to decay. The effects of various non-dimensional processing parameters on the behavior of freeze coat are also investigated.