• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.21-30
• /
• 2004

This study deals with free vibrations of laminated composite structures with a channel section using finite element method. In this paper, the mixed finite element method using Lagrangian and Hermite interpolation functions is adopted and a high-order plate theory is used to analyze laminated composite non-prismatic folded plates with a channel section more accurately for free vibration. The theory accounts for parabolic distribution of the transverse shear stress and requires no shear correction factors supposed in the first-order plate theory. An 32×32 matrix is assembled to transform the system element matrices from the local to global coordinates using a coordinate transformation matrix, in which an eighth drilling degree of freedom (DOF) per node is appended to the existing 7-DOF system. The results in this study are compared with those of available literatures for the conventional and first-order plate theory. Sample studies are carried out for various layup configurations and length-thickness ratio, and geometric shapes of plates. The significance of the high-order plate theory in analyzing complex composite structures with a channel section is enunciated in this paper.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.1A
• /
• pp.21-33
• /
• 2006

The reliability analysis can be conducted more effectively by formulating the stochastic finite element method suitable for the reliability theory about the cable stayed bridge. After conducting the initial equilibrium analysis of the cable stayed bridge, the program which can conduct the linear and nonlinear stochastic finite element analysis using the perturbation method and the reliability analysis considered to the correlation of the random variable is developed. Using the results of this program about the cable stayed bridge, the characteristic of the node displacement, element force and cable tension according to the correlation of the random variable is investigated quantitatively. Also the reliability index and the failure probability are examined by the compounding the correlation of the random variable.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.15 no.6
• /
• pp.118-126
• /
• 2011

The objective of this paper is to study the structural behavior of Composite Basement Wall (CBW) according to shear span-to-depth ratio through an experiment and predict the nonlinear behavior of CBW by using ADINA program widely has been being used for FE analysis. Especially, this study focuses on the part of CBW in which the Reinforced Concrete (RC) is under compression stress; At the region of CBW around each floor, RC part stresses by compressive force when lateral press by soil acts on the wall. The contact condition between RC wall and steel (H-Pile) including stud connector is main factor in the analysis since it governs overall structural behavior. In order to understand the structural behavior of CBW whose RC part is under compressive stress, an experimental work and finite element analysis were performed. Main parameter in the test is shear span-to-depth ratio. For simplicity in analysis, reinforcements were not modeled as a seperated element but idealized as smeared to concrete. All elements were modeled to have bi-linear relation of material properties. Three type of contact conditions such as All Generate Option (AGO), Same Element Group Option with Tie(SEGO-T) and Same Element Group Option with Not tie(SEGO-NT) were considered in the analysis. For each analysis, the stress flow and concentration were reviewed and analysis result was compared to test one. From the test result, CBW represented ductile behavior by contribution of steel member even if it had short shear span-to-depth ration which is close to "1". The global composite behavior of CBW whose concrete wall was under compressive stress could be predicted by using contact element in ADINA program. Especially, the modeling by using AGO and SEGO-T showed more close relation on comparing with test result.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

• Magazine of the Korea Concrete Institute
• /
• v.2 no.4
• /
• pp.91-98
• /
• 1990

Joints are provided in cement concrete pavements to control transverse and longitudinal cracking that occur due to restrained deformations caused by moisture and temperature variations in the slab. But the constuction of joints reduces the load-carrying capacity of the pavement at the joints, and pavements have been deteriorated by cracks at the slab edges along the joints due to traffic loads. Therefore, it is important to analyze the behavior of joints accurately in the design of cement concrete pavements. In this study, the mechanical behavior of cement concrete pavement slabs is analyzed by the plate-finite element model, and Winkler foundation model is adopted to analyze the subgrades. The load transfer mechan¬ism of joints are composed of dowel action, aggregate interlocking, and tied-key action, and the analytical pro¬gram is developed using these joint models. Using this numerical model as an analysis tool, the effects of joint parameters on the behavior of pavements are investigated.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.6
• /
• pp.19-30
• /
• 2005

A coupled three-dimensional pile group analysis method was developed by considering complex behavior of sub-structures (pile-soil-cap) which included soil nonlinearity and the behavior of super-structure (pier). As an intermediate analysis method between FBPier 3.0 and Group 0.0, it took advantages of each method. Among the components of a pile group, individual piles were modeled with stiffness matrices of pile heads and soils with nonlinear load-transfer curves (t-z, q-z and p-y curves). A pile cap was modeled with modified four-node flat shell elements and a pier with three-dimensional beam element, so that a unified analysis could be possible. A nonlinear analysis method was proposed in this study with a mixed incremental and iteration techniques. The proposed method for a pile group subjected to axial and lateral loads was compared with othe. analytical methods (i.e., Group 6.0 and FBPier 3.0). It was found that the proposed method could predict the complex behavior of a pile group well, even though piles were modelled simply in this study by using pile head stiffness matrices which were different from the method introduced in FBPier 3.0.

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

Fiber-reinforced composite materials due to their high specific strength, high stiffness and light weight are becoming increasingly used in many engineering industry, especially in the aerospace, marin and civil, etc. In this paper, the buckling load and mode shapes of composite laminates with elliptical cross-section including transverse shear deformations are analyzed. For solving this problems, a versatile flat shell element has been developed by combining a membrane element with drilling degree-of-freedom and a plate bending element. Also, an improved shell element has been established by the combined use of the addition of enhanced assumed strain and the substitute shear strain fields. The combined influence of shell geometry and elliptical cross-sectional parameter, fiber angle, and lay-up on the buckling loads of elliptical cylinder is examined. The critical buckling loads and mode shapes analyzed here may serve as a benchmark for future investigations.