• Structural Engineering and Mechanics
• /
• 제16권3호
• /
• pp.259-274
• /
• 2003

In this paper, a novel numerical solution technique, the differential cubature method is employed to study the buckling problems of thick plates with arbitrary quadrilateral planforms and non-uniform boundary constraints based on the first order shear deformation theory. By using this method, the governing differential equations at each discrete point are transformed into sets of linear homogeneous algebraic equations. Boundary conditions are implemented through discrete grid points by constraining displacements, bending moments and rotations of the plate. Detailed formulation and implementation of this method are presented. The buckling parameters are calculated through solving a standard eigenvalue problem by subspace iterative method. Convergence and comparison studies are carried out to verify the reliability and accuracy of the numerical solutions. The applicability, efficiency, and simplicity of the present method are demonstrated through solving several sample plate buckling problems with various mixed boundary constraints. It is shown that the differential cubature method yields comparable numerical solutions with 2.77-times less degrees of freedom than the differential quadrature element method and 2-times less degrees of freedom than the energy method. Due to the lack of published solutions for buckling of thick rectangular plates with mixed edge conditions, the present solutions may serve as benchmark values for further studies in the future.

• 한국소음진동공학회논문집
• /
• 제24권5호
• /
• pp.414-420
• /
• 2014

In this paper, the vibration characteristics of the trapezoidal corrugated plate with axial stiffeners and lumped masses are investigated by the analytical method. The corrugated plate can be treated as an equivalent orthotropic plate as this plate has different flexure properties in two perpendicular directions; flexible in the corrugation direction and stiff in the transverse direction. The effective extensional and flexural stiffness of the equivalent plate are considered to obtain the precise solution in the analysis. The plate is stiffened by concentric stiffeners horizontally to the corrugation direction. The discrete stiffener theory is adopted to consider the position of stiffener. To demonstrate the validity of the proposed approach, the comparison is made with the results of 3D ANSYS finite element solutions. Some numerical results are presented to check the effect of the geometric properties.

• Structural Engineering and Mechanics
• /
• 제60권2호
• /
• pp.337-350
• /
• 2016

Even though extensive researches have been performed for steam explosion due to their complex mechanisms and inherent uncertainties, establishment of severe accident management guidelines and strategies is one of state-of-the arts in nuclear industry. The goal of this research is primarily to examine effects of vessel failure modes and locations on nuclear facilities under typical steam explosion conditions. Both discrete and integrated models were employed from the viewpoint of structural integrity assessment of steel components and evaluation of the cracking and crushing in reinforced concrete structures. Thereafter, comparison of systematic analysis results was performed; despite the vessel failure modes were dominant, resulting maximum stresses at the all steel components were sufficiently lower than the corresponding yield strengths. Two failure criteria for the reinforced concrete structures such as the limiting failure ratio of concrete and the limiting strains for rebar and liner plate were satisfied under steam explosion conditions. Moreover, stresses of steel components and reinforced concrete structures were reduced with maximum difference of 12% when the integrated model was adopted comparing to those of discrete models.

• 한국농공학회논문집
• /
• 제53권6호
• /
• pp.153-163
• /
• 2011

This research describe improved algorithm that is able to decide terminal criterion of Evolutionary Structural Optimization (ESO), reducing load of calculation to search load path of concrete beam, and apply to agricultural facilities. The ESO method is that make to discrete structure, structural analyze each element stress through FEM. And repeat generation with next material condition to become for most suitable composing. Individual element introduces concept of zero stiffness, but zero stiffness decisions are gone to direction of exclusion. In this stduy, improve algorithm to be convergence by 'Rule of Alive or Die' in arrival because is most suitable. Also, existing terminal criterion lack consistency because that used depend on experience of researcher. This research procedure is fellowed. First, all modulus of elasticity assume a half of elasticity modulus of material, Second, structural analysis by FEM, Third, apply to the remove ratio and restoration ratio for the 'rule of alive or die'. Forth, reconstruct the element and material conditions. And repeat the first to forth process. The terminal time of evolutional procedure is the all elastic modulus of element changed to blank value or elasticity modulus value of original. Therefore, in this study, consist the algorithm for programming, and apply to the agricultural facilities with concrete.

• 터널과지하공간
• /
• 제20권1호
• /
• pp.15-27
• /
• 2010

국토의 약 70%가 산악지형인 우리나라의 지형적 특성상 각종 건설공사에 따른 대절토 사면의 형성으로 인해 낙석의 위험성이 급증하고 있다. 특히 해마다 반복되는 집중호우 시의 사면 붕괴는 대규모 낙석을 동반하면서 막대한 인적, 물적 피해를 초래하고 있다. 이러한 문제점들을 해결하기 위해 많은 현장조사연구 및 기초연구들이 수행되고 있으나 이들 연구의 대부분은 낙석이 발생할 경우의 보강 및 낙석 위험도 분석 등에 국한되고 있어, 근본적인 낙석의 거동에 관한 검토 및 분석은 미진한 실정이다. 따라서 본 연구에서는 수치해석적 기법을 이용하여 사면의 외적 요소에 따른 낙석의 거동 특성을 규명하였으며, 이를 위해 개별요소법의 하나인 PFC2D가 이용되었다. 사면의 외적 요소로 사면의 표면상태, 구배, 소단 개수, 소단 구배 및 낙석의 크기에 대한 소단의 너비의 비 등이 검토되었으며, 이들 요소에 대한 다양한 변화를 통해 낙석의 거동 메커니즘이 수치해석적으로 분석되었다.

• 한국터널지하공간학회 논문집
• /
• 제18권2호
• /
• pp.129-142
• /
• 2016

