A Modulus gradient elasticity model for nano-reinforced composites

Gülaşık, Hasan
In this work, nanocomposites and their numerical simulations are studied. At the beginning of the study, the properties of the polymer nanocomposites are explained based on a specific nano-inclusion and polymer matrix couple, namely, carbon nanotube (CNT) and thermoplastic polyetheretherketone (PEEK) polymer. In a literature comparison study, it is shown that the properties of the constituents, interface properties, manufacturing methods, characterization methods and therefore mechanical properties of the CNT/PEEK nanocomposites can vary significantly among different studies. Classical elasticity formulations may become inadequate for the modeling of the nanostructured materials. They do not contain any information about the size and applicable from nanometer to meter scale. Moreover, they do not properly describe stress/strain singularities and are questionable if wavelength of deformation is comparable to dominant micro-structural length scale. Therefore, some extensions of the classical elasticity formulations have been proposed in literature. Two of the widely-used extensions, the Eringen’s nonlocal elasticity and the Aifantis’s gradient elasticity formulations, are explained. It is seen that nonlocal/gradient formulations include higher order fields and boundary conditions which are not easy to understand intuitively. They also have complex formulations and are computationally expensive. In this work, a new gradient elasticity formulation, the so-called E-grad model, is proposed to overcome some of the difficulties in the nonlocal and the gradient elasticity formulations. In the new formulation, similar to the differential relation between the local strain and the gradient enhanced strain in the classical models of gradient elasticity, a differential relation is proposed for the elastic constants of linear elasticity. Analytical and finite element solutions of the proposed formulation are derived for a one-dimensional inhomogeneous rod. The results of the proposed model are compared with a classical model of gradient elasticity for a one-dimensional model problem. It is seen that the discontinuities in the modulus, displacement, strain and stress fields are removed by the proposed model. Furthermore, there are no additional higher-order fields and boundary conditions and the numerical formulations are simpler than the nonlocal/gradient elasticity models. Then, the E-grad model is extended to more general three-dimensional inhomogeneous materials with isotropic linear elastic constituents. The finite element formulation for axisymmetric problems is derived and a model problem of a soft cylindrical rod with a stiff spherical inclusion is solved. It is seen that discontinuities and/or sharp changes in the modulus, displacement, strain and stress fields that exist in local formulations are smoothed out with the proposed model. The proposed model is compared with a micromechanical model from literature and experiments conducted on polyimide/silica nanocomposites. The results obtained by the proposed approach agree well with the experimentally measured values of the nanocomposite modulus. The model is also extended to obtain anisotropic macroscopic response by choosing different length scale parameters in different directions. At the end, a CNT reinforced polymer nanocomposite problem from literature is reconsidered in which the nanocomposite is assumed to be composed of four distinct phases: CNT, interface, interphase and bulk polymer. Rather than being homogeneous, the interphase is considered to be graded by the E-grad model. By using the E-grad model and the genetic algorithm optimization, homogenized elastic constants of the transversely isotropic effective fiber are calculated. It is seen that, although the effective fiber has higher modulus in the axial direction, it has lower modulus values in transverse and shear directions compared to the polymer matrix. Then, the effect of the orientation distribution of the effective fibers in a nanocomposite is taken into account by using an orientation distribution function. It is seen that, if effective fibers are aligned in a direction, the modulus of the composite increases in that direction as expected. However, it is also seen that, isotropic distribution of the effective fibers makes the composite to have lower modulus than the matrix due to low transverse and shear moduli of the effective fiber.


Investigation of InP and SiGe nanomaterials via molecular dynamics simulations
Nayir, Nadire; Erkoç, Şakir; Taşcı, Emre; Department of Physics (2018)
In this dissertation, InP and Si/Ge nanomaterials were investigated via classical molec- ular dynamics simulations. Structural and dynamical properties of InP nanorods and nanoparticles were examined under different conditions such as heating, cool- ing, strain etc. In the second part of this dissertation, a ReaxFF reactive force field development has been targeted to provide a more generalized description for the sim- ulation of the crystallization process of amorphous Ge on Si substrate. For this aim, the...
A Density Functional Study of Ni-2 and Ni-13 Nanoclusters
Önal, Işık; Uzun, Alper; ÖZKAR, SAİM (2009-04-01)
Ni-2 and Ni-13 clusters have been studied by means of DFT calculations performed at B3LYP level using three different basis sets. Equilibrium geometry calculations resulted in a binding energy of nickel dimer of 1.078 eV/atom for the basis set m6-31G*, showing good agreement with experimental data. Ni-13 nanocluster was found to have an equilibrium structure of icosahedral for all three basis sets in accordance with experimental findings. The intercept of the line obtained by drawing a correlation line of b...
An Asymptotic-Numerical Hybrid Method for Solving Singularly Perturbed Linear Delay Differential Equations
Cengizci, Süleyman (Hindawi Limited, 2017)
In thiswork, approximations to the solutions of singularly perturbed second-order linear delay differential equations are studied. We firstly use two-term Taylor series expansion for the delayed convection term and obtain a singularly perturbed ordinary differential equation (ODE). Later, an efficient and simple asymptotic method so called Successive Complementary Expansion Method (SCEM) is employed to obtain a uniformly valid approximation to this corresponding singularly perturbed ODE. As the final step, ...
A Novel Computational Method to Calculate Nonlinear Normal Modes of Complex Structures
Samandarı, Hamed; Ciğeroğlu, Ender (2019-01-31)
In this study, a simple and efficient computational approach to obtain nonlinear normal modes (NNMs) of nonlinear structures is presented. Describing function method (DFM) is used to capture the nonlinear internal forces under periodic motion. DFM has the advantage of expressing the nonlinear internal force as a nonlinear stiffness matrix multiplied by a displacement vector, where the off-diagonal terms of the nonlinear stiffness matrix can provide a comprehensive knowledge about the coupling between the mo...
A Fully Implicit Finite Volume Lattice Boltzmann Method for Turbulent Flow
Cevik, Fatih; Albayrak, Kahraman (2017-08-01)
Almost all schemes existed in the literature to solve the Lattice Boltzmann Equation like stream & collide, finite difference, finite element, finite volume schemes are explicit. However, it is known fact that implicit methods utilizes better stability and faster convergence compared to the explicit methods. In this paper, a method named herein as Implicit Finite Volume Lattice BoltzmannMethod (IFVLBM) for incompressible laminar and turbulent flows is proposed and it is applied to some 2D benchmark test cas...
Citation Formats
H. Gülaşık, “A Modulus gradient elasticity model for nano-reinforced composites,” Ph.D. - Doctoral Program, Middle East Technical University, 2018.