Evaluation of Stress Wave Attenuation in a Polymer Matrix Composite Using Finite Element Analysis Technique

Author(s):  
Shrikant Nargund

The dynamic behavior of polymer composites is significantly affected by the properties of their micro constituents including shape and size of inclusions and inclusions/matrix adhesion properties. Wave propagation through such a composite is a complex phenomenon as it includes random scattering, absorption and transmittance of the incident wave and is dependent upon factors such as the properties, size and placement of the inclusions inside the matrix. Finite element modeling provides a viable approach for investigating the effects of micro constituent structure on the dynamic behavior of polymer composites. In this paper, we investigate the stress wave attenuation characteristics of a particulate polymer matrix composite using Finite Element (FE) analysis approach. The wave attenuation of ultrasonic sinusoidal waves of frequency ranging from 1 MHz to 4 MHz is evaluated for different FE models. The spherical inclusions are randomly distributed inside the polymer matrix with a certain minimum distance apart from each other. Inclusion-Matrix adhesion properties are studied by modeling a small region at the interface of inclusions and matrix known as interphase region. The interphase region is modeled explicitly using the cohesive zone modeling approach to study how the properties of this region will affect the wave attenuation characteristics of the polymer composite. Cohesive zone models are governed by traction separation law which helps in the measurement of the inclusion-matrix bonding strength and also allow the study of de-bonding at the interface in the critically stressed region produced due application of load. Thus the FE models consist of three phases; polymer matrix, particulate inclusions and the interphase region. Various three dimensional FE models are created using 3D tetrahedral/hexahedral elements by varying the radius of the spherical inclusions and by varying volume fraction of the inclusions. The analyses are performed using a general purpose finite element software LS-Dyna. A rate dependent viscoelastic material model with four terms in prony series expansion is used for modeling the polymer matrix. A linear elastic isotropic material model is used for modeling the inclusions. The wave attenuation is measured as reduction in the amplitude of the wave as it passes through the composite. A comparison of results for various models is done to check for general trend of attenuation coefficient as a function of size of inclusions, volume fraction of inclusions, frequency of loading and interphase region properties. Results show that volume fraction and load frequency have a maximum effect on the wave attenuation coefficient. Interphase region stiffness and interface de-bonding also plays an important role in attenuation characteristics of the polymer composite.

Author(s):  
Shank S. Kulkarni ◽  
Alireza Tabarraei ◽  
Pratik P. Ghag

The properties of the inclusions, viz. size, shape, and distribution significantly affect macroscopic properties of a polymer composite. Finite element (FE) modeling provides a viable approach for investigating the effects of the inclusions on the macroscopic properties of the polymer composite. In this paper, finite element method is used to investigate ultrasonic wave propagation in polymer matrix composite with a dispersed phase of inclusions. The finite element models are made up of three phases; viz. the polymer matrix, inclusions (micro constituent), and interphase zones between the inclusions and the polymer matrix. The analysis is performed on a three dimensional finite element model and the attenuation characteristics of ultrasonic longitudinal waves in the matrix are evaluated. The attenuation in polymer composite is investigated by changing the size, volume fraction of inclusions, and addition of interphase layer. The effect of loading frequency of the wave on the attenuation characteristics is also studied by varying the frequency in the range of 1–4 MHz. Results of the test revealed that higher volume fraction of inclusions gave higher attenuation in the polymer composite as compared to the lower volume fraction model. Smaller size of inclusions are preferred over larger size as they give higher wave attenuation. It was found that the attenuation characteristics of the polymer composite are better at higher frequencies as compared to lower frequencies. It is also concluded that the interphase later plays a significant role in the attenuation characteristics of the composite.


Author(s):  
Shank S. Kulkarni ◽  
Alireza Tabarraei ◽  
Pratik Ghag

