|
|
| (не показаны 22 промежуточные версии 4 участников) |
| Строка 1: |
Строка 1: |
| − | [[Файл:microscope_ceramic_polymer_comp.jpg|thumb|200px|Microscopic photograph the investigated ceramic-polymer-composite]] | + | [[Файл:composite.png|thumb|200px|Microscopic photograph of a TiO2-PMMA composite]] |
| − | For the research on a ceramic-polymer-composite (TiO2 in a PMMA matrix) with very high content of ceramic particles a numerical model is being investigated. To simulate the nonlinear deformation due to indentation a discrete element approach has been chosen with each ceramic particle being represented by one particle in the simulation. Several models for the interaction of the ceramic particles with one another as well as with the polymeric matrix, such as cohesion and physical bonds, are being tested using the software package EDEM and are also being compared with experimental results.
| + | [[Файл:bending2.png|thumb|200px|Bending test of simulated TiO2-PMMA specimen]] |
| | + | [[Файл:bending_results.png|thumb|200px|Load-displacement behaviour during bending using cohesive DEM]] |
| | | | |
| − | == Previous Work == | + | == Abstract == |
| − | [[Файл:old_bond_specimen.jpg|thumb|200px|Specimen of ''previous work'' with bonds representing polimeric matrix]]
| |
| − | Simulating the mechanical behaviour of the investigated material, Hertzian contacts between the particles were combined with EDEM's default bond model representing the polymeric matrix. Indentation tests revealed qualitative accordance to experimental results within the region of elastic deformation of the specimen. It was shown that a stiffer material used for the polymeric matrix, which has been modelled by bonds with larger Young's modulus, leads to a higher stiffness of the specimen.
| |
| | | | |
| − | Performing large deformations within the simulation, the created specimen deforms elastically until the stress in the bonds reaches a previously defined threshold value. Bond breakage occurs after exceeding this threshold and the interaction between the particles becomes purely Hertzian. For this region where deformation becomes plastic, the stress-strain behaviour of the model doesn't represent the real material's behaviour anymore.
| + | Interest in the simulation of material properties using numerical methods is increasing, as laboratory experiments and tests are often expensive and generally limited in flexibility. |
| | + | In this work, a DEM-approach is presented which simulates the elastic-plastic behavior of particle reinforced composites. While in the past, Discrete Element Modeling (DEM) has successfully been used to simulate elasticity and brittleness of materials, this modified approach based on the principals of Molecular Dynamics suggests a DEM-approach for modeling plastic materials. |
| | + | Simulations of a ceramic-polymer composite (TiO2-particles in a matrix PMMA) subjected to varying load conditions are presented. While Hertz-Mindlin forces are inherent in the algorithm, two variants of modeling particle attraction, particle bonds and cohesion, are applied individually. According to the results of this work, the new cohesive DEM is in good qualitative agreement with the behavior of particle-reinforced composites and therefore preferred. Ideas are presented for a modification of the computer algorithm to calibrate this cohesive DEM approach, aiming at an effective numerical method which can support experimental research on particle-reinforced composites. |
| | | | |
| − | == Representing Polymer by Cohesion of Particles == | + | == Complete Article == |
| − | [[Файл:during_indentation1.jpg|thumb|200px|Cohesion model (40,000 pts): specimen during indentation; particles colored with respect to compressive force; specimen has periodic boundaries in X- and Y-direction]] | + | [[Медиа:CohesiveDEM_DA_MSzelwis_Dez2011.pdf]] |
| − | [[Файл:after_indentation_1.jpg|thumb|200px|Cohesion model (40,000 pts): surface of specimen after indentation]]
| |
| − | [[Файл:indentationforce40k.png|thumb|200px|Cohesion model (40,000 pts): indentation-force over time]]
| |
| − | [[Файл:indentationforce10k.png|thumb|200px|Cohesion model (10,000 pts): indentation-force over time]]
| |
| − | The main problem of using bonds simulating the material's plastic behaviour is that in this case, particle rearrangement always goes along with the breakage of bonds and therefore with a drastic change of the simulated material's properties. Plastic deformation is generally associated with so-called dislocation movement where the atoms of the material rearrange during deformation due to imperfections in the crystal structure. During this process energy dissipates but the mechanical properties of the material remain more or less unchanged.
| |
| − | Therefore, it would be desirable for a model which is supposed to represent the investigated ceramic-polymer-composite, to contain consistent particle interactions before and after the relocation of particles.
| |
| | | | |
| − | The Cohesion of primary particles is an implemented contact model in EDEM in which the simulated attractive force between particles does not change with thier relocation. Also, it represents how we generally understand the constitution of substances on an atomic scale. Consequently, it might be suitable to model the mechanical properties of the investigated material including elastic as well as plastic behaviour.
| + | == См. также == |
| − | Nevertheless, it is questionable if this approach can predict the properties of the material when used on a very large scale where each particle in the simulation represents one ceramic particle embedded in a matrix of polymer.
| |
| | | | |
| − | === Details of Simulation (DRAFT) ===
| + | *[[Michael Szelwis]] |
| − | To keep the calculations as simple as possible first simulations were run using the default Hertz-Mindlin contact model as well as the model for particle cohesion which are provided within EDEM. In this case the '''Hertz-Mindlin-Force''' between two equivalent particles in contact gets calculated as follows:
| + | *[[Гамбургский проект]] |
| − | :<math>F_{HM}(\delta_n) = \frac{2}{3}\times\frac{E}{(1-\nu^2)}\times\sqrt{R/2}\times\delta_n^{\frac{3}{2}}</math>
| |
| − | with <math>E</math> and <math>\nu</math> being the Young's Modulus and the Poisson's Ratio of the particles, <math>R</math> being thier radius and <math>\delta_n</math> being the normal overlap assuming the particles remain spherical when in contact. Considering this assumption, the '''Particle Cohesion''' calculates:
| |
| − | :<math>F_{C}(\delta_n) = k\times\pi\times(R^2-(R-\delta_n/2)^2)</math>
| |
| − | with <math>k</math> being a constant parameter of the cohesive force (also referred to as energy density).
| |
| − | | |
| − | This leads to the '''Total Force''' between two interacting particles of
| |
| − | :<math>F_{Total}(\delta_n) = -\frac{1}{4}\pi k\delta_n^2+\frac{\sqrt{2R}\times E}{3(1-\nu^2)}\delta_n^{\frac{3}{2}}+\pi k R\delta_n</math>
| |
| − | If attractive and repulsive forces are of the same magnitude the particle contact reaches an equilibrium state with an equilibrium overlap <math>\delta_{n,eq}</math>. Besides the trivial solution <math>\delta_{n,eq1}=0</math> there might be found a second solution for which the particle contact becomes stable. Unfortunately, this second equilibrium distance <math>\delta_{n,eq2}</math> can not be adjusted to <math>0<\delta_{n,eq2}<R</math> which would be suitable for the intended simulations. Besides, this is the only range for which EDEM calculates the particle forces which means that with the indicated formulas, an equilibrium distance of two interacting particles can never be reached.
| |
| − | | |
| − | While a stable configuration for indentation experiments can still be reached (high friction between particles prevents them from collapsing) the simulation of a bending experiment will not be possible. To eliminate this problem a modified cohesion model with a larger radius for the calculation of cohesive forces than the particle radius will be considered in future simulations (see below).
| |
| − | | |
| − | Further characteristics of simulation:
| |
| − | | |
| − | * 10,000/40,000 particles
| |
| − | * Periodic boundaries
| |
| − | * Quasi-2D
| |
| − | * Close packing (0.63)
| |
| − | | |
| − | === Results ===
| |
| − | Simulations revealed the possibility to qualitatively predict a material's response due to indentation by approximating the structure of the investigated material using discrete particles applying Hertzian contact laws, linear cohesion as well as linear friction.
| |
| − | Furthermore, this concept is comparable to how deformation processes are generally explained on an atomic scale which emphasizes the consideration of this approach. The fact that the particles in the simulation being of magnitues larger than real atoms does not seem to have an impact on the qualitative results even though it involves a certain coarseness of the plastic area in which particle relocation occurs (see force-displacement curves of cohesion models with 10,000 & 40,000 particles).
| |
| − | | |
| − | == Next Steps within this Work ==
| |
| − | * Quantitative comparison of the results with experimental data
| |
| − | * Realizing model with varying relevant radius for calculations of cohesive- and contact forces (at the moment it is impossible to reach equilibrium state)
| |
| − | * Performing bending test (with new model of varying radii)
| |
| − | * Ensuring the results do not change quantitively by increasing the number of used particles | |
| − | * Examining the suitability using the cohesive model for the investigation of hirarchic materials | |
| | | | |
| | | | |
| | [[Category: Студенческие проекты]] | | [[Category: Студенческие проекты]] |
| | + | [[Category: Механика дискретных сред]] |
Interest in the simulation of material properties using numerical methods is increasing, as laboratory experiments and tests are often expensive and generally limited in flexibility.
In this work, a DEM-approach is presented which simulates the elastic-plastic behavior of particle reinforced composites. While in the past, Discrete Element Modeling (DEM) has successfully been used to simulate elasticity and brittleness of materials, this modified approach based on the principals of Molecular Dynamics suggests a DEM-approach for modeling plastic materials.
Simulations of a ceramic-polymer composite (TiO2-particles in a matrix PMMA) subjected to varying load conditions are presented. While Hertz-Mindlin forces are inherent in the algorithm, two variants of modeling particle attraction, particle bonds and cohesion, are applied individually. According to the results of this work, the new cohesive DEM is in good qualitative agreement with the behavior of particle-reinforced composites and therefore preferred. Ideas are presented for a modification of the computer algorithm to calibrate this cohesive DEM approach, aiming at an effective numerical method which can support experimental research on particle-reinforced composites.
