Diversos materiais de interesse tecnológico podem ser considerados como materiais heterogêneos, onde sua natureza aleatória deve ser considerada para representar corretamente o comportamento não-linear. A avaliação quantitativa do dano em materiais sujeitos a estados de tensão ou deformação tem grande importância devido ao caráter crítico desse fenômeno, que pode crescer de forma abrupta e resultar em uma falha catastrófica da estrutura. Em estudos anteriores, Carpinteri e seus colaboradores apresentaram diferentes aspectos relacionados à caracterização do dano em materiais heterogêneos. Três desses aspectos são discutidos neste trabalho: (i) a proposta do número de fragilidade como medida da fragilidade da estrutura em análise; (ii) a evolução da dimensão fractal em que o processo de dano é desenvolvido; (iii) e os índices globais obtidos pela análise de Emissão Acústica (EA). Neste trabalho, uma versão do método dos elementos discretos formado por barras é utilizada para explorar esses aspectos. Espécimes de material quase-frágil são simulados e, quando possível, os resultados numéricos são comparados com resultados experimentais. Além disso, a discussão dos resultados obtidos é realizada visando compreender o comportamento desse tipo de material, como por exemplo, o fato de que estruturas com diferentes dimensões, porém com o mesmo número de fragilidade, têm comportamentos similares. Ademais, o método numérico empregado é apresentado como uma ferramenta viável para complementar a informação obtida em ensaios experimentais na avaliação do processo de dano. Pelo sinal de EA obtido através do método numérico empregado, obtiveram-se parâmetros para a classificação do tipo de falha e para o cálculo do tensor momento, onde comparando os resultados obtidos pelo modelo numérico e pela EA se observam comportamentos coerentes. / Several materials with technological interest can be considered as heterogeneous materials, where their random nature must be accounted to correctly represent the nonlinear behavior. The quantitative evaluation of damage in materials subjected to stress or strain states have great importance due to the critical nature of this phenomenon, which abruptly can rise to catastrophic failure. In previous studies, Carpinteri and his coworkers have presented different aspects of the damage process characterization in heterogeneous materials. Three of these aspects are discussed in this work: (i) the brittleness number proposal to measure the brittleness level of the analyzed structure, (ii) the relationship of the fractal dimension in which the damage process is developed, (iii) and the global indexes obtained by the Acoustic Emission (AE) analysis. In the present work, a version of discrete element method formed by bars is used to explore these concepts. A set of quasi-brittle material specimen is simulated and, when possible, the numerical results are compared with experimental ones. Moreover, a discussion of the obtained results is carried out aiming to understand the behavior of this kind of material, for instance, the fact that structures with different dimensions, but with the same brittleness number, have similar behavior. In addition, the numerical method is presented as a viable tool to complement information from experimental test on the damage process. From the AE signal obtained by the numerical method, parameters to classify the type of crack and for calculating the moment tensor were obtained, where consistent behaviors are observed by comparing the results of the numerical model and the AE.

Lors de la première transformation du bois en scierie, les grumes sont surfacées ou équarries par des têtes de fraisages appelées slabber ou canter. Sous leur action, le copeau de bois est fragmenté en plaquettes dont la valorisation est un enjeu majeur de la filière. Débouché le plus rémunérateur de cette ressource, l’industrie de la pâte à papier impose des critères dimensionnels aux plaquettes approvisionnées notamment concernant leur épaisseur. L’objectif de l’étude est donc d’améliorer la compréhension des mécanismes mis en jeu par la coupe du bois, dans une configuration simplifiée de coupe orthogonale, afin d’optimiser la granulométrie des plaquettes produites. Une campagne expérimentale de coupe sur machine-outil à commande numérique a été réalisée sur du hêtre vert afin de déterminer les mécanismes principaux actionnés. Un modèle mécanique simplifié du comportement dynamique du bois vert est déterminé, ceci afin de développer un modèle numérique du bois vert par la Méthode des Eléments Discrets (DEM) en vue de simuler sa coupe. Une étude préliminaire pour déterminer la capacité de la méthode à modéliser à l’échelle mésoscopique des milieux fortement orthotropes a été réalisée. Sa capacité et ses limites démontrées, la démarche de calibration du modèle numérique a été élaborée et le modèle sollicité en configuration de coupe orthogonale numérique. Les premières simulations présentent des résultats encourageants. / When entering sawmills, logs are faced into cants by the mean of chipper-canters. During this machining, the ribbon produced is split into small chips whose proper valorization is a high economic stake for the industry. The paper maker industry, which is the most worthwhile chips supplier, is strongly concerned by dimensional criterions of the chips for their process optimization, especially toward their thickness. The objective of this work is to improve the comprehension of cutting and fragmentation mechanisms, in a simplified orthogonal cutting configuration, to provide cutting rules to optimize the produced chips geometry. An experimental campaign dedicated to green beech cutting on a computer numerical command machining center is done to study those mechanisms. A simple mechanical model is derived, in order to develop a Discrete Element Method (DEM) model of the material to simulate cutting operations. A preliminary study dedicated to prove the capability of DEM to model a wood-like orthotropic continuous media is presented. The numerical model is then calibrated and cutting simulations are designed to copy the experimental conditions. First results are encouraging.

