scholarly journals Explicit Dynamic Finite Element Method for Failure with Smooth Fracture Energy Dissipations

2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Jeong-Hoon Song ◽  
Thomas Menouillard ◽  
Alireza Tabarraei

A numerical method for dynamic failure analysis through the phantom node method is further developed. A distinct feature of this method is the use of the phantom nodes with a newly developed correction force scheme. Through this improved approach, fracture energy can be smoothly dissipated during dynamic failure processes without emanating noisy artifact stress waves. This method is implemented to the standard 4-node quadrilateral finite element; a single quadrature rule is employed with an hourglass control scheme in order to decrease computational cost and circumvent difficulties associated with the subdomain integration schemes for cracked elements. The effectiveness and robustness of this method are demonstrated with several numerical examples. In these examples, we showed the effectiveness of the described correction force scheme along with the applicability of this method to an interesting class of structural dynamic failure problems.

1996 ◽  
Vol 3 (4) ◽  
pp. 259-268 ◽  
Author(s):  
M.S. Yao

The large number of unknown variables in a finite element idealization for dynamic structural analysis is represented by a very small number of generalized variables, each associating with a generalized Ritz vector known as a basis vector. The large system of equations of motion is thereby reduced to a very small set by this transformation and computational cost of the analysis can be greatly reduced. In this article nonlinear equations of motion and their transformation are formulated in detail. A convenient way of selection of the generalized basis vector and its limitations are described. Some illustrative examples are given to demonstrate the speed and validity of the method. The method, within its limitations, may be applied to dynamic problems where the response is global in nature with finite amplitude.


2008 ◽  
Vol 400-402 ◽  
pp. 613-619
Author(s):  
Hui Xiong ◽  
Shou Ping Shang ◽  
Liang Huang

Combined with the respective advantages in S-R(Sway-Rocking) impedance concept and finite-element method, a simplified 3D structural dynamic FEM considering composite pile-group-soil effects is presented. The structural members including piles are modeled by spacial beam or shell elements, and raft-base is divided into thick-shell elements with its spring-dashpot boundary coefficient obtained by impedance backcalculated. The mass-spring elements for soil between piles are set to simulate vertical, horizontal pile-group effects by strata-equivalent approach. The soil beside composite body is separated into near-field and far-field parts. The former is modeled by nonlinear spring-dashpot elements based on Winkler’s hypothesis, while the latter is modeled by a series of linear mass-spring-dashpots. With the effects of boundary track forces and energy radiation, the presented model enables researchers to conduct the time-domain nonlinear analysis in a relatively simple manner which avoids sophisticated boundary method and solid-element mesh bringing with tremendous computational cost. The seismic effect on dynamic interaction of pile-soil-complicated structures would be efficiently annotated from two structural engineering and geotechnical engineering aspects and the numerical calculation effort would be drastically decreased too. The complete procedure is mainly performed using the parametric design language assembled in the Finite Element Code Ansys. With the dynamic analysis of foundation and superstructure for a pile-supported 15-storey building, the influence of the participant effect on structural dynamic response will be depicted by various dynamic parameters of pile-soil-raft foundation in detail. Not only do the results have an agreement with some conclusions drawn by the general interaction theory, but also certain of phenomena which would be disagree with that by general analysis is involved. Even with the finite-element meshes for 68 piles, the time-history analysis procedure for PGSS (Pile-Group-Soil-Superstructure) system and the qualitative evaluation with various SSI parameters can be also fulfilled efficiently and rapidly by presented means. These results may be of help to the designers to quickly assess the significance of interaction effect for the high-rise buildings resting on any type or layout of pile-group foundation.


Author(s):  
Kai Zhou ◽  
Pei Cao ◽  
Jiong Tang

Uncertainty quantification is an important aspect in structural dynamic analysis. Since practical structures are complex and oftentimes need to be characterized by large-scale finite element models, component mode synthesis (CMS) method is widely adopted for order-reduced modeling. Even with the model order-reduction, the computational cost for uncertainty quantification can still be prohibitive. In this research, we utilize a two-level Gaussian process emulation to achieve rapid sampling and response prediction under uncertainty, in which the low- and high-fidelity data extracted from CMS and full-scale finite element model are incorporated in an integral manner. The possible bias of low-fidelity data is then corrected through high-fidelity data. For the purpose of reducing the emulation runs, we further employ Bayesian inference approach to calibrate the order-reduced model in a probabilistic manner conditioned on multiple predicted response distributions of concern. Case studies are carried out to validate the effectiveness of proposed methodology.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Lihong Zhang ◽  
Tianyun Liu ◽  
Qingbin Li

This paper presents a new robust and efficient time integration algorithm suitable for various complex nonlinear structural dynamic finite element problems. Based on the idea of composition, the three-point backward difference formula and a generalized central difference formula are combined to constitute the implicit algorithm. Theoretical analysis indicates that the composite algorithm is a single-solver algorithm with satisfactory accuracy, unconditional stability, and second-order convergence rate. Moreover, without any additional parameters, the composite algorithm maintains a symmetric effective stiffness matrix and the computational cost is the same as that of the trapezoidal rule. And more merits of the proposed algorithm are revealed through several representative finite element examples by comparing with analytical solutions or solutions provided by other numerical techniques. Results show that not only the linear stiff problem but also the nonlinear problems involving nonlinearities of geometry, contact, and material can be solved efficiently and successfully by this composite algorithm. Thus the prospect of its implementation in existing finite element codes can be foreseen.


