• Structural Engineering and Mechanics
• /
• v.31 no.5
• /
• pp.481-506
• /
• 2009

The present study is concerned with the stochastic linear free vibration study of laminated composite plate embedded with piezoelectric layers with random material properties. The system equations are derived using higher order shear deformation theory. The lamina material properties of the laminate are modeled as basic random variables for accurate prediction of the system behavior. A $C^0$ finite element is used for spatial descretization of the laminate. First order Taylor series based mean centered perturbation technique in conjunction with finite element method is outlined for the problem. The outlined probabilistic approach is used to obtain typical numerical results, i.e., the mean and standard deviation of natural frequency. Different combinations of simply supported, clamped and free boundary conditions are considered. The effect of side to thickness ratio, aspect ratio, lamination scheme on scattering of natural frequency is studied. The results are compared with those available in literature and an independent Monte Carlo simulation.

• International Journal of Safety
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2003

The present study deals with an approximate integral equation approach to finite deflection of elastic plates with arbitrary plane form. An integral formulation leads to a system of boundary integral equations involving values of deflection, slope, bending moment and transverse shear force along the edge. The basic principles of the development of boundary element technique are reviewed. A computer program for solving for stresses and deflections in a isotropic, homogeneous, linear and elastic bending plate is developed. The fundamental solution of deflection and moment is employed in this program. The deflections and moments are assumed constant within the quadrilateral element. Numerical solutions for sample problems, obtained by the direct boundary element method, are presented and results are compared with known solutions.

• Structural Engineering and Mechanics
• /
• v.35 no.3
• /
• pp.283-299
• /
• 2010

This paper addresses the numerical simulation of fatigue crack growth in arbitrary 2D geometries under constant amplitude loading by the using a new finite element software. The purpose of this software is on the determination of 2D crack paths and surfaces as well as on the evaluation of components Lifetimes as a part of the damage tolerant assessment. Throughout the simulation of fatigue crack propagation an automatic adaptive mesh is carried out in the vicinity of the crack front nodes and in the elements which represent the higher stresses distribution. The fatigue crack direction and the corresponding stress-intensity factors are estimated at each small crack increment by employing the displacement extrapolation technique under facilitation of singular crack tip elements. The propagation is modeled by successive linear extensions, which are determined by the stress intensity factors under linear elastic fracture mechanics (LEFM) assumption. The stress intensity factors range history must be recorded along the small crack increments. Upon completion of the stress intensity factors range history recording, fatigue crack propagation life of the examined specimen is predicted. A consistent transfer algorithm and a crack relaxation method are proposed and implemented for this purpose. Verification of the predicted fatigue life is validated with relevant experimental data and numerical results obtained by other researchers. The comparisons show that the program is capable of demonstrating the fatigue life prediction results as well as the fatigue crack path satisfactorily.

• Advances in nano research
• /
• v.12 no.5
• /
• pp.441-455
• /
• 2022

This study explores the linear and nonlinear solutions of sigmoid functionally graded material (S-FGM) nanoplate with porous effects. A size-dependent numerical solution is established using the strain gradient theory and isogeometric finite element formulation. The nonlinear nonlocal strain gradient is developed based on the Reissner-Mindlin plate theory and the Von-Karman strain assumption. The sigmoid function is utilized to modify the classical functionally graded material to ensure the constituent volume distribution. Two different patterns of porosity distribution are investigated, viz. pattern A and pattern B, in which the porosities are symmetric and asymmetric varied across the plate's thickness, respectively. The nonlinear finite element governing equations are established for bending analysis of S-FGM nanoplates, and the Newton-Raphson iteration technique is derived from the nonlinear responses. The isogeometric finite element method is the most suitable numerical method because it can satisfy a higher-order derivative requirement of the nonlocal strain gradient theory. Several numerical results are presented to investigate the influences of porosity distributions, power indexes, aspect ratios, nonlocal and strain gradient parameters on the porous S-FGM nanoplate's linear and nonlinear bending responses.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.387-394
• /
• 2004

In the structural analysis procedure using finite element technique, the performance of a linear equation solver is critical because the linear equation solving part spends most of the computing time for finite element analysis codes. However, most of researchers are still using inefficient profile-based direct solvers such as the band solver or the skyline solver. In this research, we introduce the multifrontal solution method as an efficient direct solution method for structural analysis, and show the efficiency and performance of the multifrontal solution method by comparing the performance of our own implementation of the multifrontal method with the band solver or the skyline solver. In addition, we also compare the performance of our solver with other implementations of the multifrontal method such as WSMP and MUMPS as well as commercial structural analysis packages such as ABAQUS and NASTRAN. Through the performance test results, the usefulness and efficiency of our domain-wise multifrontal solver for structural analysis is shown.

• Proceedings of the KIEE Conference
• /
• 1997.07a
• /
• pp.168-170
• /
• 1997