Abstract This work studies the damping property of Nanocomposites through simulating wave propagation using the Finite Element Method (FEM). For this purpose Representative Volume Element (RVE) of the composite material is created using Random Sequential Absorption (RSA) algorithm. Damping property is represented using the wave attenuation coefficient. The matrix material is assumed to be isotropic visco-elastic in nature with randomly dispersed stiff elastic spherical fillers. In order to model mechanical imperfections at the boundary of fillers and matrix, the interphase layer is modeled surrounding the spherical fillers. Determining the thickness and stiffness of this interphase region is a challenging task. Therefore this study aims at investigating the effect of variation in thickness and stiffness values of the interphase region on damping property of whole composite using sensitivity analysis. Two specific cases with a volume fraction of 5 % and 8.6 % are selected for sensitivity analysis. It has been found that both the thickness and stiffness of the interphase region plays an important role in deciding the damping properties of the polymer composite. Value of attenuation coefficient is more sensitive to the thickness of interphase than stiffness and hence it is important to choose the value of thickness correctly for accurate predictions.


2021 ◽  
Vol 3 (1) ◽  
Author(s):  
Anas Kanan ◽  
Michael Kaliske

AbstractDielectric elastomer actuators (DEA) have been demonstrated to exhibit a quasi-immediate electro-mechanical actuation response with relatively large deformation capability. The properties of DEA make them suitable to be used in the form of major active components within soft robotics and biomimetic artificial muscles. However, some of the electro-active material properties impose limitations on its applications. Therefore, researchers attempt to modify the structure of the homogeneous DEA material by the incorporation of fillers that possess distinct electro-mechanical properties. This modification of the material’s structure leads to a fabricated inhomogeneous composite. From the point of mathematical material modelling and numerical simulation, we propose a material model and a computational framework using the finite element method, which is capable of emulating nonlinear electro-elastic interactions. We consider a coupled electro-mechanical description with the electric and the electro-mechanical properties of the material assumed to be nonlinearly dependent on the deformation. Furthermore, we demonstrate a coupled ansatz that expresses the electric response as dielectrically quasi-linear with only density-dependent electric permittivity. We couple the electro-mechanical models to the extended tube model, which is a suitable approach for the realistic emulation of the hyperelastic response of rubber-like materials. Thereafter, we demonstrate analytical and numerical solutions of a homogeneous electro-elastic body with the Neo-Hookean material model and the extended tube model to express the hyperelastic response. Finally, we use the finite element method to investigate several heterogeneous configurations consisting of soft DEA matrix filled with spherical stiff inclusions with changing volume fraction and ellipsoidal inclusions with varying aspect ratio.


Author(s):  
D. S. Pritchard

The effect of varying the strain rate loading conditions in compression on a copper single crystal dispersion-hardened with SiO2 particles has been examined. These particles appear as small spherical inclusions in the copper lattice and have a volume fraction of 0.6%. The structure of representative crystals was examined prior to any testing on a transmission electron microscope (TEM) to determine the nature of the dislocations initially present in the tested crystals. Only a few scattered edge and screw dislocations were viewed in those specimens.


TAPPI Journal ◽  
2012 ◽  
Vol 11 (1) ◽  
pp. 61-66 ◽  
Author(s):  
DOEUNG D. CHOI ◽  
SERGIY A. LAVRYKOV ◽  
BANDARU V. RAMARAO

Delamination between layers occurs during the creasing and subsequent folding of paperboard. Delamination is necessary to provide some stiffness properties, but excessive or uncontrolled delamination can weaken the fold, and therefore needs to be controlled. An understanding of the mechanics of delamination is predicated upon the availability of reliable and properly calibrated simulation tools to predict experimental observations. This paper describes a finite element simulation of paper mechanics applied to the scoring and folding of multi-ply carton board. Our goal was to provide an understanding of the mechanics of these operations and the proper models of elastic and plastic behavior of the material that enable us to simulate the deformation and delamination behavior. Our material model accounted for plasticity and sheet anisotropy in the in-plane and z-direction (ZD) dimensions. We used different ZD stress-strain curves during loading and unloading. Material parameters for in-plane deformation were obtained by fitting uniaxial stress-strain data to Ramberg-Osgood plasticity models and the ZD deformation was modeled using a modified power law. Two-dimensional strain fields resulting from loading board typical of a scoring operation were calculated. The strain field was symmetric in the initial stages, but increasing deformation led to asymmetry and heterogeneity. These regions were precursors to delamination and failure. Delamination of the layers occurred in regions of significant shear strain and resulted primarily from the development of large plastic strains. The model predictions were confirmed by experimental observation of the local strain fields using visual microscopy and linear image strain analysis. The finite element model predicted sheet delamination matching the patterns and effects that were observed in experiments.