2015 ◽  
Vol 138 (1) ◽  
Author(s):  
Charanjeet Kaur Malhi ◽  
Rudra Pratap

In this work, we address the issue of modeling squeeze film damping in nontrivial geometries that are not amenable to analytical solutions. The design and analysis of microelectromechanical systems (MEMS) resonators, especially those that use platelike two-dimensional structures, require structural dynamic response over the entire range of frequencies of interest. This response calculation typically involves the analysis of squeeze film effects and acoustic radiation losses. The acoustic analysis of vibrating plates is a very well understood problem that is routinely carried out using the equivalent electrical circuits that employ lumped parameters (LP) for acoustic impedance. Here, we present a method to use the same circuit with the same elements to account for the squeeze film effects as well by establishing an equivalence between the parameters of the two domains through a rescaled equivalent relationship between the acoustic impedance and the squeeze film impedance. Our analysis is based on a simple observation that the squeeze film impedance rescaled by a factor of jω, where ω is the frequency of oscillation, qualitatively mimics the acoustic impedance over a large frequency range. We present a method to curvefit the numerically simulated stiffness and damping coefficients which are obtained using finite element analysis (FEA) analysis. A significant advantage of the proposed method is that it is applicable to any trivial/nontrivial geometry. It requires very limited finite element method (FEM) runs within the frequency range of interest, hence reducing the computational cost, yet modeling the behavior in the entire range accurately. We demonstrate the method using one trivial and one nontrivial geometry.


Author(s):  
Aabhas Singh ◽  
Matthew S. Allen ◽  
Robert J. Kuether

Abstract Structural dynamic finite element models typically use multipoint constraints (MPC) to condense the degrees of freedom (DOF) near bolted joints down to a single node, which can then be joined to neighboring structures with linear springs or nonlinear elements. Scalability becomes an issue when multiple joints are present in a system, because each requires its own model to capture the nonlinear behavior. While this increases the computational cost, the larger problem is that the parameters of the joint models are not known, and so one must solve a nonlinear model updating problem with potentially hundreds of unknown variables to fit the model to measurements. Furthermore, traditional MPC approaches are limited in how the flexibility of the interface is treated (i.e. with rigid bar elements the interface has no flexibility). To resolve this shortcoming, this work presents an alternative approach where the contact interface is reduced to a set of modal DOF which retain the flexibility of the interface and are capable of modeling multiple joints simultaneously. Specifically, system-level characteristic constraint (S-CC) reduction is used to reduce the motion at the contact interface to a small number of shapes. To capture the hysteresis and energy dissipation that is present during microslip of joints, a hysteretic element is applied to a small number of the S-CC Shapes. This method is compared against a traditional MPC method (using rigid bar elements) on a two-dimensional finite element model of a cantilever beam with a single joint near the free end. For all methods, a four-parameter Iwan element is applied to the interface DOF to capture how the amplitude dependent modal frequency and damping change with vibration amplitude.


2015 ◽  
Vol 8 (4) ◽  
pp. 582-604
Author(s):  
Zhengqin Yu ◽  
Xiaoping Xie

AbstractThis paper proposes and analyzes semi-discrete and fully discrete hybrid stress finite element methods for elastodynamic problems. A hybrid stress quadrilateral finite element approximation is used in the space directions. A second-order center difference is adopted in the time direction for the fully discrete scheme. Error estimates of the two schemes, as well as a stability result for the fully discrete scheme, are derived. Numerical experiments are done to verify the theoretical results.


2012 ◽  
Vol 204-208 ◽  
pp. 4455-4459 ◽  
Author(s):  
Liu Hong Chang ◽  
Chang Bo Jiang ◽  
Man Jun Liao ◽  
Xiong Xiao

The explicit dynamic finite element theory is applied on the collision of ships with buoys for computer simulation. Using ANSYS/LS-DYNA finite element analysis software, the numerical simulation of the collision between the ton ship and the buoy with different structures and impact points. The collision force, deformation, displacement parameters and the weak impact points of a buoy are obtained. Based on the numerical simulation results, analysis of buoys and structural collision damages in anti-collision features are discussed, and several theoretical sugestions in anti-collision for the design of buoy are provided.


2014 ◽  
Vol 931-932 ◽  
pp. 994-998
Author(s):  
Rangsan Wannapop ◽  
Thira Jearsiripongkul ◽  
Thawatchai Boonluang

This research represents a design and analysis of Automatic loading copper wire machine for the actuator arm (ALCM). The process of copper wire placement on a single actuator arm type compensates human workers. In this research, copper wire placement set is made as a 3D model by computer program before undergoes arrangement analysis via explicit dynamic finite element method to study a suitable speed for copper wire placing. It is considered by characteristics of copper wire after placed and failures occurred during the process that will define suitable speed of motor rotation. The suitable speed is corresponding to copper wire characteristic as preferred, prevent copper wire fracture and time reduction compare to human work.


Sign in / Sign up

Export Citation Format

Share Document