• Steel and Composite Structures
• /
• v.24 no.1
• /
• pp.65-77
• /
• 2017

The purpose of this research is to study the nonlinear free vibration and post-buckling analysis of functionally graded carbon nanotube reinforced composite (FG-CNTRC) beams resting on a nonlinear elastic foundation. Uniformly and functionally graded distributions of single walled carbon nanotubes as reinforcing phase are considered in the polymeric matrix. The modified form of rule of mixture is used to estimate the material properties of CNTRC beams. The governing equations are derived employing Euler-Bernoulli beam theory along with energy method and Hamilton's principle. Applying von $K\acute{a}rm\acute{a}n's$ strain-displacement assumptions, the geometric nonlinearity is taken into consideration. The developed governing equations with quadratic and cubic nonlinearities are solved using variational iteration method (VIM) and the analytical expressions and numerical results are obtained for vibration and stability analysis of nanocomposite beams. The presented comparative results are indicative for the reliability, accuracy and fast convergence rate of the solution. Eventually, the effects of different parameters, such as foundation stiffness, volume fraction and distributions of carbon nanotubes, slenderness ratio, vibration amplitude, coefficients of elastic foundation and boundary conditions on the nonlinear frequencies, vibration response and post-buckling loads of FG-CNTRC beams are examined. The developed analytical solution provides direct insight into parametric studies of particular parameters of the problem.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.485-498
• /
• 2020

A diffusion synthetic acceleration (DSA) technique for the S_N transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

• Composites Research
• /
• v.25 no.6
• /
• pp.217-223
• /
• 2012

This study investigates a geometrical nonlinear dynamic behaviors of laminated skew plates made of advanced composite materials (ACM). Based on the first-order shear deformation plate theory (FSDT), the Newmark method and Newton-Raphson iteration are used for the nonlinear dynamic solution. The effects of cutout sizes, skew angles and lay up sequences on the nonlinear dynamic response for various parameters are studied using a nonlinear dynamic finite element program developed for this study. The several numerical results were in good agreement with those reported by other investigators for square composite plates with or without central cutouts, and the new results reported in this paper show the significant interactions between the cutout, skew angles and layup sequence in the laminate. Key observation points are discussed and a brief design guideline of skew laminates is given.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• Computational Structural Engineering : An International Journal
• /
• v.1 no.1
• /
• pp.31-38
• /
• 2001

In the case of the non-proportionally damped system such as the soil-structure interaction system, the structural control system and composite structures, the eigenproblem with the damping matrix should be necessarily performed to obtain the exact dynamic response. However, most of the eigenvalue analysis methods such as the subspace iteration method and the Lanczos method may miss some eigenvalues in the required ones. Therefore, the eigenvalue analysis method must include a technique to check the missed eigenvalues to become the practical tools. In the case of the undamped or proportionally damped system the missed eigenvalues can easily be checked by using the well-known Sturm sequence property, while in the case of the non-proportionally damped system a checking technique has not been developed yet. In this paper, a technique of checking the missed eigenvalues for the eigenproblem with the damping matrix is proposed by applying the argument principle. To verify the effectiveness of the proposed method, two numerical examples are considered.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.385-393
• /
• 2005

This paper established the dynamic model of a flexible Timoshenko beam capable of geometrical nonlinearities subject to large overall motions by using the finite element method. Equations of motion are derived by using Hamilton principle and are formulated in terms of finite elements in which the nonlinear constraint equations are adjoined to the system using Lagrange multipliers. The Newmark direct integration method and the Newton-Raphson iteration are employed here for the numerical study which is to demonstrate the efficiency of the proposed formulation.

• Journal of Korean Society of Steel Construction
• /
• v.8 no.3 s.28
• /
• pp.21-34
• /
• 1996

A numerical method is presented for computation of eigenvector derivatives used an iterative procedure with guaranteed convergence. An approach for treating the singularity in calculating the eigenvector derivatives is presented, in which a shift in each eigenvalue is introduced to avoid the singularity. If the shift is selected properly, the proposed method can give very satisfactory results after only one iteration. A criterion for choosing an adequate shift, dependent on computer hardware is suggested ; it is directly dependent on the eigenvalue magnitudes and the number of bits per numeral of the computer. Another merit of this method is that eigenvector derivatives with repeated eigenvalues can be easily obtained if the new eigenvectors are calculated. These new eigenvectors lie "adjacent" to the m (number of repeated eigenvalues) distinct eigenvectors, which appear when the design parameter varies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The results are compared with those of Nelson's method which can find the exact eigenvector derivatives. For the case of repeated eigenvalues, a cantilever beam is considered. The results are compared with those of Dailey's method which also can find the exact eigenvector derivatives. The design parameter of the cantilever plate is its thickness, and that of the cantilever beam its height.

• Structural Engineering and Mechanics
• /
• v.46 no.5
• /
• pp.735-753
• /
• 2013

The wind load is always the dominant load of cooling tower due to its large size, complex geometry and thin-wall structure. At present, when computing the wind-induced response of the large-scale cooling tower, the wind pressure distribution is obtained based on code regulations, wind tunnel test or computational fluid dynamic (CFD) analysis, and then is imposed on the tower structure. However, such method fails to consider the change of the wind load with the deformation of cooling tower, which may result in error of the wind load. In this paper, the analysis of the large cooling tower based on the iterative method for wind pressure is studied, in which the advantages of CFD and finite element method (FEM) are combined in order to improve the accuracy. The comparative study of the results obtained from the code regulations and iterative method is conducted. The results show that with the increase of the mean wind speed, the difference between the methods becomes bigger. On the other hand, based on the design of experiment (DOE), an approximate model is built for the optimal design of the large-scale cooling tower by a two-level optimization strategy, which makes use of code-based design method and the proposed iterative method. The results of the numerical example demonstrate the feasibility and efficiency of the proposed method.

• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.165-176
• /
• 2008

In this paper, a new iterative method for solving vehicle-bridge interaction problems is proposed. Iterative methods have advantages over the non-iterative methods in that it is not necessary to update the system matrix for a given wheel location, and the method can be applied for a new type of car or bridge with few or no modifications. In the proposed method, the necessity of system matrices update is eliminated using the equivalent interaction force acting on the bridge, which is obtained iteratively. Ballast stiffness is included in the interaction forces and the geometric compatibility at the contact points are used as convergence criteria. The bridge is considered as an elastic Bernoulli-Euler beam with surface irregularity and ballast stiffness. The moving vehicle is modeled as a multi-axle mass-spring-damper system having many degrees of freedom depending on the number of axles. The pitching effect, which is the interaction effect between the rear and front wheels when a vehicle begins to enter or leave the bridge, is also considered in the formulation including extended ground boundaries having surface irregularity and ballast stiffness. The applicability of the proposed method is illustrated in the numerical studies.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.