• Aerospace Engineering and Technology
• /
• v.13 no.2
• /
• pp.60-65
• /
• 2014

In this paper, numerical integration method using operational matrix of integration is studied. Using the operational matrix of integration, modified fixed point iteration method is introduced in order to solve rapidly an initial value problem for non-linear equation of motion. As an example, an initial value problem for orbital motion is considered. Through the numerical example, it is shown that the algorithm is efficient from the computational time point of view.

• Proceedings of the Korea Electromagnetic Engineering Society Conference
• /
• 2000.11a
• /
• pp.117-120
• /
• 2000

When analyzing the scatting problem of multi-layered structures using closed-form Green's function, one of the main difficulties is that the numerical integrations for the evaluation of diagonal matrix elements converge slowly and are not so stable. Accordingly, even when the integration for the singularity of type e$\^$-jkr//${\gamma}$/, corresponding to the source dipole itself, is performed using such a mathod, this difficulty persists in the integration corresponding to the finite number of complex images. In order to resolve this difficulty, a new technique based upon the Gaussian quadrature in polar coordinates for the evaluation of the two-dimensional generalized exponential integral is presented. Stability of the algorithm and convergence is discussed. Performance is demonstrated for the example of a microstrip patch antenna.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.139-148
• /
• 1997

A finite element method is programmed to analyse the nonlinear behavior of axisymmetric structures. The lst order Mindlin shell theory which takes into account the transversal shear deformation is used to formulate a conical two node element with six degrees of freedom. To evade the shear locking phenomenon which arises in Mindlin type element when the effect of shear deformation tends to zero, the reduced integration of one point Gauss Quadrature at the center of element is employed. This method is the Updated Lagrangian formulation which refers the variables to the state of the most recent iteration. The solution is searched by Newton-Raphson iteration method. The tangent matrix of this method is obtained by a finite difference method by perturbating the degrees of freedom with small values. For the moment this program is limited to the analyses of non-linear elastic problems. For structures which could have elastic stability problem, the calculation is controled by displacement.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.38 no.12
• /
• pp.56-65
• /
• 2001

This paper proposes an algorithm for the capacitance extraction of general 3-dimensional conductors in an ideal uniform dielectric that uses a high-order quadrature approximation method combined with the typical first-order collocation method to enhance the accuracy and adopts an efficient matrix-vector product algorithm for the model-order reduction to achieve efficiency. The proposed method enhances the accuracy using the quadrature method for interconnects containing corners and vias that concentrate the charge density. It also achieves the efficiency by reducing the model order using the fact that large parts of system matrices are of numerically low rank. This technique combines an SVD-based algorithm for the compression of rank-deficient matrices and Gram-Schmidt algorithm of a Krylov-subspace iterative technique for the rapid multiplication of matrices. It is shown through the performance evaluation procedure that the combination of these two techniques leads to a more efficient algorithm than Gaussian elimination or other standard iterative schemes within a given error tolerance.

• Structural Engineering and Mechanics
• /
• v.27 no.6
• /
• pp.713-739
• /
• 2007

Nowadays computers can perform symbolic computations in addition to mere number crunching operations for which they were originally designed. Symbolic computation opens up exciting possibilities in Structural Mechanics and engineering. Classical areas have been increasingly neglected due to the advent of computers as well as general purpose finite element software. But now, classical analysis has reemerged as an attractive computer option due to the capabilities of symbolic computation. The repetitive cycles of simultaneous - equation sets required by the finite element technique can be eliminated by solving a single set in symbolic form, thus generating a truly closed-form solution. This consequently saves in data preparation, storage and execution time. The power of Symbolic computation is demonstrated by six examples by applying symbolic computation 1) to solve coupled shear wall 2) to generate beam element matrices 3) to find the natural frequency of a shear frame using transfer matrix method 4) to find the stresses of a plate subjected to in-plane loading using Levy's approach 5) to draw the influence surface for deflection of an isotropic plate simply supported on all sides 6) to get dynamic equilibrium equations from Lagrange equation. This paper also presents yet another computationally efficient and accurate numerical method which is based on the concept of derivative of a function expressed as a weighted linear sum of the function values at all the mesh points. Again this method is applied to solve the problems of 1) coupled shear wall 2) lateral buckling of thin-walled beams due to moment gradient 3) buckling of a column and 4) static and buckling analysis of circular plates of uniform or non-uniform thickness. The numerical results obtained are compared with those available in existing literature in order to verify their accuracy.

