Prediction of mode I delamination resistance of z-pinned laminates using the embedded finite element technique

A new finite modelling approach is presented to analyse the mode I delamination fracture toughness of z-pinned laminates using the computationally efficient embedded element technique. In the FE model,each z-pin is represented by a single one-dimensional truss element that is embedded within the laminate. Each truss is given the material, geometric and spatial properties associated with the global crackbridging traction response of a z-pin in the laminate; this simplification provides a computationally efficient and flexible model where pin elements are independent of the underlying structural mesh for thelaminate. The accuracy of the FE modelling approach is assessed using mode I interlaminar fracture toughness data for a carbon-epoxy laminate reinforced with z-pins made of copper, titanium or stainless steel. The model is able to predict with good accuracy the crack growth resistance curves and fracture toughness properties for the different types of z-pinned laminate.

Porohyperelastic analysis to explore mechanical properties of chondrocytes using numerical modeling and experiments : a finite element study

There are many continuum mechanical models have been developed such as liquid drop models, solid models, and so on for single living cell biomechanics studies. However, these models do not give a fully approach to exhibit a clear understanding of the behaviour of single living cells such as swelling behaviour, drag effect, etc. Hence, the porohyperelastic (PHE) model which can capture those aspects would be a good candidature to study cells behaviour (e.g. chondrocytes in this study). In this research, an FEM model of single chondrocyte cell will be developed by using this PHE model to simulate Atomic Force Microscopy (AFM) experimental results with the variation of strain rate. This material model will be compared with viscoelastic model to demonstrate the advantages of PHE model. The results have shown that the maximum value of force applied of PHE model is lower at lower strain rates. This is because the mobile fluid does not have enough time to exude in case of very high strain rate and also due to the lower permeability of the membrane than that of the protoplasm of chondrocyte. This behavior is barely observed in viscoelastic model. Thus, PHE model is the better model for cell biomechanics studies.

Stiffness analysis of polymeric composites using the finite element method

The ability to predict the mechanical behavior of polymer composites is crucial for their design and manufacture. Extensive studies based on both macro- and micromechanical analyses are used to develop new insights into the behavior of composites. In this respect, finite element modeling has proved to be a particularly powerful tool. In this article, we present a Galerkin scheme in conjunction with the penalty method for elasticity analyses of different types of polymer composites. In this scheme, the application of Green's theorem to the model equation results in the appearance of interfacial flux terms along the boundary between the filler and polymer matrix. It is shown that for some types of composites these terms significantly affect the stress transfer between polymer and fillers. Thus, inclusion of these terms in the working equations of the scheme preserves the accuracy of the model predictions. The model is used to predict the most important bulk property of different types of composites. Composites filled with rigid or soft particles, and composites reinforced with short or continuous fibers are investigated. For each case, the results are compared with the available experimental results and data obtained from other models reported in the literature. Effects of assumptions made in the development of the model and the selection of the prescribed boundary conditions are discussed.

Statistical analysis of finite element modeling in sheet metal forming and springback analysis

This work focuses on development of a method to statistically study forming and springback problems of TRansformation Induced Plasticity (TRIP) through an industrial case study. A Design of Experiments (DOE) approach was used to study the sensitivity of predictions to four user input parameters in implicit and explicit sheet metal forming codes. Numerical results were compared to experimental measurements of parts stamped in an industrial production line. The accuracy of forming strain predictions for TRIP steel were comparable with conventional steel, but the springback predictions of TRIP steel were far less accurate. The statistical importance of selected parameters for forming and springback prediction is also discussed. Changes of up to ±10% in Young's modulus and coefficient of friction were found to be insignificant in improving or deteriorating the statistical correlation of springback accuracies.

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.

Dynamic characterisation of four nine-story large-panel R.C. buildings in Bishkek (Kyrgyzstan): A comparison between experimental ambient vibration analysis and numerical finite element modeling

With the outlook of improving seismic vulnerability assessment for the city of Bishkek (Kyrgyzstan), the global dynamic behaviour of four nine-storey r.c. large-panel buildings in elastic regime is studied. The four buildings were built during the Soviet era within a serial production system. Since they all belong to the same series, they have very similar geometries both in plan and in height. Firstly, ambient vibration measurements are performed in the four buildings. The data analysis composed of discrete Fourier transform, modal analysis (frequency domain decomposition) and deconvolution interferometry, yields the modal characteristics and an estimate of the linear impulse response function for the structures of the four buildings. Then, finite element models are set up for all four buildings and the results of the numerical modal analysis are compared with the experimental ones. The numerical models are finally calibrated considering the first three global modes and their results match the experimental ones with an error of less then 20%.

