Hybrid special finite element method for bi-material crack

In the paper, two novel 2-D hybrid special finite elements each containing an interfacial edge crack, which lies along or vertical to the interface between two materials, are developed. These proposed elements can assure the high precision especially in the vicinity of crack tip and provide a better description of its singularity with only one hybrid element surrounding one interfacial crack, thus, the numerical modeling of fracture analysis on bi-material crack can be greatly simplified. Numerical examples are provided to demonstrate the validity and versatility of the proposed method.

Nonlinear dynamic modeling of ionic polymer conductive network composite actuators using rigid finite element method

Ionic polymer conductive network composite (IPCNC) actuators are a class of electroactive polymer composites that exhibit some interesting electromechanical characteristics such as low voltage actuation, large displacements, and benefit from low density and elastic modulus. Thus, these emerging materials have potential applications in biomimetic and biomedical devices. Whereas significant efforts have been directed toward the development of IPMC actuators, the establishment of a proper mathematical model that could effectively predict the actuators' dynamic behavior is still a key challenge. This paper presents development of an effective modeling strategy for dynamic analysis of IPCNC actuators undergoing large bending deformations. The proposed model is composed of two parts, namely electrical and mechanical dynamic models. The electrical model describes the actuator as a resistive-capacitive (RC) transmission line, whereas the mechanical model describes the actuator as a system of rigid links connected by spring-damping elements. The proposed modeling approach is validated by experimental data, and the results are discussed. © 2014 Elsevier B.V. All rights reserved.

A study of guided wave propagation in timber pole using spectral finite element method

Timber is one of the most widely used structural material all over the world. Round timbers can be seen as a structural component in historical buildings, jetties, short span bridges and also as piles for foundation and poles for electrical and power distribution. To evaluate the current condition of these cylindrical type timber structures, guided wave has a great potential. However, the difficulties associated with the guided wave propagation in timber materials includes orthotropic behaviour of wood, moisture contents, temperature, grain direction, etc. In addition, the effect of fully or partially filled surrounding media, such as soil, water, etc. causes attenuation on the generated stress wave. In order to investigate the effects of these parameters on guided wave propagation, extensive numerical simulation is required to conduct parametric studies. Moreover, due to the presence of multi modes in guided wave propagation, dispersion curves are of great importance. Even though conventional finite element method (FEM) can determine dispersion curves along with wave propagation in time domain, it is highly computationally expensive. Furthermore, incorporating orthotropic behaviour and surrounding media to model a thick cylindrical wave (large diameter cylindrical structures) make conventional FEM inefficient for this purpose. In contrast, spectral finite element method (SFEM) is a semi analytical method to model the guided wave propagation which does not need fine meshes compared to the other methods, such as FEM or finite difference method (FDM). Also, even distribution of mass and stiffness of structures can be obtained with very few elements using SFEM. In this paper, the suitability of SFEM is investigated to model guided wave propagation through an orthotropic cylindrical waveguide with the presence of surrounding soil. Both the frequency domain analysis (dispersion curves) and time domain reconstruction for a multi-mode generated input signal are presented under different loading location. The dispersion curves obtained from SFEM are compared against analytical solution to verify its accuracy. Lastly, different numerical issues to solve for the dispersion curves and time domain results using SFEM are also discussed.

A finite element simulation of micro-mechanical frictional behaviour in metal forming

Friction is a critical factor for sheet metal forming (SMF). The Coulomb friction model is usually used in most finite element (FE) simulation for SMF. However, friction is a function of the local contact deformation conditions, such as local pressure, roughness and relative velocity. Frictional behaviour between contact surfaces can be based on three cases: boundary, hydrodynamic and mixed lubrication. In our microscopic friction model based on the finite element method (FEM), the case of dry contact between sheet and tool has been considered. In the view of microscopic geometry, roughness depends upon amplitude and wavelength of surface asperities of sheet and tool. The mean pressure applied on the surface differs from the pressure over the actual contact area. The effect of roughness (microscopic geometric condition) and relative speed of contact surfaces on friction coefficient was examined in the FE model for the microscopic friction behaviour. The analysis was performed using an explicit FE formulation. In this study, it was found that the roughness of deformable sheet decreases during sliding and the coefficient of friction increases with increasing roughness of contact surfaces. Also, the coefficient of friction increases with the increase of relative velocity and adhesive friction coefficient between contact surfaces.

