Composites Science and Technology
|
|
- Naomi Norman
- 6 years ago
- Views:
Transcription
1 Composites Science and Technology 69 (2009) Contents lists available at ScienceDirect Composites Science and Technology journal homepage: Simulation of interphase percolation and gradients in polymer nanocomposites Rui Qiao a, L. Catherine Brinson b, * a Department of Mechanical Engineering, Northwestern University, Evanston, IL 60208, United States b Department of Materials Science and Engineering, Northwestern University, Evanston, IL 60208, United States article info abstract Article history: Received 4 August 2008 Received in revised form 25 October 2008 Accepted 20 November 2008 Available online 6 December 2008 Keywords: A. Polymer-matrix composites (PMCs) A. Nanocomposites B. Interphase C. Finite element analysis Interphase percolation Experimental data suggests that well dispersed nanoparticles within a polymer matrix induce a significant interphase zone of altered polymer mobility surrounding each nanoparticle, which can lead to a percolating interphase network inside of the composite. To investigate this concept and the nature of the interphase, a two-dimensional finite element model is developed to study the impact of interphase zones on the overall properties of the composite. Thirty non-overlapping identical circular inclusions are randomly distributed in the matrix with layers of interphase surrounding the inclusions. The simulation results clearly show that the loss moduli of composites are either broadened or shifted corresponding to the absence or presence of a geometrically percolating interphase network. Our numerical study correlates well with experimental data showing broadening of loss peaks for unfunctionalized composites and a large shift of the loss modulus for functionalized nanotube polymer composites. Further, our results indicate the existence of a gradient in properties of the interphase layer and that incorporating this gradient into modeling is critical to reflect the behavior of polymer nanocomposites. Ó 2008 Elsevier Ltd. All rights reserved. 1. Introduction * Corresponding author. Tel.: ; fax: address: cbrinson@northwestern.edu (L. Catherine Brinson). Polymer nanocomposites have attracted intense attention in the past decade, due to their unique properties and great potential as future materials. Recent experimental work demonstrated that by incorporating nanoscale inclusions, polymer nanocomposites can exhibit significant improvement in mechanical, electrical, thermal, and other physical properties in comparison to their parent polymer systems. Effects of nanoreinforcement include dramatically enhanced strength, stiffness, fracture toughness, increased thermal stability and heat distortion temperature, increased chemical resistance, and increased electrical conductivity while maintaining optical clarity [1 11]. Most important is that these superior properties are achieved at very low loading levels of inclusions, so the parent polymer does not sacrifice the advantages of low density and high processibility. These extraordinary behaviors make polymer nanocomposites a promising multifunctional material in many fields, including the aerospace, automotive, and medical device industries. A variety of nanoparticles morphologies have been considered, including spherical particles (e.g. silica), platelets (e.g. clay and graphite) and nanotubes. In addition to the nanoinclusions themselves, the interphase, a special region of polymer chains in the vicinity of the nanofillers, also plays an important role in the improvement of polymer composite properties. The existence of nanofiller surfaces in the polymer alters the mobility of polymer chains surrounding them. Such perturbations in polymer molecular mobility extend several radii of gyration and create regions of polymer, the interphase, with properties and response different from that of the host bulk polymer. In conventional fiber or particulate composites, the interphase is only present an amount considerably smaller than the matrix and the reinforcements, thus the interphase is mainly considered for its contribution to the load transfer. However, due to the large surface to volume ratio of nanofillers (several orders of magnitude larger than conventional fillers), the amount of interphase polymer generated in nanocomposites can be substantial. Although direct measurements of interphase region in nanocomposites are lacking, recent work [12] has correlated thin film and nanocomposite data, providing quantitative evidence that local polymer properties are altered at substantial distances from nanoparticle surfaces and scale with the results of thin films confined by substrates of the same surface chemistry. Moreover, recent experimental work by Rittigstein et al. [13] observed T g -confinement effects hundreds of nanometers away from confining surfaces in doubly supported thin films. These results on model nanocomposites suggest similar extents for interphase regions could exist in polymer nanocomposites. Thus, in nanocomposites with a large extent of interphase and large number of inclusions, the interphase region has the potential to penetrate through the nanocomposite to form a percolating altered polymer network that may fundamentally change the response of the bulk polymer. The remarkable change of viscoelastic properties is one example that cannot be well explained without considering the contribution /$ - see front matter Ó 2008 Elsevier Ltd. All rights reserved. doi: /j.compscitech
2 492 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) of interphase material. Recent experimental studies on the nanocomposites show that the transition zones of storage/loss modulus and tand curves may be broadened or shifted in the time/temperature domain [8,14 16]. Since the properties of nanoinclusions are not time/temperature dependent within the experimental measurement scale, those enhancements cannot be simply attributed to the existence of nanoparticles. Moreover, some other observations indicate that the morphology of the interphase region should be taken into consideration to better explain the property changes in such composites [14,15,17,18]. For example, a remarkable shift in relaxation characteristics by 30 C relative to the pristine polymer matrix is observed in the functionalized single-wall nanotube composite, while, in contrast, the unmodified nanotube composite retains the same relaxation characteristics as the bulk polymer [15]. Khaled et al. [14] also reported that the strong interfacial interactions between functionalized T i O 2 nanofibers and a PMMA matrix leads to an increase of T g. In these cases, a percolated network of interphase polymer surrounding the nanoparticles may be significantly impacting the T g of the composites. Thus, a good understanding of the interphase effect on the thermomechanical properties is important for the accurate prediction of overall behavior. Micromechanical methods such as Halpin Tsai, self-consistent scheme, and Mori-Tanaka have been applied to multiphase polymeric materials (polymer blends and composites) to understand and predict the thermomechanical properties. Although these models were developed for elasticity, it is straightforward, using correspondence principle, to invoke an analogy between elasticity solutions and viscoelasticity solutions. For example, Li and Weng studied the overall viscoelastic properties of a two-phase composite containing randomly oriented ellipsoidal inclusions using Mori-Tanaka method [19]. Colombini et al used a self-consistent scheme to analyze the mechanical properties of polymer blends [20]. An unavoidable difficulty for most of micromechanical models is that in order to obtain closed-form analytical solutions, complex shapes of inclusions are prohibited. Moreover, the assumption of a continuous matrix phase is usually fundamental in those models. For nanocomposites, both of these features of homogenization schemes are problematic because of the interphase: the interphase region surrounding the inclusion complicates the inclusion morphology and large amounts of interphase could form an irregular, connective network rendering the micromechanics approach inadmissible. In this work, we use a finite element approach to study the impact of interphase on the viscoelastic properties as well as thermal response of polymeric nanocomposites. We focus on analysis of the influence of a connective interphase network in the matrix to provide a qualitative explanation of T g enhancement in the functionalized systems. The finite element method is employed in this study because it allows us to consider a very complex microstructure and is able to predict the macroscopic behaviors in a length/time scale close to the experimental tests. For example, the elastic properties of polymer nanocomposites have been predicted by finite element (FE) analysis and good agreement with experimental results have been reported [21,22]. The FE approach has been also adopted to study the viscoelastic properties in traditional polymer composites [23,24]. It is also noteworthy that, unlike some finite element simulations for polymer blends [25,26] in which each phase is discrete because the immiscibility of the component polymers leads to a heterogeneous multiphase system, a gradient of material properties is explicitly considered in this study. Our approach is designed to capture the evolution of the interphase in nanocomposites from discrete regions to a connective network and allows us to better understand its influence on the overall viscoelastic properties at different stages. This paper is organized as follows: The finite element model will be introduced in next section. Results and discussion will follow in Section 3. In Section 4, the conclusion will be presented. 2. Finite element model To examine the influence of the interphase on matrix dominated properties, a viscoelastic, 2D plain strain model was adopted. This configuration simplifies the computational requirements and, although it represents explicitly the case of aligned nanotubes, the results on time dependent properties are generalizable to other fiber arrangements. Since such composites can be idealized as transversely isotropic material, a representative volume element (RVE) with periodic structure is created to obtain the effective mechanical behavior of the composite. The results provide a qualitative understanding for the impact of interphase interactions on the viscoelastic response of nanocomposite Multiparticles unit cell generation A square unit cell containing a random dispersion of thirty nonoverlapping identical circles is generated to simulate the in situ configurations of aligned continuous nanotube embedded in polymer matrix. Each circle represents a nanotube in a transverse section of nanocomposite. To represent the continuum properties of the composite with many inclusions, periodic boundary conditions are applied to replicate an infinite but repetitive configuration so that the results are comparable to the experimental results. The periodic boundary conditions in two dimensions can be expressed as a function of the displacement u: uðx 1 ; 0ÞþU 2 ¼ uðx 1 ; LÞ uð0; X 2 ÞþU 1 ¼ uðl; X 2 Þ where L is the length of square edge. U 1 and U 2 depend on the particular loading applied on the cell. For example, if the cell is loaded uniformly along the X 1 direction with the deformation of d, the vectors U are written as U 1 =(d,0), U 2 = (0,u 2 ), where u 2 is computed from the traction free condition: Z T 2 dc ¼ 0 on X 2 ¼ L ð2þ C where T 2 stands for the normal traction acting on the boundary X 2 = L. In order to properly enforce the periodic boundary conditions, the geometry of the unit cell must be periodic. This condition has particular impact on the inclusions close to or intersecting the edge of the unit cell. Moreover, the centers of inclusions must be generated randomly so that the particle arrangement is statistically isotropic. Finally, the microstructure should be suitable for finite element discretization. To fulfill all the requirements, the random sequential adsorption (RSA) algorithm with additional conditions is adopted to generate the coordinates of the inclusion centers [27,28]. According this method, the inclusions are generated randomly and sequentially inside of the cell such that each inclusion is accepted if it does not overlap any of the already existing inclusions. Given the radius of inclusions as r, the new particle i is accepted if the following conditions are satisfied: The distance between the centroid of particle i and the centroids of all existing particles j =1,2,...,i 1 have to exceed a minimum value s, where s > 2r. The reason to impose this condition is not only to avoid distorted elements between circles, but also to save sufficient space for interphase layers. By controlling the magnitude of s, we can generate distributions of particles that range from uniform to random dispersions (see below for details). Mathematically, we can express this condition as: jx i x j þ lj P s ð1þ ð3þ
3 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) where x is the centroid of each particle and l =(l 1,l 2 ). l 1 and l 2 are zero except for particles near the boundaries, in which cases l 1 and/or l 2 take on the values L or L to account for minimum distances to particles near the opposite edges (because of periodic boundary conditions). The proximity of particle edges to the boundary of the unit cell is controlled to prevent the presence of distorted elements during meshing. Mathematically, this condition can be expressed as: jx i k rj > t; k ¼ 1; 2 jx i þ r Lj > t; k ¼ 1; 2 ð4þ k where t is a small value depended on the quality of mesh. The modified RSA algorithm, as expressed by the conditions of Eqs. (3) and (4), was used to generate centers of particles in the unit cell. Herein, we assigned a consistent radius, one of unity, to all of inclusions since we want to isolate the effects of interphase percolation in the composite. In this study, the lateral dimension of the cell was chosen as 70 times the radius of the inclusions so that the volume fraction of nanoparticles is about 1.9% for all samples which is among the fractions frequently used in nanocomposites. Two models with different dispersion were generated as shown in Fig. 1. For the first sample (referred to as the uniform model in the remainder of this paper), the value of s in Eq. (3) was set to be 10r, causing a large excluding volume of each particle, so the dispersion is very uniform. Meanwhile, the value of s was set to be equal to 2r for the second sample (called random in the following) which produces randomly dispersed particles. The interphase region in the polymeric composites was modeled as an annulus surrounding the particle with the thickness varying from half to several times diameter of the nanotube. Subsequently, several cells with different volume fractions of interphase were created. In this study, perfect bonding between matrix and inclusions was assumed, thus no interface contact models need be considered in the simulations. Fig. 2 shows two cells with uniformly dispersed particles surrounded by interphase layers: in the left plot, the thickness of interphase layers is set to be 1.5 times the nanotube diameter which leads to a volume fraction of interphase (V t ) as 28.8%, while, in the right plot, the thickness is 3 times the nanotube diameter and V t is 81.4%. Herein, all volume fractions are calculated based on the mesh so there is no error induced by finite element discretization. In the model construction, the overlapped annuli were merged to form a union and the overlapped interphase area was only counted once. From Fig. 2, we can clearly see a percolated network of the interphase region appears in the matrix as the V t exceeds the percolation threshold, which Fig. 1. Schematic of the unit cell for (a) uniformly (b) randomly dispersed nanocomposite where nanotubes are represented by black disks. The volume fractions of nanotubes in both configurations are 1.9%. See text for details. Fig. 2. Schematics of composite configurations with uniformly distributed inclusions (from Fig. 1a) and different volume fractions of interphase layers. The black interior inclusions represent nanotubes, while the blue annuli surrounding the interior nanotubes represent the interphase layers and rest of area is the bulk matrix. (a) Before interphase percolation, V t = 28.8%. (b) After interphase percolation, V t = 81.4%.
4 494 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) is 70.6% for this uniform dispersion (as shown in Fig. 1a). It should be noted that the percolation threshold in this paper refers to the minimum volume fraction of the interphase to form a connective network throughout the matrix under given locations of each particle. Particularly, this threshold is calculated numerically by detecting geometrical percolation in the composite given the dispersion showing Fig. 1a. Therefore, the percolation threshold defined here strongly depends on the geometry of the RVE and varies as the distribution of inclusions changes, differing from classic geometrical percolation theory. The percolation threshold for the random case shown in Fig. 1b is also calculated and is much lower, only 43.0% Material properties The frequency domain response of polycarbonate (PC) measured by dynamic mechanical analysis [29] was used as the material property of matrix in this model. The master curves of bulk PC at a reference temperature of 150 C was subsequently trimmed and approximated with a 28-term Prony Series using the DYNAFIT program previously developed in this lab [30]. The comparison between the experimental results and the Prony Series approximation is shown in Fig. 3a. In addition, the Williams Landel Ferry (WLF) equation was utilized to interpolate shift factors for temperatures greater than T g of PC in order to smooth the simulation results in the temperature domain (Fig. 3b, see text in section on Impact on glass transition temperature for details). To meet the requirement of viscoelastic material input of ABA- QUS, the complex Young s modulus of PC must be converted to complex shear and bulk moduli and the values of the instantaneous properties (the glassy equilibrium values at zero time) must be provided. The instantaneous Poisson s ratio of PC is 0.4 [31] and the instantaneous bulk modulus K * was calculated from this value and the instantaneous Young s modulus using elasticity conversion expressions. Since the variation of bulk modulus is relatively small in the frequency domain, while the changes in the Poisson ratio are known to be important, K * was assumed to be constant. Therefore, the complex shear modulus can be calculated at each frequency point via the correspondence principle [31] using the data for Young s modulus, E *, and the constant K * : G ¼ 3K E 9K E ð5þ The axial modulus of nanotubes is generally two orders of magnitude higher than the transverse modulus [32] and the absolute values of the moduli depend significantly on the type of nanotube and its processing method. However, our simulations (results not shown here) illustrate that neither the anisotropy nor the magnitude of the nanoparticle modulus has measurable impact on the viscoelastic properties of the nanocomposites. Thus, since our results exclusively focus on the viscoelastic properties of the nanocomposites, for simplification, the material properties of nanotubes were assigned to be linear, elastic, and isotropic with Young s modulus of 1 TPa and Poisson s ratio of 0.3. Although the properties of nanoparticles and bulk polymer are readily available, it is more of challenge to determine the interphase properties. First, it is not yet possible to measure those properties using current experimental techniques without biasing the results with significant assumptions. On the other hand, the numerical calculation of molecular dynamics (MD) can capture some structural and dynamical details of interphase only at the molecular time and length scale, while the continuum modeling requires material properties in many orders of magnitudes longer. Nevertheless, both numerical and experimental work provides guidance for the trend in continuum properties of interphase compared to the host bulk polymer. It has been noted that the attractive interactions between nanoinclusions and matrix can restrict mobility of polymer chains [33 38]. For example, Wei et al. [36] performed MD study that demonstrated an increase of 20 C to 60 C int g for a polyethylene cell containing nanotubes. Other MD simulations also showed T g to increase near the nanoparticles for attractive interaction [33 35]. Similarly, at the other end of the spectrum, experimental observations show that the T g of a nanocomposite can be on the order of ten degrees higher than the pure matrix materials [15,39]. Therefore, based on such simulation and experimental results, we hypothesize that the interphase properties are related to those of the bulk polymer matrix by a shift in the relaxation times. This approximation captures one of the most important characteristics of interphase, that of altered polymer mobility, by introducing the simplest possible assumption on the properties. In this paper, the interphase properties were determined by shifting the PC master curve two decades lower in the frequency domain, representative of strong positive interactions between the nanoparticle and the host matrix material. 3. Results and discussion Finite element models were created with different interphase volume fractions and the calculation was performed using the commercial FEA package, ABAQUS. The impact of interphase volume and structure were analyzed separately in this study Frequency response Based on the particle distribution from the uniform case (Fig. 1a), seven models were created with different thicknesses Fig. 3. (a) Experimentally measured complex modulus of PC and its Prony Series approximation. (b) The shift factor versus temperature where the red line is the data fit by WLF equation (original data is from reference [23]).
5 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) of interphase layers resulting in volume fractions of V t = 5.7%, 28.8%, 46.1%, 57.7%, 65.2%, 81.4%, and 87.4% (the corresponding interphase thicknesses are 0.5d, 1.5d, 2d, 2.25d, 2.5d, 3d, and 3.25d, where d is the diameter of inclusion). The comparisons of the transverse storage and loss moduli for those models are shown in Fig. 4. The results clearly show that the transition zones of nanocomposites shift towards lower frequency as compared with pure PC. The loss modulus peak, shown in Fig. 4b, is broadened compared to pure matrix, and the broadening becomes more apparent as the volume fraction of interphase increases. In addition, the loss peaks are shifted to lower frequency corresponding to the transition zone of interphase. Note that the relaxation characteristics of the nanocomposite transitions from being dominated by the bulk PC to the interphase when the volume fraction of interphase exceeds the percolation threshold (70.6% in this configuration). This result is more evident by comparing the tand curves of samples with different interphase volume fraction (Fig. 4c). For comparison, we can also consider the configuration of the random sample (Fig. 1b). In this case, we can expect a lower percolation threshold of interphase compared to the uniform dispersion. Fig. 5 shows the tand results for these simulations with several different volume fractions of interphase. Consistent with the geometric percolation threshold of 43.0% in this configuration, we also observe that the simulation results demonstrate a transition from PC to interphase dominated behavior once this percolation threshold is reached. Taken together, these results indicate that not only is the amount of interphase important, but critically the connectivity of the interphase (its percolation) has significant impact on the overall viscoelastic properties of nanocomposites. In order to further verify the impact of the interphase percolation, a comparison between the finite element simulation and the simple rule of mixture (ROM) results has been conducted. For the ROM, the transverse formulation is used as appropriate for the transverse modulus of unidirectional composites, which is: 1 E c ¼ m f E f þ m i E i þ m m E m ð6þ where subscripts, c, f, i, and m represent the composite, nanotube, interphase, and matrix, respectively. Since both moduli of interphase and polymer matrix are complex values, the modulus of composite calculated by Eq. (6) is also a complex number. Fig. 6 shows the comparison of tand in frequency space that is computed by the FE models and ROM method respectively. Note that the FE models used here are exactly same as those shown before, thus the percolation thresholds are 70.6% and 43.0% with respect to uniform case and random case. Fig. 6a shows that in the case with low interphase volume fraction no percolation occurs in either model, from FE or ROM predictions. Fig. 6b displays the high volume case in which the interphase is percolated in both of the models. Since the transverse ROM treats the composite in a sandwich mode, the overall properties are dominated by the weakest phase. Thus, because interphase modulus exceeds that of pure PC, the ROM predictions for tand of the composite are still dominated by the characteristics of the pure PC and do not capture the shift of tand peak. This comparison to the ROM provides additional evidence that the percolation of the interphase significantly affects the properties of polymeric nanocomposites, in such a way that a simple ROM prediction can no longer capture the essence Impact on glass transition temperature (T g ) As the vast majority of thermomechanical experimental data on nanocomposites is in the temperature domain, we also use the FEM approach to predict the temperature dependent response of the nanocomposites with different interphase fractions. The finite Fig. 4. Predicted moduli of nanocomposites of uniformly distributed inclusions (Fig. 1a) with different volume fraction of interphase (V t ). Properties for bulk PC and the interphase are also shown. (a) Comparison of storage moduli. (b) Comparison of loss moduli. (c) Comparison of tand. element model and material properties used are identical to the uniform dispersion case described earlier. The material behavior in temperature domain was obtained by calculating the complex modulus of nanocomposites at a fixed frequency, in this case 1 Hz dynamic loading, at different temperatures. The material properties for each phase at each temperature were obtained by shifting the master curves at the reference temperature by the appropriate shift factors. The shift factors were obtained in the
6 496 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) Fig. 5. Tand in frequency domain of nanocomposites with randomly distributed inclusions (Fig. 1b). Fig. 7. Prediction of glass transition temperature of nanocomposites with uniform nanotube distribution (Fig. 1a). The data at each temperature point is calculated at frequency of 1 Hz and the tand curves have been normalized. time temperature superposition (TTSP) process for frequency domain response master curves of the pure PC material (Fig. 3) and the WLF equation was used to fit the experimental curve of shift factors for T > T g of PC. The temperature domain responses for two nanocomposites with different volume fraction of interphase, V t = 15% and 81.4%, were plotted in Fig. 7. The glass transition temperature (T g ) of the nanocomposite with larger interphase, indicated by the tand peak, clearly shifts to higher temperature compared to the pure PC and is dominated by the interphase properties. This result is again consistent with the concept of a percolated interphase dominating the effective viscoelastic properties of the composite. Experimental data for composites with surface modified nanotubes indicates that the larger region of reducedmobility polymer (i.e. interphase) in those systems leads to a higher T g than unfunctionalized nanocomposites [15,39]. The broadening of the transition zone is also observed for the composite with 15% interphase, but is less pronounced, which could be attributed to the low volume fraction of interphase in that model. Note that a broadened tand peak is not always observed [15] and further work will focus on the origin of the broadening in the temperature domain Gradient of Interphase properties One concern that arises from the previous results compared to experimental data is the dual-peak appearing on the tand curves for nanocomposites with moderate volume fraction of interphase (e.g. 50% interphase for the uniform distribution case in Fig. 4c). One possible reason for such a sharp distribution is that only two relaxation modes were used in our models, one for bulk PC and the other for an interphase with two decades less mobility. In the geometric configuration of the numerical model, the interphase properties were assigned to a single discrete domain around each nanoparticle with a sudden discontinuous change to the bulk matrix properties at its boundary. Considering experimental data on ultrathin polymer films on an attractive surfaces [37], the polymer mobility changes in a gradient fashion away from the polymer surface interface. Thus, to determine if the gradient in properties could lead to a change in response which would diminish the artificial dual-peak at certain volume fractions, we considered the simplest possible gradient: two concentric layers of interphase surround the nanotubes with modulated properties (Fig. 8). Maintaining the basic geometry of previous models, the interphase annulus was divided Fig. 6. Comparison between the FE simulations and rule of mixtures (ROM) results. (a) Shows the case for a low interphase volume fraction which is below the percolation thresholds for both uniform and random distributions. (b) Shows the high volume fraction case in which the interphase is percolated in both configurations.
7 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) Fig. 8. Schematic of unit cell for a uniform nanotube distribution (Fig. 1a) with double-layer of interphase surrounding the nanotubes. The ratios of thickness (inner:outer) are 1:1 for all models. The total interphase volume fractions (V t ) are 46.1%, 65.2%, and 81.6%, respectively. into two concentric annuli of the same thickness, such that the sum of volume fractions of the two interphase layers remains the same as the interphase volume fraction of the corresponding single-layer interphase models showed earlier. The properties of the first interphase, that closest to the nanotube and denoted as Interphase-1, remained the same as before (shifted 2 decades toward lower frequency). The properties of the second interphase, closest to the PC and denoted as Interphase-2, were also derived from bulk PC but shifted only one decade in frequency, thereby simulating the decay of the properties toward that of the matrix material. Three examples are studied here for the uniformly distributed inclusion configuration (Fig. 1a). In the first case (Fig. 8a), the inner layer and outer layer thicknesses are identical to the diameter of the nanotube, and the total volume fraction of two layers is 46.1%, the same as the V t in one sample studied earlier. In the second sample (Fig. 8b), the thicknesses of inner layer and outer layer are 1.25 times the nanotube diameter, resulting in a total volume fraction for two interphase layers of 65.2%. In the third model, both thicknesses of the two interphase layers are equal to 1.5 times the nanotube diameter for a total volume fraction of 81.6%. The dynamic responses of new models were calculated and tand curves of these new models compared to those of the previous singlelayer models are plotted in Fig. 9. The results clearly show that, for the double-layer cases, the tand curves are much smoother and double peaks are substituted by a single peak which moves towards to the bulk PC region. These results are important, as they agree fundamentally with experimental data on nanocomposites in which gradual shifting of the tand peaks are observed, but never double peaks. The implication then is that the gradient interphase ultimately leads to a unique relaxation characteristic for the nanocomposites. Comparing the double-layer model tand curves with those of previous single-layer models, we see that the peaks for the gradient interphase shift to higher frequency. This phenomenon is more pronounced for the samples where the interphase dominates the overall viscoelastic properties (i.e. Fig. 9c). This result occurs because the outer layer interphase, which plays important role in the double-layer model, has a transition zone occurring at higher frequency than the inner layer interphase. The inner layer interphase, with a lower frequency loss peak, is used only for half of the total interphase. In contrast, the single-layer model uses the inner layer interphase properties for the full extent of the interphase zone. It is also observed that the magnitudes of the tand peaks differ for the single-layer and double-layer interphase models. Since the properties of interphase were obtained by shifting the master curve of bulk PC horizontally to lower frequency, the mechanical damping capabilities are identical for all those materials. Therefore, we should expect the same areas under the tand curves for both single-layer and double-layer models, as shown in Fig. 9. For the cases where two maxima appeared on the tand curve of the single-layer model (a and b in Fig. 9), the main peak is necessarily lower than the peak of the corresponding double-layer model to maintain the same damping effect. In the case of Fig. 9c, both single- and double-layer models are fully percolated with only one peak, suppressing the relaxation signature of bulk PC. Thus in this case, both models have the same intensity. Comparing the movement of the tand for each of the double-layer cases, we observe that the tand peaks gradually shift to lower frequency as the volume fraction of interphase increases. In particular, when a percolated interphase network forms inside of the composite (Fig. 8c), the predicted tand peak moves to a position between two interphase peaks (Fig. 9c), indicating that the interphase dominates the overall relaxation properties of the composite. 4. Conclusion The influence of the interphase and its structure on the viscoelastic response of polymer nanocomposites was studied via a two-dimensional finite element analysis. A unit cell with thirty dispersed particles was created to represent the transverse section of unidirectional nanotube reinforced polymer composites. 1 Each particle was coated by a layer of interphase, whose properties were obtained by shifting the bulk matrix properties in the frequency domain, corresponding to altered mobility of the polymer in the vicinity of nanotubes. Given the experimentally measured frequency domain response of bulk polymer, the viscoelastic properties of the nanocomposites in both frequency and temperature domains were calculated. Furthermore, varying the thickness of interphase layers and modifying the dispersion of nanotubes caused changes in the interphase connectivity in the polymer. The simulation results demonstrated that relaxation characteristics of nanocomposites are greatly dependent upon the volume of interphase and formation of an interphase network. Quantitative comparison between the experimental data and the model prediction is not possible at current stage, due to many uncertainties, including nanotube and interphase properties, in situ nanotube morphology, and the 2D nature of the simulation. Nevertheless, 1 We note that unit cells with different numbers of particles were also studied, all resulted in the same conclusion.
8 498 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) For the case of low volume fraction of interphase such that the interphase domains remain discrete, the predicted transition zones of the nanocomposites are broadened but no significant shift occurs in frequency or temperature domains compared to the properties of the bulk polymer. A discrete interphase situation is consistent with relatively thin interphase layers, very low volume fractions of nanotubes, or the nanotubes aggregated to clusters. Results correspond to data on unmodified surface of MWNT nanocomposites and unfunctionalized SWNT nanocomposites. For the case of high volume fraction of interphase such that the interphase is percolated through the composite, a clear shift of the predicted relaxation peaks of the nanocomposites occur compared to the peak of the pure polymer. A percolated microstructure corresponds to the results of functionalized SWNT nanocomposites and MWNT nanocomposites with surface modification which have better dispersion of nanotubes and a thicker region of reduced-mobility polymer around the nanotubes. It is important to use a gradient in the properties of the interphase to avoid an artificial dual-peak phenomenon. While an artificial dual-peak was observed at some intermediate volume fractions for the single-layer interphase cases, implementation of a simple double-layer interphase with modulated properties completely removed this anomaly. The results from this work provide powerful insight into the importance of percolating interphase in nanocomposites for viscoelastic properties. However, a significant limitation of the work is the 2D nature of the models. The two-dimensionality does not permit us to study the influence of aspect ratio or shape of the nanoparticles, which is known to be important. The percolation thresholds for 2D models are also significantly higher than for 3D network structures and the resulting connectivity of the interphase domains is more limited. Thus, a 3D model is necessary to better study the properties of nanocomposite with different inclusions. Moreover, a quantitative comparison between the numerical and experimental results would be enhanced by the extension to a 3D system. Acknowledgment The authors wish to thank the support of the National Science Foundation NIRT program under Grant No References Fig. 9. Comparison of tand in frequency domain for single-layer model and doublelayer model. The interphase volume fractions (V t ) of these models are: (a) 46.1%, (b) 65.2%, and (c) 81.4%, respectively, refer to text for details. the predicted pattern of the impact of the interphase on the overall performance of the nanocomposites is consistent with experimental observations [15,29]. Significant observations are: [1] Biercuk MJ et al. Carbon nanotube composites for thermal management. Appl Phys Lett 2002;80(15): [2] Cadek M et al. Morphological and mechanical properties of carbon-nanotubereinforced semicrystalline and amorphous polymer composites. Appl Phys Lett 2002;81(27): [3] Coleman JN et al. High-performance nanotube-reinforced plastics: understanding the mechanism of strength increase. Adv Funct Mater 2004;14(8): [4] Coleman JN, Khan U, Gun ko YK. Mechanical reinforcement of polymers using carbon nanotubes. Adv Mater 2006;18(6): [5] Favier V et al. Mechanical percolation in cellulose whisker nanocomposites. Polym Eng Sci 1997;37(10): [6] Hajji P et al. Tensile behavior of nanocomposites from latex and cellulose whiskers. Polym Compos 1996;17(4): [7] Messersmith PB, Giannelis EP. Synthesis and barrier properties of poly(epsiloncaprolactone)-layered silicate nanocomposites. J Polym Sci Part A: Polym Chem 1995;33(7): [8] Putz KW et al. Elastic modulus of single-walled carbon nanotube/poly(methyl methacrylate) nanocomposites. J Polym Sci Part B: Polym Phys 2004;42(12): [9] Rong MZ et al. Improvement of tensile properties of nano-sio 2 /PP composites in relation to percolation mechanism. Polymer 2001;42(7):
9 R. Qiao, L. Catherine Brinson / Composites Science and Technology 69 (2009) [10] Sandler J et al. Development of a dispersion process for carbon nanotubes in an epoxy matrix and the resulting electrical properties. Polymer 1999;40(21): [11] Sandler JKW et al. Ultra-low electrical percolation threshold in carbonnanotube-epoxy composites. Polymer 2003;44(19): [12] Bansal A et al. Quantitative equivalence between polymer nanocomposites and thin polymer films. Nat Mater 2005;4(9): [13] Rittigstein P et al. Model polymer nanocomposites provide an understanding of confinement effects in real nanocomposites. Nat Mater 2007;6(4): [14] Khaled SM et al. Synthesis of TiO 2 PMMA nanocomposite: using methacrylic acid as a coupling agent. Langmuir 2007;23(7): [15] Ramanathan T, Liu H, Brinson LC. Functionalized SWNT/polymer nanocomposites for dramatic property improvement. J Polym Sci Part B: Polym Phys 2005;43(17): [16] Ramanathan T et al. Graphitic nanofillers in PMMA nanocomposites an investigation of particle size and dispersion and their influence on nanocomposite properties. J Polym Sci Part B: Polym Phys 2007;45(15): [17] Zhang H et al. Property improvements of in situ epoxy nanocomposites with reduced interparticle distance at high nanosilica content. Acta Mater 2006;54(7): [18] Schadler L, Brinson L, Sawyer W. Polymer nanocomposites: a small part of the story. JOM J Miner Metals Mater Soc 2007;59(3): [19] Li J, Weng GJ. Strain-rate sensitivity, relaxation behavior, and complex moduli of a class of isotropic viscoelastic composites. J Eng Mater Technol Trans ASME 1994;116(4): [20] Colombini D, Merle G, Alberola ND. Use of mechanical modeling to study multiphase polymeric materials. Macromolecules 2001;34(17): [21] Selmi A et al. Prediction of the elastic properties of single walled carbon nanotube reinforced polymers: a comparative study of several micromechanical models. Compos Sci Technol 2007;67(10):2071. [22] Linjie Zhu KAN. Numerical simulation of the tensile modulus of nanoclay-filled polymer composites. J Polym Sci Part B: Polym Phys 2004;42(12): [23] Liu H, Brinson LC. A hybrid numerical-analytical method for modeling the viscoelastic properties of polymer nanocomposites. J Appl Mech Trans ASME 2006;73(5): [24] Fisher FT, Brinson LC. Viscoelastic interphases in polymer-matrix composites: theoretical models and finite-element analysis. Compos Sci Technol 2001;61(5): [25] Read DJ et al. Theoretical and finite-element investigation of the mechanical response of spinodal structures. Eur Phys J E Soft Matter 2002;8(1): [26] Takaaki Matsuoka SY. Computer simulation of phase structure and mechanical properties of polymer mixtures. J Appl Polym Sci 1998;68(5): [27] Rintoul MD, Torquato S. Reconstruction of the structure of dispersions. J Colloids Interf Sci 1997;186(2): [28] Segurado J, Llorca J. A numerical approximation to the elastic properties of sphere-reinforced composites. J Mech Phys Solids 2002;50(10): [29] Fisher FT et al. Spectral response and effective viscoelastic properties of MWNT-reinforced polycarbonate. Adv Compos Lett 2004;13(2): [30] Bradshaw RD, Brinson LC. A sign control method for fitting and interconverting material functions for linearly viscoelastic solids. Mech Time-Depend Mater 1997;1:85. [31] Brinson Hal F, Brinson LC. Polymer engineering science and viscoelasticity: an introduction. New York, USA: Springer; [32] Shen L, Li J. Transversely isotropic elastic properties of single-walled carbon nanotubes. Phys Rev B 2004;69(4): [33] Smith GD et al. A molecular dynamics simulation study of the viscoelastic properties of polymer nanocomposites. J Chem Phys 2002;117(20): [34] Smith JS, Bedrov D, Smith GD. A molecular dynamics simulation study of nanoparticle interactions in a model polymer-nanoparticle composite. Compos Sci Technol 2003;63(11): [35] Starr FW, Schroder TB, Glotzer SC. Molecular dynamics simulation of a polymer melt with a nanoscopic particle. Macromolecules 2002;35(11): [36] Wei CY, Srivastava D, Cho KJ. Thermal expansion and diffusion coefficients of carbon nanotube-polymer composites. Nano Letters 2002;2(6): [37] Ellison CJ, Torkelson JM. The distribution of glass-transition temperatures in nanoscopically confined glass formers. Nat Mater 2003;2(10): [38] Park CH et al. Thickness and composition dependence of the glass transition temperature in thin random copolymer films. Polymer 2004;45(13): [39] Eitan A et al. Reinforcement mechanisms in MWCNT-filled polycarbonate. Compos Sci Technol 2006;66(9):
MACROSCALE EXPERIMENTAL EVIDENCE OF A REDUCED- MOBILITY NON-BULK POLYMER PHASE IN NANOTUBE- REINFORCED POLYMERS
MACROSCALE EXPERIMETAL EVIDECE OF A REDUCED- MOBILITY O-BULK POLYMER PHASE I AOTUBE- REIFORCED POLYMERS F.T. Fisher and L.C. Brinson Department of Mechanical Engineering orthwestern University Evanston,
More informationModule 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method
Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without
More informationEffect of different crosslink densities on the thermomechanical properties of polymer nanocomposites
Effect of different crosslink densities on the thermomechanical properties of polymer nanocomposites *Byungjo Kim 1), Joonmyung Choi 2), Suyoung Yu 3), Seunghwa Yang 4) and Maenghyo Cho 5) 1), 2), 3),
More informationDIGIMAT for NANO-COMPOSITES
DIGIMAT for NANO-COMPOSITES Document Version 2.0.1, February 2010 info@e-xstream.com www.e-xstream.com Materials: Engineering Plastics, Reinforced Plastics, Mineral and Clay fillers. e-xstream Technology:
More informationTwo-step multiscale homogenization of polymer nanocomposites for large size RVEs embodying many nanoparticles
Two-step multiscale homogenization of polymer nanocomposites for large size RVEs embodying many nanoparticles *Kyungmin Baek 1), Hyunseong Shin 2), Jin-Gyu Han 3) and Maenghyo Cho 4) 1), 2), 3), 4) Department
More informationCHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES
CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES 4. Introduction Fillers added to polymer matrices as additives are generally intended for decreasing the cost (by increase in bulk) of
More informationMechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach
Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Sandhyarani Biswas Department of Mechanical Engineering, N.I.T Rourkela, INDIA Abstract: Now-a-days,
More informationNumerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites
Computational Materials Science 39 (2007) 198 204 www.elsevier.com/locate/commatsci Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites S. Kari
More informationPolymer Nanocomposites: Introduction
Chapter 1 Polymer Nanocomposites: Introduction Downloaded via 148.251.232.83 on April 5, 2019 at 22:23:00 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published
More informationEVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD
U.P.B. Sci. Bull., Series D, Vol. 79, Iss. 3, 2017 ISSN 1454-2358 EVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD Mohamed said BOUTAANI 1, Salah
More informationThe effect of interfacial bonding on the damping behavior of nanocomposites
The effect of interfacial bonding on the damping behavior of nanocomposites Literature review The rapid and continuous growth in aerospace, automotive, and military applications requires special materials
More informationMechanical Properties of Fiber Reinforced Composites Using Buckminster Fullerene Reinforcement
IJRMET Vo l. 4, Is s u e Sp l - 1, No v 2013- Ap r i l 2014 ISSN : 2249-5762 (Online ISSN : 2249-5770 (Print Mechanical Properties of Fiber Reinforced Composites Using Buckminster Fullerene Reinforcement
More informationMicromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading
Indian Journal of Engineering & Materials Sciences Vol. 15, October 2008, pp. 382-390 Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading K Sivaji Babu a *, K Mohana
More informationModule 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains
Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model
More informationAnalysis of high loss viscoelastic composites
Analysis of high loss viscoelastic composites by C. P. Chen, Ph.D. and R. S. Lakes, Ph.D. Department of Engineering Physics Engineering Mechanics Program; Biomedical Engineering Department Materials Science
More informationThe stress transfer efficiency of a single-walled carbon nanotube in epoxy matrix
JOURNAL OF MATERIALS SCIENCE 39 (2 004)4481 4486 The stress transfer efficiency of a single-walled carbon nanotube in epoxy matrix K. Q. XIAO, L. C. ZHANG School of Aerospace, Mechanical and Mechatronic
More informationModeling percolation in high-aspect-ratio fiber systems. I. Soft-core versus hard-core models
PHYSICAL REVIEW E 75, 041120 2007 Modeling percolation in high-aspect-ratio fiber systems. I. Soft-core versus hard-core models L. Berhan 1, * and A. M. Sastry 2 1 Department of Mechanical, Industrial
More informationContinuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites
International OPEN ACCESS Journal Of Modern Engineering Research (IJMER) Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites Sonali Gholap 1, Dr. Dhananjay R. Panchagade
More informationBuckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate
Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate Outline Introduction Representative Volume Element (RVE) Periodic Boundary Conditions on RVE Homogenization Method Analytical
More informationCHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles
CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles Mark Hughes 14 th March 2017 Today s learning outcomes To understand the role of reinforcement, matrix and
More informationAM11: Diagnostics for Measuring and Modelling Dispersion in Nanoparticulate Reinforced Polymers. Polymers: Multiscale Properties.
AM11: Diagnostics for Measuring and Modelling Dispersion in Nanoparticulate Reinforced Polymers Polymers: Multiscale Properties 8 November 2007 Aims Provide diagnostic tools for quantitative measurement
More informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP014259 TITLE: Mechanical Behavior of Polymer Nanocomposites: A Discrete Simulation Approach DISTRIBUTION: Approved for public
More informationA Hybrid Numerical-Analytical Method for Modeling the Viscoelastic Properties of Polymer Nanocomposites
Hua Liu Department of Mechanical Engineering, orthwestern University, Evanston, IL 60208 L. Catherine Brinson 1 Department of Materials Science and Engineering, orthwestern University, Evanston, IL 60208
More informationPredicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics
Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Neelam Mishra 1, and Kaushik Das *1 1 School of Minerals Metallurgical and Materials Engineering, Indian
More informationCHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement
CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement Mark Hughes 19 th September 2016 Outline Fibre architecture Volume fraction and the rule of mixtures Principle
More informationMODELLING INTERACTION EFFECT OF NANOSILICA PARTICLES ON NANOSILICA/EPOXY COMPOSITE STIFFNESS
MODELLING INTERACTION EFFECT OF NANOSILICA PARTICLES ON NANOSILICA/EPOXY COMPOSITE STIFFNESS Mulyadi a*, I. Gitman a, C. Pinna a, C. Soutis b a Department of Mechanical Engineering, The University of Sheffield,
More informationModeling percolation in high-aspect-ratio fiber systems. II. The effect of waviness on the percolation onset
PHYSICAL REVIEW E 75, 041121 2007 Modeling percolation in high-aspect-ratio fiber systems. II. The effect of waviness on the percolation onset L. Berhan 1, * and A. M. Sastry 2 1 Department of Mechanical,
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationDevelopment of a code to generate randomly distributed short fiber composites to estimate mechanical properties using FEM
International Journal of Theoretical and Applied Mechanics. ISSN 0973-6085 Volume 12, Number 4 (2017) pp. 863-872 Research India Publications http://www.ripublication.com Development of a code to generate
More informationFinite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements
Proceedings o the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. Finite Element Modeling o Residual Thermal Stresses in Fiber-Reinorced Composites Using Dierent Representative
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
More informationAspect Ratio Requirements for Nanotube-Reinforced, Polymer-Matrix Composites
Aspect Ratio Requirements for Nanotube-Reinforced, Polymer-Matrix Composites J.A.Nairn Wood Science and Engineering, Oregon State University, Corvallis, OR 97330, USA Abstract A fiber s efficiency in a
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationModification of Electrically Conductive Adhesives for Better Mechanical and Electrical Properties
D3 Advanced packaging, new packages & materials, optoelectronics, optical interconnection techn. Modification of Electrically Conductive Adhesives for Better Mechanical and Electrical Properties Pavel
More informationApparent thermal conductivity of periodic two-dimensional composites
Computational Materials Science 25 (2002) 329 338 www.elsevier.com/locate/commatsci Apparent thermal conductivity of periodic two-dimensional composites M. Jiang a, I. Jasiuk b, *, M. Ostoja-Starzewski
More informationSupporting Information. Controlled Structure Evolution of Graphene Networks in Polymer Composites
Supporting Information Controlled Structure Evolution of Graphene Networks in Polymer Composites Stephen C. Boothroyd,,# David W. Johnson,,# Michael P. Weir, Carl D. Reynolds, James M. Hart, Andrew J.
More informationEFFECT OF VACANCY DEFECTS ON THE MECHANICAL PROPERTIES OF CARBON NANOTUBE REINFORCED POLYPROPYLENE
International Journal of Mechanical Engineering and Technology (IJMET) Volume 8, Issue 7, July 2017, pp. 1370 1375, Article ID: IJMET_08_07_148 Available online at http://www.iaeme.com/ijmet/issues.asp?jtype=ijmet&vtype=8&itype=7
More informationSize Effects In the Crushing of Honeycomb Structures
45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 19-22 April 2004, Palm Springs, California AIAA 2004-1640 Size Effects In the Crushing of Honeycomb Structures Erik C.
More informationEffect of fibre shape on transverse thermal conductivity of unidirectional composites
Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and
More informationInternational Journal of Engineering Science
International Journal of Engineering Science 63 (2013) 1 9 Contents lists available at SciVerse ScienceDirect International Journal of Engineering Science journal homepage: www.elsevier.com/locate/ijengsci
More informationFunctionalized Carbon Nanotubes a key to nanotechnology?
1 27th Max Born Symposium Multiscale Modeling of Real Materials Wroclaw, Sep 19, 2010 Functionalized Carbon Nanotubes a key to nanotechnology? Karolina Milowska, Magda Birowska & Jacek A. Majewski Faculty
More informationA Continuum Model for Carbon Nanotube-Infused Polyimides
A Continuum Model for Carbon Nanotube-Infused Polyimides Heather Wilson, Sumanth Banda 2, Ralph C. Smith 3, and Zoubeida Ounaies 4 Department of Mathematics,3 Center for Research in Scientific Computation
More informationMECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy
MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS Vinu Unnithan and J. N. Reddy US-South American Workshop: Mechanics and Advanced Materials Research and Education
More informationINFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS
Djoković, J. M., et.al.: Influence of Temperature on Behavior of the Interfacial THERMAL SCIENCE: Year 010, Vol. 14, Suppl., pp. S59-S68 S59 INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK
More informationMechanics of wafer bonding: Effect of clamping
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 1 1 JANUARY 2004 Mechanics of wafer bonding: Effect of clamping K. T. Turner a) Massachusetts Institute of Technology, Cambridge, Massachusetts 0219 M. D. Thouless
More informationMulti-scale mechanics and structure of semi-hard cheese
Multi-scale mechanics and structure of semi-hard cheese T.J. Faber a,b, P.J.Schreurs b, J.M.J.G. Luyten a, H.E.H.Meijer b a FrieslandCampina Research, Deventer, The Netherlands (timo.faber@frieslandcampina.com))
More informationModeling and Simulation of Artificial Core-Shell Based Nanodielectrics for Electrostatic Capacitors Applications
Modeling and Simulation of Artificial Core-Shell Based Nanodielectrics for Electrostatic Capacitors Applications -D. Musuwathi Ekanath 1, N. Badi* 1,2, A. Bensaoula 2 1 Center for Advanced Materials, University
More informationMICROMECHANICAL DEFORMATIONS IN PARTICULATE FILLED POLYMERS: THE EFFECT OF ADHESION
MICROMECHANICAL DEFORMATIONS IN PARTICULATE FILLED POLYMERS: THE EFFECT OF ADHESION K. Renner, J. Móczó, B. Pukánszky Laboratory of Plastics and Rubber Technology, Department of Physical Chemistry and
More informationAuthors: Correspondence: ABSTRACT:
Multi-Scale Modeling of the Impact and Failure of Fiber Reinforced Polymer Structures using DIGIMAT to LS-DYNA Interface Authors: L. Adam, R. Assaker & R. Ramaya e-xstream engineering S.A. 7, Rue du Bosquet.
More information3D Compression Molding
Autodesk Simulation Moldflow Insight 2014 3D Compression Molding Executive summary In this work, the simulation results from a program developed for the three-dimensional analysis of compression molding
More informationMechanical and Thermoviscoelastic Behavior of Clay/Epoxy Nanocomposites
Mat. Res. Soc. Symp. Proc. Vol. 7 23 Materials Research Society I6.5.1 Mechanical and Thermoviscoelastic Behavior of Clay/Epoxy Nanocomposites Jandro L. Abot, Asma Yasmin and Isaac M. Daniel Robert McCormick
More informationCHAPTER 4: VISCOELASTIC BEHAVIOR OF CARBON NANOTUBE-REINFORCED POLYMERS
CHAPTER 4: VISCOELASTIC BEHAVIOR OF CARBON NANOTUBE-REINFORCED POLYMERS Recent experimental results demonstrate that substantial improvements in the elastic properties of a polymer can be attained by using
More informationA micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds
International Journal of Solids and Structures 42 (2005) 4179 4191 www.elsevier.com/locate/ijsolstr A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial
More informationEffect of Interlayers on Mechanical Properties and Interfacial Stress Transfer of 2D Layered Graphene- Polymer Nanocompsites
University of Kentucky UKnowledge Theses and Dissertations--Mechanical Engineering Mechanical Engineering 2017 Effect of Interlayers on Mechanical Properties and Interfacial Stress Transfer of 2D Layered
More informationMechanical Behavior of Fullerene Reinforced Fiber Composites with Interface Defects through Homogenization Approach and Finite Element Method
, pp.67-82 http://dx.doi.org/1.14257/ijast.215.78.6 Mechanical Behavior of Fullerene Reinforced Fiber Composites with Interface Defects through Homogenization Approach and Finite Element Method P. Prasanthi
More informationA Molecular Modeling Approach to Predicting Thermo-Mechanical Properties of Thermosetting Polymers
A Molecular Modeling Approach to Predicting Thermo-Mechanical Properties of Thermosetting Polymers Natalia Shenogina, Wright State University Mesfin Tsige, University of Akron Soumya Patnaik, AFRL Sharmila
More informationConstitutive Modeling of Nanotube-Reinforced Polymer Composite Systems. G. M. Odegard, V. M. Harik, K. E. Wise, and T. S. Gates
Constitutive Modeling of Nanotube-Reinforced Polymer Composite Systems G. M. Odegard, V. M. Harik, K. E. Wise, and T. S. Gates ABSTRACT In this study, a technique has been proposed for developing constitutive
More informationNITRILE RUBBER (NBR) NANOCOMPOSITES BASED ON DIFFERENT FILLER GEOMETRIES (Nanocalcium carbonate, Carbon nanotube and Nanoclay)
CHAPTER 5 NITRILE RUBBER (NBR) NANOCOMPOSITES BASED ON DIFFERENT FILLER GEOMETRIES (Nanocalcium carbonate, Carbon nanotube and Nanoclay) 5.1 Introduction Nanocalcium carbonate (NCC) is a particulate nanofiller
More informationMICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING
MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING N. Krishna Vihari 1, P. Phani Prasanthi 1, V. Bala Krishna Murthy 2* and A. Srihari Prasad 3 1 Mech. Engg. Dept., P. V. P. Siddhartha
More informationThe role of interparticle contact in conductive properties of random particulate materials
Available online at www.sciencedirect.com Acta Materialia 56 (2008) 2810 2818 www.elsevier.com/locate/actamat The role of interparticle contact in conductive properties of random particulate materials
More informationEFFECT OF INTERPHASE CHARACTERISTIC AND PROPERTY ON AXIAL MODULUS OF CARBON NANOTUBE BASED COMPOSITES
15 EFFECT OF INTERPHASE CHARACTERISTIC AND PROPERTY ON AXIAL MODULUS OF CARBON NANOTUBE BASED COMPOSITES Md. Abdulla Al Masud and A.K.M. Masud * Department of Industrial and Production Engineering, Bangladesh
More informationIntensity (a.u.) Intensity (a.u.) Raman Shift (cm -1 ) Oxygen plasma. 6 cm. 9 cm. 1mm. Single-layer graphene sheet. 10mm. 14 cm
Intensity (a.u.) Intensity (a.u.) a Oxygen plasma b 6 cm 1mm 10mm Single-layer graphene sheet 14 cm 9 cm Flipped Si/SiO 2 Patterned chip Plasma-cleaned glass slides c d After 1 sec normal Oxygen plasma
More informationSupplementary Figures
Supplementary Figures Supplementary Figure S1. a, the cross-sectional and b, top view SEM images of a PC/SWNT bilayer (SWNT film thickness of ~ 1µm). S1 Supplementary Figure S2. The obtained SWNT film
More informationChapter 5. Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice
Chapter 5 Effects of Photonic Crystal Band Gap on Rotation and Deformation of Hollow Te Rods in Triangular Lattice In chapter 3 and 4, we have demonstrated that the deformed rods, rotational rods and perturbation
More informationDynamic viscoelastic properties of polyvinyl chloride with physical aging
Korea-Australia Rheology Journal, 27(4), 259-266 (November 2015) DOI: 10.1007/s13367-015-0026-8 www.springer.com/13367 Dynamic viscoelastic properties of polyvinyl chloride with physical aging Fang Tian
More informationVIBRATION CHARACTERISTICS OF EMBEDDED DOUBLE WALLED CARBON NANOTUBES SUBJECTED TO AN AXIAL PRESSURE
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS VIBRATION CHARACTERISTICS OF EMBEDDED DOUBLE WALLED CARBON NANOTUBES SUBJECTED TO AN AXIAL PRESSURE X. W. Lei 1, T. Natsuki 2, J. X. Shi 1, Q. Q. Ni
More informationSTUDIES ON NANO-INDENTATION OF POLYMERIC THIN FILMS USING FINITE ELEMENT METHODS
STUDIES ON NANO-INDENTATION OF POLYMERIC THIN FILMS USING FINITE ELEMENT METHODS Shen Xiaojun, Yi Sung, Lallit Anand Singapore-MIT Alliance E4-04-0, 4 Engineering Drive 3, Singapore 7576 Zeng Kaiyang Institute
More informationThermal load-induced notch stress intensity factors derived from averaged strain energy density
Available online at www.sciencedirect.com Draft ScienceDirect Draft Draft Structural Integrity Procedia 00 (2016) 000 000 www.elsevier.com/locate/procedia 21st European Conference on Fracture, ECF21, 20-24
More informationApplication of Time-Temperature-Stress Superposition Principle to Nonlinear Creep of Poly(methyl methacrylate)
Key Engineering Materials Vols. 34-341 (27) pp. 191-196 online at http://www.scientific.net (27) rans ech Publications, Switzerland Application of ime-emperature-stress Superposition Principle to Nonlinear
More informationKINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)
KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP) 1 University of Science & Technology Beijing, China, niukm@ustb.edu.cn 2 Tsinghua University, Department of Engineering Mechanics, Beijing, China,
More informationThree-dimensional flexible and conductive interconnected graphene networks grown by chemical vapour deposition
SUPPLEMENTARY INFORMATION Three-dimensional flexible and conductive interconnected graphene networks grown by chemical vapour deposition S1. Characterization of the graphene foam (GF) and GF/PDMS composites
More informationNanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation. Mark Hardiman
Nanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation Mark Hardiman Materials and Surface Science Institute (MSSI), Department of Mechanical and Aeronautical
More informationEFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM
Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE
More informationA MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES
A MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES Anastassia N. Rissanou b,c*, Vagelis Harmandaris a,b,c* a Department of Applied Mathematics, University of Crete, GR-79, Heraklion, Crete,
More informationCARBON NANOTUBE-POLYMER COMPOSITES: AN OVERVIEW Brian Grady University of Oklahoma
CARBON NANOTUBE-POLYMER COMPOSITES: AN OVERVIEW Brian Grady University of Oklahoma Abstract Carbon nanotubes are in many ways similar to polymers. Both molecules have contour lengths typically on the order
More informationOn the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar
NDT&E International 33 (2000) 401 407 www.elsevier.com/locate/ndteint On the study of elastic wave scattering and Rayleigh wave velocity measurement of concrete with steel bar T.-T. Wu*, J.-H. Sun, J.-H.
More informationDirection sensitive deformation measurement with epoxy/cnt nanocomposites
Direction sensitive deformation measurement with epoxy/cnt nanocomposites S.T. Buschhorn*, M.H.G. Wichmann, J. Gehrmann, L. Böger, K. Schulte Technische Universität Hamburg-Harburg, Institute of Polymers
More informationDYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,
More informationDEVELOPMENT OF SOPHISTICATED MATERIALS USING THE COMPUTER SIMULATION
7 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE Baia Mare, Romania, May 17-18, 2007 ISSN-1224-3264 DEVELOPMENT OF SOPHISTICATED MATERIALS USING THE COMPUTER SIMULATION Zuzana MURČINKOVÁ Ing. Zuzana MURČINKOVÁ,
More informationPP/MWCNT/OC Nanocomposites Rheological and Mechanical Properties
International Workshop, Action COST FA0904 Characterization, Mechanics and Performance of Innovative Polymer Nanomaterials for Food Packaging Application, September 24-25, 2013, Varna, Bulgaria PP/MWCNT/OC
More informationDEVELOP WEAR-RESISTANT POLYMERIC COMPOSITES BY USING NANOPARTICLES
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS DEVELOP WEAR-RESISTANT POLYMERIC COMPOSITES BY USING NANOPARTICLES 1 Abstract L. Chang 1 *, K. Friedrich 2 1 School of Aerospace, Mechanical & Mechatronic
More informationPrediction of Micromechanical Behaviour of Elliptical Frp Composites
Prediction of Micromechanical Behaviour of Elliptical Frp Composites Kiranmayee.Nerusu Dept. of Mechanical Engg. P. V. P. Siddhartha Institute of Technology, Vijayawada 520 007, A.P, India. P. Phani Prasanthi
More informationMicromechanics modeling for the stiffness and strength properties of glass fibers/cnts/epoxy composites
High Performance Structures and Materials V 279 Micromechanics modeling for the stiffness and strength properties of glass fibers/cnts/epoxy composites M. Kim, F. A. Mirza & J. I. Song School of Mechatronics
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationModeling the Dynamic Propagation of Shear Bands in Bulk Metallic Glasses
Modeling the Dynamic Propagation of Shear Bands in Bulk Metallic Glasses B.J. Edwards, K. Feigl, M.L. Morrison*, B. Yang*, P.K. Liaw*, and R.A. Buchanan* Dept. of Chemical Engineering, The University of
More informationElectrical connection network within an electrically conductive adhesive
Electrical connection network within an electrically conductive adhesive D.Busek, P. Mach Department of Electrotechnology, Faculty of Electrical Engineering Technická 2, 166 27 Prague, Czech Republic,
More informationMechanical Properties of Polymers. Scope. MSE 383, Unit 3-1. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Mechanical Properties of Polymers Scope MSE 383, Unit 3-1 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Structure - mechanical properties relations Time-dependent mechanical
More informationXI. NANOMECHANICS OF GRAPHENE
XI. NANOMECHANICS OF GRAPHENE Carbon is an element of extraordinary properties. The carbon-carbon bond possesses large magnitude cohesive strength through its covalent bonds. Elemental carbon appears in
More informationPARTICLE SIZE EFFECT FROM MICRO TO NANO ON THE THERMAL AND MECHANICAL PROPERTIES OF POLYMER COMPOSITES
Indian Institute of Technology Delhi (IITD), New Delhi, 2013 PARTICLE SIZE EFFECT FROM MICRO TO NANO ON THE THERMAL AND MECHANICAL PROPERTIES OF POLYMER COMPOSITES by: M.HOSSEIN ALAEI Department of Applied
More informationA METHOD FOR CALCULATING SURFACE STRESS AND ELASTIC CONSTANTS BY MOLECULAR DYNAMICS
A METHOD FOR CALCULATING SURFACE STRESS AND ELASTIC CONSTANTS BY MOLECULAR DYNAMICS Satoshi Izumi Dept. of Mechanical Engineering, Univ. of Tokyo., JAPAN Shotaro Hara Dept. of Mechanical Engineering, Univ.
More informationDYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX
66 6th RILEM Symposium PTEBM'03, Zurich, 2003 DYNAMIC AND TRANSIENT TESTING OF ASPHALT BINDER AND PAVING MIX L. Zanzotto, O.J. Vacin and J. Stastna University of Calgary, Canada Abstract: A commercially
More informationModule 4: Behaviour of a Laminae-II. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites
Module 4: Behaviour of a Laminae-II Learning Unit 1: M1 M4.1 Mechanics of Composites M4.1.1 Introduction to Mechanics of Composites The relation between ply uniaxial strengths and constituent properties
More informationDependence of equivalent thermal conductivity coefficients of single-wall carbon nanotubes on their chirality
Journal of Physics: Conference Series PAPER OPEN ACCESS Dependence of equivalent thermal conductivity coefficients of single-wall carbon nanotubes on their chirality To cite this article: V S Zarubin and
More informationProbabilistic Micromechanics Analysis of CNT Nanocomposites with Three-Dimensional Karhunen-Loève Expansion
Probabilistic Micromechanics Analysis of CNT Nanocomposites with Three-Dimensional Karhunen-Loève Expansion Fei-Yan Zhu 1, Sungwoo Jeong 2 and Gunjin Yun 3 1), 2, 3) Department of Mechanical & Aerospace
More informationOverview of Probabilistic Modeling of Woven Ceramic Matrix Composites ABSTRACT
ABSTRACT An overview of current research efforts in understanding the cause of variability in the thermo-mechanical properties of woven ceramic matrix composites is presented. Statistical data describing
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationSIMULATION OF CNT COMPOSITES USING FAST MULTIPOLE BEM
194 Journal of Marine Science and Technology, Vol. 17, No. 3, pp. 194-22 (29) SIMULATION OF CNT COMPOSITES USING FAST MULTIPOLE BEM hen-han ao*, Jun-Dong u*, Hai-Tao Wang**, and iao-ping heng* Key words:
More informationTHE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF ORIENTED STRANDBOARD: A FINITE ELEMENT SIMULATION
Proceedings of McMat5: 5 Joint ASME/ASCE/SES Conference on Mechanics and Materials June 1-3, 5, Baton Rouge, Louisiana, USA 255 THE INFLUENCE OF IN-PLANE DENSITY VARIATION ON ENGINEERING PROPERTIES OF
More information