• International Journal of Highway Engineering
• /
• v.14 no.5
• /
• pp.1-9
• /
• 2012

PURPOSES: A viscoelastic axisymmetric finite element analysis code has been developed for stress analysis of asphalt pavement structures. METHODS: Generalized Maxwell Model (GMM) and 4-node isoparametric element were employed for finite element formulation. The code was developed using $C^{+}^{+}$ computer program language and named as KICTPAVE. For the verification of the developed code, a structural model of a pavement system was constructed. The structural model was composed of three layers: asphalt layer, crushed stone layer, and soil subgrade. Two types of analysis were considered for the verification: (1)elastic static analysis, (2)viscoelastic time-dependent analysis. For the elastic static analysis, linear elastic material model was assigned to all the layers, and a static load was applied to the structural model. For the viscoelastic time-dependent analysis, GMM and linear elastic material model were assigned to the asphalt layer and all the other layers respectively, and a cyclic loading condition was applied to the structural model. RESULTS: The stresses and deformations from KICTPAVE were compared with those from ABAQUS. The analysis results obtained from the two codes showed good agreement in time-dependent response of the element under the loading area as well as the surface deformation of asphalt layer, and horizontal and vertical stresses along the axisymmetric axis. CONCLUSIONS: The validity of KICTPAVE was confirmed by showing the agreement of the analysis results from the two codes.

• Structural Engineering and Mechanics
• /
• v.88 no.1
• /
• pp.83-94
• /
• 2023

An accurate and highly efficient inverse element labelled iPCB is developed based on the inverse finite element method (iFEM) for real-time shape estimation of plane-curved structures (such as arch bridges) utilizing onboard strain data. This inverse problem, named shape sensing, is vital for the design of smart structures and structural health monitoring (SHM) procedures. The iPCB formulation is defined based on a least-squares variational principle that employs curved Timoshenko beam theory as its baseline. The accurate strain-displacement relationship considering tension-bending coupling is used to establish theoretical and measured section strains. The displacement fields of the isoparametric element iPCB are interpolated utilizing nonuniform rational B-spline (NURBS) basis functions, enabling exact geometric modelling even with a very coarse mesh density. The present formulation is completely free from membrane and shear locking. Numerical validation examples for different curved structures subjected to different loading conditions have been performed and have demonstrated the excellent prediction capability of iPCBs. The present formulation has also been shown to be practical and robust since relatively accurate predictions can be obtained even omitting the shear deformation contributions and considering polluted strain measures. The current element offers a promising tool for real-time shape estimation of plane-curved structures.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.151-159
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory) concept such as the effect of compression softening in cracked concrete, and macroscopic and rotating crack models etc. was presented for the nonlinear behaviour of structural concrete. In this category, tangential or secant material stiffnesses for cracked concrete were also defined and discussed in view of the iterative solution schemes for nonlinear equations. Considering the computational efficiency and the ability of modelling the post-ultimate behaviour as major concerns, the incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Nonlinear Analysis of structrual Concrete by FEM : Monotonic Loading) developed baed on the CFT constitutive relationships and the incremetal solution strategy described enables the predictions of strength and deformation capacities in a full range. crack patterns and their corresponding widths, and yield extents of reinforcement. As the verfication purpose of NASCOM, the prediction of Cervenka's panel test results including the load resistance and the deformation history was made. A limited number of predictions indicate a good correlation in a general sense.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.455-483
• /
• 2018

The study deals with physical modeling of a typical three storeyed building frame supported by a pile group of four piles ($2{\times}2$) embedded in cohesive soil mass using three dimensional finite element analysis. For the purpose of modeling, the elements such as beams, slabs and columns, of the superstructure frame; and that of the pile foundation such as pile and pile cap are descretized using twenty noded isoparametric continuum elements. The interface between the pile and the soil is idealized using sixteen node isoparametric surface element. The soil elements are modeled using eight nodes, nine nodes and twelve node continuum elements. The present study considers the linear elastic behaviour of the elements of superstructure and substructure (i.e., foundation). The soil is assumed to behave non-linear. The parametric study is carried out for studying the effect of soil- structure interaction on response of the frame on the premise of sub-structure approach. The frame is analyzed initially without considering the effect of the foundation (non-interaction analysis) and then, the pile foundation is evaluated independently to obtain the equivalent stiffness; and these values are used in the interaction analysis. The spacing between the piles in a group is varied to evaluate its effect on the interactive behaviour of frame in the context of two embedment depth ratios. The response of the frame included the horizontal displacement at the level of each storey, shear force in beams, axial force in columns along with the bending moments in beams and columns. The effect of the soil- structure interaction is observed to be significant for the configuration of the pile groups and in the context of non-linear behaviour of soil.

• Coupled systems mechanics
• /
• v.3 no.3
• /
• pp.305-318
• /
• 2014

