• 한국전산구조공학회논문집
• /
• 제20권2호
• /
• pp.113-125
• /
• 2007

탄성지반 위의 비대칭 개/폐단면의 박벽보에 대한 탄성해석 및 안정성해석을 수행하기 위해 엄밀한 강성행렬을 계산하기 위한 개선된 수치해석 기법을 새롭게 제시한다. 본 연구에서 제시한 수치해석기법은 박벽보의 안정성 해석을 위한 엄밀한 강성행렬을 산정하는 선행된 수치해석기법의 결점을 보완하고 있다. 본 연구에서 제시한 기법은 일반화된 고유치 문제에 관한 해를 얻는 것으로서 일반화된 14개의 변위에 대한 고유치 문제를 평형방정식에 관한 1차의 연립상미분 방정식으로 변환함으로써 얻어진다. '0'의 고유치에 대응되는 변위파라미터에 대해 다항식이 가정되며 항등조건으로부터 '0'의 고유치의 수와 동일한 미결정된 파라미터를 포함하는 고유 모우드가 결정되고 이로부터 'non-zero'의 고유치와 다항식의 해를 조합함으로써 엄밀한 변위함수가 결정된다. 이후 부재력-변위의 관계를 이용하여 엄밀한 강성행렬을 산정하게 된다. 본 연구에서 개발한 수치해석 기법의 타당성을 검증하기 위해서 본 연구에서 제시한 이론에 의한 해를 제시하고 보요소 및 쉘요소을 사용한 유한요소해와 비교 검토한다.

• 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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 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.

• Structural Engineering and Mechanics
• /
• 제13권4호
• /
• pp.437-453
• /
• 2002

In this paper, dynamic stiffness and flexibility for circular membranes are analytically derived using an efficient mixed-part dual boundary element method (BEM). We employ three approaches, the complex-valued BEM, the real-part and imaginary-part BEM, to determine the dynamic stiffness and flexibility. In the analytical formulation, the continuous system for a circular membrane is transformed into a discrete system with a circulant matrix. Based on the properties of the circulant, the analytical solutions for the dynamic stiffness and flexibility are derived. In deriving the stiffness and flexibility, the spurious resonance is cancelled out. Numerical aspects are discussed and emphasized. The problem of numerical instability due to division by zero is avoided by choosing additional constraints from the information of real and imaginary parts in the dual formulation. For the overdetermined system, the least squares method is considered to determine the dynamic stiffness and flexibility. A general purpose program has been developed to test several examples including circular and square cases.

• Structural Engineering and Mechanics
• /
• 제65권1호
• /
• pp.91-98
• /
• 2018

Three structure-dependent integration methods with no numerical dissipation have been successfully developed for time integration. Although these three integration methods generally have the same numerical properties, such as unconditional stability, second-order accuracy, explicit formulation, no overshoot and no numerical damping, there still exist some different numerical properties. It is found that TLM can only have unconditional stability for linear elastic and stiffness softening systems for zero viscous damping while for nonzero viscous damping it only has unconditional stability for linear elastic systems. Whereas, both CEM and CRM can have unconditional stability for linear elastic and stiffness softening systems for both zero and nonzero viscous damping. However, the most significantly different property among the three integration methods is a weak instability. In fact, both CRM and TLM have a weak instability, which will lead to an adverse overshoot or even a numerical instability in the high frequency responses to nonzero initial conditions. Whereas, CEM possesses no such an adverse weak instability. As a result, the performance of CEM is much better than for CRM and TLM. Notice that a weak instability property of CRM and TLM might severely limit its practical applications.

• 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.

• 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.

• 한국해양공학회지
• /
• 제26권5호
• /
• pp.11-17
• /
• 2012

In this study, a damage detection method that uses sensitivity-based finite (FE) element model updating with the natural frequency and zero frequency was proposed. The stiffness matrix for a structure was modified using the sensitivity-based FE model updating method. A sensitivity analysis was used to update the FE model, and the natural frequencies and zero frequencies were considered as target parameters to supplement the information on the vibration characteristics. The locations and values of the damages were estimated from the modified stiffness matrix. Several numerical examples were considered to verify the performance of the proposed method.

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

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

• 한국농공학회지
• /
• 제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.