2008 ◽  
Vol 36 (1) ◽  
pp. 63-79 ◽  
Author(s):  
L. Nasdala ◽  
Y. Wei ◽  
H. Rothert ◽  
M. Kaliske

Abstract It is a challenging task in the design of automobile tires to predict lifetime and performance on the basis of numerical simulations. Several factors have to be taken into account to correctly estimate the aging behavior. This paper focuses on oxygen reaction processes which, apart from mechanical and thermal aspects, effect the tire durability. The material parameters needed to describe the temperature-dependent oxygen diffusion and reaction processes are derived by means of the time–temperature–superposition principle from modulus profiling tests. These experiments are designed to examine the diffusion-limited oxidation (DLO) effect which occurs when accelerated aging tests are performed. For the cord-reinforced rubber composites, homogenization techniques are adopted to obtain effective material parameters (diffusivities and reaction constants). The selection and arrangement of rubber components influence the temperature distribution and the oxygen penetration depth which impact tire durability. The goal of this paper is to establish a finite element analysis based criterion to predict lifetime with respect to oxidative aging. The finite element analysis is carried out in three stages. First the heat generation rate distribution is calculated using a viscoelastic material model. Then the temperature distribution can be determined. In the third step we evaluate the oxygen distribution or rather the oxygen consumption rate, which is a measure for the tire lifetime. Thus, the aging behavior of different kinds of tires can be compared. Numerical examples show how diffusivities, reaction coefficients, and temperature influence the durability of different tire parts. It is found that due to the DLO effect, some interior parts may age slower even if the temperature is increased.


2001 ◽  
Vol 29 (3) ◽  
pp. 186-196 ◽  
Author(s):  
X. Yan

Abstract A method is described to predict relative body turn up endurance of radial truck tires using the finite element method. The elastomers in the tire were simulated by incompressible elements for which the nonlinear mechanical properties were described by the Mooney-Rivlin model. The belt, carcass, and bead were modeled by an equivalent orthotropic material model. The contact constraint of a radial tire structure with a flat foundation and rigid rim was treated using the variable constraint method. Three groups of tires with different body turn up heights under inflation and static footprint loading were analyzed by using the finite element method. Based on the detail analysis for stress analysis parameters in the critical regions in the tires, the relative body turn up edge endurance was predicted.


Mathematics ◽  
2020 ◽  
Vol 9 (1) ◽  
pp. 25
Author(s):  
Bagh Ali ◽  
Rizwan Ali Naqvi ◽  
Amna Mariam ◽  
Liaqat Ali ◽  
Omar M. Aldossary

The below work comprises the unsteady flow and enhanced thermal transportation for Carreau nanofluids across a stretching wedge. In addition, heat source, magnetic field, thermal radiation, activation energy, and convective boundary conditions are considered. Suitable similarity functions use to transmuted partial differential formulation into the ordinary differential form, which is solved numerically by the finite element method and coded in Matlab script. Parametric computations are made for faster stretch and slowly stretch to the surface of the wedge. The progressing value of parameter A (unsteadiness), material law index ϵ, and wedge angle reduce the flow velocity. The temperature in the boundary layer region rises directly with exceeding values of thermophoresis parameter Nt, Hartman number, Brownian motion parameter Nb, ϵ, Biot number Bi and radiation parameter Rd. The volume fraction of nanoparticles rises with activation energy parameter EE, but it receded against chemical reaction parameter Ω, and Lewis number Le. The reliability and validity of the current numerical solution are ascertained by establishing convergence criteria and agreement with existing specific solutions.


Sign in / Sign up

Export Citation Format

Share Document