Asphalt mixture compaction is an important procedure of asphalt mixture construction and can significantly affect the performance of asphalt pavement. Many laboratory compaction methods (or devices), have been developed to study the asphalt mixture compaction. Nevertheless, the whole process from the selection of aggregate to laboratory compaction is still time-consuming and requires significant human and material resources. In order to better understand asphalt mixture compaction, some researchers began to use finite element method (FEM) to study and analyze mixture compaction. However, FEM is a continuum approach and lacks the ability to take into account the slippage and interlocking of aggregates during compaction. Discrete Element Method (DEM) is a discontinuum analysis method, which can simulate the deformation process of joint systems or discrete particle assembly under quasi-static and dynamic condition. Therefore, it can overcome the shortcomings of FEM and is a more effective tool than FEM to simulate asphalt mixture compaction. In this study, an open source 3D DEM code implemented with the C++ programming language was modified and applied to simulate the compaction of hot-mix asphalt (HMA). A viscoelastic contact model was developed in the DEM code and was verified through comparison with well established analytical solutions. The input parameters of the newly developed contact model were obtained through nonlinear regression analysis of dynamic modulus test results. Two commonly used compaction methods (Superpave gyratory compaction and asphalt vibratory compaction) and one linear kneading compaction based on APA machine were simulated using the DEM code, and the DEM compaction models were verified through the comparison between the DEM predicted results and the laboratory measured test results. The air voids distribution within the asphalt specimens was also analyzed by post processing virtual DEM compaction digital specimens and the level of heterogeneity of the air void distribution within the specimens in the vertical and lateral directions was studied. The DEM simulation results in this study were in a relatively good agreement with the experimental data and previous research results, which demonstrates that the DEM is a feasible method to simulate asphalt mixture compaction under different loading conditions and, with further research, it could be a potentially helpful tool for asphalt mix design by reducing the number of physical compactions in the laboratory.

In this thesis, the effects of size and stress on permeability, deformability and strength of fractured rock masses are investigated. A comparison study was carried out to examine the effects of considering, or not considering, the correlation between distributions of fracture apertures and fracture trace lengths on the hydro-mechanical behavior of fractured rocks. The basic concepts used are the fundamental principles of the general theory of elasticity, Representative Elementary Volume (REV), the tensor of equivalent permeability, and the strength criteria of the fractured rocks. Due to the size and stress dependence of the hydro-mechanical properties of rock fractures, the overall effective (or equivalent) hydro-mechanical properties of the fractured rocks are also size and stress-dependent. However, such dependence cannot be readily investigated in laboratory using small samples, and so numerical modeling becomes a necessary tool for estimating their impacts. In this study, a closed-form relation is established for representing the correlation between a truncated lognormal distribution of fracture apertures and a truncated power law distribution of trace lengths, as obtained from field mapping. Furthermore, a new nonlinear algorithm is developed for predicting the relationship between normal stress and normal displacement of fractures, based on the Bandis model and the correlation between aperture and length. A large number of stochastic Discrete Fracture Network (DFN) models of varying sizes were extracted from some generated large-sized parent realizations based on a realistic fracture system description from a site investigation programme at Sellafield, UK, for calculating the REV of hydro-mechanical properties of fractured rocks. Rotated DFN models were also generated and used for evaluation of the distributions of directional permeabilities, such that tensors of equivalent permeability could be established based on stochastically established REVs. The stress-dependence of the permeability and the stress-displacement behaviour were then investigated using models of REV sizes. The Discrete Element Method (DEM) was used for numerical simulation of the fluid flow, deformability properties and mechanical strength behavior of fractured rocks. The results show significant scale-dependency of rock permeability, deformability and strength, and its variation when the correlation between aperture and trace length of fractures are concerned, with the overall permeability and deformability more controlled by dominating fractures with larger apertures and higher transmissivity and deformability, compared with fracture network models having uniform aperture. As the second moment of aperture distribution increases, a fractured rock mass shows more discrete behavior and an REV is established in smaller value of second moment with much larger model size, compared with the models with uniform fracture aperture. When the fracture aperture pattern is more scattered, the overall permeability, Young’s modulus and mechanical strength change significantly. The effect of stress on permeability and fluid flow patterns in fractured rock is significant and can lead to the existence or non-existence of a permeability tensor. Stress changes the fluid flow patterns and can cause significant channeling and the permeability tensor, and REV may be destroyed or re-established at different applied stress conditions. With an increase in the confining stress on the DEM models, the strength is increased. Compared with the Hoek-Brown criterion, the Mohr-Coulomb strength envelope provides a better fit to the results of numerical biaxial compression tests, with significant changes of the strength characteristic parameters occurring when the second moment of the aperture distribution is increased. / QC 20100702

