• Structural Engineering and Mechanics
• /
• v.22 no.5
• /
• pp.613-629
• /
• 2006

In this paper a recently developed scaled boundary finite element method (SBFEM) is applied to simulate stress concentration for two-dimensional structures. In addition, a simple and independent formulation for evaluating the coefficients, not only of the singular term but also higher order non-singular terms, of the stress fields near crack-tip is presented. The formulation is formed by comparing the displacement along the radial points ahead of the crack-tip with that of standard Williams' eigenfunction solution for the crack-tip. The validity of the formulation is examined by numerical examples with different geometries for a range of crack sizes. The results show good agreement with available solutions in literatures. Based on the results of the study, it is conformed that the proposed numerical method can be applied to simulate stress concentrations in both cracked and uncracked structure components more easily with relatively coarse and simple model than other computational methods.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.2
• /
• pp.199-209
• /
• 2012

An advanced aeroelastic formulation for flutter analyses is presented in this paper. Refined 1D structural models were coupled with the doublet lattice method, and the g-method was used for flutter analyses. Structural models were developed in the framework of the Carrera Unified Formulation (CUF). Higher-order 1D structural models were obtained by using Taylor-like expansions of the cross-section displacement field of the structure. The order (N) of the expansion was considered as a free parameter since it can be arbitrarily chosen as an input of the analysis. Convergence studies on the order of the structural model can be straightforwardly conducted in order to establish the proper 1D structural model for a given problem. Flutter analyses were conducted on several wing configurations and the results were compared to those from literature. Results show the enhanced capabilities of CUF 1D in dealing with the flutter analysis of typical wing structures with high accuracy and low computational costs.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Structural Engineering and Mechanics
• /
• v.69 no.4
• /
• pp.427-437
• /
• 2019

In this study, an alternative solution procedure presented by using variational methods for analysis of shear deformable functionally graded material (FGM) beams with mixed formulation. By using the advantages of $G{\hat{a}}teaux$ differential approaches, a refined complex general functional and boundary conditions which comprises seven independent variables such as displacement, rotation, bending moment and higher-order bending moment, shear force and higher-order shear force, is derived for general thick-thin FGM beams via shear deformation beam theories. The mixed-finite element method (FEM) is employed to obtain a beam element which have a 2-nodes and total fourteen degrees-of-freedoms. A computer program is written to execute the analyses for the present study. The numerical results of analyses obtained for different boundary conditions are presented and compared with results available in the literature.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.659-683
• /
• 2014

The static analysis of structures with arbitrary cross-section geometry and material lamination via a refined one-dimensional (1D) approach is presented in this paper. Higher-order 1D models with a variable order of expansion for the displacement field are developed on the basis of Carrera Unified Formulation (CUF). Classical Euler-Bernoulli and Timoshenko beam theories are obtained as particular cases of the first-order model. Numerical results of displacement, strain and stress are provided by using the finite element method (FEM) along the longitudinal direction for different configurations in excellent agreement with three-dimensional (3D) finite element solutions. In particular, a layered thin-walled cylinder is considered as first assessment with a laminated conventional cross-section. An atherosclerotic plaque is introduced as a typical structure with arbitrary cross-section geometry and studied for both the homogeneous and nonhomogeneous material cases through the 1D variable kinematic models. The analyses highlight limitations of classical beam theories and the importance of higher-order terms in accurately detecting in-plane cross-section deformation without introducing additional numerical problems. Comparisons with 3D finite element solutions prove that 1D CUF provides remarkable three-dimensional accuracy in the analysis of even short and nonhomogeneous structures with arbitrary geometry through a significant reduction in computational cost.

• Structural Engineering and Mechanics
• /
• v.89 no.2
• /
• pp.199-211
• /
• 2024

The loading capacity of engineering structures/components reduces after the initiation and propagation of crack eventually leads to the final failure. Hence, it becomes essential to deal with the crack and its effects at the design and simulation stages itself, by detecting the prone area of the fracture. The phase-field (PF) method has been accepted widely in simulating fracture problems in complex geometries. However, most of the PF methods are formulated with second order continuity theoryinvolving C⁰ continuity. In the present study, PF method based on fourth-order (i.e., higher order) theory, maintaining C¹ continuity has been proposed for ductile fracture simulation. The formulation includes fourth-order derivative terms of phase field variable, varying between 0 and 1. Applications of fourth-order PF theory to ductile fracture simulation resulted in novelty in this area. The proposed formulation is numerically solved using a two-dimensional finite element (FE) framework in 3-layered manner system. The solutions thus obtained from the proposed fourth order theory for different benchmark problems portray the improvement in the accuracy of the numerical results and are well matched with experimental results available in the literature. These results are also compared with second-order PF theory and a comparison study demonstrated the robustness of the proposed model in capturing ductile behaviour close to experimental observations.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.663-668
• /
• 2007

