• Interaction and multiscale mechanics
• /
• v.2 no.1
• /
• pp.69-89
• /
• 2009

Reinforced and prestressed concrete (RC and PC) thin walls are crucial to the safety and serviceability of structures subjected to shear. The shear strengths of elements in walls depend strongly on the softening of concrete struts in the principal compression direction due to the principal tension in the perpendicular direction. The past three decades have seen a rapid development of knowledge in shear of reinforced concrete structures. Various rational models have been proposed that are based on the smeared-crack concept and can satisfy Navier's three principles of mechanics of materials (i.e., stress equilibrium, strain compatibility and constitutive laws). The Cyclic Softened Membrane Model (CSMM) is one such rational model developed at the University of Houston, which is being efficiently used to predict the behavior of RC/PC structures critical in shear. CSMM for RC has already been implemented into finite element framework of OpenSees (Fenves 2005) to come up with a finite element program called Simulation of Reinforced Concrete Structures (SRCS) (Zhong 2005, Mo et al. 2008). CSMM for PC is being currently implemented into SRCS to make the program applicable to reinforced as well as prestressed concrete. The generalized program is called Simulation of Concrete Structures (SCS). In this paper, the CSMM for RC/PC in material scale is first introduced. Basically, the constitutive relationships of the materials, including uniaxial constitutive relationship of concrete, uniaxial constitutive relationships of reinforcements embedded in concrete and constitutive relationship of concrete in shear, are determined by testing RC/PC full-scale panels in a Universal Panel Tester available at the University of Houston. The formulation in element scale is then derived, including equilibrium and compatibility equations, relationship between biaxial strains and uniaxial strains, material stiffness matrix and RC plane stress element. Finally the formulated results with RC/PC plane stress elements are implemented in structure scale into a finite element program based on the framework of OpenSees to predict the structural behavior of RC/PC thin-walled structures subjected to earthquake-type loading. The accuracy of the multiscale modeling technique is validated by comparing the simulated responses of RC shear walls subjected to reversed cyclic loading and shake table excitations with test data. The response of a post tensioned precast column under reversed cyclic loads has also been simulated to check the accuracy of SCS which is currently under development. This multiscale modeling technique greatly improves the simulation capability of RC thin-walled structures available to researchers and engineers.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.1
• /
• pp.89-94
• /
• 2014

The high Young's modulus and tensile strength of carbon nanotubes has attracted great attention from the research community given the potential for developing super-strong, super-stiff composites with carbon nanotube reinforcements. Over the decades, the strength and stiffness of carbon nanotube-reinforced polymer nanocomposites have been researched extensively. However, unfortunately, such strong composite materials have not been developed yet. It has been reported that the efficiency of load transfer in such systems is critically dependent on the quality of adhesion between the nanotubes and the polymer chains. In addition, the waviness and orientation of the nanotubes embedded in a matrix reduce the reinforcement effectiveness. In this study, we carried out performed micromechanics-based numerical modeling and analysis by varying the geometry of carbon nanotubes including their aspect ratio, orientation, and waviness. The results of this analysis allow for a better understanding of the load transfer capabilities of carbon nanotube-reinforced polymer composites.

• Journal of the Korean Ceramic Society
• /
• v.40 no.8
• /
• pp.784-790
• /
• 2003

Present investigation shows the analysis results for ceramic reinforced metal matrix composite under uniaxial transverse tensile loading. The resulting deformation, the projected damage type, and stress-strain behavior were computed depending on microstructure details such as the type of periodic reinforcement array, and the type of interface bonding. A two-dimensional finite element analysis was conducted based on the unit-cell of square, hexagonal, or diagonal periodic away For composite with strong interface bonding, the transverse stress vs. strain curve was generally increased with the increase of the ceramic volume fraction. For the composite with weakly bonded interface, however, the transverse stress vs. strain curve was reduced against the ceramic volume fraction. The decrease was caused by the interface debonding-induced stiffness reduction of the composite. For the composite of weakly bonded interface, the relative reduction rate in the final limit stress for hexagonal array was larger than that for square array. Outcome of the present study was compared favorably with the published literature data.

