• Journal of the Korea Concrete Institute
• /
• v.18 no.2 s.92
• /
• pp.189-198
• /
• 2006

The composition of most engineering materials is heterogeneous at some degree. It is simply a question of scale at which the level of heterogeneity becomes apparent. In the case of cementitious granular materials such as concrete the heterogeneity appears at the mesoscale where it is comprised of aggregate particles, a hardened cement paste and voids. Since it is difficult to consider each separate particle in the topological description explicitly, numerical models of the meso-structure are normally confined to two-phase matrix particle composites in which only the larger inclusions are accounted for. 2-D and 3-D concrete blocks(Representative Volume Element, RVE) are used to simulating heterogeneous concrete meso-structures in the form of aggregates in the hardened mortar with nearly zero-thickness linear or planar interfaces. The numerical sensitivity of these meso-structures are Investigated with respect to the different morphologies of heterogeneity and the different level of coupling constant among fracture mode I, II and III. In addition, a numerically homogenized concrete block in 3-D using Hashin-Shtrikman variational bounds provides an evidence of the effective cracking paths which are quite different with those of heterogenous concrete block. However, their average force-displacement relationship show a pretty close match each other.

• Korean Journal of Materials Research
• /
• v.28 no.11
• /
• pp.671-679
• /
• 2018

Finite element analyses are carried out to understand the piezoelectric behaviors of ZnO nanowires. Three different types of ZnO nanowires, with aspect ratios of 1:2. 1:31, and 1:57, are analyzed for uniaxial compression, pure bending, and buckling. Under the uniaxial compression with a strain of $1.0{\times}10^{-4}$ as the reference state, it is predicted that all three types of nanowires develop the same magnitude of the piezoelectric fields, which suggests that longer nanowires exhibit higher piezoelectric potential. However, this prediction is not in agreement with the experimental results previously reported in the literature. Such discrepancy is understood when the piezoelectric behaviors under bending and buckling are considered. When only the strain field due to bending is present in bending or buckling, the antisymmetric nature of the through-thickness stain distribution indicates that two piezoelectric fields, the same in magnitude and opposite in sign, develop along the thickness direction, which cancels each other out, resulting in a zero net piezoelectric field. Once additional strain contribution due to axial deformation is superposed on the bending, such field cancelling is compensated for due to the axial component of the piezoelectric field. Such numerical predictions seem to explain the reported experimental results while providing a guideline for the design of nanowire-based piezoelectric devices.

• Structural Engineering and Mechanics
• /
• v.71 no.5
• /
• pp.469-483
• /
• 2019

In this paper, an improved mathematical model is presented for the bending analysis of doubly curved functionally graded material (FGM) sandwich rhombic conoids. The mathematical model includes expansion of Taylor's series up to the third degree in thickness coordinate and normal curvatures in in-plane displacement fields. The condition of zero-transverse shear strain at upper and lower surface of rhombic conoids is implemented in the present model. The newly introduced feature in the present mathematical model is the simultaneous inclusion of normal curvatures in deformation field and twist curvature in strain-displacement equations. This unique introduction permits the new 2D mathematical model to solve problems of moderately thick and deep doubly curved FGM sandwich rhombic conoids. The distinguishing feature of present shell from the other shells is that maximum transverse deflection does not occur at its center. The proposed new mathematical model is implemented in finite element code written in FORTRAN. The obtained numerical results are compared with the results available in the literature. Once validated, the current model was employed to solve numerous bending problems by varying different parameters like volume fraction indices, skew angles, boundary conditions, thickness scheme, and several geometric parameters.

• Structural Engineering and Mechanics
• /
• v.46 no.2
• /
• pp.269-294
• /
• 2013

This paper is concerned with the determination of exact buckling loads and vibration frequencies of variable thickness isotropic plates using well known finite difference technique. The plates are subjected to uni, biaxial compression and shear loadings and various combinations of boundary conditions are considered. The buckling load is found out as the in plane load that makes the determinant of the stiffness matrix equal to zero and the natural frequencies are found out by carrying out eigenvalue analysis of stiffness and mass matrices. New and exact results are given for many cases and the results are in close agreement with the published results. In this paper, like finite element method, finite difference method is applied in a very simple manner and the application of boundary conditions is also automatic.

• Computers and Concrete
• /
• v.24 no.1
• /
• pp.51-61
• /
• 2019

