Summary of doctoral thesis in Mechanical engineering and Engineering mechanics: An approach to approximate and fem-model the conductivity and elasticity of multicomponent material

Chia sẻ: _ _ | Ngày: | Loại File: PDF | Số trang:28

Thêm vào BST

Báo xấu

26
lượt xem 4
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Equivalent inclusion approach is then developed to account for possible diversions, such as non-idealistic geometric forms of the inhomogeneities, imperfect matrix-inclusion contacts, filler dispersions, and when the particular values of the fillers’ properties are unspecified, using available numerical or experimental reference conductivity data for particular composites.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Summary of doctoral thesis in Mechanical engineering and Engineering mechanics: An approach to approximate and fem-model the conductivity and elasticity of multicomponent material

1 MINISTRY OF EDUCATION AND VIETNAM ACADEMY OF TRAINING SCIENCE AND TECHNOLOGY GRADUATE UNIVERSITY SCIENCE AND TECHNOLOGY ----------------------------- DO QUOC HOANG AN APPROACH TO APPROXIMATE AND FEM-MODEL THE CONDUCTIVITY AND ELASTICITY OF MULTI- COMPONENT MATERIAL Major: Engineering mechanics Code: 9 52 01 01 SUMMARY OF DOCTORAL THESIS IN MECHANICAL ENGINEERING AND ENGINEERING MECHANICS Hanoi – 2019
2 The thesis has been completed at: Graduate University Science and Technology – Vietnam Academy of Science and Technology. Supervisors: 1. Assoc. Prof. DrSc. Pham Duc Chinh 2. Assoc. Prof. Dr. Tran Anh Binh Reviewer 1: Reviewer 2: Reviewer 3: Thesis is defended at Graduate University Science and Technology-Vietnam Academy of Science and Technology at ….., on ……….. Hardcopy of the thesis be found at : - Library of Graduate University Science and Technology - Vietnam national library
3 INTRODUCTION 1. Relevance of the thesis topic Most effective medium approximations for isotropic inhomogeneous materials are based on dilute solutions of some typical inclusions in an infinite matrix medium, while the simplest approximations are those for the composites with spherical and circular inclusions. Practical particulate composites often involve inhomogeneities of more complicated geometry than that of the spherical (or circular) one. In our approach, those inhomogeneities are supposed to be substituted by simple equivalent spherical (circular) inclusions from a comparison of their respective dilute solution results. Then the available simple approximations for the equivalent spherical (circular) inclusion material can be used to estimate the effective conductivity of the original composite. Numerical illustrations of the approach are performed on some 2D and 3D geometries involving elliptical and ellipsoidal inclusions. 2. Thesis objective Develop near interaction approximations for the conductivity and elasticity of multi-component materials with spherical (circular) form inclusions. Equivalent inclusion approach is then developed to account for possible diversions, such as non-idealistic geometric forms of the inhomogeneities, imperfect matrix-inclusion contacts, filler dispersions, and when the particular values of the fillers’ properties are unspecified, using available numerical or experimental reference conductivity data for particular composites. We use the eXtended Finite Elements Method (XFEM) to estimate the effective conductivity of 2D macroscopically-isotropic composites containing elliptic inclusions and the equivalent ones with circular inclusions for comparisons with the approximations. 3. Scope The thesis focuses on conductivity and elasticity of multi- component materials, the Finite Element Method (FEM) and approximation schems
4 4. Research methods  Near interaction approximations has been constructed from the minimum energy for the macroscopic conductivity and elasticity of the multi-component matrix composites with spherical (circular) inclusions. Equivalent replacement of complex-geometry inclusions by the equivalent spherical, circular, disk and needle ones with equivalent properties using polarization approximation, dilute solutions, and experimental referemce results.  Numerical method: use Matlab program to homogenize some periodic material models in the framework of FEM method (XFEM). The results of FEM are considered as the accurate reference results for comparisons with the approximation ones. 5. The contributions of the thesis Beside Introduction section, the thesis contains 3 Chapters, a Conclusion section and a list of publications relevant to the thesis. References cited in the thesis are listed at the end of the thesis. CHAPTER 1. OVERVIEW 1.1. Opening Multi-component materials have complex structures, different individual mechanical properties. Many authors offered different evaluation methods, including the effective medium approximations and the variational ones. Geometric parameters have bên added to improve the étimates. In this chapter, the author presents the concept of hômgenization and an overview of the constructions of approximation methods for complex multi-component materials. The stress field  (x) is related to the strain field  ( x) by Hook’s law:  (x)  C(x) :  (x), (1.1) The average values of the stress and strain on V is defined as: 1 1    dx ,    dx. V V V V (1.2)
5 Assume homogeneous boundary conditions for displacements: u(x)   0  x. (1.3) Or the respective ones for the tractions  n   0 n (1.4) With the solutions σ, ε on V, the relationship between the averaged stress and strain on V is presented through the effective elastic tensor Ceff:   Ceff :  , Ceff  T(k eff ,  eff ). (1.5) k eff and  eff are effective elastic bulk and shear moduli. Another approach is to determine the effective elasticity coefficients by finding the infimum of the energy function on V (the fields  need to be compatible):  0 : Ceff :  0  inf   : C :  dx , (1.6)     0 V Or through the dual principle (the fields  need to be equilibrated):  0 : (Ceff )1 :  0  inf   : (C) 1 :  dx. (1.7)   0 V Similarly, the equations for the conductivity problem: The flux J must satisfy the equilibrium condition: ·J (x)  0 With the solutions J, E  T on V, the thermal conductivity coefficient (effective) ceff is determined as: J  c eff E  c eff T . (1.8) The minimum energy principles are also the main tools to find the macroscopic conductivity: c eff E0 ·E0  inf 0  cE·Edx, (1.9)  E  E V and:
6 (c eff ) 1 J 0 ·J 0  inf 0  c1 J·Jdx, (1.10)  J  J V 1.2. Overview of approximation methods for multi-component materials 1.2.1. Dilute solutions The effective conductivity ceff of the dilute solution of ellipsoidal inclusions with axes ratio a: b: c, randomly oriented in a continuous matrix is expressed in the form: c eff  cM  vI (cI  cM ) Dc (cI , cM ) , vI  1 , (1.11) cM 1 1 1 Dc (cI , cM )  [   ], 3 cI A  cM (1  A) cI B  cM (1  B ) cI C  cM (1  C ) The general fomula of Dc (cI , cM ) for spherical (d=3) and circular (d=2) inclusions is: dcM Dc (cI , cM )  . cI  (d  1)cM 1.2.2. Maxwell Approximation Maxwell approximation is built for 2-phase material from the matrix + spherical inclusions with any volume ratios, not limited by dilute distribution case (M - matrix symbol, I – inclusion symbol). 1 eff  vI v  c    M    d  1 cM , c  I   d  1 c M dcM  1 (1.12) eff  vI vM  2(d  1) K      K*M ; K*M  M ,  I  K  d  1  *M K K M  K *M  d vI vM d 2 K M  2(d  1)(d  2) M  eff   MA  (  ) 1  *M ; *M  M .  I  * M M  * M 2dK M  4d  M 1.2.3. Differential Approximation - DA we obtain the following differential equations for the effective conductivity ceff = c(1) of the composite
7 dc 1 n   vI (cI  c) Dc (cI , c), dt 1  vI t  1 (1.16a) n c(0)  cM , 0  t  1 , vI   vI ,  1 For elastic coefficient dK 1 n   vI ( K I  K ) DK ( K I , I , K ,  ), dt 1  vI t  1 d 1 n   vI ( I   ) DK ( K I ,  I , K ,  ), dt 1  vI t  1 (1.16b) n K (0)  K M ,  (0)   M 0  t  1 , vI   vI  ,  1 1.2.4. Self-consistent approximations - SA The Self-consistent approximation method (SA) for composite materials n components, is cSA=c solution of the following equation: n   v  (c   c ) D  (c  , c )  0 . 1 I I c I (1.17) SA for the moduli of elasticity are the solutions KSA=K and  SA=  of a system of two equations n   v  ( K   K ) D  ( K  ,   , K ,  )  0, 1 I I K I I n (1.18)   v  (    )D 1 I I M ( K I ,  I , K ,  )  0 . 1.2.5. Mori-Tanaka Approximation (MTA) Mori-Tanaka type approximation (MTA), based on the assumption that the fields in an circular inclusion are determined as if it is embedded in the matrix with remote average strain of the matrix, yields [Le Quang] for two-phase composites c MTA  c M  vI (c I  c M )·{vM [I  p·c M1·(c I  c M )]  vI I}1 . (1.19)
8 While for the multi-phase ones (matrix + n inclusions) n c MTA  {vM c M   vI c I ·[I  p ·c M1·(c I  c M )]1}  1 n (1.20) 1 1 1 ·{vM I   vI [I  p ·c ·(c I  c M )] } . M  1 MTA for the effective conductivity of d-dimentional multi- component isotropic materials with spherical inclusions (circular) has the following form n   v  (c   c 1 I I M )dcM / [cI  (d  1)cM ] cMTA  cM  n . (1.21) vM   vI dcM / [cI  (d  1)cM ]  1 1.3. Three-point correlation estimates of Phạm ĐC Three-point correlation estimates of Phạm ĐC are for the effective thermal conductivity of the multi-component materials. The bounds have been built from the minimum energy principles. The general expression of the upper bound for ceff is c eff  Pc (2c0 )  c** , (1.22) where c0 is a positive parameter, 1  v  Pc (c** )       c** , (1.23)   c  c*  n n v c**  3( ) 2  (c  c0 )A X  X  . (1.24)  1 c  2c0  ,  , 1 Similarly, the respective expression of the dual principle is written in the form: n n v c**  3c02 (1  2c0  )2  (c1  c01 )A X  X  . (1.25)  1 c  2c0  ,  , 1
9 Where n v 1 X     1 c  2c0 c  2c0 n (1.26) v 1 X     1 c  2c0 c  2c0 We choose the value c0 to eliminate the component c** , c** to make the inequality stronger to get the respective bounds. 1.3. Finite element method for homogenization solution The finite element method (FEM), is a numerical method for solving problems of engineering and mathematical physics. Typical problem areas of interest include structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential. The analytical solution of these problems generally require the solution to boundary value problems for partial differential equations. The finite element method formulation of the problem results in a system of algebraic equations. The method approximates the unknown function over the domain. To solve the problem, it subdivides a large system into smaller, simpler parts that are called finite elements. The simple equations that model these finite elements are then assembled into a larger system of equations that models the entire problem. FEM then uses variational methods from the calculus of variations to approximate a solution by minimizing an associated error function. A typical work out of the method involves dividing the domain of the problem into a collection of subdomains, with each subdomain represented by a set of element equations to the original problem, followed by systematically recombining all sets of element equations into a global system of equations for the final calculation. The global system of equations has known solution techniques, and can be calculated from the initial values of the original problem to obtain a numerical answer.
10 CHƯƠNG 2. FINITE ELEMENT METHOD 2.1. Introduce 2.1.1. FEM for thermal solution The finite element method with fine meshes is now used for reference comparisons. Let us consider a periodic cell Ω with the external boundary ∂Ω. The strong form of the conductivity problem is written as: q  x  0 in , q  x  c  x E  x in , E  x   T  x  in , (2.1) T x periodicin , q  x  .n antiperiodicin , The weak form associated with the above equations is given by:  cN .N T d   i i i (2.2) Using linear form functions for triangular elements with 3 joint as: N  x, y   ax  by  c, (2.3) Equation in matrix form:  E e ( x)    B e  T e  . (2.4) The effective thermal conductivity in x1 direction can be computed as: L1 c1eff  q1d . V T V The domain V is the periodic cell, so the effective thermal conductivity of RVE is given by cFeffEM  c1eff
11 We solve the temperature at each node element position, and find the effective conductivity coefficient according to the heat transfer equation: T q avg   k eff . (2.5) X1 2.1.2. Model of thermal solution We consider the three-dimensional cubic periodic microstructure as given in Fig.2.1. The finite element mesh for a periodic cell is given in Fig. 2.2. Fig 2.1: Mesing up for Body center cubic Some results about thermal solution: Fig 2.2: Thermal distribution in the model
12 2.1.2. Model for elastic solution The displacement field according to the degrees of freedom at the joint element qe ue   N .qe . (2.6) The deformation state of the joint elements will be:  e   ue     N qe   B qe . (2.7) The equation for stress of elements:  e   D  e . (2.8) The full potential of the element:  u   U e e  Ae . (2.9)  e q   1 q  K  q  q  P  . e T e e e T e e (2.10) 2 As a result, we get the equation to solve by finite element method:   Ne T  Ne T    Le  K e  Le  q     Le  P e  0. (2.11)  q   e 1  e 1 Or  K  q    P . (2.12) 2.2. Extended – Finite Element Method (XFEM) 2.2.1. LevelSet function The XFEM displacement approximation can be expressed by u h ( x)   N i ( x)ui   N *j ( x) ( x)a j . (2.13) i i e One important example of such a function would be the signed distance function:
13  ( x)  s( x) x  x , (2.14) Fig 2.3: LevelSet function An example of level set function for one and four inclusions is provided in Figure 3: Fig 2.4 Level set function defining a circular interface 2.2.2. Using eXtend – Finite element method for thermal solution The field equations of thermo-static problem are given by   q ( x)  r ( x)  0 in  (2.15) q ( x)  C ( x)T ( x) Above, q(x) denotes the heat flux, r(x) is a heat source term and c(x) the conductivity tensor. The periodic boundary conditions are: q.n is antiperiodic on  , and T is periodic on  . (C) is the conductivity matrix. More precisely, C(x) = CI for anisotropic inclusions is defined as CI   RT C  R, (2.16) where C∗ is conductivity matrix in local coordinate system, The weak form associated is given by finding such that
14  q(T )d    rTd  0,   By substituting the temperature field defined in (8) into the weak form (16) we obtain the discrete system of linear ordinary equations T Kd  Q, d  T a where d are nodal unknowns and K and Q are the global stiffness matrix and external flux, respectively. More precisely, the matrix K and vector Q are defined by K   BT c( x) Bd, Q   N T rd.   where B and N are the matrices of shape function derivatives and shape functions associated with the approximation scheme (8). The periodic boundary condition is finally introduced to (17) by mean of multiplicator Lagrange. 2.2.2. Results of thermal solution, using XFEM The calculation results indicate the temperature change in the calculation model Hình 2.5: Temperature distribution in the model CHAPTER 3. POLARIZATION APPROXIMATION (PA) 3.1. Introdution Consider a representative volume element (RVE) of an isotropic n-icomponent material that occupies spherical region V of Euclidean space. The center of the sphere is also the origin of the Cartesian system of coordinates {x}. The RVE consists of n components occupying regions Va V of volumes va and having conductivities (thermal, electrical, etc.) ca(a = 1,. . .,n; the volume of V is assumed
15 to be the unity). Starting from the minimum energy principles and using Hashin–Shtrikman polarization trial fields, one derives the following three-point correlation bounds on the effective conductivity ceff of the composite (Le & Pham, 1991; Pham, 1993) Pc (2c0 )  c**  c eff  [ Pc1 (2c0 )  c** ]1 . (3.1) where n v 1 Pc (c* )  ( )  c* . (3.2)  1 c  c* v A   ,ij,ij dx    , (3.3) V 3 where conventional summation on repeating Latin indices (but not on the Greek indices) is assumed; Latin indices after comma designate differentiation with respective Cartesian coordinates; the arbitrary positive constant c0 is often referred to as the conductivity of a comparison material; the harmonic potentials ua(x) appear in the expressions of Hashin– Shtrikman polarization trial fields; the three- point correlation parameters Ab ac relate the microgeometries of the three phases Va, Vb, Vc. Simple property (polarization) functions P, being monotonously increasing functions of their arguments and sharing the same structure, shall take a central place in our bounds and estimates for an easy qualitative comparisons between them. If one takes c0 = cmax = max{c1,. . .,cn} (or c0 = cmin = min{c1,. . .,cn}), then c**  0 (or c**  0 ) and can be neglected to strengthen the inequalities in (1), and subsequently one obtains Hashin–Shtrikman bounds Pc (2cmax )  c eff  Pc (2cmin ) . (3.4) With (3.5) in hands, one substitutes c0 = cM into (3.7) to obtain c**  c**  0 . Then, the bounds (1) converge to the unique value of the effective conductivity of the model c eff  Pc (2cM ) . (3.5)
16 Then taking c0 = cM, one finds that ** c c 0 ** and deduces the polarization approximation (PA) for the effective conductivity of our n-component matrix-inclusion composite (generally in d dimensions) expressed through a property (polarization) function P: c eff  cPA  Pc ((d  1)cM ). (3.6) 3.2. Result Using Ansys softwate mesh model for FEM, and author build a program to calculate by Matlab, show results as graphs. 3.2.1. Two-dimensional periodic three-component composites formed We examine two-dimensional periodic three-component composites formed from a continuous matrix phase and two inclusion phases, which have isotropic effective isotropic properties. The first one is bodycentered square periodic cell as given in the figure 1 and the second one is body-centered hexagonal microstructure as shown in the figure 2. The diameter of inclusions of each phases is taken such that vI2=vI3. The effective thermal conductivity ceff is computed with theparameters shown in the table 1. The obtained results are reported in the figure 3.a), 3.b), 3.c), 3.d) corresponding the data in the table 1.a) 1.b), 1.c), 1.d) respectively. Figure 3.1: body-centered square and hexagonal periodic three- component microstructure The thermal conductivity of inclusions and matrix: cM c1 c2 1 10 3
17 Figure 3.2: approximations and finite element results for the effective conductivity of the three-phase matrix mixtures. 3.2.2. Three-dimensional periodic three-component composites formed Figure 3.3: 3D cubic periodic three-component BCC Figure 3.4: 3D cubic periodic three-component FCC Example of calculation according to the data in the table (a) cM = 1 c1 = 3 c2 = 10 (b) cM = 3 c1 = 1 c2 = 10 (c) cM = 3 c1 = 10 c2 = 1
18 (d) cM = 10 c1 = 1 c2 = 3 Figure 3.5: Grahp of results for 3D solution 3.2.1. Effective medium approximations for the elastic moduli Building models for elastic problems according to Body-Center Cubic Example (a) KM = 4  M=2 KI2 = 1  I2=0.4 KI3 = 20  I3=12 (b) KM = 4  M=2 KI2 = 20  I2=12 KI3 = 1  I3=0.4 (c) KM = 1  M=0.4 KI2 = 4  I2=2 KI3 = 20  I3=12 (d) KM = 20  M=12 KI2 = 4  I2=2 KI3 = 10  I3=0.4
19 Fig 3.6: Graph of effective elastic modulu results CHAPTER 4. EQUIVALENT APPROXIMATION 4.1 Equivalent inclusion approach 4.1.1 Dilute solution for equivalent circle inclusions Presume one has particles of certain shapes from a particular component material, and the effective conductivity of a dilute suspension of those randomly oriented particles, having conductivity cα and volume proportion vα (α = 2,...,n) in a matrix of conductivity c1 = cM, has the form c eff  cM  v (c  cM ) D(c , cM ), v  1. (4.1)
20 In the meantime, the dilute suspension of d-dimensional spherical particles having conductivity c¯α and volume proportion vα in the matrix of the same conductivity cM has the particular expression dcM c eff  cM  v (c  cM ) ,v  1. (4.2) c  (d  1)cM Equalizing (1) and (2), one finds dcM2  (d  1)cM (c  cM ) D(c , cM ) c  . (4.3) dcM  (c  cM ) D(c , cM ) In special case the anisotropic inclusions have the elipse shape, has the particular expression cM (c  cM )(1  r ) 2 D(c , cM )  . (4.4) 2(c  r cM )(r c  cM ) and inclusions have ellipsoid (3D) shape cM  1 1 1  (4.5) D(c , cM )     . 3 c A     M c (1  A ) c B    cM (1  B ) c C    cM (1  C )   aˆ bˆ cˆ  dt aˆ bˆ cˆ  dt A  , 0 (aˆ 2  t )(t ) B   (bˆ 2 , 2 2 0   t )(t ) aˆ bˆ cˆ  dt C      (cˆ 2 ,(t )  (aˆ 2  t )(b 2  t )(c 2  t ). 2 0  t )(t ) 4.1.2 Materials with anisotropic inclusions We consider the two-component 2D square-periodic suspension of anisotropic inclusions having conductivity cI1 and cI2 in a matrix of conductivity cM D  cI 1 ,..., cId , cM  (d  1)  2 cI  cM . (4.6) 2  D  cI 1 ,..., cId , cM  4.1.3 Equivalent inclusion approach spherical inclusions (platelet, fibrous)