• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.615-626
• /
• 2010

The discrete-layer model is proposed to analyze the edge-effect problem of laminates under extension and flexure. Based on three-dimensional elasticity theory, the displacement fields of each layer in a laminate have been treated discretely in terms of three displacement components across the thickness. The displacement fields at bottom and top surfaces within a layer are approximated by two-dimensional shape functions. Then two surfaces are connected by one-dimensional high order shape functions. Thus the p-convergent refinement on approximated one- and two-dimensional shape functions can be implemented independently of each other. The quality of present model is mostly determined by polynomial degrees of shape functions for given displacement fields. For nodal modes with physical meaning, the linear Lagrangian polynomials are considered. Additional modes without physical meaning, which are created by increasing nodeless degrees of shape functions, are derived from integrals of Legendre polynomials which have an orthogonality property. Also, it is assumed that mapping functions are linear in the light of shape of laminated plates. The results obtained by this proposed model are compared with those available in literatures. Especially, three-dimensional out-of-plane stresses in the interior and near the free edges are evaluated and convergence performance of the present model is established with the stress results.

• International Journal of Automotive Technology
• /
• v.7 no.3
• /
• pp.353-360
• /
• 2006

The nominal structural stress calculation method proposed by Radaj has included some problems as follows: (a) How the value of the diameter D is decided in the method; (b) It is not possible to estimate nominal structural stress of the spot welded joints with the balanced sheet in-plane load that no general loads are obtained by FE shell analysis. In this paper, the new method for calculating nominal structural stress was proposed to solve above-mentioned problems. The proposed method calculates the nominal structural stress through the circular plate theory in theory of elasticity. This theoretical analysis uses not only general loads but also nodal displacements around spot welding provided by FE shell analysis as boundary condition. Fatigue test data of various spot-welded joints could be arranged in a narrow bandwidth on S-N chart using the nominal structural stresses calculated by proposed method. The fatigue life prediction method using the proposed method for calculating nominal structural stress is useful for the prior evaluation technique that can predict the fatigue life of spot welding by CAE.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.3
• /
• pp.344-352
• /
• 2002

In the analysis of the 3 dimensional plate structures by the finite element method, the triangular elements are generally used for the global stiffness matrix of the analyzed system. But the triangular elements of the plates have some problems in the process of formulation and in the precision of analysis. The formulation of the finite element method to analyze 3 dimensional plate structures using quadrilateral elements is presented in this paper. The degree of freedom off nodal point is 6, that is, the displacements in the direction off-y-z is and the rotations about x-y-z axis and then the degree of freedom off element is 24. For the comparison of the analysis using triangular elements and quadrilateral elements, the rectangular plates subjected to the uniform load and a concentrated load on the centroid of the plate, for which the theoretical solutions have been obtained, are analyzed. The calculated deflections of the rectangular plates using the finite element method by the triangular elements and the quadrilateral elements are also compared with the deflections of the plates calculated by theoretical solutions. The defections of the rectangular plates calculated by the finite element method using the quadrilateral elements are closer to the theoretical solutions than the defections calculated by the finite element method using the triangular elements. The deflection of the centroid of plate, calculated by the finite element method, converges to that of theoretical solution as the number of elements is increased. This convergence is much more rapid for the case of using the quakrilateral elements than fir the case of using triangular elements.

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

In this study we proposed a simplified finite element model of anchor bolt system inserted through concrete structures considering clamping forces. The three different finite element types using LS-DYNA are applied for numerical efficiency of the anchor bolt modeling. Combined beam and solid elements are used to reflect the tension state at internal part of anchor bolt due to torques. The clamping forces due to torques are considered by introducing a compression for a nut plane modeled by beam elements. The numerical examples show good agreement with different element types. Parametric studies are focused on the various effects of different element types on the induced axial and shear forces of anchor bolts considering clamping forces.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2008.11a
• /
• pp.985-988
• /
• 2008

A linear analysis method using reduced secant stiffness was developed for inelastic earthquake design of reinforced concrete structures. In the proposed method, the beam-column element and plane element, which are the same as used in conventional elastic analysis, are used for structural modeling. Based on the structural plastic mechanism intended by engineer, the distribution of inelastic members is determined. The secant stiffness of the inelastic members is determined based on the target ductility of the structure. Inelastic strengths of the members are calculated by using linear analysis on the structure modeled with secant stiffness. Plastic rotations in the inelastic members are calculated with the nodal rotations resulting from the secant stiffness analysis. For verification, the proposed method was applied to the inelastic earthquake designs of a moment-resisting frame and a dual system of two dimensions, and also a dual system of three dimensions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.121-131
• /
• 1986