An enhanced first shear deformation theory for composite plate is developed. The detailed process is as follows. Firstly, the theory is formulated by modifying higher order zigzag theory. That is, the higher order theory is separated into the warping function representing the higher order terms and lower order terms. Secondly, the relationships between higher order zig-zag field and averaged first shear deformation field based on the Reissner-Mindlin's plate theory are derived. Lastly, the effective shear modulus is calculated by minimizing error between higher order energy and first order energy. Then the governing equation of FSDT is solved by substituting shear modulus into effective shear modulus. The recovery processing with the nodal unknown obtained from governing equation is performed. The accuracy of the present proposed theory is demonstrated through numerical examples. The proposed method will serve as a powerful tool in the prediction of laminated composite plate.

• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.297-320
• /
• 2005

In recent years the pure displacement formulation for plate elements has not been as popular as other formulations. We revisit the pure displacement formulation for shear-deformable plate elements and propose a family of N-node, displacement-compatible, fully-integrated, pure-displacement, triangular, Mindlin plate elements, MIN-N. The development has been motivated by the relative simplicity of the pure displacement formulation and by the success of the existing 3-node plate element, MIN3. The formulation of MIN3 is generalized to obtain the MIN-N family, which possesses complete, fully compatible kinematic fields, in which the interpolation functions for transverse displacement are one degree higher than those for rotations. General element-level formulas for the thin-limit Kirchhoff constraints are developed. The 6-node, 18 degree-of-freedom element MIN6, with cubic displacement and quadratic rotations, is implemented and tested extensively. Numerical results show that MIN6 exhibits good performance for both static and dynamic analyses in the linear, elastic regime. The results illustrate that the fully-integrated MIN6 element has excellent performance in the thin limit, even for coarse meshes, and that it does not require shear relaxation.

• Structural Engineering and Mechanics
• /
• v.60 no.1
• /
• pp.1-19
• /
• 2016

In this study, the supersonic panel flutter of doubly curved composite sandwich panels with variable thickness is considered under aerothermoelastic loading. Considering different radii of curvatures of the face sheets in this paper, the thickness of the core is a function of plane coordinates (x,y), which is unique. For the first time in the current model, the continuity conditions of the transverse shear stress, transverse normal stress and transverse normal stress gradient at the layer interfaces, as well as the conditions of zero transverse shear stresses on the upper and lower surfaces of the sandwich panel are satisfied. The formulation is based on an enhanced higher order sandwich panel theory and the vertical displacement component of the face sheets is assumed as a quadratic one, while a cubic pattern is used for the in-plane displacement components of the face sheets and the all displacement components of the core. The formulation is based on the von $K{\acute{a}}rm{\acute{a}}n$ nonlinear approximation, the one-dimensional Fourier equation of the heat conduction along the thickness direction, and the first-order piston theory. The equations of motion and boundary conditions are derived using the Hamilton principle and the results are validated by the latest results published in the literature.

• Steel and Composite Structures
• /
• v.43 no.1
• /
• pp.1-17
• /
• 2022

The free and forced nonlinear dynamic behaviors of Porous Functionally Graded Material (PFGM) plates are examined by means of a High-Order Implicit Algorithm (HOIA). The formulation is developed using the Third-order Shear Deformation Theory (TSDT). Unlike previous works, the formulation is written without resorting to any homogenization technique neither rule of mixture nor considering FGM as a laminated composite, and the distribution of the porosity is assumed to be gradually variable through the thickness of the PFGM plates. Using the Hamilton principle, we establish the governing equations of motion. The Finite Element Method (FEM) is used to compute approximations of the resulting equations; FEM is adopted using a four-node quadrilateral finite element with seven Degrees Of Freedom (DOF) per node. Nonlinear equations are solved by a HOIA. The accuracy and the performance of the proposed approach are verified by presenting comparisons with literature results for vibration natural frequencies and dynamic response of PFGM plates under external loading. The influences of porosity volume fraction, porosity distribution, slenderness ratio and other parameters on the vibrations of PFGM plate are explored. The results demonstrate the significant impact of different physical and geometrical parameters on the vibration behavior of the PFGM plate.