• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2009년도 추계학술대회 논문집
• /
• pp.321-321
• /
• 2009

This paper presents quasi-zero-stiffness (QZS) characteristic of a passive isolator using flexures under compression force. The passive isolator consists of a positive stiffness element (a vertical coil spring) and a negative stiffness element (flexures under compression force), and their proper combination of the positive and negative stiffness elements can produce both substantial static and zero dynamic stiffness, so called QZS. Firstly, a nonlinear dimensionless expression of a flexure under compression force is derived. A dynamic model of the passive isolator is developed and numerical simulations of its time and frequency response are performed. Then, undesirable nonlinear vibration is quantified using a period doubling bifurcation diagram and a Poincare's map of the isolator under forced excitation. Finally, experiments are performed to validate the QZS characteristic of the passive isolator.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior 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. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. 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 of materials. 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 hierarchical void 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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.589-596
• /
• 2006

An improved numerical method to obtain the exact element stiffness matrix is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric thin-walled beam-columns with two-types of elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column. This numerical technique is firstly accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Then exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions.

• 한국강구조학회 논문집
• /
• 제14권4호
• /
• pp.509-518
• /
• 2002

비대칭 단면을 갖는 박벽보의 3차원 휨-배틂 좌굴해석 및 정적해석을 위하여, 평형방정식과 힘-변위 관계식을 이용하여 엄밀한 정적요소강성행렬을 수치적으로 산정하는 개선된 기법을 제시한다. 먼저 14개의 변위피라미터를 도입하여 고차의 연립미분방정식을 1차 연립미분방정식으로 변환하고, 복소수 영역에서 선형고유치문제를 해를 구한다. 이 경우 동적강성행렬을 산정하는 경우와는 달리 복수개의 '영'의 고유치가 발생한다. 이에 대응하는 변위피라미터의 다항식을 항등식 조거능로부터 구하고, 이를 고유치와 결합하여 박벽보 요소의 엄밀한 처짐함수를 구한다. 이렇게 구한 엄밀한 처짐함수에 재단력-변위 관계식을 적용하여 세가지 초기단면력 조건에 대응하는 엄밀한 정적요소강성행렬을 산정한다. 본 방법의 타당성을 보이기 위하여 비대칭 박벽보의 좌굴하중과 처짐값을 계산하고 해석해나 ABAQUS 쉘요소를 이용한 해석결과 및 직선보요소를 사용한 유한요소해의 결과와 비교, 검증한다.

• Structural Engineering and Mechanics
• /
• 제8권6호
• /
• pp.591-605
• /
• 1999

An accurate and efficient shell element is presented. The stiffness of the shell element is decomposed into two parts with one part corresponding to stretching and bending deformation and the other part corresponding to shear deformation of the shell. Both parts of the stiffness are calculated with reduced integration rules, thereby improving computational efficiency. Shear strains are averaged on the reference surface such that neither locking phenomena nor any zero energy mode can occur. The satisfactory behaviour of the element is demonstrated in several numerical examples.

• 한국농공학회지
• /
• 제27권2호
• /
• pp.38-46
• /
• 1985

The reinforced concrete farm silos on the elastic foundatin are widely used in agricultural engineering because of their superior structural performance, economy and attractive appearance. Various methods for the analysis and design of farm silo, such as the analytical method, the finite difference method, and the finite element methods, can be used. But the analytical procedure can not be applied for the intricate conditions in practice. Therefore lately the finite element method has been become in the structural mechanics. In this paper, a method of finite element analysis for the cylindrical farm silo on ffness matrix for the elastic foundation governed by winkler's assumption. A complete computer programs have been developed in this paper can be applicable not only to the shell structures on elastic foundation but also to the arbitrary three dimensional structures. Assuming the small deflection theory, the membrane and plate bending behaviours of flat plate element can be assumed mutually uncoupled. In this case, the element has 5 degrees of freedom per node when defined in the local coordinate system. However, when the element properties are transformed to the global coordinates for assembly, the 6th degree of freedom should be considered. A problem arises in this procedure the resultant stiffness in the 6th degree of freedom at this node will be zero. But this singularity of the stiffness matrix can be eliminated easily by merely replacing the zero diagonal by dummy stiffness.

• Structural Engineering and Mechanics
• /
• 제21권5호
• /
• pp.539-551
• /
• 2005

Finite elements based on isoparametric formulation are known to suffer spurious stiffness properties and corresponding stress oscillations, even when care is taken to ensure that completeness and continuity requirements are enforced. This occurs frequently when the physics of the problem requires multiple strain components to be defined. This kind of error, commonly known as locking, can be circumvented by using reduced integration techniques to evaluate the element stiffness matrices instead of the full integration that is mathematically prescribed. However, the reduced integration technique itself can have a further drawback - rank deficiency, which physically implies that spurious energy modes (e.g., hourglass modes) are introduced because of reduced integration. Such instability in an existing stiffness matrix is generally detected by means of an eigenvalue test. In this paper we show that a knowledge of the dimension of the solution space spanned by the column vectors of the strain-displacement matrix can be used to identify the instabilities arising in an element due to reduced/selective integration techniques a priori, without having to complete the element stiffness matrix formulation and then test for zero eigenvalues.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1996년도 봄 학술발표회 논문집
• /
• pp.182-189
• /
• 1996

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear Programming formulation of the topology Problem is developed and examined. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.

• 대한기계학회논문집
• /
• 제13권3호
• /
• pp.317-322
• /
• 1989

본 연구에서는 Liu의 매트릭스를 강소성 유한요소법에 도입하여 통상의 소성가공 공정중에 있는 피가공물의 3차원 변형을 실용적인 수준에서 해석 가능케 하는 강소성 유한요소법을 도입하여 통상의 소성가공 공정중에 있는 피가공물의 3차원 변형을 실용적인 수준에서 해석 가능케 하는 강소성 유한요소법을 제안하고 실례를 통하여 제안한수법에 의하여 얻어진 해의 안정성과 계산효율을 검토한다.

• Structural Engineering and Mechanics
• /
• 제39권4호
• /
• pp.541-558
• /
• 2011

This paper presents an analysis of stochastic eigenvalue problem of plane bar structures. Particular attention is paid to the effect of spatial variations of the flexural properties of the structure on the first four eigenvalues. The problem of spatial variations of the structure properties and their effect on the first four eigenvalues is analyzed in detail. The stochastic eigenvalue problem was solved independently by stochastic finite element method (stochastic FEM) and Monte Carlo techniques. It was revealed that the spatial variations of the structural parameters along the structure may substantially affect the eigenvalues with quite wide gap between the two extreme cases of zero- and full-correlation. This is particularly evident for the multi-segment structures for which technology may dictate natural bounds of zero- and full-correlation cases.