• Structural Engineering and Mechanics
• /
• 제68권4호
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.

• 수산해양기술연구
• /
• 제35권1호
• /
• pp.54-59
• /
• 1999

High-speed electronic digital computers have enabled engineers to employ various numerical discretization techniques for solutions of complex problems. The Finite Element Method is one of the such technique. The Finite Element Method is one of the numerical analysis based on the concepts of fundamental mathematical approximation. Three dimensional plate structures used often in partition of ship, box girder and frame are analyzed by Finite Element Method. In design of structures, the static deflections, stress concentrations and dynamic deflections must be considered. However, these problem belong to geometrically nonlinear mechanical structure analysis. The analysis of each element is independent, but coupling occurs in assembly process of elements. So, to overcome such a difficulty the shell theory which includes transformation matrix and a fictitious rotational stiffness is taken into account. Also, the Mindlin's theory which is considered the effect of shear deformation is used. The Mindlin's theory is based on assumption that the normal to the midsurface before deformation is "not necessarily normal to the midsurface after deformation", and is more powerful than Kirchoff's theory in thick plate analysis. To ensure that a small number of element can represent a relatively complex form of the type which is liable to occur in real, rather than in academic problem, eight-node quadratic isoparametric elements are used. are used.

• 한국전산구조공학회논문집
• /
• 제20권3호
• /
• pp.239-246
• /
• 2007

본 논문에서는 한두 부재의 순간적 결손에 따른 동적 거동을 정적해석을 통하여 합리적이고 효율적으로 해석할 수 있는 등가정적 연쇄붕괴 해석기법을 제시한다. 제시된 기법은 부재 결손에 따른 구조물 강성 변화 및 순간적 결손에 따른 동적거동 확대 효과를 등가의 하중으로 치환한 강성등가하중을 초기 구조물에 적용하여 해석하는 방법으로서 기둥을 하나씩 제거해 가며 반복해석을 수행해야 하는 연쇄붕괴해석 특성에 매우 효율적이면서도 신뢰성이 높은 장점을 갖는다. 제시한 강성등가하중에 의한 해석결과를 시간이력해석결과 및 GSA에 의한 해석결과와 비교한 결과, 휨모멘트, 축력, 및 수직변위 등의 측면에서 GSA에 의한 해석결과에 비해 시간이력해석결과에 상당히 근접하는 결과를 나타냈다. 이를 통해 강성등가 하중에 의한 해석기법이 GSA에 의한 정적해석방법을 대체하는 새로운 정적해석기법으로서 효용성이 있음을 확인하였다.

• Structural Engineering and Mechanics
• /
• 제27권6호
• /
• pp.713-739
• /
• 2007

Nowadays computers can perform symbolic computations in addition to mere number crunching operations for which they were originally designed. Symbolic computation opens up exciting possibilities in Structural Mechanics and engineering. Classical areas have been increasingly neglected due to the advent of computers as well as general purpose finite element software. But now, classical analysis has reemerged as an attractive computer option due to the capabilities of symbolic computation. The repetitive cycles of simultaneous - equation sets required by the finite element technique can be eliminated by solving a single set in symbolic form, thus generating a truly closed-form solution. This consequently saves in data preparation, storage and execution time. The power of Symbolic computation is demonstrated by six examples by applying symbolic computation 1) to solve coupled shear wall 2) to generate beam element matrices 3) to find the natural frequency of a shear frame using transfer matrix method 4) to find the stresses of a plate subjected to in-plane loading using Levy's approach 5) to draw the influence surface for deflection of an isotropic plate simply supported on all sides 6) to get dynamic equilibrium equations from Lagrange equation. This paper also presents yet another computationally efficient and accurate numerical method which is based on the concept of derivative of a function expressed as a weighted linear sum of the function values at all the mesh points. Again this method is applied to solve the problems of 1) coupled shear wall 2) lateral buckling of thin-walled beams due to moment gradient 3) buckling of a column and 4) static and buckling analysis of circular plates of uniform or non-uniform thickness. The numerical results obtained are compared with those available in existing literature in order to verify their accuracy.

• 대한토목학회논문집
• /
• 제13권1호
• /
• pp.55-63
• /
• 1993

탄성계수의 확률분포를 고려한 트러스 구조물의 유한요소 해석을 수행하였다. 탄생계수의 분포함수를 이용한 형상함수를 유도하여 강성행열을 구성하고 Monte Carlo Simulation 기법을 이용하여 수치해석을 수행하였다. 공간 트러스 구조물 해석을 통해서 해석과정을 검증하고 사장교에서 케이블 장력의 변동폭을 예측하였다.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2003년도 추계 학술발표회논문집
• /
• pp.191-198
• /
• 2003

Flat slab system has been adopted in many buildings constructed recently because of the advantage of reduced floor heights to meet the economical and architectural demands. Structural engineers commonly use the equivalent frame method(EFM) with equivalent beams proposed by Jacob S. Grossman in practical engineering for the analysis of flat slab structures. However, in many cases, when it is difficult to use the EFM, it is necessary to use a refined finite element model for an accurate analysis. But it would take significant amount of computational time and memory if the entire building structure were subdivided into a finer mesh. An efficient analytical method is proposed in this study to obtain accurate results in significantly reduced computational time. The proposed method employs super elements developed using the matrix condensation technique and fictitious beams are used in the development of super elements to enforce the compatibility at the interfaces of super elements. The stiffness degradation of flat slab system considered in the EFM was taken into account by reducing the elastic modulus of floor slabs in this study. Static and dynamic analyses of example structures were peformed and the efficiency and accuracy of the proposed method were verified by comparing the results with those of the refined finite element model and the EFM.