Partial delamination modeling in composite beams using a finite element method

A new method of modeling partial delamination in composite beams is proposed and implemented using the finite element method. Homogenized cross-sectional stiffness of the delaminated beam is obtained by the proposed analytical technique, including extension-bending, extension-twist and torsion-bending coupling terms, and hence can be used with an existing finite element method. A two noded C1 type Timoshenko beam element with 4 degrees of freedom per node for dynamic analysis of beams is implemented. The results for different delamination scenarios and beams subjected to different boundary conditions are validated with available experimental results in the literature and/or with the 3D finite element simulation using COMSOL. Results of the first torsional mode frequency for the partially delaminated beam are validated with the COMSOL results. The key point of the proposed model is that partial delamination in beams can be analyzed using a beam model, rather than using 3D or plate models. (c) 2013 Elsevier B.V. All rights reserved.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Fracture analysis of spot-welds with edge cracks using 2-D hybrid special finite element

This paper employed a systematic analysis using a 2-D hybrid special finite element containing an edge crack in order to describe the fracture behavior of spot-welds in automotive structures. The 2-D hybrid special finite element is derived form a mixed formulation with a complex potential function with the description of the singularity of a stress field. The hybrid special finite element containing an edge crack can give a better description of its singularity with only one hybrid element surrounding one crack. The advantage of this special element is that it can greatly simplify the numerical modeling of the spot welds. Some numerical examples demonstrate the validity and versatility of the present analysis method. The lap-shear, lap-tension and angle-clip specimens are analyzed and some useful fracture parameters such as the stress intensity factor and the initial direction of crack growth are obtained simultaneously.

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.

Complex modal analysis of a vertical rotor through finite element method

Natural frequencies were analyzed (axial, torsional and flexural) and frequency response of a vertical rotor with a hard disk at the edge through the classical and complex modal analysis. The mathematical modeling was based on the theory of Euler-Bernoulli beam. The equation that rules the movement was obtained through the Lagrangian formulation. The model considered the effects of bending, torsion and axial deformation of the shaft, besides the gravitational and gyroscopic effects. The finite element method was used to discretize the structure into hollow cylindrical elements with 12 degrees of freedom. Mass, stiffness and gyroscopic matrices were explained consistently. This type of tool, based on the use of complex coordinates to describe the dynamic behavior of rotating shaft, allows the decomposition of the system in two submodes, backward and forward. Thus, it is possible to clearly visualize that the orbit and direction of the precessional motion around the line of the rotating shaft is not deformed. A finite element program was developed using Matlab ®, and numerical simulations were performed to validate this model.

The Galerkin finite element approximation of the fractional cable equation

The cable equation is one of the most fundamental equations for modeling neuronal dynamics. Cable equations with a fractional order temporal derivative have been introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, the fractional cable equation involving two integro-differential operators is considered. The Galerkin finite element approximations of the fractional cable equation are proposed. The main contribution of this work is outlined as follow: • A semi-discrete finite difference approximation in time is proposed. We prove that the scheme is unconditionally stable, and the numerical solution converges to the exact solution with order O(Δt). • A semi-discrete difference scheme for improving the order of convergence for solving the fractional cable equation is proposed, and the numerical solution converges to the exact solution with order O((Δt)2). • Based on the above semi-discrete difference approximations, Galerkin finite element approximations in space for a full discretization are also investigated. • Finally, some numerical results are given to demonstrate the theoretical analysis.

Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data

The finite element (FE) analysis is an effective method to study the strength and predict the fracture risk of endodontically-treated teeth. This paper presents a rapid method developed to generate a comprehensive tooth FE model using data retrieved from micro-computed tomography (μCT). With this method, the inhomogeneity of material properties of teeth was included into the model without dividing the tooth model into different regions. The material properties of the tooth were assumed to be related to the mineral density. The fracture risk at different tooth portions was assessed for root canal treatments. The micro-CT images of a tooth were processed by a Matlab software programme and the CT numbers were retrieved. The tooth contours were obtained with thresholding segmentation using Amira. The inner and outer surfaces of the tooth were imported into Solidworks and a three-dimensional (3D) tooth model was constructed. An assembly of the tooth model with the periodontal ligament (PDL) layer and surrounding bone was imported into ABAQUS. The material properties of the tooth were calculated from the retrieved CT numbers via ABAQUS user's subroutines. Three root canal geometries (original and two enlargements) were investigated. The proposed method in this study can generate detailed 3D finite element models of a tooth with different root canal enlargements and filling materials, and would be very useful for the assessment of the fracture risk at different tooth portions after root canal treatments.