Browsing by Author "Javili, Ali"
Now showing 1 - 20 of 49
Results Per Page
Sort Options
Item Embargo A novel constitutive model for surface elasticity at finite strains suitable across compressibility spectrum(Elsevier Masson, 2023-03-24) Javili, Ali; Dörtdivanlioğlu, BerkinThe surface elasticity theory of Gurtin–Murdoch has proven to be remarkably successful in predicting the behavior of materials at the nano scale, which can be attributed to the fact that the surface-to-volume ratio increases as the problem dimension decreases. On the other hand, surface tension can deform soft elastic solids even at the macro scale resulting e.g. in elastocapillary instabilities in soft filaments reminiscent of Plateau–Rayleigh instabilities in fluids. Due to the increasing number of applications involving nanoscale structures and soft solids such as gels, the surface elasticity theory has experienced a prolific growth in the past two decades. Despite the large body of literature on the subject, the constitutive models of surface elasticity theory at large deformations are not suitable to capture the surface behavior from fully compressible to nearly incompressible elasticity, especially from a computational perspective. A physically meaningful and proper decomposition of the surface free energy density in terms of area-preserving and area-varying contributions remains yet to be established. We show that an immediate and intuitive generalization of the small-deformation surface constitutive models does not pass the simple extension test at large deformations and results in unphysical behavior at lower Poisson’s ratios. Thus, the first contribution of the manuscript is to introduce a novel decomposed surface free energy density that recovers surface elasticity across the compressibility spectrum. The second objective of this paper is to formulate an axisymmetric counterpart of the elastocapillary theory methodically derived from its three-dimensional format based on meaningful measures relevant to the proposed surface elasticity model. Various aspects of the problem are elucidated and discussed through numerical examples using the finite element method enhanced with surface elasticity.Item Open Access A versatile implicit computational framework for continuum-kinematics-inspired peridynamics(Springer, 2023-11-13) Firooz, S.; Javili, Ali; Steinmann, P.Continuum-kinematics-inspired peridynamics (CPD) has been recently proposed as a novel reformulation of peridynamics that is characterized by one-, two- and three-neighbor interactions. CPD is geometrically exact and thermodynamically consistent and does not suffer from zero-energy modes, displacement oscillations or material interpenetration. In this manuscript, for the first time, we develop a computational framework furnished with automatic differentiation for the implementation of CPD. Thereby, otherwise tedious analytical differentiation is automatized by employing hyper-dual numbers (HDN). This differentiation method does not suffer from round-off errors, subtractive cancellation errors or truncation errors and is thereby highly stable with superb accuracy being insensitive to perturbation values. The computational framework provided here is compact and model-independent, thus once the framework is implemented, any other material model can be incorporated via modifying the potential energy solely. Finally, to illustrate the versatility of our proposed framework, various potential energies are considered and the corresponding material response is examined for different scenarios.Item Open Access Aspects of computational homogenization in magneto-mechanics: Boundary conditions, RVE size and microstructure composition(Elsevier, 2018-01) Zabihyan, R.; Mergheim, J.; Javili, Ali; Steinmann, P.In the present work, the behavior of heterogeneous magnetorheological composites subjected to large deformations and external magnetic fields is studied. Computational homogenization is used to derive the macroscopic material response from the averaged response of the underlying microstructure. The microstructure consists of two materials and is far smaller than the characteristic length of the macroscopic problem. Different types of boundary conditions based on the primary variables of the magneto-elastic enthalpy and internal energy functionals are applied to solve the problem at the micro-scale. The overall responses of the RVEs with different sizes and particle distributions are studied under different loads and magnetic fields. The results indicate that the application of each set of boundary conditions presents different macroscopic responses. However, increasing the size of the RVE, solutions from different boundary conditions get closer to each other and converge to the response obtained from periodic boundary conditions.Item Open Access Aspects of interface elasticity theory(Sage Publications Ltd., 2018) Javili, Ali; Ottosen, N. S.; Ristinmaa, M.; Mosler, J.Interfaces significantly influence the overall material response especially when the area-to-volume ratio is large, for instance in nanocrystalline solids. A well-established and frequently applied framework suitable for modeling interfaces dates back to the pioneering work by Gurtin and Murdoch on surface elasticity theory and its generalization to interface elasticity theory. In this contribution, interface elasticity theory is revisited and different aspects of this theory are carefully examined. Two alternative formulations based on stress vectors and stress tensors are given to unify various existing approaches in this context. Focus is on the hyper-elastic mechanical behavior of such interfaces. Interface elasticity theory at finite deformation is critically reanalyzed and several subtle conclusions are highlighted. Finally, a consistent linearized interface elasticity theory is established. We propose an energetically consistent interface linear elasticity theory together with its appropriate stress measures.Item Open Access Atomistic two-, three- and four-body potentials. Spatial and material settings(Elsevier Ltd, 2021-09) Steinmann, P.; Smith, A.; Birang, E.; McBridge, A.; Javili, AliIn molecular dynamics or molecular statics (MD/MS) multi-body potentials empirically capture the energetic interactions in atomistic systems enabling the computation of the corresponding atomistic forces as energetic conjugates to the atomistic positions. We distinguish here between spatial and material atomistic positions and consequently between the corresponding spatial and material atomistic forces. In quasi-statics, i.e. MS, the former, also denoted as deformational atomistic forces, contribute to the classical deformational mechanics (i.e., equilibrium) problem that seeks to minimise the total potential energy of an atomistic system with respect to the atomistic positions relative to the ambient space. The latter, also denoted as configurational atomistic forces, contribute to the configurational mechanics (i.e., non-equilibrium) problem that determines the release of total potential energy of an atomistic system upon variation of the atomistic positions relative to the ambient material, i.e., due to perturbations of the material (initial) atomistic configuration. The importance of material atomistic forces is that they drive energetically favourable re-organisations of the material atomistic configuration, thereby characterising the tendency of generic atomistic defects to propagate. In this contribution we focus on two-, three-, and four-body potentials, whereby we distinguish between novel stretch- and classical angle-based potentials for the two latter cases. Taken together, as the main contribution, we derive expressions for the corresponding spatial and, for the first time, material atomistic forces and highlight their striking formal similarity. The derivations are detailed but the final expression compact and well-suited for numerical implementation.Item Open Access Bifurcation behavior of compressible elastic half-space under plane deformations(Elsevier, 2020) Bakiler, A. Derya; Javili, AliA finitely deformed elastic half-space subject to compressive stresses will experience a geometric instability at a critical level and exhibit bifurcation. While the bifurcation of an incompressible elastic half-space is commonly studied, the bifurcation behavior of a compressible elastic half-space remains elusive and poorly understood to date. The main objective of this manuscript is to study the bifurcation of a neo-Hookean compressible elastic half-space against the well-established incompressible case. The formulation of the problem requires a novel description for a non-linear Poisson’s ratio, since the commonly accepted definitions prove insufficient for the current analysis. To investigate the stability of the domain and the possibility of bifurcation, an incremental analysis is carried out. The incremental analysis describes a small departure from an equilibrium configuration at a finite deformation. It is shown that at the incompressibility limit, our results obtained for a compressible elastic half-space recover their incompressible counterparts. Another key feature of this contribution is that we verify the analytical solution of this problem with computational simulations using the finite element method via an eigenvalue analysis. The main outcome of this work is an analytical expression for the critical stretch where bifurcation arises. We demonstrate the utility of our model and its excellent agreement with the numerical results ranging from fully compressible to incompressible elasticity. Moving forward, this approach can be used to comprehend and harness the instabilities in bilayer systems, particularly for compressible ones.Item Open Access Boundary viscoelasticity theory at finite deformations and computational implementation using isogeometric analysis(Elsevier BV, 2021-02-01) Dortdivanlioglu, B.; Javili, AliUse of surface elasticity theory has experienced a prolific growth recently due to its utility in understanding the mechanics of nanomaterials and soft solids at small scales. Various extensions of surface elasticity theory have been proposed. The main objective of this contribution is to formulate a finite deformation theory for boundary viscoelasticity in principal stretches by accounting for strain-dependent boundary stresses. We present a model that utilizes a nonlinear evolution law and thus is not restricted to the states that are close to the thermodynamic equilibrium. Boundary contributions include both surface and curve effects wherein boundary elasticity as well as boundary tension are accounted for. The boundary constitutive models are formulated such that fluid-like and solid-like viscoelastic behavior of boundaries are considered. A geometrically exact computational framework using isogeometric analysis inherently suited to account for boundaries is developed. Equipped with the theoretical and computational framework, the influence of boundary viscoelasticity on the material response is illustrated. Non-equilibrium counterpart of surface tension is introduced and its effects are elucidated via examples. Through numerical examples, various applications of the bulk–boundary coupled formulation which require further investigation are highlighted.Item Open Access Bounds on size effects in composites via homogenization accounting for general interfaces(Springer, 2020-01) Firooz, Soheil; Javili, Ali; Chatzigeorgiou, G.This manuscript provides novel bounds and estimates, for the first time, on size-dependent properties of composites accounting for generalized interfaces in their microstructure, via analytical homogenization verified by computational analysis. We extend both the composite cylinder assemblage and Mori–Tanaka approaches to account for the general interface model. Our proposed strategy does not only determine the overall response of composites, but also it provides information about the local fields for each phase of the medium including the interface. We present a comprehensive study on a broad range of interface parameters, stiffness ratios and sizes. Our analytical solutions are in excellent agreement with the computational results using the finite element method. Based on the observations throughout our investigations, two notions of size-dependent bounds and ultimate bounds on the effective response of composites are introduced which yield a significant insight into the size effects, particularly important for the design of nano-composites.Item Open Access Bounds on size-dependent behaviour of composites(Taylor & Francis, 2018) Saeb, S.; Steinmann, P.; Javili, AliComputational homogenisation is a powerful strategy to predict the effective behaviour of heterogeneous materials. While computational homogenisation cannot exactly compute the effective parameters, it can provide bounds on the overall material response. Thus, central to computational homogenisation is the existence of bounds. Classical firstorder computational homogenisation cannot capture size effects. Recently, it has been shown that size effects can be retrieved via accounting for elastic coherent interfaces in the microstructure. The primary objective of this contribution is to present a systematic study to attain computational bounds on the sizedependent response of composites. We show rigorously that interface-enhanced computational homogenisation introduces two relative length scales into the problem and investigate the interplay between them. To enforce the equivalence of the virtual power between the scales, a generalised version of the Hill–Mandel condition is employed, and accordingly, suitable boundary conditions are derived. Macroscopic quantities are related to their microscopic counterparts via extended average theorems. Periodic boundary conditions provide an effective behaviour bounded by traction and displacement boundary conditions. Apart from the bounds due to boundary conditions for a given size, the size-dependent response of a composite is bounded, too. The lower bound coincides with that of a composite with no interface. Surprisingly, there also exists an upper bound on the size-dependent response beyond which the expected ‘smaller is stronger’ trend is no longer observed. Finally, we show an excellent agreement between our numerical results and the corresponding analytical solution for linear isotropic materials which highlights the accuracy and broad applicability of the presented scheme.Item Open Access The computational framework for continuum-kinematics-inspired peridynamics(Springer Science and Business Media B.V., 2020) Javili, Ali; Firooz, Soheil; McBride, A. T.; Steinmann, P.Peridynamics (PD) is a non-local continuum formulation. The original version of PD was restricted to bond-based interactions. Bond-based PD is geometrically exact and its kinematics are similar to classical continuum mechanics (CCM). However, it cannot capture the Poisson effect correctly. This shortcoming was addressed via state-based PD, but the kinematics are not accurately preserved. Continuum-kinematics-inspired peridynamics (CPD) provides a geometrically exact framework whose underlying kinematics coincide with that of CCM and captures the Poisson effect correctly. In CPD, one distinguishes between one-, two- and three-neighbour interactions. One-neighbour interactions are equivalent to the bond-based interactions of the original PD formalism. However, two- and three-neighbour interactions are fundamentally different from state-based interactions as the basic elements of continuum kinematics are preserved precisely. The objective of this contribution is to elaborate on computational aspects of CPD and present detailed derivations that are essential for its implementation. Key features of the resulting computational CPD are elucidated via a series of numerical examples. These include three-dimensional problems at large deformations. The proposed strategy is robust and the quadratic rate of convergence associated with the Newton–Raphson scheme is observed.Item Open Access Configurational peridynamics(Elsevier B.V., 2023-07-31) Steinmann, P. ; de Villiers, A.M. ; McBride, A.T. ; Javili, AliConfigurational forces that drive the evolution of material structures such as defects are introduced into a geometrically-exact peridynamics framework. The concept of bond-number double-density facilitates the definition of a peridynamic potential energy functional that inherits the key features of its conventional (local) continuum and discrete counterparts. The spatial and material variations of the peridynamic potential energy functional give rise to familiar Piola- and Cauchy-type bond-wise interaction forces that enter the pointwise force balance in the spatial and material setting, respectively. It is shown that the point-wise material body force density is a result of a non-local pull-back of the bond-wise spatial interaction force, and thereby captures non-local contributions. Several key features of configurational peridynamics are demonstrated via a computational example and a comparison to conventional configurational continuum mechanics.Item Open Access Continuum-kinematics-inspired peridynamics. Mechanical problems(Elsevier, 2019) Javili, Ali; McBride, A. T.; Steinmann, P.The main objective of this contribution is to develop a novel continuum-kinematicsinspired approach for peridynamics (PD), and to revisit PD’s thermodynamic foundations. We distinguish between three types of interactions, namely, one-neighbour interactions, two-neighbour interactions and three-neighbour interactions. While one-neighbour interactions are equivalent to the bond-based interactions of the original PD formalism, twoand three-neighbour interactions are fundamentally different to state-based interactions in that the basic elements of continuum kinematics are preserved exactly. In addition, we propose that an externally prescribed traction on the boundary of the continuum body emerges naturally and need not vanish. This is in contrast to, but does not necessarily violate, standard PD. We investigate the consequences of the angular momentum balance and provide a set of appropriate arguments for the interactions accordingly. Furthermore, we elaborate on thermodynamic restrictions on the interaction energies and derive thermodynamically-consistent constitutive laws through a Coleman–Noll-like procedure.Item Restricted Continuum-kinematics-inspired peridynamics: thermo-mechanical problems(Springer, 2021-03-31) Javili, Ali; Ekiz, Ekim; McBride, A. T.; Steinmann, P.The recently proposed continuum-kinematics-inspired peridynamics (CPD) is extended to account for thermo-mechanical coupling at large deformations. The key features of CPD are that it is geometrically exact and is built upon multi-neighbour interactions. The bond-based interactions of the original PD formalism are equivalent to one-neighbour interactions of CPD. Two- and three-neighbour interactions, however, are fundamentally different from state-based PD in that the basic elements of continuum kinematics are preserved exactly. We elaborate on thermodynamic restrictions on the interaction energies and derive thermodynamically consistent constitutive laws through a Coleman–Noll-like procedure. Notably, we show that various choices for temperature, or coldness, satisfy the dissipation inequality and provide meaningful temperature, or coldness, evolution equations together with Fourier-like conduction relationsItem Open Access Correction to: The computational framework for continuum-kinematics-inspired peridynamics(Springer Science and Business Media Deutschland GmbH, 2020-09-03) Javili, Ali; Firooz, Soheil; McBride, A. T.; Steinmann, P.Item Open Access Describing droplet motion on surface-textured ratchet tracks with an inverted double pendulum model(American Chemical Society, 2021-04-27) Naji, Mayssam; Yelekli Kirici, Ecem; Javili, Ali; Erdem, Emine YeganWe describe the motion of a droplet on a textured ratchet track using a nonlinear resonator model. A textured ratchet track is composed of a semicircular pillar array that induces a net surface tension local gradient on a droplet placed on it. When a vertical vibration is applied, hysteresis is overcome, and the droplet moves toward the local lower energy barrier; however, due to the repetitive structure of texture, it keeps moving until the end of the track. The droplet motion depends on the amplitude and frequency of the vertical oscillation, and this dependence is nonlinear. Therefore, finding a fully analytic solution to represent this motion is not trivial. Consequently, the droplet motion remains poorly understood. In this study, we elaborate on the utility of a double pendulum as a basis for modeling the droplet motion on surfaces inducing asymmetric force. Similar to the droplet motion, resonators, such as a double pendulum, are simple, yet nonlinear systems. Moreover, an inverted double pendulum motion has key characteristics such as the two-phase motion and the double peak motion, which are also observed in the droplet motion. We use various data-processing methods to highlight the similarity between these two systems both qualitatively and quantitatively. After establishing this comparison, we propose a model that utilizes an inverted double pendulum mounted on a moving cart to successfully simulate the motion of a droplet on a ratchet track. This methodology will lead to the development of an accurate droplet-motion modeling approach, and we believe that it will be useful to understand droplet dynamics more deeply.Item Open Access Designing tunable composites with general interfaces(Elsevier, 2019) Saeb, S.; Steinmann, P.; Javili, AliIn this manuscript, we employ interface enhanced computational homogenization to explore and detail on a number of unfamiliar characteristics that composites can exhibit at different length scales. Here, the interface between the constituents is general in the sense that both displacement and traction jumps across the interface are admissible. We carry out numerous computational investigations using the finite element method for a broad range of various material parameters. Our numerical results reveal that the effective response of a microstructure embedding general interfaces is intuitively unpredictable and highly complex. In particular, for certain ranges of material parameters the overall response shows insensitivity with respect to either microstructure size or stiffness-ratio between inclusion and matrix. This unique behavior is observed likewise for two- and three-dimensional unit-cells. Our findings provide a valuable guideline to design tunable composites utilizing interfaces.Item Open Access A different catch for Poisson(Springer, 2022-05-04) Bakiler, A. Derya; Javili, Ali; Giorgio, Ivan; Placidi, Luca; Barchiesi, Emilio; Emek Abali, Bilen; Altenbach, HolmPoisson’s ratio, similar to other material parameters of isotropic elasticity, is determined via experiments corresponding to small strains. Yet at small-strain linear elasticity, Poisson’s ratio has a dual nature; although commonly understood as a geometrical parameter, Poisson’s ratio is also a material parameter. From a geometrical perspective only, the concept of Poisson’s ratio has been extended to large deformations by Beatty and Stalnaker. Here, through a variational analysis, we firstly propose an alternative relationship between the Poisson ratio and stretches at finite deformations such that the nature of Poisson’s ratio as a material parameter is retained. In doing so, we introduce relationships between the Poisson ratio and stretches at large deformations different than those established by Beatty and Stal naker. We show that all the nonlinear definitions of Poisson’s ratio coincide at the reference configuration and thus, material and geometrical descriptions too coincide, at small-strains linear elasticity. Secondly, we employ this variational approach to bring in the notion of nonlinear Poisson’s ratio in peridynamics, for the first time. In particular, we focus on bond-based peridynamics. The nonlinear Poisson’s ratio of bond-based peridynamics coincides with 1/3 for two-dimensional and 1/4 for three-dimensional problems, at the reference configuration.Item Open Access A displacement-based approach to geometric instabilities of a film on a substrate(Sage Publications, 2019-02) Javili, Ali; Bakiler, A. DeryaWhen a thin film adhered to a compliant substrate is growing, it will eventually buckle in order to release the compressive stresses accumulated within the film due to growth. Such geometric instabilities caused by compressive stresses prevail among all living systems in nature and their outcomes range from highly beneficial to destructive. Therefore, understanding compression induced instabilities is of crucial importance. Note that the origin of the “compression” need not necessarily be differential growth, as it may be due to pre-stretch or thermal expansion. A commonly accepted solution strategy for instabilities in bilayer structures dates back to the seminal work of Allen and employs the Airy stress functions. Owing to its reliance on a stress-based approach, the Allen solution is limited to linear two-dimensional problems and its success depends entirely on choosing an appropriate Airy function. The main objective of this contribution is to circumvent these limitations via a displacement-based approach formally suitable for three-dimensional problems, anisotropic materials, and even applicable to finite deformations. Furthermore, the Allen solution in its original form is valid for the plane-stress condition but often it is mistakenly compared with the numerical simulations corresponding to the plane-strain condition. We analyze the subtle difference between the solutions associated with the plane-strain and plane-stress conditions. Next, the analytical solution is compared against the computational results using the finite element method via eigenvalue analysis. Finally, it is briefly explained how the current approach can be utilized beyond the classical bilayer systems.Item Embargo Extended general interfaces: Mori–Tanaka homogenization and average fields(Elsevier Ltd, 2022-08-24) Firooz, S.; Chatzigeorgiou, G.; Steinmann, P.; Javili, AliA well-established methodology to capture interphases in heterogeneous materials is to replace them by a zero-thickness interface model. Commonly accepted interface models intuitively assume that to satisfy the angular momentum balance, interfaces must coincide with the mid-layer of their corresponding interphases. Recently, via adopting weighted averages, an extended general interface model has been developed that allows for arbitrary interface positions while fulfilling the angular momentum balance. This manuscript incorporates this novel interface model into the Mori–Tanaka method within the framework of homogenization. Analytical solutions are developed to determine effective properties as well as average local fields for fiber-reinforced and particle-reinforced composites. Computational simulations using the finite element method (FEM) are carried out to compare with the analytical solutions. Through a set of numerical examples, the significance of the interface position on the overall response of heterogeneous materials is highlighted. Our extended framework clarifies various ambiguous observations originating from the trivial assumption of restricting the interface position to the mid-plane. One advantage of the current interface model is that it covers both the elastic and cohesive interface models at its limits and therefore the analytical solutions are widely applicable regardless of the interface type.Item Open Access From beams to bilayers: A unifying approach towards instabilities of compressible domains under plane deformations(Elsevier Ltd, 2021-10) Bakiler, A. Derya; Dörtdivanoğlu, B.; Javili, AliInstabilities that form when a domain of compliant elastic material goes under compressive forces are prevalent in nature and have found many applications. Even though instabilities are observed in a myriad of fields and materials, the large deformation bifurcation analysis of compressible domains, may it be beams, half-spaces, or bilayers, remains understudied compared to the incompressible case. In this work, we present a unifying approach for the instability analysis of a compressible elastic domain under plane deformations, wherein the unifying approach is then particularized for beams, half-spaces, and bilayers. First, the large-deformation incremental analysis for a rectangular, compressible, hyperelastic domain under plane deformations is developed, which serves as a generic and all-encompassing framework for other geometries. Subsequently, this generic framework is applied to the specific domains of beam, half-space, and lastly as the superimposition of the two; bilayer. Obtained analytical results for the onset of wrinkling in the beam, half-space and bilayer geometries are explored in the full range of compressibility and for various geometrical parameters, including their comparison with computational simulations using the finite element method, cultivating excellent agreements between analytical and numerical results all across the material and geometrical parameter spectrum. The analytical framework presented here provides grounds for further works on other modes of instabilities and more complex geometries.
- «
- 1 (current)
- 2
- 3
- »