• Proceedings of the KSR Conference
• /
• 2001.10a
• /
• pp.500-507
• /
• 2001

In this study, Interaction among each element was construed by the unit of new stiffness matrix to analyze the finite element about Japanese precast concrete slab track and improved slab track. Dynamic analysis which is assumed a static analysis and a trainload on the transverse and the longitudinal load of the train into a series periodic function was performed by using the common program. And then, the difference of the movement between an improved section and an existing structure type was realized. Longitudinal static analysis indicated that the stress of the improved section is smaller than that of the protrusion of the existing slab track. And static and dynamic analysis on transverse load showed a little decrease of the displacement on new slab track. But the dynamic analysis result showed that new track system was considerably decreased by 30% compared with the existing Japanese slab track.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.67-78
• /
• 2012

A three-dimensional approximate computer-based method, YSPR (Yonsei Piled Raft), was developed for analysis of behavior of piled raft foundations. The raft was modeled as a flat shell element having 6 degrees of freedom at each node and the pile was modeled as a beam-column element. The behaviors of pile head and soil were controlled by using $6{\times}6$ stiffness matrix. To model the non-linear behavior, the soil-structure interaction between soil and pile was modeled by using nonlinear load-transfer curves (t-z, q-z and p-y curves). Comparison with previous model and FEM analysis showed that YSPR gave similar load-displacement behaviors. Comparison with field measurement also indicated that YSPR gave a reasonable result. It was concluded that YSPR could be effectively used in analysis and design of piled raft foundations.

• Structural Engineering and Mechanics
• /
• v.13 no.5
• /
• pp.557-568
• /
• 2002

For the evaluation of the capability of a tubular member of an offshore structure to absorb the collision energy, a simple method can be employed for the collision analysis without performing the detailed analysis. The most common simple method is the rigid-plastic method. However, in this method any characteristics for horizontal movement and rotation at the ends of the corresponding tubular member are not included. In a real structural system of an offshore structure, tubular members sustain a certain degree of elastic support from the adjacent structure. End fixity has influences in the behaviors of a tubular member. Three-dimensional FEM analysis can include the effect of end fixity fully, however in viewpoints of the inherent computational complexities of the 3-D approach, this is not the recommendable analysis at the initial design stage. In this paper, influence of end fixity on the behaviors of a tubular member is investigated, through a new approach and other approaches. A new analysis approach that includes the flexibility of the boundary points of the member is developed here. The flexibility at the ends of a tubular element is extracted using the rational reduction of the modeling characteristics. The property reduction is based on the static condensation of the related global stiffness matrix of a model to end nodal points of the tubular element. The load-displacement relation at the collision point of the tubular member with and without the end flexibility is obtained and compared. The new method lies between the rigid-plastic method and the 3-demensional analysis. It is self-evident that the rigid-plastic method gives high strengthening membrane effect of the member during global deformation, resulting in a steeper slope than the present method. On the while, full 3-D analysis gives less strengthening membrane effect on the member, resulting in a slow going load-displacement curve. Comparison of the load-displacement curves by the new approach with those by conventional methods gives the figures of the influence of end fixity on post-yielding behaviors of the relevant tubular member. One of the main contributions of this investigation is the development of an analytical rational procedure to figure out the post-yielding behaviors of a tubular member in offshore structures.

• Structural Engineering and Mechanics
• /
• v.4 no.5
• /
• pp.553-568
• /
• 1996