Under loading conditions found in many geotechnical structures, it is common to observe failure in zones of high localized strain called shear bands. Existing models predict these localizations, but provide little insight into the micromechanics within the shear bands. This research captures the variation in microstructure inside and outside of shear bands that were formed in laboratory plane strain and two-dimensional discrete element method (DEM) biaxial compression experiments. Plane strain compression tests were conducted on dry specimens of Ottawa 20-30 sand to calibrate the device, assess global response repeatability, and develop a procedure to quantitatively define the onset of localization. A new methodology was employed to quantify and correct for the additional stresses imparted by the confining membrane in the vicinity of the shear band. Unsheared and sheared specimens of varying dilatancy were solidified using a two-stage resin impregnation procedure. DEM tests were performed using an innovative servo-controlled flexible lateral confinement algorithm to provide additional insights into laboratory results. The solidified specimens were sectioned and the resulting surfaces prepared for microstructure observation using bright field microscopy and morphological analysis. Local void ratio distributions and their statistical properties were determined and compared. Microstructural parameters for subregions in a grid pattern and along predefined inclined zones were also calculated. Virtual surfaces parallel to the shear band were identified and their roughnesses assessed. Similar calculations were performed on the DEM simulations at varying strain levels to characterize the evolution of microstructure with increasing strain. The various observations showed that the mean, standard deviation, and entropy of the local void ratio distributions all increased with increasing strain levels, particularly within regions of high local strains. These results indicate that disorder increases within a shear band and that the soil within the shear band does not adhere to the classical concept of critical state, but reaches a terminal void ratio that is largely a function of initial void ratio. Furthermore, there appears to be a transition zone between the far field and the fully formed shear block, as opposed to an abrupt delineation as traditionally inferred.

Integrated theoretical and experimental analysis is carried out in this work to investigate the fundamental failure mechanisms and flow patterns involved in the process of fluid injection into dense granular media. The experimental work is conducted with aqueous glycerin solutions, utilizing a novel setup based on a Hele-Shaw cell filled with dense dry sand. The two dimensional nature of the setup allows direct visualization and imaging analysis of the real-time fluid and grain kinematics. The experimental results reveal that the fluid flow patterns show a transition from simple radial flow to a ramified morphology while the granular media behaviors change from that of rigid porous media to localized failure that lead to development of fluid channels. Based on the failure/flow patterns, four distinct failure/flow regimes can be identified, namely, (i) a simple radial flow regime, (ii) an infiltration-dominated regime, (iii) a grain displacement-dominated regime, and (iv) a viscous fingering-dominated regime. These distinct failure/flow regimes emerge as a result of competition among various energy dissipation mechanisms, namely, viscous dissipation through infiltration, dissipation due to grain displacements, and viscous dissipation through flow in thin channels and can be classified based on the characteristic times associated with fluid injection, hydromechanical coupling and viscoelastoplasticity. The injection process is also analyzed numerically using the discrete element method (DEM) coupled with two fluid flow scheme, a fixed coarse grid scheme based on computational fluid dynamics (CFD) and a pore network modeling scheme. The numerical results from the two complementary methods reproduce phenomena consistent with the experimental observations and justify the concept of associating the displacement regimes with the partition among energy dissipation mechanisms. The research in this work, though fundamental in nature, will have direct impacts on many engineering problems in civil, environmental and petroleum engineering such as ground improvement, environmental remediation and reservoir stimulation.

