• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.4
• /
• pp.335-341
• /
• 2017

In this paper, beam finite elements with higher-order derivatives' continuity are formulated and evaluated for various boundary conditions. All the beam elements are based on Euler-Bernoulli beam theory. These higher-order beam elements are often required to analyze structures by using newly developed higher-order beam theories and/or non-classical beam theories based on nonlocal elasticity. It is however rare to assess the performance of such elements in terms of boundary and loading conditions. To this end, two higher-order beam elements are formulated, in which $C^2$ and $C^3$ continuities of the deflection are enforced, respectively. Three different boundary conditions are then applied to solve beam structures, such as cantilever, simply-support and clamped-hinge conditions. In addition to conventional Euler-Bernoulli beam boundary conditions, the effect of higher-order boundary conditions is investigated. Depending on the boundary conditions, the oscillatory behavior of deflections is observed. Especially the geometric boundary conditions are problematic, which trigger unstable solutions when higher-order deflections are prescribed. It is expected that the results obtained herein serve as a guideline for higher-order derivatives' continuous finite elements.

• Steel and Composite Structures
• /
• v.23 no.1
• /
• pp.41-51
• /
• 2017

A $C^0$ FE model developed based on an efficient higher order zigzag theory is used for hygrothermal analysis of laminated composite plates. The $C^0$ FE model satisfies the inter-laminar shear stress continuity at the interfaces and zero transverse shear stress conditions at plate top and bottom. In this model the first derivatives of transverse displacement have been treated as independent variables to circumvent the problem of $C^1$ continuity associated with the above plate theory. In the present theory the above mentioned $C^0$ continuity of the present element is compensated in the stiffness matrix formulation by using penalty parameter approach. In order to avoid stress oscillations observed in the displacement based finite element, the stress field derived from temperature/moisture fields (initial strains) must be consistent with total strain field. Special steps are introduced by field consistent approach (e.g., sampling at gauss points) to compensate this problem. A nine noded $C^0$ continuous isoparametric element is used in the proposed FE model. Comparison of present numerical results with other existing solutions shows that the proposed FE model is efficient, accurate and free of locking.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.849-856
• /
• 2006

This paper is to examine the applicability of ordinary Kriging interpolation(OK) to the p-adaptivity of the finite element analysis that is based on variogram. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. In case of non-uniform p-distribution, the continuity between elements with different polynomial orders is achieved by assigning zero higher-order derivatives associated with the edge in common with the lower-order derivatives. It is demonstrated that the validity of the proposed approach by analyzing results for stress singularity problem.

• Steel and Composite Structures
• /
• v.38 no.3
• /
• pp.281-291
• /
• 2021

An equivalent single-layer theory (EST) is put forward for analyzing free vibrations of steel-concrete composite beams (SCCB) based on a higher-order beam theory. In the EST, the effect of partial interaction between sub-beams and the transverse shear deformation are taken into account. After using the interlaminar shear force continuity condition and the shear stress free conditions at the top and bottom surface, the displacement function of the EST does not contain the first derivatives of transverse displacement. Therefore, the C0 interpolation functions are just demanded during its finite element implementation. Finally, the EST is validated by comparing the results of two simply-supported steel-concrete composite beams which are tested in laboratory and calculated by ANSYS software. Then, the influencing factors for free vibrations of SCCB are analyzed, such as, different boundary conditions, depth to span ratio, high-order shear terms, and interfacial shear connector stiffness.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.147-161
• /
• 2023

Based on the ideas of variational differential quadrature (VDQ) and finite element method (FEM), a numerical approach named as VDQFEM is applied herein to study the large deformations of plate-type structures under static loading with arbitrary shape hole made of functionally graded graphene platelet-reinforced composite (FG-GPLRC) in the context of higher-order shear deformation theory (HSDT). The material properties of composite are approximated based upon the modified Halpin-Tsai model and rule of mixture. Furthermore, various FG distribution patterns are considered along the thickness direction of plate for GPLs. Using novel vector/matrix relations, the governing equations are derived through a variational approach. The matricized formulation can be efficiently employed in the coding process of numerical methods. In VDQFEM, the space domain of structure is first transformed into a number of finite elements. Then, the VDQ discretization technique is implemented within each element. As the last step, the assemblage procedure is performed to derive the set of governing equations which is solved via the pseudo arc-length continuation algorithm. Also, since HSDT is used herein, the mixed formulation approach is proposed to accommodate the continuity of first-order derivatives on the common boundaries of elements. Rectangular and circular plates under various boundary conditions with circular/rectangular/elliptical cutout are selected to generate the numerical results. In the numerical examples, the effects of geometrical properties and reinforcement with GPL on the nonlinear maximum deflection-transverse load amplitude curve are studied.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.216-221
• /
• 2008

Shape design optimization for linear elasticity problem is performed using isogeometric analysis method. In many design optimization problems for real engineering models, initial raw data usually comes from CAD modeler. Then designer should convert this CAD data into finite element mesh data because conventional design optimization tools are generally based on finite element analysis. During this conversion there is some numerical error due to a geometry approximation, which causes accuracy problems in not only response analysis but also design sensitivity analysis. As a remedy of this phenomenon, the isogeometric analysis method is one of the promising approaches of shape design optimization. The main idea of isogeometric analysis is that the basis functions used in analysis is exactly same as ones which represent the geometry, and this geometrically exact model can be used shape sensitivity analysis and design optimization as well. In shape design sensitivity point of view, precise shape sensitivity is very essential for gradient-based optimization. In conventional finite element based optimization, higher order information such as normal vector and curvature term is inaccurate or even missing due to the use of linear interpolation functions. On the other hands, B-spline basis functions have sufficient continuity and their derivatives are smooth enough. Therefore normal vector and curvature terms can be exactly evaluated, which eventually yields precise optimal shapes. In this article, isogeometric analysis method is utilized for the shape design optimization. By virtue of B-spline basis function, an exact geometry can be handled without finite element meshes. Moreover, initial CAD data are used throughout the optimization process, including response analysis, shape sensitivity analysis, design parameterization and shape optimization, without subsequent communication with CAD description.