Finite element stiffness matrix methods are presented for finding natural frequencies (or buckling loads) and modes of repetitive structures. The usual approximate finite element formulations are included, but more relevantly they also permit the use of 'exact finite elements', which account for distributed mass exactly by solving appropriate differential equations. A transcendental eigenvalue problem results, for which all the natural frequencies are found with certainty. The calculations are performed for a single repeating portion of a rotationally or linearly (in one, two or three directions) repetitive structure. The emphasis is on rotational periodicity, for which principal advantages include: any repeating portions can be connected together, not just adjacent ones; nodes can lie on, and members along, the axis of rotational periodicity; complex arithmetic is used for brevity of presentation and speed of computation; two types of rotationally periodic substructures can be used in a multi-level manner; multi-level non-periodic substructuring is permitted within the repeating portions of parent rotationally periodic structures or substructures and; all the substructuring is exact, i.e., the same answers are obtained whether or not substructuring is used. Numerical results are given for a rotationally periodic structure by using exact finite elements and two levels of rotationally periodic substructures. The solution time is about 500 times faster than if none of the rotational periodicity had been used. The solution time would have been about ten times faster still if the software used had included all the substructuring features presented.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.13 no.3
• /
• pp.21-28
• /
• 2009

The precise analysis of soil-pile-structure interaction requires a proper description of soil layer, pile, and structure. In commonly used finite element simulations, mesh boundaries should match the material discontinuity line. However, in practice, the geometry of soil profiles and piles may be so complex that mesh alignment becomes a wasteful and difficult task. To overcome these difficulties, a different integration method is adopted in this paper, which enables easy integration over a regular element with material discontinuity regardless of the location of the discontinuity line. By applying this integration method, the mesh can be generated rapidly and in a highly structured manner, leading to a very regular stiffness matrix. The influence of the shape of the soil profile and piles on the response is examined, and the validity of the proposed soil-pile structure interaction analysis method is demonstrated through several examples. It is seen that the proposed analysis method can be easily used on soil-pile-structure interaction problems with complex interfaces between materials to produce reliable results regardless of the material discontinuity line.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.307-314
• /
• 2009

The present paper briefly describes the space frame element and the fundamental strategies in computational elastic bifurcation theory of geometrically nonlinear, single load parameter conservative elastic spatial structures. A method for large deformation(rotation) analysis of space frame is based on an eulerian formulation, which takes into consideration the effects of large joint translations and rotations with finite deformation(rotation). The local member force-deformation relationships are based on the beam-column approach, and the change in member chord lengths caused by axial strain and flexural bowing are taken into account. and the derived geometric stiffness matrix is unsymmetric because of the fact that finite rotations are not commutative under addition. To detect the singular point such as bifurcation point, an iterative pin-pointing algorithm is proposed. And the path switching mode for bifurcation path is based on the non-negative eigen-value and it's corresponding eigen-vector. Some numerical examples for bifurcation analysis are carried out for a plane frame, plane circular arch and space dome structures are described.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.1
• /
• pp.21-30
• /
• 1985

This paper deals whth the buckling as well as postbuckling analysis of axisymmertric shells taking the initial deflection effects into account. Incremental equilibrium equations, based on the principle of virtual work, were derived by the finite element method, the successive step-by-step Newton-Raphson iterative technique was adopted. To define the transition pattern of postbuckling behavior from the prebuckling state more accurately, a simple solution method was developed, i.e. the critical load was calculated by the load extrapolation method with the determinant of tangent stiffness matrix and the equilibrium configuration in the immediate postbuckling stage was obtained by perturbation scheme and eigenvalue analysis. Degenerated isoparametric shell elements were used to analyse the axisymmetric shell of revolution. And by the method developed in this paper, the computer program applicable to the nonlinear analysis of both thin and moderately thick shells was constructed. To verify the capabilities and accuracies of the present solution method, the computed results were compared with the results of analytical solutions. These results coincided fairly well in both the small deflection and large deflection ranges. Various numerical analyses were done to show the effect of initial deflection and shape of shells on buckling load and postbuckling behavior. Futhermore, corrected directions of applied loads at every increment steps were used to determine the actual effects of large deflection in non-conservative load systems such as hydrostatic pressure load. The following conclusions can be obtained. (1) The method described in this paper was found to be both economic and effective in calculating buckling load and postbuckling behavior of shell structure. (2) Buckling and postbuckling behavior of spherical caps is critically dependent upon their geometric configuration, i.e. the shape of spherical cap and quantities of the initial deflection. (3) In the analysis of large deflection problems of shells by the incremental method, corrections of the applied load directions are needed at every incremental step to compensate the follower force effects.