Compressive ability is one of the most important mechanical properties of concrete material. The compressive failure process of concrete is pretty complex with internal tension, shear damage and friction between cracks. To simulate the complex fracture process of concrete at meso level, methodology for meso-structural analysis of concrete specimens is developed; the zero thickness cohesive elements are pre-inserted to simulate the crack initiation and propagation; the constitutive applied in cohesive element is established to describe the mechanism of crack separation, closure and friction behavior between the fracture surfaces. A series of simulations were carried out based on the model proposed in this paper. The results reproduced the main fracture and mechanical feature of concrete under compression condition. The effect of key material parameters, structure size, and aggregate content on the concrete fracture pattern and loading carrying capacities was investigated. It is found that the inner friction coefficient has a significant influence on the compression character of concrete, the compression strength raises linearly with the increase of the inner friction coefficient, and the fracture pattern is sensitive to the mesostructure of concrete.

• Geotechnical Engineering
• /
• v.12 no.3
• /
• pp.35-48
• /
• 1996

Reinforced soil structures may experience large local movements between soil and reinforcement. The failure modes of a reinforced structure depend on several factors which are governed by deformation and slipping of the reinforcement. In some cases, pulling out of the reinforcement may occur instead of rupturing, The growing use of geosynthetic liner system for storage of solid and liquid wastes has led to a number of slope instability problems where the synthetic liner may undergo a large amount of stretching and slipping as a result of the loading. The conventional finite element model for the soil-reinforcement interface uses a zero thickness joint element with normal and shear stiffnesses and can only accommodate a small amount of deformation. When a large slippage occurs, the model provides an i ncorrect mechanism for deformation. This paper presents a new interface finite element model which is able to simulate a large amount of slippage between soil and reinforcement. The formulation of the model is presented and the capability of the model is demonstrated using illustrative examples.

• Proceedings of the Korea Society for Simulation Conference
• /
• 2001.10a
• /
• pp.215-219
• /
• 2001

The resonators made of piezoelectric crystals such as a quartz crystal are widely used. Their frequency-temperature characteristics are of primary importance for their applications to the frequency control devices. The characteristics estimation is useful for determining their design parameters. In the present paper, several types of resonators are numerically analyzed. The numerical solutions are made using 3-D Finite Element Modeling, and the results are compared with the theoretical values whenever they are available. To demonstrate the validity of the present numerical approach, the application is made to the analysis of the plates with some well-established cutting angles. For the resonator stable with temperature change, the cutting angle is important in which the temperature coefficient of the first order is chosen to be zero. The rotated Y-cut plates in thickness-shear mode are considered. The equivalent circuit representation is often used fur describing the characteristics at the electrical terminals which enables the circuit analysis including the effect of temperature change by using the circuit simulators. The equivalent circuit parameters are obtained by fitting the admittance-frequency curve from the finite element analysis.

• Earthquakes and Structures
• /
• v.17 no.3
• /
• pp.257-270
• /
• 2019

The purpose of this work is to analyze the bending behavior of laminated composite plates using the refined fourvariable theory and the finite element method approach using an ANSYS 12 computational code. The analytical model is based on the multilayer plate theory of shear deformation of the nth-order proposed by Xiang et al 2011 using the theory principle developed by Shimpi and Patel 2006. Unlike other theories, the number of unknown functions in the present theory is only four, while five or more in the case of other theories of shear deformation. The formulation of the present theory is based on the principle of virtual works, it has a strong similarity with the classical theory of plates in many aspects, it does not require shear correction factor and gives a parabolic description of the shear stress across the thickness while filling the condition of zero shear stress on the free edges. The analysis is validated by comparing results with those in the literature.

• International Journal of Highway Engineering
• /
• v.11 no.3
• /
• pp.151-161
• /
• 2009

This study was conducted to investigate the characteristics of the curling behavior of long-span pavement slabs under environmental loads. By developing and using finite element models of the long-span pavement slabs, the stress distribution and the effects of slab length, slab thickness, stiffness of underlying layers, and the restraints of the slab ends on the curling behavior were analyzed. In addition, the field experiments were performed with the actual long-span pavement slab to obtain the curling behavior of the real structure under environmental loads. As a result of this study, it was found that the vertical displacements of the long-span pavement slab along the centerline due to the curling behavior were zero except for the areas near the slab ends, and the curling stresses were maximum and constant where the displacements were zero. The slab length and the stiffness of underlying layers did not affect the maximum curling stresses. The restraints at the slab ends made the curling stresses occur near the slab ends, but did not much affect the maximum curling stresses.

• Structural Engineering and Mechanics
• /
• v.33 no.5
• /
• pp.621-633
• /
• 2009

In this paper we compare the solution of the one-dimensional beam model and the numerical solution of the two-dimensional linearized elasticity problem for rectangular domain of the beam-like form. We first derive the beam model starting from the two-dimensional linearized elasticity, the same way it is derived from the three-dimensional linearized elasticity. Then we present the numerical solution of the two-dimensional problem by finite element method. As expected the difference of two approximations becomes smaller as the thickness of the beam tends to zero. We then analyze the applicability of the one-dimensional model and verify the main properties of the beam modeling for thin beams.