This paper presents the dynamic analysis of a reciprocating actuation system based on moving magnet actuator for gas compression. For the analysis of the linear actuator, an axisymmetric finite element method (FEM) considering the saturation effect of the magnetic material is used, and electrical circuit equation, mechanical dynamic equation and pressure dynamics are coupled. In the FE analysis, we adopt a moving line technique. The pressure dynamics of the gas in the compressor is modeled by using the law of thermodynamics. The analysis results are compared fairly well with experimental ones.

• Wind and Structures
• /
• v.18 no.3
• /
• pp.235-252
• /
• 2014

An effective numerical technique to calculate the reactions of a multi-spanned transmission line conductor system, under arbitrary loads varying along the spans, is developed. Such variable loads are generated by High Intensity Wind (HIW) events in the form of tornadoes and downburst. First, a semi-closed form solution is derived to obtain the displacements and the reactions at the ends of each conductor span. The solution accounts for the nonlinearity of the system and the flexibility of the insulators. Second, a numerical scheme to solve the derived closed-form solution is proposed. Two conductor systems are analyzed under loads resulting from HIW events for validation of the proposed technique. Non-linear Finite Element Analyses (FEA) are also conducted for the same two systems. The responses resulting from the technique are shown to be in a very good agreement with those resulting from the FEA, which confirms the technique accuracy. Meanwhile, the semi-closed form technique shows superior efficiency in terms of the required computational time. The saving in computational time has a great advantage in predicting the response of the conductors under HIW events, since this requires a large number of analyses to cover different potential locations and sizes of those localized events.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.49-54
• /
• 1997

Rice's monotone streamline upwind finite element method, which was proposed to treat convection-dominated flows, is applied to the linear triangular element. An alignment technique of unstructured grids with given velocity fields is used to prevent the interpolation error produced in evaluating the convection term in the upwind method. The alignment of grids is accomplished by optimizing a target function defined with the inner-product of a properly chosen side vector in the element with the velocity field. Two pure advection problems are considered to demonstrate the superiorities of the present approach in solving the convection-dominated flow on the unstructured grid. Solutions obtained with aligned grids are much closer to the exact solutions than those with initial regular grids. The capability of the present approach in predicting the appearance of the secondary vortex in the laminar confined jet impingement is shown by comparing streamlines to those produced by SIMPLE on a highly stretched grid toward the impingement plate.

• Structural Engineering and Mechanics
• /
• v.46 no.2
• /
• pp.213-229
• /
• 2013

Although performance based assessment procedures are mainly developed for reinforced concrete and steel buildings, URM (Unreinforced Masonry) buildings occupy significant portion of buildings in earthquake prone areas of the world as well as in IRAN. Variability of material properties, non-engineered nature of the construction and difficulties in structural analysis of masonry walls make analysis of URM buildings challenging. Despite sophisticated finite element models satisfy the modeling requirements, extensive experimental data for definition of material behavior and high computational resources are needed. Recently, nonlinear equivalent frame models which are developed assigning lumped plastic hinges to isotropic and homogenous equivalent frame elements are used for nonlinear modeling of URM buildings. The equivalent frame models are not novel for the analysis of masonry structures, but the actual potentialities have not yet been completely studied, particularly for non-linear applications. In the present paper an effective tool for the non-linear static analysis of 2D masonry walls is presented. The work presented in this study is about performance assessment of unreinforced brick masonry buildings through nonlinear equivalent frame modeling technique. Reliability of the proposed models is tested with a reversed cyclic experiment conducted on a full scale, two-story URM building at the University of Pavia. The pushover curves were found to provide good agreement with the experimental backbone curves. Furthermore, the results of analysis show that EFM (Equivalent Frame Model) with Dolce RO (rigid offset zone) and shell element have good agreement with finite element software and experimental results.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.369-376
• /
• 2000

For the general loading condition and boundary condition, it is very difficult to obtain closed-form solutions for buckling loads and natural frequencies of thin-walled structures because its behaviour is very complex due to the coupling effect of bending and torsional behaviour. Consequently most of previous finite element formulations introduced approximate displacement fields using shape functions as Hermitian polynomials, isoparametric interpoation function, and so on. The purpose of this study is to calculate the exact displacement field of a thin-walled straight beam element with the non-symmetric cross section and present a consistent derivation of the exact dynamic stiffness matrix. An exact dynamic element stiffness matrix is established from Vlasov's coupled differential equations for a uniform beam element of non-symmetric thin-walled cross section. This numerical technique is accomplished via a generalized linear eigenvalue problem by introducing 14 displacement parameters and a system of linear algebraic equations with complex matrices. The natural frequencies are evaluated for the non-symmetric thin-walled straight beam structure, and the results are compared with available solutions in order to verify validity and accuracy of the proposed procedures.