지표침하를 발생시키는 다양한 원인 중에 하나로 알려진 지하에서 발생되는 공동은 지반의 불연속적 특징이다. 그러므로 기존의 연속체해석으로 분석하기에는 한계가 있을 것으로 판단된다. 본 연구에서는 지표침하를 발생시키는 구형지하공동에 대해서 개별요소법을 활용하여 분석하였다. 지표침하에 영향을 미치는 구형지하공동의 특성인자로서 지반의 물성치, 구형지하공동의 심도 및 크기를 선택하였고, 각 특성인자들의 값의 변화가 지표침하에 미치는 영향에 대해 분석하였다. 또한 지하공동의 직경에 대한 상대적 깊이와 지표침하의 관계를 분석하였고, 이를 통해 지하공동의 붕괴 예측 및 보강유무를 결정하기 위한 근간을 제시하였다.

• 한국재료학회지
• /
• 제26권11호
• /
• pp.611-622
• /
• 2016

Particle morphology change and different experimental condition analysis during composite fabrication process by traditional ball milling with discrete element method (DEM) simulation were investigated. A simulation of the three dimensional motion of balls in a traditional ball mill for research on the grinding mechanism was carried out by DEM simulation. We studied the motion of the balls, the ball behavior energy and velocity; the forces acting on the balls were calculated using traditional ball milling as simulated by DEM. The effect of the operational variables such as the rotational speed, ball material and size on the flow velocity, collision force and total impact energy were analyzed. The results showed that increased rotation speed with interaction impact energy between balls and balls, balls and pots and walls and balls. The rotation speed increases with an increase of the impact energy. Experiments were conducted to quantify the grinding performance under the same conditions. Furthermore, the results showed that ball motion affects the particle morphology, which changed from irregular type to plate type with increasing rotation speed. The evolution was also found to depend on the impact energy increase of the grinding media. These findings are useful to understand and optimize the particle motion and grinding behavior of traditional ball mills.

• 한국지진공학회논문집
• /
• 제9권2호통권42호
• /
• pp.47-58
• /
• 2005

이산폰툰형 부유식교량의 시간영역 지진응답해석을 지진파의 공간분포를 고려하여 수행하였다. 지반운동의 공간변화는 파의 통과 영향, 엇결성 뎡향 및 국부지반 영향을 반영하는 결맞음 함수 모델로 고려하였다. 부유식 교량의 상부구조는 공간뼈대요소와 탄성현수케이블요소를 이용하여 모델링 하였고, FHWA 가이드라인의 스프링 모델을 사용하여 교대의 지반 구조물 상호작용에 대한 영향을 고려하엿으며, 경계요소법으로 산정한 폰툰 동수력계수의 주파수 의종성을 고려하기 위해서 시간지연함수를 사용하였다. 지반운동의 공간변화를 고려한 다중지점 지진입력을 교량의 양단에 도입하였고 응답의 시간이력을 동시 가진 시와 비교하였다. 장주기의 동적특성을 갖는 이산폰툰형 부유식 교량의 지진응답 해석에서 지반운동의 공간분포를 고려할 경우에 주몰할만한 응답의 증폭이 나타날 수 있음을 확인하였다.

• 한국지반신소재학회논문집
• /
• 제17권4호
• /
• pp.129-139
• /
• 2018

개별요소법은 입자형태의 재료를 다루는 분야에서 널리 사용되는 방법이다. 특히, 지반에서와 같이 변형이 크게 발생하는 특성을 지닌 재료의 경우에도 유용하게 사용할 수 있는 해석방법 중 하나이다. 본 연구에서는 개별요소법을 사용하여 지반을 생성하고 지반에 TBM 형상을 관입시킴으로써 굴진을 모사하고자 하였다. 해석에 사용된 TBM은 7.73m 직경의 토압식 쉴드 TBM을 사용하였으며 해석을 통해 굴진 성능을 예측할 수 있는 수치해석 모델을 검토하고자 하였다. 대상 모델의 스포크 개수는 8개이며, 커터헤드의 개구율은 약 21.31%이다. 또한, 커터헤드에는 52개 디스크 커터와 167개 커터비트가 굴착 도구로 장착되어 있으며, 커터헤드는 게이지 커터가 장착된 커터헤드 외주면이 굴곡진 세미돔 타입이다. 해석을 통해 커터헤드와 각각의 절삭 도구에 작용하는 반력과 저항 토크를 검토할 수 있었다. 또한 커터헤드 중심에서부터 거리에 따른 비교를 통해 커터헤드의 절삭 도구 위치별로 발생하는 반력과 토크를 검토하였다.

• Structural Engineering and Mechanics
• /
• 제84권4호
• /
• pp.547-560
• /
• 2022

By presetting various reinforcement diameters in topology optimization with the discrete model finite element analysis, an algorithm of bidirectional evolutionary structural optimization of multi-level reinforcement diameter is presented to obtain the optimal reinforcement topologies which describe the degree of stress of different parts. The results of a comparative study on different reinforcement feasible domain demonstrate that the more angle types of reinforcement are arranged in the initial domain, the higher utilization rate of reinforcement of the optimal topology becomes. According to the nonlinear finite element analysis of some deep beam examples, the ones designed with the optimization results have a certain advantage in ultimate bearing capacity, although their failure modes are greatly affected by the reinforcement feasible domain. Furthermore, the bearing capacity can be improved when constructional reinforcements are added in the subsequent design. However the adding would change the relative magnitude of the bearing capacity between the normal and inclined section, or the relative magnitude between the flexural and shear capacity within the inclined section, which affects the failure modes of components. Meanwhile, the adding would reduce the deformation capacity of the components as well. It is suggested that the inclined reinforcement and the constructional reinforcement should be added properly to ensure a desired ductile failure mode for components.