The study deals with the physical modeling of a typical building frame resting on pile foundation and embedded in cohesive soil mass using complete three-dimensional finite element analysis. Two different pile groups comprising four piles ($2{\times}2$) and nine piles ($3{\times}3$) are considered. Further, three different pile diameters along with the various pile spacings are considered. The elements of the superstructure frame and those of the pile foundation are descretized using twenty-node isoparametric continuum elements. The interface between the pile and pile and soil is idealized using sixteen-node isoparametric surface elements. The current study is an improved version of finite element modeling for the soil elements compared to the one reported in the literature (Chore and Ingle 2008). The soil elements are discretized using eight-, nine- and twelve-node continuum elements. Both the elements of superstructure and substructure (i.e., foundation) including soil are assumed to remain in the elastic state at all the time. The interaction analysis is carried out using sub-structure approach in the parametric study. The total stress analysis is carried out considering the immediate behaviour of the soil. The effect of various parameters of the pile foundation such as spacing in a group and number piles in a group, along with pile diameter, is evaluated on the response of superstructure. The response includes the displacement at the top of the frame and bending moment in columns. The soil-structure interaction effect is found to increase displacement in the range of 58 -152% and increase the absolute maximum positive and negative moments in the column in the range of 14-15% and 26-28%, respectively. The effect of the soil- structure interaction is observed to be significant for the configuration of the pile groups and the soil considered in the present study.

• Computational Structural Engineering : An International Journal
• /
• v.1 no.2
• /
• pp.107-116
• /
• 2001

The present paper investigates the nonlinear behavior of reinforced concrete shear walls with a crank based on a finite element analysis. The loading type is a horizontal cyclic one such as earthquake loads. Experiments of the shear walls with and without cranks, performed previously to see flow the behavior changes depending on the crank, are compared with the results obtained from the finite element analysis. The finite element analysis is based on an isoparametric degenerated shell formulation. The nonlinear constitutive equations fur concrete are modeled adopting the formulation based on a concept of Ring Typed-Lattice Model. The experiments indicate that the shear walls with a crank have low stiffness and relatively low carrying capacity compared with an ordinary plane shear wall without cranks and that they are more ductile, and the tendency is a1so confirmed based on the finite element analysis. Moreover, a good agreement between the experiments and analyses is obtained, accordingly, it is confined that the present numerical analysis scheme based on the Lattice Model is a powerful one to evaluate the behavior of reinforced concrete shear walls with cranks and without cranks.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.445-454
• /
• 1991

In this thesis, Mindlin plate element with nine nodes and three degrees-of-freedom at each node is formulated and is employed in eigen-analysis of a rectangular plates in order to alleviate locking phenomenon of eigenvalues. Eigenvalues and their modes may be locked if conventional $C_{0}$-isoparametric element is used. In order to reduce stiffness locking phenomenon, two methods (1, the general reduced and selective integration, 2, the new element that use of modified shape function) are studied. Additionally in order to reduce the error due to mass matrix, two mass matrixes (1, Gauss-Legendre mass matrix, 2, Gauss-Lobatto mass matrix) are considered. The results of eigen-analysis for two models (the square plate with all edges simply-supported and all edges built-in), computed by two methods for stiffness matrix and by two mass matrixes are compared with theoretical solutions and conventional numerical solutions. These comparisons show that the performance of the two methods with Gauss-Lobatto mass matrix is better than that of the conventional plate element. But, by considering the spurious rigid body motions, the element which employs modified shape function with full integration and Gauss-Lobatto mass matrix can elevate the accuracy and convergence of numerical solutions.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2067-2077
• /
• 1996

In adaptive finite element analysis, r- and h-methods are generally used on the basis of a discretization error estimator. In this paper, an rh-method is proposed as a new adaptive method which can improve the adaptivity performance by using both of them. This suggested rh-method moves nodal coordinates of initially given model to adjust element discretization errors and thereafter performes the h-method tdo obtain the specified accuracy of finite element solutions. Numerical experiments for various plane problems were performed using 4-noded isoparametric quadrilateral elements. As a result, the rh-method has been shown to be an accurate and efficient adaptive analysis method to obtain as improved solution.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.923-937
• /
• 2014

The rational finite element method is different from the standard finite element method, which is constructed using basic solutions of the governing differential equations as interpolation functions in the elements. Therefore, it is superior to the isoparametric approach because of its obvious physical meaning and accuracy; it has successfully been applied to the isotropic elasticity problem. In this paper, the formulation of rational finite elements for plane orthotropic elasticity problems is deduced. This method is formulated directly in the physical domain with full consideration of the requirements of the patch test. Based on the number of element nodes and the interpolation functions, different approaches are applied with complete polynomial interpolation functions. Then, two special stiffness matrixes of elements with four and five nodes are deduced as a representative application. In addition, some typical numerical examples are considered to evaluate the performance of the elements. The numerical results demonstrate that the present method has a high level of accuracy and is an effective technique for solving plane orthotropic elasticity problems.

• Structural Engineering and Mechanics
• /
• v.67 no.3
• /
• pp.301-310
• /
• 2018

In this paper, a multi-objective multiparameter optimization procedure is developed by combining rigorously developed metamodels with an evolutionary search algorithm-Genetic Algorithm (GA). Response surface methodology (RSM) is used for developing the metamodels to replace the tedious finite element analyses. A nine-node isoparametric plate bending element is used for conducting the finite element simulations. Highly accurate numerical data from an author compiled FORTRAN finite element program is first used by the RSM to develop second-order mathematical relations. Four material parameters-${\frac{E_1}{E_2}}$, ${\frac{G_{12}}{E_2}}$, ${\frac{G_{23}}{E_2}}$ and ${\upsilon}_{12}$ are considered as the independent variables while simultaneously maximizing fundamental frequency, ${\lambda}_1$ and frequency separation between the $1^{st}$ two natural modes, ${\lambda}_{21}$. The optimal material combination for maximizing ${\lambda}_1$ and ${\lambda}_{21}$ is predicted by using a multi-objective GA. A general sensitivity analysis is conducted to understand the effect of each parameter on the desired response parameters.