• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.4
• /
• pp.293-300
• /
• 2012

In this study, the delamination analysis has been implemented to investigate the initiation and propagation of crack in composite laminates composed of orthotropic materials. A simple modeling was achieved by moving nodal technique without re-meshing work when crack propagation occurred. This paper aims at achieving two specific objectives. The first is to suggest a very simple modeling scheme compared with those applied to conventional h-FEM based models. To verify the performance of the proposed model, analysis of double cantilever beams with composite materials was implemented and then the results were compared with reference values in literatures. The second one is to investigate the behavior of interior delamination problems using the proposed model. To complete these objectives, the full-discrete-layer model based on Lobatto shape functions was considered and energy release rates were calculated using three-dimensional VCCT(virtual crack closure technique) based on linear elastic fracture mechanics.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.481-494
• /
• 2001

An enhanced degenerated shell finite element (FE), which has been developed for inelastic analysis of reinforced concrete structures is described in this paper. Generally, Reissner-Mindlin (RM) assumptions are adopted to develop the degenerated shell FE so that transverse shear deformation effects is considered. However, it is found that there are serious defects such as locking phenomena in RM degenerated shell FE since the stiffness matrix has been overestimated in some situations. As remedies of locking phenomena, reduced integration, incompatible mode and assumed strain method have been used. Especially, the assumed strain method has been successfully used in many FEs. But contrarily, there is a few investigation on the performance of the assumed strains in the inelastic analysis of concrete structures. Therefore, shell formulation is provided in this paper with emphasis on the terms related to the stiffness matrix based on assumed strain method and microscopic concrete material model. Finally, the performance of the present shell element is tested and demonstrated with several numerical examples. From the numerical tests, the present result shows a good agreement with experimental data or other numerical results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3A
• /
• pp.447-455
• /
• 2006

A general purpose of 4-node co-rotational resultant shell element is developed for the solution of nonlinear problems of reinforced concrete, steel and fiber-reinforced composite structures. The formulation of the geometrical stiffness presented here is defined on the mid-surface by using the second order kinematic relations and is efficient for analyzing thick plates and shells by incorporating bending moment and transverse shear resultant forces. The present element is free of shear locking behavior by using the ANS (Assumed Natural Strain) method such that the element performs very well as thin shells. Inelastic behaviour of concrete material is based on the plasticity with strain hardening and elasto-plastic fracture model. The plasticity of steel is based on Von-Mises Yield and Ivanov Yield criteria with strain hardening. The transverse shear stiffness of laminate composite is defined by an equilibrium approach instead of using the shear correction factor. The proposed formulation is computationally efficient and versitile for most civil engineering application and the test results showed good agreement.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.711-719
• /
• 1994

This paper deals with the development of an efficient algorithm which can consider the bond-slip effect in numerical analysis of reinforced concrete structures. Unlike the classical bond-link or bond-zone element using double nodes, the proposed model is considering the bond-slip effect without taking double nodes by incorporation of the equivalent steel stiffness. Moreover after calculation of nodal displacements, the deformation of steel at each node can be found out through the back-substitution technique from first up to final steel element using a governing equation constructed based on the equilibrium at each node of steel and the compatibility condition between steel and concrete.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.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.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.31 no.4
• /
• pp.423-434
• /
• 1995

It was used boundary element method(BEM) and analysed axisymmetric problem to solve hydraulic cylinder for large vessel acting uniform internal pressure(25N/m super(2)) within elastic limit. This paper was utilized the carbon steel tubes for machine structural purposed model, inner radius was 150mm and outer radius was 250mm, axial length was semi-infinite and the isoparametric element was used. The important results obtained in this study were summarized as follows. Radial, tangential and shearing stress occured the maximum stresses(48, -20 and 34MPa) at the inner radius and the minimum stresses(32, -4 and 18MPa) at the outer radius of the hydraulic cylinder for large vessel. But negative signs have meaning compressive stress and stress diminution ratio was about 0.15MPa/mm. The use of isoparametric element raised accuracy and the increment of input data lessened the error in internal point but computer run-time was increased. The double node was improved the internal solutions to settle discontinuity at corner and the double exponential formula lessened error of stress value at boundary neighborhood. And then coincidence between the analytical and exact results is found to be fairly good, showing that the proposed analytical by BEM is reliable.

• Computational Structural Engineering
• /
• v.9 no.4
• /
• pp.141-154
• /
• 1996

Delaunay triangulation is a very powerful method of mesh generation for its versatility such as handling complex geometries, element density control, and local/global remeshing capability, The limit of generating simplex elements(3-node elements in 2-D) only is resolved by adding generation module of 6-node quadratic elements. Since proposed adjacency does not change from 3-node element mesh to 6-node mesh, generation module can utilize the original simplex element generator. Therefore, versatility of the Delaunay triangulation is preserved. A simple upsetting problem is employed to show the possibility of the algorithm.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.4 no.1
• /
• pp.35-50
• /
• 2000

본 논문에서는 반복하중을 받는 철근 콘크리트 쉘구조물의 해석을 위한 비선형 유한요소 해법을 제시하였다 유한 요소로서는 충상화기법을 이용한 부재회전강성도를 갖는 4절점 평면 쉘요소가 개발되었다 두께 방향에 대한 철근과 콘크리트의 재료성질을 고려하기 위하여 충상화기법이 도입되었다. 재료적 비선형성에 대해서는 균열콘크리트에 대한 인장, 압축, 전단모델과 콘크리트중에 있는 철근모델을 조합하여 고려하였다. 이에 대한 콘크리트의 균열모델로서는 분산균열모델을 사용하였으며 철근에 대해서는 1축 응력상태로가정하여 등가의 분산 분포된 철근량으로 모델화하였다 구성모델은 재하, 제하 그리고 재재하과정을 포함하여 요소는 반복하중하에서 철근콘크리트 쉘의 거동을 파악할 수 있다 신뢰성 있는 실험결과와 비교를 통하여 본 논문의 해석방법이 반복하중을 받는 철근콘크리트 쉘구조의 비선형 해석에 적합한 방법임을 입증하고자 한다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.307-315
• /
• 2008

Post-tensioned precast concrete system (PPS) consists of U-shaped precast wide beams and concrete column. The continuity of beam-column joint is provided with the topping concrete on the PC shell beam and post-tensioning. Nonlinear analysis was conducted, using ANSYS, a finite-element analysis program, to obtain data for determining the characteristics of the structure and to allow various parametric analyses for post-tensioned wide beam-column connections. In this analysis, the Solid 65 element was used, in which concrete element had 8 nodes and each node had 3 degrees of freedomIn. Solid 65, the shear-transfer factor reflects a decrease of shear strength for the positions with cracks, as an impact factor to make the analysis value approximate the experiment value. In this study, the behavior of test specineus were most closely predicted to the experimental results, when the shear-transfer coefficient 0.85 was used for a closed crack, and 0.2 was used for an open crack.