• Steel and Composite Structures
• /
• v.37 no.6
• /
• pp.711-731
• /
• 2020

This paper deals with free vibration of FG sandwich annular sector plates on Pasternak elastic foundation with different boundary conditions, based on the three-dimensional theory of elasticity. The plates with simply supported radial edges and arbitrary boundary conditions on their circular edges are considered. The influence of carbon nanotubes (CNTs) waviness, aspect ratio, internal pores and graphene platelets (GPLs) on the vibrational behavior of functionally graded nanocomposite sandwich plates is investigated in this research work. The distributions of CNTs are considered functionally graded (FG) or uniform along the thickness of upper and bottom layers of the sandwich sectorial plates and their mechanical properties are estimated by an extended rule of mixture. In this study, the classical theory concerning the mechanical efficiency of a matrix embedding finite length fibers has been modified by introducing the tube-to-tube random contact, which explicitly accounts for the progressive reduction of the tubes' effective aspect ratio as the filler content increases. The core of structure is porous and the internal pores and graphene platelets (GPLs) are distributed in the matrix of core either uniformly or non-uniformly according to three different patterns. The elastic properties of the nanocomposite are obtained by employing Halpin-Tsai micromechanics model. A semi-analytic approach composed of 2D-Generalized Differential Quadrature Method (2D-GDQM) and series solution is adopted to solve the equations of motion. The fast rate of convergence and accuracy of the method are investigated through the different solved examples. Some new results for the natural frequencies of the plate are prepared, which include the effects of elastic coefficients of foundation, boundary conditions, material and geometrical parameters. The new results can be used as benchmark solutions for future researches.

• Steel and Composite Structures
• /
• v.48 no.3
• /
• pp.275-291
• /
• 2023

This paper has focused on presenting vibration analysis of trapezoidal sandwich plates with 3D-graphene foam reinforced polymer matrix composites (GrF-PMC) core and FG wavy CNT-reinforced face sheets. The porous graphene foam possessing 3D scaffold structures has been introduced into polymers for enhancing the overall stiffness of the composite structure. Also, 3D graphene foams can distribute uniformly or non-uniformly in the plate thickness direction. The effective Young's modulus, mass density and Poisson's ratio are predicted by the rule of mixture. In this study, the classical theory concerning the mechanical efficiency of a matrix embedding finite length fibers has been modified by introducing the tube-to-tube random contact, which explicitly accounts for the progressive reduction of the tubes' effective aspect ratio as the filler content increases. The First-order shear deformation theory of plate is utilized to establish governing partial differential equations and boundary conditions for trapezoidal plate. The governing equations together with related boundary conditions are discretized using a mapping-generalized differential quadrature (GDQ) method in spatial domain. Then natural frequencies of the trapezoidal sandwich plates are obtained using GDQ method. Validity of the current study is evaluated by comparing its numerical results with those available in the literature. It is explicated that 3D-GrF skeleton type and weight fraction, carbon nanotubes (CNTs) waviness and CNT aspect ratio can significantly affect the vibrational behavior of the sandwich structure. The plate's normalized natural frequency decreased and the straight carbon nanotube (w=0) reached the highest frequency by increasing the values of the waviness index (w).

• Steel and Composite Structures
• /
• v.34 no.2
• /
• pp.215-226
• /
• 2020