Modeling and simulation of granular matter has important applications in both natural science and industry. One widely used method is the discrete element method (DEM). It can be used for simulating granular matter in the gaseous, liquid as well as solid regime whereas alternative methods are in general applicable to only one. Discrete element analysis of large systems is, however, limited by long computational time. A number of solutions to radically improve the computational efficiency of DEM simulations are developed and analysed. These include treating the material as a nonsmooth dynamical system and methods for reducing the computational effort for solving the complementarity problem that arise from implicit treatment of the contact laws. This allow for large time-step integration and ultimately more and faster simulation studies or analysis of more complex systems. Acceleration methods that can reduce the computational complexity and degrees of freedom have been invented. These solutions are investigated in numerical experiments, validated using experimental data and applied for design exploration of iron ore pelletising systems. / <p>This work has been generously supported by Algoryx Simulation, LKAB (dnr 223-</p><p>2442-09), Umeå University and VINNOVA (2014-01901).</p>

Tyrimų sritis ir darbo aktualumas. Kuriant modernias įvairios paskirties mechanines sistemas, technologijas ir įrangą, svarbiomis tampa jas sudarančios medžiagos. Savaime suprantama, kad žinomos ir naujai kuriamos medžiagos dabar kur kas išsamiau nagrinėjamos daugelyje mokslo šakų, įskaitant ir me-džiagų mechaniką. Visos medžiagos mezo- ir mikrostruktūros požiūriu yra ne-vienalytės. Jų mikroskopinės savybės skirtingos, lyginant su įprastu kontinuu-mu. Medžiagų savybėms tirti dažniausiai taikomi eksperimentiniai metodai. Eksperimentiniais metodais ištirti medžiagos struktūras ir jose vykstančius procesus ir įvertinti tam tikras jų savybes labai brangu. Tai viena priežasčių, kodėl skaitinis modeliavimas tampa realia tyrimų alternatyva. Skaitinį eksperi-mentą galima kartoti daug kartų, valdant bandinio parametrus, išlaikant tas pa-čias sąlygas, ir stebėti reiškiniui būdingus rodiklius visame tūryje. Šiuolaikiniai modeliavimo metodai yra kompleksiniai. Jie jungia fenome-nologines ir statistines idėjas, o matematiniai modeliai sudaromi taikant konti-nuumo mechanikos ir jų diskrečiųjų modelių bei molekulinės dinamikos pri-klausomybes. Diskrečiųjų elementų metodas (DEM) taip pat priskiriamas šiuo-laikinių metodų kategorijai. Jis skirtas kontaktuojančių dalelių sistemų dinami-niam modeliavimui. Kintanti dalelių sistemos topologija – būdingas metodo požymis. Pastaruoju metu DEM jau taikomas kontinuumui modeliuoti ir praktikoje aktualiems irimo uždaviniams spręsti. Reikia pastebėti... [toliau žr. visą tekstą] / Research area and topicality of the work. Mechanical properties and their evolution under loading are the most significant factors for the development of various mechanical structures, technologies and equipment. It seems to be natu-ral that deeper understanding of the behaviour of existing and design of new materials presents a challenge in different research areas. It should be noted, that all the materials are heterogeneous in meso- and micro- scales. They exhibit essential differences, compared to the macroscopic continuum behaviour. Basically, both experimental and numerical simulation methods are extensively applied for investigation purposes. Experimental techniques, capable of giving a realistic view of the inside of the material and extracting the real data, are very expensive. Therefore, the nu-merical simulation tools are extensively used as an alternative for investigation purposes. They have considerable advantages allowing the reproduction of multiple experiments and providing comprehensive data about ongoing phe-nomena. Recently, numerical technologies have become highly multidisciplinary subjects. They comprise phenomenological and statistical ideas, while mathe-matical models employ the relations of continuum mechanics, classical discre-tization methods and molecular dynamics. The Discrete Element Method (DEM) is one of new methods. It is aimed at simulating the dynamic behaviour of the contacting particles. Variable topology of the system of particles is an... [to full text]

