• 한국해양공학회지
• /
• 제12권2호통권28호
• /
• pp.1-21
• /
• 1998

In this paper, a non-conforming 3-D 8-node solid element(MQM10) has beets applied to the analyses of static and dynamic bending problems of laminated composite plates The QM10 element exhibits stiffer bending stiffness which is caused by the reduction of degree of freedom from Q11 element. As an effective way to correct the relative stiffness stiffening phenomenon the modification of Gauss sampling points for composite plates is proposed. The quantity of modification is a function of material properties. Also, another two modified equations are obtained, one is modification for stress, and the other is modification of coefficient of shear modulus in free vibration. It is noted that MQM10 element can analyse the static and free vibration problems of various 3-dimensional composite plates composed of unidirectional laminae, woven laminae or braided laminae. The results of MQM10 element are in good agreement with those of 20-node element.

• 전산구조공학
• /
• 제10권4호
• /
• pp.117-130
• /
• 1997

본 논문에서는 등방성판의 인장이나 전단변형은 물론 굽힘문제에도 적용할 수 있는 3차원 고체요소들 중에서 최소의 자유도를 갖는 수정 10절점 상당요소를 제안하였다. 제안된 수정 10절점 상당요소는 Q11요소나 20절점요소로부터 자유도가 줄어듬에 기인한 과대한 굽힘강성을 나타낸다. 이러한 상대적 강성과잉 현상을 수정하기 위한 효과적인 방법으로 가우스 적분점 수정 방법을 제안하였다. 수정량은 포아송 비의 함수이다. 수정 10절점 상당요소의 효과를 여러 가지 예에 적용하여 검증하였다. 제안된 수정 10절점 상당요소에 의한 등방성판의 정적해석과 자유진동 해석의 결과들은 20절점요소를 사용한 결과들과 잘 일치하였다.

• Structural Engineering and Mechanics
• /
• 제34권6호
• /
• pp.779-799
• /
• 2010

In order to consider high-order effects on the actual limit state function, a new response surface method is proposed for structural reliability analysis by the use of high-order approximation concept in this study. Hermite polynomials are used to determine the highest orders of input random variables, and the sampling points for the determination of highest orders are located on Gaussian points of Gauss-Hermite integration. The cross terms between two random variables, only in case that their corresponding percent contributions to the total variation of limit state function are significant, will be added to the response surface function to improve the approximation accuracy. As a result, significant reduction in computational cost is achieved with this strategy. Due to the addition of cross terms, the additional sampling points, laid on two-dimensional Gaussian points off axis on the plane of two significant variables, are required to determine the coefficients of the approximated limit state function. All available sampling points are employed to construct the final response surface function. Then, Monte Carlo Simulation is carried out on the final approximation response surface function to estimate the failure probability. Due to the use of high order polynomial, the proposed method is more accurate than the traditional second-order or linear response surface method. It also provides much more efficient solutions than the available high-order response surface method with less loss in accuracy. The efficiency and the accuracy of the proposed method compared with those of various response surface methods available are illustrated by five numerical examples.

• 대한기계학회논문집
• /
• 제18권10호
• /
• pp.2640-2652
• /
• 1994

The stiffness of 6-node isotropic element is stiffer than that of 8-node isotropic element of same configuration. This phenomenon was called 'Relative Stiffness Stiffening Phenomenon'. In this paper, an equation of sampling point modification which correct this phenomenon was derived for the composite plate, as well as an equation for an isotropic plate. The relative stiffness stiffening phenomena of an isotropic plate element could be corrected by modifying Gauss sampling points in the numerical integration of stiffness matrix. This technique could also be successfully applied to the static analyses of composite plate modeled by the 3-dimensional 16-node elements. We predicted theoretical errors of stiffness versus the number of layers that result from the reduction of numerical integration order. These errors coincide very well with the actual errors of stiffness. Therefore, we can choose full integration of reduced integration based upon the permissible error criterion and the number of layers by using the thoretically predicted error.

• Steel and Composite Structures
• /
• 제23권1호
• /
• pp.41-51
• /
• 2017

A $C^0$ FE model developed based on an efficient higher order zigzag theory is used for hygrothermal analysis of laminated composite plates. The $C^0$ FE model satisfies the inter-laminar shear stress continuity at the interfaces and zero transverse shear stress conditions at plate top and bottom. In this model the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of $C^1$ continuity associated with the above plate theory. In the present theory the above mentioned $C^0$ continuity of the present element is compensated in the stiffness matrix formulation by using penalty parameter approach. In order to avoid stress oscillations observed in the displacement based finite element, the stress field derived from temperature/moisture fields (initial strains) must be consistent with total strain field. Special steps are introduced by field consistent approach (e.g., sampling at gauss points) to compensate this problem. A nine noded $C^0$ continuous isoparametric element is used in the proposed FE model. Comparison of present numerical results with other existing solutions shows that the proposed FE model is efficient, accurate and free of locking.

• 대한기계학회논문집A
• /
• 제20권4호
• /
• pp.1320-1332
• /
• 1996

In the present study, we consider modified 5-node equivalent solid element which has smallest degree of freedom among 2-dimensional solid elements accounting bending deformation as well as extensional and shear deformations, We shall investigate static and dynamic characteristics of this element, which is very effective in thin beam, thick beam, large displacement problems, beam of variable thickness, and asymmetrically stepped beam, etc., as well as relatively simple problems of beam. The degree of freedom of this element is 10, which is smaller than 18 of 9-node element, 16 of 8-node elemtns, 12 of modified 6-node element and Q6 element. Therefore, this element is expected to broaden the effective range of application of the solid elements in the bending problems further.

• 전산구조공학
• /
• 제8권4호
• /
• pp.87-97
• /
• 1995

2차원 유한요소 모델의 동일한 형상과 하중 조건에 있어서 6절점 요소의 굽힘 강성은 8절점 요소의 굽힘 강성보다 더 크게 나타난다. 이와 같은 현상은 3차원 16절점 요소와 20절점 요소에서도 나타나며, 완전 요소의 중간 절점들을 제거하므로 인하여 나타난다. 따라서 이 현상을 상대적 강성강화 현상이라 할 수 있다. 강성강화 현상을 보정하기 위한 매우 효과적인 방법으로 가우스 적분점 수정법을 도출하였으며, 이 방법은 확장적인 강성과 같이 다른 종류의 강성을 변화시키지 않으며, 또한 패취시험을 통과하였다. 적분점 수정량은 재료의 포아송비의 함수로 나타나며, 2차원 평면응력 상태와 평면변형율 상태에 대한 두개의 수정식을 구하였고, 또한 3차원 고체요소에 대하여 확장하였다. 가우스 적분점 수정법의 효과를 검증하기 위하여 보와 판의 자유 및 강제운동 문제를 해석하였으며, 등방성 적층 보와 판에 대해서도 단층보와 단층판과 같은 방법으로 적용하여 그 효율성을 입증하였다.