• 한국구조물진단유지관리공학회 논문집
• /
• 제20권3호
• /
• pp.10-17
• /
• 2016

본 연구는 단경간 교량의 정적하중입력/변위출력관계를 이용한 새로운 교량 유한요소모델 개선 방법을 제안하였고, 실내 모형교량 실험을 통해 검증하였다. 기존의 유한요소모델개선기법은 실험으로부터 얻어진 모드계수와 유한요소모델로부터 예측된 모드계수가 유사해지도록 유한요소모델을 개선하는데, 이 과정에서 구조계의 질량행렬에 대한 가정을 필요로 한다. 제안된 기법은 질량행렬을 가정하지 않고, 오히려 질량행렬 추정을 가능하게 하는 장점을 가진다. 제안된 기법은 두 단계로 구성된다. 첫째, 정적 하중입력-변위응답으로부터 강성행렬을 개선하고, 둘째, 실측된 고유진동수를 이용하여 질량행렬을 개선한다. 실험검증을 위하여 실내 모형교량을 제작하였고, 제안된 기법을 이용하여 모형교량의 탄성계수를 추정하였으며, Universal Testing Machine으로 부터 얻어진 탄성계수와 비교하였다. 또한 기존의 유한요소모델개선기법으로 추정된 탄성계수와 비교하였다. 실험의 결과들로부터 제안된 기법이 합리적으로 탄성계수와 질량밀도를 추정하는 것이 관찰되었고, 기존의 유한요소모델개선기법은 고차모드를 사용했을 때 상대적으로 큰 오차를 주는 것이 관찰되었다. 추가적으로 유한요소모델링 오차에 대하여 토의하였다.

• Structural Engineering and Mechanics
• /
• 제26권3호
• /
• pp.343-361
• /
• 2007

This article presents a numerical tool for dimensioning two-threaded fasteners connecting prismatic parts subjected to fatigue tension loads that are coplanar with the screw axis. A simplified numerical model is developed from unidirectional finite elements, modeling the connected parts and screws with bent elements and the elastic contact layer between the parts with springs. An algorithm updating the contact stiffness matrix, calculating forces and displacements at each node of the structure and thus normal stresses in the screws in both static and fatigue is further developed using C language. An experimental study is also conducted in parallel with the numerical approach to validate the developed model assumptions, the numerical model and the 3D finite element results. Since stiffness values for the compressive zones in the parts are analytically difficult to determine, a statistical software method is used, from which a tuning factor is derived for identifying these stiffness values. The method is also applied to set out the influence of each parameter on the fatigue behaviour of each screw. Finally, the developed model will be used to establish a new, sophisticated, fast and accurate tool for dimensioning bolted mechanical structures.

• 한국정밀공학회지
• /
• 제12권3호
• /
• pp.101-109
• /
• 1995

Axisymmetric and quasi-static finite element analysis of an inflated tire rotating with constant angular velocity and contact to road has been performed. Centrifugal force effect was added to load stiffness matrix and equation of effective material properties were calculated by the Halpin-Tsai formulation. In this report, radial truck/bus tire was analyzed. It was inflated and rotated at speeds up to 140 km/h. Then, contact problem was performed to calculate stress-strain field of tire wiht flat rigid road under the load due to the self-weight of a vehicle. Significant changes of stress-strain field of tire were observed in the finite element analysis. Shear stress, strain and strain energy density were rapidly increased at the dege of #2 belt at freely rotating state. This concentrated stress and strain made belt edge sparation. Under the condition of flat riged road contact, strain energy density of #2 belt, carcass turn-up part were concentrated and bigger values than only freely rotation state. Therefore, dynamic behaivor of tire has to considered as design factors which are affected to belt edge separation and bead breakage.

• Interaction and multiscale mechanics
• /
• 제2권3호
• /
• pp.321-332
• /
• 2009

Recent development in composites containing phase-transforming particles, such as vanadium dioxide or barium titanate, reveals the overall stiffness and viscoelastic damping of the composites may be unbounded (Lakes et al. 2001, Jaglinski et al. 2007). Negative stiffness is induced from phase transformation predicted by the Landau phase transformation theory. Although this unbounded phenomenon is theoretically supported with the composite homogenization theory, detailed stress analyses of the composites are still lacking. In this work, we analyze the stress distribution of the Hashin-Shtrikman (HS) composite and its two-dimensional variant, namely a circular inclusion in a square plate, under the assumption that the Young's modulus of the inclusion is negative. Assumption of negative stiffness is a priori in the present analysis. For stress analysis, a closed form solution for the HS model and finite element solutions for the 2D composite are presented. A static loading condition is adopted to estimate the effective modulus of the composites by the ratio of stress to average strain on the loading edges. It is found that the interfacial stresses between the circular inclusion and matrix increase dramatically when the negative stiffness is so tuned that overall stiffness is unbounded. Furthermore, it is found that stress distributions in the inclusion are not uniform, contrary to Eshelby's theorem, which states, for two-phase, infinite composites, the inclusion's stress distribution is uniform when the shape of the inclusion has higher symmetry than an ellipse. The stability of the composites is discussed from the viewpoint of deterioration of perfect interface conditions due to excessive interfacial stresses.