Baigiamajame magistro darbe nagrinėjama dalelių mišinio dinaminė elgsena kvazistatiniame pusiausvyros būvyje. Pagrindinis dėmesys skiriamas mišinio makro- ir mikrobūsenai charakterizuoti. Darbo tikslas – ištirti betono užpildams būdingomis fizinėmis-mechaninėmis savybėmis pasižyminčio sferinių dalelių mišinio struktūrą ir elgseną fliuktuojančioje aplinkoje kvazistatinio pusiausvyros būvio sąlygomis. Dalelių mišinio elgsenos dinaminis modeliavimas atliktas diskrečiųjų elementų metodu (DEM). Modeliuoti viendispersis ir daugiadispersis mišiniai. Pagal gautus rezultatus, mišinio makrobūsena charakterizuojama struktūros užimama tūrio dalimi (tankumu), dalelių kontaktinėmis jėgomis bei koordinacijos skaičiumi. Atskirų mišinio dalelių mikrobūsena charakterizuojama taikant unikalų geometrinį metodą, sukurtą pagal Hamiltono mechanikos, bendrosios reliatyvumo teorijos ir kvantinės gravitacijos teorijų principus. Metodui realizuoti parengtas kompiuterinis algoritmas patvirtino iškeltas hipotezes. / In the Master Thesis an investigation of dynamical behavior of the mixture of particles under quasistatic equilibrium conditions was performed. It focuses on a characterization of the macro- and microstate of the mixture. The purpose of the research is to explore structure and behavior of the mixture of spherical particles, characterized by the typical physical-mechanical properties of concrete aggregates, in a fluctuating environment under quasistatic equilibrium conditions. Dynamic simulation of the behavior of the particles was performed via the discrete element method (DEM). According to the results obtained, the macrostate of the mixture was investigated in view of a volume fraction, contact forces of particles and a coordination number. A characterization of the microstate of distinct particle was done using the unique geometrical method which had been created according to Hamiltonian mechanics, the theory of relativity and the principles of quantum gravity. The developed computer algorithm allowed to confirm the hypotheses of the method.

Recently, powders of the size d (0.1 μm < d < 10 μm) have been referred to ultrafine particles. The particle shape considered is assumed to be a sphere of the diameter d. The handling of powders is of great importance for processing of pharmaceuticals, cement, chemicals and other products. Most of these technological processes involve powder compaction, storage, transportation, mixing, etc, therefore, understanding of the fundamentals of particles interaction behaviour is very essential in the design of machines and equipment as well as in powder technology, cleaning of environment and other areas. The dynamic behaviour of particulate systems is very complicated due to the complex interactions between individual particles and their interaction with the surroundings. Understanding the underlying mechanisms can be effectively achieved via particle scale research. The problem of a normal contact may be resolved in a number of ways. In spite of huge progress in experimental techniques, direct lab tests with individual particles are still rather time-consuming and expensive. The interaction of particles as solid bodies is actually a classical problem of contact mechanics. In the case of ultrafine particles, the reduction of the particle size shifts the contact zones into the nanoscale or subnanoscale. Thus, steadily increasing contribution of adhesion has to be considered in the development of the physically correct constitutive models and numerical tools. Consequently, it may... [to full text] / Ultrasmulkios dalelės yra šiuolaikinės chemijos, farmacijos, maisto ir kitų pramonės šakų produktų sudėtinė dalis. Tiriant pramoninius technologinius procesus, neišvengiamai reikalingos teorinės žinios apie ultrasmulkių dalelių elgseną. Išsamus supratimas įmanomas tik atlikus įvairius tyrimus. Pastaruoju metu milteliai, klasifikuojami kaip ultrasmulkios (0,1 < d < 10 μm) dalelės, imti plačiai naudoti pramoniniuose procesuose, todėl suprasti ultrasmulkių dalelių elgsenos fundamentalumą miltelių technologijoje yra labai svarbu. Ultrasmulki dalelė yra itin maža, todėl su ja atlikti fizinį eksperimentą, kuris reikalauja specialios įrangos bei žinių, labai sunku. Tokiu atveju dažniausiai naudojamas skaitinis eksperimentas, kurį galima atlikti virtualiai. Skaitinio eksperimento metu yra tiriamos dinaminės ultrasmulkios dalelės savybės bei sprendžiamas dinaminis uždavinys. Taikant skaitinius modelius bei dalelės judėjimą aprašančias jėgų lygtis, naudojami sąveikos modeliai, apimantys adhezinę, klampią, tamprią bei tampriai plastinę sąveikas. Mikroskopinis adhezinės sąveikos modeliavimas – aktualus mechanikos mokslo uždavinys. Taikant sąveikos modelius, svarbu pritaikyti ir diskrečiųjų elementų metodą, kadangi, norint aprašyti dalelių elgseną, visų pirma reikia su-vokti ir aprašyti dalelės modelį. Dalelės elgsenos skaitiniam modeliavimui siūlomi teoriniai modeliai leidžia tirti dalelės sąveiką su dalele ar tampria puserdve bei sąveikos dinamiką. Šie modeliai galėtų būti pritaikyti... [toliau žr. visą tekstą]