Object oriented finite element analysis: a distributed approach to mesh generation

The object-oriented finite element method (OOFEM) has attracted the attention of many researchers. Compared with the traditional finite element method, OOFEM software has the advantages of maintenance and reuse. Moreover, it is easier to expand the architecture to a distributed one. In this paper, we introduce a distributed architecture of a object-oriented finite element preprocessor. A comparison between the distributed system and the centralised system shows that the former, presented in the paper, greatly improves the performance of mesh generation. Other finite element analysis modules could be expanded according to this architecture.

The design and development of an integrated system for object-oriented finite element computing

Finite Element Method (FEM) is widely used in Science and Engineering since 1960’s. The vast majority of FEM software is procedure-oriented. However, this conventional style of designing FEM software encounters problems in maintenance, reuse, and expansion of the software. Recently the object-oriented finite element method attracts the attention of lots of researchers, and now there is a growing interest in this method. In this paper, the object-oriented finite element (OOFE) is briefly introduced. Then the design and development of an integrated OOFE system is described. A comparison of the integrated OOFE system and a procedure-oriented system shows that our OOFE system has many advantages.

Finite element simulation of deformation processes of micro-components

The aim of this paper is to improve the understanding of deformation of micro medical needle and thread during assembly and then to develop an economical and flexible deformation method. Therefore, the swaging process is computationally simulated with the finite element method in this paper. A commercially available explicit nonlinear finite element analysis code, LS-Dyna, is used to model the 3-D deformation and contact problem. As the firmness of the assembly on the needle depends on the contact force and friction, the contact and the slide between the needle and thread are taken into account in the simulation. The general surface-to-surface contact algorithm (STS) is used to simulate the contact. The paper provides an insight into the deformation of the micro products.

Three-dimensional Finite Element modelling of laser shock peening process

Laser shock peening (LSP) is an innovative surface treatment technique for metal alloys, with the great improvement of their fatigue, corrosion and wear resistance performance. Finite element method has been widely applied to simulate the LSP to provide the theoretically predictive assessment and optimally parametric design. In the current work, 3-D numerical modelling approaches, combining the explicit dynamic analysis, static equilibrium analysis algorithms and different plasticity models for the high strain rate exceeding 106s-1, are further developed. To verify the proposed methods, 3-D static and dynamic FEA of AA7075-T7351 rods subject to two-sided laser shock peening are performed using the FEA package–ABAQUS. The dynamic and residual stress fields, shock wave propagation and surface deformation of the treated metal from different material modelling approaches have a good agreement.

Finite element simulations of the cyclic elastoplastic behaviour of copper thin films

A large-scale computational and statistical strategy is presented to investigate the development of plastic strain heterogeneities and plasticity induced roughness at the free surface in multicrystalline films subjected to cyclic loading conditions, based on continuum crystal plasticity theory. The distribution of plastic strain in the grains and its evolution during cyclic straining are computed using the finite element method in films with different ratios of in-plane grain size and thickness, and as a function of grain orientation (grains with a {1 1 1} or a {0 0 1} plane parallel to the free surface and random orientations). Computations are made for 10 different realizations of aggregates containing 50 grains and one large aggregate with 225 grains. It is shown that overall cyclic hardening is accompanied by a significant increase in strain dispersion. The case of free-standing films is also addressed for comparison. The overall surface roughness is shown to saturate within 10 to 15 cycles. Plasticity induced roughness is due to the higher deformation of {0 0 1} and random grains and due to the sinking or rising at some grain boundaries.

Finite element modeling and analysis of thrombosis in middle cerebral artery

Thrombotic stroke, which is caused by blood clot in the cerebral artery, is a major source of increased mortality and morbidity. Considering as efficient and fastest methods, mathematical approaches have gained significant importance for analyzing and understanding the biological events like thrombosis. This paper presents a computational model to analyze the effects of thrombosis using the theory of coupled fluid dynamics-structure interaction. The finite element method is used for the modeling of thrombosis (blood clot) of different stages in the middle cerebral artery with physiological compliance. The developed model is used to investigate the consequences that occur due to the various sizes of clots in the artery in the form of blood flow velocity, blood pressure, and artery wall stress. Such numerical assessment will facilitate better understanding of the biophysical process in case of thrombosis and thus would support medical practitioners to take faster curing steps.