• Steel and Composite Structures
• /
• v.35 no.1
• /
• pp.77-92
• /
• 2020

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

• Smart Structures and Systems
• /
• v.25 no.5
• /
• pp.605-617
• /
• 2020

The existence of damages in structures causes changes in the physical properties by reducing the modal parameters. In this paper, we develop a two-stages approach based on normalized Modal Strain Energy Damage Indicator (nMSEDI) for quick applications to predict the location of damage. A two-dimensional IsoGeometric Analysis (2D-IGA), Machine Learning Algorithm (MLA) and optimization techniques are combined to create a new tool. In the first stage, we introduce a modified damage identification technique based on frequencies using nMSEDI to locate the potential of damaged elements. In the second stage, after eliminating the healthy elements, the damage index values from nMSEDI are considered as input in the damage quantification algorithm. The hybrid of Teaching-Learning-Based Optimization (TLBO) with Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) are used along with nMSEDI. The objective of TLBO is to estimate the parameters of PSO-ANN to find a good training based on actual damage and estimated damage. The IGA model is updated using experimental results based on stiffness and mass matrix using the difference between calculated and measured frequencies as objective function. The feasibility and efficiency of nMSEDI-PSO-ANN after finding the best parameters by TLBO are demonstrated through the comparison with nMSEDI-IGA for different scenarios. The result of the analyses indicates that the proposed approach can be used to determine correctly the severity of damage in beam structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.159-168
• /
• 2005

An improved stability design method for beam-columns of plane frames is proposed based on system buckling analysis and second-order elastic analysis. For this, the tangent stiffness matrix of beam-column elements is first derived using stability functions and a procedure for evaluating effective buckling lengths is reviewed using elastic system buckling analysis. And then the second-order analysis procedure is presented considering $P-\Delta$ effects and is compared with the closed-form solution through numerical examples. Design examples showing the validity of the proposed method we presented and their numerical results are compared with those obtained from the conventional stability design methods. Finally some useful conclusions are drawn.

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

This paper presents a reliable numerical procedure for nonlinear time-history analysis of space steel frames subjected to dynamic loads. Geometric nonlinearities of member (P-$\delta$) and frame (P-$\Delta$) are taken into account by the use of stability functions in framed stiffness matrix formulation. The gradual yielding along the member length and over the cross section is included by using a tangent modulus concept and a softening plastic hinge model based on the New-Orbison yield surface. A computer program utilizing the average acceleration method for the integration scheme is developed to numerically solve the equation of motion of framed structure formulated in an incremental form. The results of several numerical examples are compared with those derived from using beam element model of ABAQUS program to illustrate the accuracy and the computational efficiency of the proposed procedure.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.31-37
• /
• 2017

In this paper, computation results of reducible modeling, stress recovery and energy release rate were compared with the results of VABS, Virtual Crack Closure Technique. The result of stress recovery analysis for 1-D model including the stiffness matrix is compared with stress results of three-dimensional 3-D FEM. Energy release rate of composite beam with longitudinal cracks is calculated and compare verifications of numerical analysis results of 3-D FEM and VABS. The procedure of calculating energy release rate through dimensional reduction and stress recovery is intended to be efficient and be utilized in the life-cycle of high-altitude uav's wing, wind blades and tilt rotor blade.

• Proceedings of the KSR Conference
• /
• 2001.10a
• /
• pp.500-507
• /
• 2001

In this study, Interaction among each element was construed by the unit of new stiffness matrix to analyze the finite element about Japanese precast concrete slab track and improved slab track. Dynamic analysis which is assumed a static analysis and a trainload on the transverse and the longitudinal load of the train into a series periodic function was performed by using the common program. And then, the difference of the movement between an improved section and an existing structure type was realized. Longitudinal static analysis indicated that the stress of the improved section is smaller than that of the protrusion of the existing slab track. And static and dynamic analysis on transverse load showed a little decrease of the displacement on new slab track. But the dynamic analysis result showed that new track system was considerably decreased by 30% compared with the existing Japanese slab track.

• Advances in Computational Design
• /
• v.2 no.3
• /
• pp.123-141
• /
• 2017

Solving a system of linear or non-linear equations is required to analyze any kind of structures. There are many ways to solve a system of equations, and they can be classified as implicit and explicit techniques. The explicit methods eliminate round-off errors and use less memory. The dynamic relaxation method (DR) is one of the powerful and simple explicit processes. The important point is that the DR does not require to store the global stiffness matrix, for which it just uses the residual loads vector. In this paper, a new approach to the DR method is expressed. In this approach, the damping, mass and time steps are similar to those of the traditional method of dynamic relaxation. The difference of this proposed method is focused on the method of calculating the damping. The proposed method is expressed such that the time step is constant, damping is equal to zero except in steps with maximum energy and the concentrated damping can be applied to minimize the energy of system in this step. In this condition, the calculation of damping in all steps is not required. Then the volume of computation is reduced. The DR method for form-finding of membrane structures is employed in this paper. The form-finding of the three plans related to the membrane structures with different loading is considered to investigate the efficiency of the proposed method. The numerical results show that the convergence rate based on the proposed method increases in all cases than other methods.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.121-131
• /
• 1986

The stiffness matrix of circularly curved frame elements including the warping effects is formulated by the solutions of vlasov's differential equations, and the procedure for the elastic static analysis of curved girder systems by the displacement method is presented. The validity of this method has been demonstrated by comparing the analysis results with other solutions. And if the tangential lines of the two frame element axes connected at any nodal point coincide, the transformation to the global coordinate system can be omitted when we analyze the structures consisting of circularly curved elements. The theory introduced in this thesis can be applied with sufficient accuracy to the structures built up with horizontally circular curved frame elements which have closed or open cross sections and are symmetric to the axis perpendicular to the plane of the curvature, such as prestressed concrete box girder bridges.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.9-21
• /
• 1992

The effectiveness of the isolation system installed at the base of the primary structure and at the support of the substructure mounted on the primary structure is evaluated for reducing of structural responses under different earthquakes in this paper. The structural responses are analyzed to identify its behavior due to the input motion characteristics such as various peak acceleration and frequency content. Three analytical models are used to evaluate the effectiveness of the isolation system in this study as follows: fixed-base primary structure with support-fixed substructure, base-isolated primary structure with support-fixed substructure, and fixed-base primary structure with support-isolated substruciure. A computer code (KBISAP) is used for numerical integration of equation of motion considering the interaction between the primary structure and the secondary structure. The matrix condensation technique and constant average acceleration method are utilized in this program. And also, the effective stiffness of the base-isolator on reducing the structural response are evaluated for various earthquakes through the relationship of the acceleration - displacement.