The aim of this study is to investigate free vibration of functionally graded porous nanocomposite rectangular plates where the internal pores and graphene platelets (GPLs) are distributed in the matrix either uniformly or non-uniformly according to three different patterns. The elastic properties of the nanocomposite are obtained by employing Halpin-Tsai micromechanics model. The GPL-reinforced plate is modeled using a semi-analytic approach composed of generalized differential quadrature method (GDQM) and series solution adopted to solve the equations of motion. The proposed rectangular plates have two opposite edges simply supported, while all possible combinations of free, simply supported and clamped boundary conditions are applied to the other two edges. The 2-D differential quadrature method as an efficient and accurate numerical tool is used to discretize the governing equations and to implement the boundary conditions. The convergence of the method is demonstrated and to validate the results, comparisons are made between the present results and those reported by well-known references for special cases treated before, have confirmed accuracy and efficiency of the present approach. New results reveal the importance of porosity coefficient, porosity distribution, graphene platelets (GPLs) distribution, geometrical and boundary conditions on vibration behavior of porous nanocomposite plates. It is observed that the maximum vibration frequency obtained in the case of symmetric porosity and GPL distribution, while the minimum vibration frequency is obtained using uniform porosity distribution.

• Steel and Composite Structures
• /
• v.44 no.1
• /
• pp.65-79
• /
• 2022

The main goal of this paper is to study the vibration of damaged core laminated annular plates with FG face sheets based on a three-dimensional theory of elasticity. The structures are made of a damaged isotropic core and two external face sheets. These skins are strengthened at the nanoscale level by randomly oriented Carbon nanotubes (CNTs) and are reinforced at the microscale stage by oriented straight fibers. These reinforcing phases are included in a polymer matrix and a three-phase approach based on the Eshelby-Mori-Tanaka scheme and on the Halpin-Tsai approach, which is developed to compute the overall mechanical properties of the composite material. In this study the effect of microcracks on the vibrational characteristic of the sandwich plate is considered. In particular, the structures are made by an isotropic core that undergoes a progressive uniform damage, which is modeled as a decay of the mechanical properties expressed in terms of engineering constants. These defects are uniformly distributed and affect the central layer of the plates independently from the direction, this phenomenon is known as "isotropic damage" and it is fully described by a scalar parameter. Three complicated equations of motion for the sectorial plates under consideration are semi-analytically solved by using 2-D differential quadrature method. Using the 2-D differential quadrature method in the r- and z-directions, allows one to deal with sandwich annular plate with arbitrary thickness distribution of material properties and also to implement the effects of different boundary conditions of the structure efficiently and in an exact manner. The fast rate of convergence and accuracy of the method are investigated through the different solved examples. The sandwich annular plate is assumed to have any arbitrary boundary conditions at the circular edges including simply supported, clamped and, free. Several parametric analyses are carried out to investigate the mechanical behavior of these multi-layered structures depending on the damage features, through-the-thickness distribution, and boundary conditions.

• Steel and Composite Structures
• /
• v.48 no.1
• /
• pp.1-17
• /
• 2023

The main goal of this paper is to study vibration of damaged core laminated sectorial plates with Functionally graded (FG) face sheets based on three-dimensional theory of elasticity. The structures are made of a damaged isotropic core and two external face sheets. These skins are strengthened at the nanoscale level by randomly oriented Carbon nanotubes (CNTs) and are reinforced at the microscale stage by oriented straight fibers. These reinforcing phases are included in a polymer matrix and a three-phase approach based on the Eshelby-Mori-Tanaka scheme and on the Halpin-Tsai approach, which is developed to compute the overall mechanical properties of the composite material. Three complicated equations of motion for the sectorial plates under consideration are semi-analytically solved by using 2-D differential quadrature method. Using the 2-D differential quadrature method in the r- and z-directions, allows one to deal with sandwich annular sector plate with arbitrary thickness distribution of material properties and also to implement the effects of different boundary conditions of the structure efficiently and in an exact manner. The fast rate of convergence and accuracy of the method are investigated through the different solved examples. The sandwich annular sector plate is assumed to be simply supported in the radial edges while any arbitrary boundary conditions are applied to the other two circular edges including simply supported, clamped and free. Several parametric analyses are carried out to investigate the mechanical behavior of these multi-layered structures depending on the damage features, through-the-thickness distribution and boundary conditions.