• Steel and Composite Structures
• /
• 제20권1호
• /
• pp.205-225
• /
• 2016

In this study the finite element method is utilized to predict the deflection and vibration characteristics of rectangular plates made of saturated porous functionally graded materials (PFGM) within the framework of the third order shear deformation plate theory. Material properties of PFGM plate are supposed to vary continuously along the thickness direction according to the power-law form and the porous plate is assumed of the form where pores are saturated with fluid. Various edge conditions of the plate are analyzed. The governing equations of motion are derived through energy method, using calculus of variations while the finite element model is derived based on the constitutive equation of the porous material. According to the numerical results, it is revealed that the proposed modeling and finite element approach can provide accurate deflection and frequency results of the PFGM plates as compared to the previously published results in literature. The detailed mathematical derivations are presented and numerical investigations are performed while the emphasis is placed on investigating the effect of the several parameters such as porosity volume fraction, material distribution profile, mode number and boundary conditions on the natural frequencies and deflection of the PFGM plates in detail. It is explicitly shown that the deflection and vibration behaviour of porous FGM plates are significantly influenced by these effects. Numerical results are presented to serve as benchmarks for future analyses of FGM plates with porosity phases.

• Computers and Concrete
• /
• 제16권6호
• /
• pp.881-896
• /
• 2015

A simultaneous analytical, experimental and numerical analysis of crack initiation, propagation and breaking process of the Central Straight through Crack Brazilian Disk (CSCBD) specimens under diametrical compression is carried out. Brazilian disc tests are being accomplished to evaluate the fracturing process based on stress intensity factors (SIFs). The effects of crack inclination angle and crack length on the fracturing processes have been investigated. The same experimental specimens have been numerically modeled by a higher order indirect boundary element method (HDDM). These numerical results are compared with the existing experimental results proving the accuracy and validity of the proposed numerical method.

• Coupled systems mechanics
• /
• 제8권1호
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2001년도 추계학술발표대회 논문집
• /
• pp.191-194
• /
• 2001

A partially coupled thermo-piezoelectric-mechanical triangular finite element model of composite laminates with surface bonded piezoelectric actuators, subjected to externally applied mechanical load, temperature change load, electric field load is developed. The governing differential equations are obtained by applying the principle of free energy and variational techniques. A higher order zigzag theory displacement field is employed to accurately capture the transverse shear and normal effects in laminated composite plates of arbitrary thickness. Nonconforming shape functions by Specht are employed in the transverse displacement variables. Numerical examples demonstrate the accuracy and efficiency of the proposed triangular plate element.

• 대한기계학회논문집A
• /
• 제30권1호
• /
• pp.18-25
• /
• 2006

The purpose of this research work is to demonstrate a successful application of hybrid-mixed formulation and nodeless degrees of freedom in developing a very accurate in-plane curved beam element for free vibration analysis. To resolve the numerical difficulties due to the spurious constraints, the present element, based on the Hellinger-Reissner variational principle and considering the effect of shear deformation, employed consistent stress parameters corresponding to cubic displacement polynomials with additional nodeless degrees. The stress parameters were eliminated by the stationary condition, and the nodeless degrees were condensed by Guyan Reduction. Several numerical examples indicated that the property of the mass matrix as well as that of the stiffness matrix have a great effect on the numerical performance. The element with consistent mass matrix produced best results on convergence and accuracy in the numerical analysis of Eigenvalue problems. Also, the higher-order mixed curved beam element showed a superior numerical behavior for the free vibration analyses.

• Structural Engineering and Mechanics
• /
• 제1권1호
• /
• pp.87-102
• /
• 1993

A $C^{\circ}$ continuous finite element formulation of a higher order displacement theory is presented for predicting linear and geometrically non-linear in the sense of von Karman transient responses of composite and sandwich plates. The displacement model accounts for non-linear cubic variation of tangential displacement components through the thickness of the laminate and the theory requires no shear correction coefficients. In the time domain, the explicit central difference integrator is used in conjunction with the special mass matrix diagonalization scheme which conserves the total mass of the element and included effects due to rotary inertia terms. The parametric effects of the time step, finite element mesh, lamination scheme and orthotropy on the linear and geometrically non-linear responses are investigated. Numerical results for central transverse deflection, stresses and stress resultants are presented for square/rectangular composite and sandwich plates under various boundary conditions and loadings and these are compared with the results from other sources. Some new results are also tabulated for future reference.

• 대한기계학회논문집A
• /
• 제21권12호
• /
• pp.1973-1983
• /
• 1997

Higher order isoparametric elements are usually used in the finite element analysis because they can represent easily the geometric shape of a complex structure ad can improve the solution quality. When these elements are used, the position of internal nodes affects greatly on the solution accuracy. Decreasing of the accuracy related to the position of internal nodes is due to non-conformal mapping often results in an unstable Jacobian value. This paper, in order to remove this difficulty, suggests a modified shape function which can establish conformal mapping between two coordinate systems. Numerical experiments with the proposed shape function show that a stable solution can be obtained without respect to the position of internal nodes, and offer convenience that one can take arbitrarily the position of internal nodes considering only the geometric shape of element boundaries.

• 대한기계학회논문집B
• /
• 제26권9호
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Structural Engineering and Mechanics
• /
• 제72권2호
• /
• pp.203-216
• /
• 2019

Reliability analysis of composite structures considering random variation of involved parameters is quite important as composite materials revealed large statistical variations in their mechanical properties. The reliability analysis of such structures by the first order reliability method (FORM) and Monte Carlo Simulation (MCS) based approach involves repetitive evaluations of performance function. The response surface method (RSM) based metamodeling technique has emerged as an effective solution to such problems. In the application of metamodeling for uncertainty quantification and reliability analysis of composite structures; the finite element model is usually formulated by either classical laminate theory or first order shear deformation theory. But such theories show significant error in calculating the structural responses of composite structures. The present study attempted to apply the RSM based MCS for reliability analysis of composite shell structures where the surrogate model is constructed using higher order shear deformation theory (HSDT) of composite structures considering the uncertainties in the material properties, load, ply thickness and radius of curvature of the shell structure. The sensitivity of responses of the shell is also obtained by RSM and finite element method based direct approach to elucidate the advantages of RSM for response sensitivity analysis. The reliability results obtained by the proposed RSM based MCS and FORM are compared with the accurate reliability analysis results obtained by the direct MCS by considering two numerical examples.

• Structural Engineering and Mechanics
• /
• 제58권6호
• /
• pp.949-966
• /
• 2016

The second order analysis taking place due to non-linear behavior of the structures under the mechanical and geometric factors through implementing exact and approximate methods is an indispensible issue in the analysis of such structures. Among the exact methods is the slope-deflection method that due to its simplicity and efficiency of its relationships has always been in consideration. By solving the differential equations of the modified slope-deflection method in which the effect of axial compressive force is considered, the stiffness matrix including trigonometric entries would be obtained. The complexity of computations with trigonometric functions causes replacement with their Maclaurin expansion. In most cases only the first two terms of this expansion are used but to obtain more accurate results, more elements are needed. In this paper, the effect of utilizing higher order terms of Maclaurin expansion on reducing the number of required elements and attaining more rapid convergence with less error is investigated for the Bernoulli beam with various boundary conditions. The results indicate that when using only one element along the beam length, utilizing higher order terms in Maclaurin expansion would reduce the relative error in determining the critical buckling load and kinematic parameters in the second order analysis.