The stiffness matrix of circularly curved frame elements including the warping effects is formulated by the solutions of vlasov's differential equations, and the procedure for the elastic static analysis of curved girder systems by the displacement method is presented. The validity of this method has been demonstrated by comparing the analysis results with other solutions. And if the tangential lines of the two frame element axes connected at any nodal point coincide, the transformation to the global coordinate system can be omitted when we analyze the structures consisting of circularly curved elements. The theory introduced in this thesis can be applied with sufficient accuracy to the structures built up with horizontally circular curved frame elements which have closed or open cross sections and are symmetric to the axis perpendicular to the plane of the curvature, such as prestressed concrete box girder bridges.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.73-84
• /
• 1993

In this research, configuration design optimization of plane truss structure has been tested by using decomposition technique. In the first level, the problem of transferring the nonlinear programming problem to linear programming problem has been effectively solved and the number of the structural analysis necessary for doing the sensitivity analysis can be decreased by developing stress constraint into member stress approximation according to the design space approach which has been proved to be efficient to the sensitivity analysis. And the weight function has been adopted as cost function in order to minimize structures. For the design constraint, allowable stress, buckling stress, displacement constraint under multi-condition and upper and lower constraints of the design variable are considered. In the second level, the nodal point coordinates of the truss structure are used as coordinating variable and the objective function has been taken as the weight function. By treating the nodal point coordinates as design variable, unconstrained optimal design problems are easy to solve. The decomposition method which optimize the section areas in the first level and optimize configuration variables in the second level was applied to the plane truss structures. The numerical comparisons with results which are obtained from numerical test for several truss structures with various shapes and any design criteria show that convergence rate is very fast regardless of constraint types and configuration of truss structures. And the optimal configuration of the truss structures obtained in this study is almost the identical one from other results. The total weight couldbe decreased by 5.4% - 15.4% when optimal configuration was accomplished, though there is some difference.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.109-122
• /
• 1983

An accurate thick beam element (TB4) which includes the effects of the shear deformation and rotatory inertia has been degenerated from the three dimensional continuum by employing the Timoshenko beam assumptions. The proposed TB4 element has four nodes and two degrees of freedom at each node, totally eight degrees of freedom. The transverse deflection W and plane rotation ${\theta}$ with the cubic interpolation functions are selected as nodal variables. The element characteristics are formulated by discretizing the beam equations of motion, using the Galerkin weighted residual method, and are numerically integrated by the reduced shear integration technique, using the three-point Gauss quadrature with the various shear coefficients. Several numerical examples are analyzed to demonstrate the accuracy and the monotonic convergence behavior of the proposed TB4 beam element. The result indicates that the TB4 element shows the more excellent performance and the monotonic convergence behavior than the other existing Timoshenko beam type elements for the whole range of the beam aspect ratios, in both static and free vibration analyses.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.463-474
• /
• 2009

Double symmetric patch repair of existing structures always causes membrane action only, however, in many cases this technique is not practical. On the other hand, the bending stiffness of the patch and the skin increases as tensile loading is increased and affects the bending deformation significantly in the case of single-sided patch repair. In this study, the p-convergent full layerwise model has been proposed to determine the stress concentration factor in the vicinity of a circular hole as well as across the thickness of plates with single-sided patch repair. In assumed displacement field, the strain-displacement relations and 3-D constitutive equations of a layer are obtained by the combination of 2-D and 3-D hierarchical shape functions. The transfinite mapping technique has been used to represent a circular boundary and Gauss-Lobatto numerical integration is implemented in order to directly obtain stresses occurred at the nodal points of each layer without other extrapolation techniques. The accuracy and simplicity of the present model are verified with comparison of the previous results in literatures using experiment and conventional 3-D finite element. Also, the bending effect has been investigated with various patch types like square, circular and annular shape.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.9
• /
• pp.1041-1050
• /
• 2011

An elasto-plastic finite element method using the theory of strain gradient plasticity is proposed to evaluate the size dependency of structural plasticity that occurs when the configuration size decreases to micron scale. For this method, we suggest a low-order plane and three-dimensional displacement-based elements, eliminating the need for a high order, many degrees of freedom, a mixed element, or super elements, which have been considered necessary in previous researches. The proposed method can be performed in the framework of nonlinear incremental analysis in which plastic strains are calculated and averaged at nodes. These strains are then interpolated and differentiated for gradient calculation. We adopted a strain-gradient-hardening constitutive equation from the Taylor dislocation model, which requires the plastic strain gradient. The developed finite elements are tested numerically on the basis of typical size-effect problems such as micro-bending, micro-torsion, and micro-voids. With respect to the strain gradient plasticity, i.e., the size effects, the results obtained by using the proposed method, which are simple in their calculation, are in good agreement with the experimental results cited in previously published papers.