A DEM study of silo discharge of a cohesive solid

Bulk handling of powders and granular solids is common in many industries and often gives rise to handling difficulties especially when the material exhibits complex cohesive behaviour. For example, high storage stresses in a silo can lead to high cohesive strength of the stored solid, which may in turn cause blockages such as ratholing or arching near the outlet during discharge. This paper presents a Discrete Element Method study of discharge of a granular solid with varying levels of cohesion from a flat-bottomed silo. The DEM simulations were conducted using the commercial EDEM code with a recently developed DEM contact model for cohesive solids implemented through an API. The contact model is based on an elasto-plastic contact with adhesion and uses hysteretic non-linear loading and unloading paths to model the elastic-plastic contact deformation. The adhesion parameter is a function of the maximum contact overlap. The model has been shown to be able to predict the stress history dependent behaviour depicted by a flow function of the material. The effects of cohesion on the discharge rate and flow pattern in the silo are investigated. The predicted discharge rates are compared for the varying levels of cohesion and the effect of adhesion is evaluated. The ability of the contact model to qualitatively predict the phenomena that are present in the discharge of a silo has been shown with the salient feature of mixed flow from a flat bottomed hopper identified in the simulation.

Effect of roller leveling on the bending behaviour of aged steel strip

The common grades of steel used in roll forming are: hot rolled carbon steel, high strength low alloy and recovery annealed cold rolled sheet. These steels are prone to ageing and are often skin passed and/or roller leveled to eliminate ageing as it can lead to problems in forming. In roll forming, shape defects such as bow, twist and camber are considered to be related to very small plastic strains in the longitudinal direction and hence knowledge of the material properties in the elastic plastic transition range is necessary if the process is to be modelled accurately. Previous studies with aluminium have indicated that skin pass rolling can lead to residual stresses in the strip. In this work, the study was extended to aged carbon steel and to the effect of roller leveling on both aged material and strip that had been given a light cold rolling to simulate a skin pass treatment. The results suggest that roller leveling reduced the magnitude of residual stresses resulting from skin pass rolling.

The significant differences observed between tensile and bending test results, at and near, the elastic plastic transition reinforces the need to consider bending properties when assessing the effect of prior processing on strip for roll forming.

Effect of tensile pre-strain on bending and unloading of automotive steels

In recent years, advanced high strength steels (AHSS) have been used in a wide range of automotive applications; they may have property variations through the thickness and the properties may also be dependent of prior processing including pre-straining. In order to model forming processes precisely using, for example, finite element analysis, it is important that material input data should adequately reflect these effects. It is known that shape defects in roll forming are related to small strains in material that has undergone prior deformation in a different strain path. Much research has already been performed on the change in the Young’s Modulus once a steel sheet has been plastically deformed,however many of these tests have only been conducted using tensile testing, and therefore may not take into account differences in compressive and tensile unloading. This research investigates the effect of tensile pre-straining on bending behaviour for various types of material;in bending, one half of the sheet will load and unload in compression and hence experience deformation under a reversed stress. Four different materials were pre-strained in tension with 1%, 3%, 7%, 11% and 25% elongation. Using a free bending test, moment curvature diagrams were obtained for bending and unloading. The results showed that the characteristics of the moment curvature diagram depended on the degree of pre-straining; more highly strained samples showed an earlier elastic-plastic transformation and a decreased Young's Modulus during unloading. This was compared to previous literature results using only tensile tests. Our results could influence the modeling of springback in low tension sheet operations, such as roll forming.

Análise não linear física e geométrica de pórticos espaciais

Apresenta-se uma formulação do tipo incrementaliterativa destinada a análise não linear de pórticos espaciais. Considera-se os efeitos não lineares introduzidos pelas mudanças de configuração geométrica da estrutura e também pela combinação destes efeitos com aqueles inerentes ao comportamento plástico exibido pelo material. As relações cinemáticas empregadas permitem a consideração de deslocamentos arbitrariamente grandes, acompanhadas de pequenas deformações . A modelagem do comportamento plástico do material é efetuada através do conceito de rótula plástica, estabelecido a partir de um critério de plastificação generalizado. Adota-se uma matriz de rigidez geométrica de barra baseada em momentos semitangenciais. Para elementos com extremos plastificados, é deduzida uma matriz de rigidez elasto-plástica. Emprega-se um método numérico do tipo incremental-iterativo, que utiliza como condição básica de controle da análise a constância do trabalho realizado pelos incrementos de cargas, em cada passo incremental (Método de Controle por Trabalho).A formulação permite uma descricão completa do desempenho mecânico da estrutura, inclusive em estágio de deformação pós-crítico em que ocorre regressão do carregamento com aumento de deslocamentos, ou vice-versa. A formulação foi implementada em um programa computacional elaborado em linguagem FORTRAN. Vários exemplos numéricos são apresentados para mostrar a eficiência das procedimentos propostos.

Ultrasonic material characterization using diffraction-free PVDF receivers

This work describes the use of a large aperture PVDF receiver in the measurement of density of liquids and elastic constants of composite materials. The density measurement of several liquids is obtained with the accuracy of less than 0.2% using a conventional NDT emitter transducer and a 70-mm diameter, 52-μm P(VDF-TrFE) membrane with gold electrodes. The determination of the elastic constants of composite materials is based in the measurement of phase velocity. It is shown that the diffraction can lead to errors around 1% in the velocity measurement when using a pair of ultrasonic transducers (1MHz and 19mm diameter) operating in transmission-reception mode separated by a distance of 100 mm. This effect is negligible when using a pair of 10-MHz transducers. On the other hand, the dispersion at 10 MHz can result in errors of about 0.5%, measuring the velocity in composite materials. The use of an 80-mm diameter, 52-μm thick PVDF membrane receiver allows measuring the phase velocity without the diffraction effects.

Fractografia quantitativa: relações entre tenacidade e o comportamento fractal em ligas metálicas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Estudo do comportamento dos aços ferramentas soldados em fadiga

Sabendo-se que a fadiga é uma redução gradual da capacidade de carga do componente pela ruptura lenta e gradual do material. E que este defeito decorre do avanço infinitesimal de microtrincas, que se formam no interior do material, imperceptível a olho nu, como também é notório que a presença dos elementos de liga nos aços, propicia alterações nas propriedades metalurgias e mecânicas no material, aplicado à obra. Por outro lado, ao ser submetido a processamentos dos mais diversos, os mesmos, deformam nas zonas elásticas, plásticas e ao fim rompem-se. Ressaltando-se o objetivo deste estudo, destina se abordar as uniões soldadas, nos estágios em que ficam sujeitas a altas temperaturas e resfriam sem controle, até a temperatura ambiente. Ao fim, são solicitados por carregamentos cíclicos constantes ou alternados. Nesta particularidade, espera-se detectar mudanças estruturais profundas na Zona Termicamente Afetada – ZTA, em razão do superaquecimento sofrido, tanto na proximidade do ponto central da poça de fusão como na região localizada na vizinhança. Como não dispomos de parâmetros suficientes e necessários ao controle destas alterações, propomos analisar a ZTA da junta soldada, com a finalidade de avaliar o comportamento metalúrgico e suas implicações causadas pelas interações gás-metal. Analisar os efeitos resultantes do triangulo formado pelo material base aço AISI/SAE 4340, SAE 1020 chanfrados em V e o eletrodo revestido E-6013. A análise da estrutura será realizada pelas técnicas de ensaio metalográfico usando o método da microscopia óptica – MO o qual é de ampla difusão nas comunicações técnicas e cientificas, através das quais, distinguirão as transformações multifásicas. Transformações estas, distintas pelas transformações das austenitas diretas em martensitas sem passar pelas ferritas e perlitas. Com estas identidades, o metalógrafo pode caracterizar e predizer as reações futuras das estruturas mediante a análise das solicitações, a que possam estar submetidas. Para analisar o comportamento em fadiga dos materiais em questão, adota-se como parâmetro auxiliar, o ensaio do pêndulo de Charpy. Portanto, pelos resultados obtidos, conclui-se que as estruturas são comprometidas por tensões internas em conseqüência fragilizam, mas o procedimento de soldagem pode ser adotado, seguido de recozimento para alivio de tensões, se não fizer trincas.

Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Modelos de elementos finitos explícitos para explosiones en estructuras reticuladas de hormigón armado : aplicaciones al estudio del colapso de edificios

La frecuencia con la que se producen explosiones sobre edificios, ya sean accidentales o intencionadas, es reducida, pero sus efectos pueden ser catastróficos. Es deseable poder predecir de forma suficientemente precisa las consecuencias de estas acciones dinámicas sobre edificaciones civiles, entre las cuales las estructuras reticuladas de hormigón armado son una tipología habitual. En esta tesis doctoral se exploran distintas opciones prácticas para el modelado y cálculo numérico por ordenador de estructuras de hormigón armado sometidas a explosiones. Se emplean modelos numéricos de elementos finitos con integración explícita en el tiempo, que demuestran su capacidad efectiva para simular los fenómenos físicos y estructurales de dinámica rápida y altamente no lineales que suceden, pudiendo predecir los daños ocasionados tanto por la propia explosión como por el posible colapso progresivo de la estructura. El trabajo se ha llevado a cabo empleando el código comercial de elementos finitos LS-DYNA (Hallquist, 2006), desarrollando en el mismo distintos tipos de modelos de cálculo que se pueden clasificar en dos tipos principales: 1) modelos basados en elementos finitos de continuo, en los que se discretiza directamente el medio continuo mediante grados de libertad nodales de desplazamientos; 2) modelos basados en elementos finitos estructurales, mediante vigas y láminas, que incluyen hipótesis cinemáticas para elementos lineales o superficiales. Estos modelos se desarrollan y discuten a varios niveles distintos: 1) a nivel del comportamiento de los materiales, 2) a nivel de la respuesta de elementos estructurales tales como columnas, vigas o losas, y 3) a nivel de la respuesta de edificios completos o de partes significativas de los mismos. Se desarrollan modelos de elementos finitos de continuo 3D muy detallados que modelizan el hormigón en masa y el acero de armado de forma segregada. El hormigón se representa con un modelo constitutivo del hormigón CSCM (Murray et al., 2007), que tiene un comportamiento inelástico, con diferente respuesta a tracción y compresión, endurecimiento, daño por fisuración y compresión, y rotura. El acero se representa con un modelo constitutivo elastoplástico bilineal con rotura. Se modeliza la geometría precisa del hormigón mediante elementos finitos de continuo 3D y cada una de las barras de armado mediante elementos finitos tipo viga, con su posición exacta dentro de la masa de hormigón. La malla del modelo se construye mediante la superposición de los elementos de continuo de hormigón y los elementos tipo viga de las armaduras segregadas, que son obligadas a seguir la deformación del sólido en cada punto mediante un algoritmo de penalización, simulando así el comportamiento del hormigón armado. En este trabajo se denominarán a estos modelos simplificadamente como modelos de EF de continuo. Con estos modelos de EF de continuo se analiza la respuesta estructural de elementos constructivos (columnas, losas y pórticos) frente a acciones explosivas. Asimismo se han comparado con resultados experimentales, de ensayos sobre vigas y losas con distintas cargas de explosivo, verificándose una coincidencia aceptable y permitiendo una calibración de los parámetros de cálculo. Sin embargo estos modelos tan detallados no son recomendables para analizar edificios completos, ya que el elevado número de elementos finitos que serían necesarios eleva su coste computacional hasta hacerlos inviables para los recursos de cálculo actuales. Adicionalmente, se desarrollan modelos de elementos finitos estructurales (vigas y láminas) que, con un coste computacional reducido, son capaces de reproducir el comportamiento global de la estructura con una precisión similar. Se modelizan igualmente el hormigón en masa y el acero de armado de forma segregada. El hormigón se representa con el modelo constitutivo del hormigón EC2 (Hallquist et al., 2013), que también presenta un comportamiento inelástico, con diferente respuesta a tracción y compresión, endurecimiento, daño por fisuración y compresión, y rotura, y se usa en elementos finitos tipo lámina. El acero se representa de nuevo con un modelo constitutivo elastoplástico bilineal con rotura, usando elementos finitos tipo viga. Se modeliza una geometría equivalente del hormigón y del armado, y se tiene en cuenta la posición relativa del acero dentro de la masa de hormigón. Las mallas de ambos se unen mediante nodos comunes, produciendo una respuesta conjunta. En este trabajo se denominarán a estos modelos simplificadamente como modelos de EF estructurales. Con estos modelos de EF estructurales se simulan los mismos elementos constructivos que con los modelos de EF de continuo, y comparando sus respuestas estructurales frente a explosión se realiza la calibración de los primeros, de forma que se obtiene un comportamiento estructural similar con un coste computacional reducido. Se comprueba que estos mismos modelos, tanto los modelos de EF de continuo como los modelos de EF estructurales, son precisos también para el análisis del fenómeno de colapso progresivo en una estructura, y que se pueden utilizar para el estudio simultáneo de los daños de una explosión y el posterior colapso. Para ello se incluyen formulaciones que permiten considerar las fuerzas debidas al peso propio, sobrecargas y los contactos de unas partes de la estructura sobre otras. Se validan ambos modelos con un ensayo a escala real en el que un módulo con seis columnas y dos plantas colapsa al eliminar una de sus columnas. El coste computacional del modelo de EF de continuo para la simulación de este ensayo es mucho mayor que el del modelo de EF estructurales, lo cual hace inviable su aplicación en edificios completos, mientras que el modelo de EF estructurales presenta una respuesta global suficientemente precisa con un coste asumible. Por último se utilizan los modelos de EF estructurales para analizar explosiones sobre edificios de varias plantas, y se simulan dos escenarios con cargas explosivas para un edificio completo, con un coste computacional moderado. The frequency of explosions on buildings whether they are intended or accidental is small, but they can have catastrophic effects. Being able to predict in a accurate enough manner the consequences of these dynamic actions on civil buildings, among which frame-type reinforced concrete buildings are a frequent typology is desirable. In this doctoral thesis different practical options for the modeling and computer assisted numerical calculation of reinforced concrete structures submitted to explosions are explored. Numerical finite elements models with explicit time-based integration are employed, demonstrating their effective capacity in the simulation of the occurring fast dynamic and highly nonlinear physical and structural phenomena, allowing to predict the damage caused by the explosion itself as well as by the possible progressive collapse of the structure. The work has been carried out with the commercial finite elements code LS-DYNA (Hallquist, 2006), developing several types of calculation model classified in two main types: 1) Models based in continuum finite elements in which the continuous medium is discretized directly by means of nodal displacement degrees of freedom; 2) Models based on structural finite elements, with beams and shells, including kinematic hypothesis for linear and superficial elements. These models are developed and discussed at different levels: 1) material behaviour, 2) response of structural elements such as columns, beams and slabs, and 3) response of complete buildings or significative parts of them. Very detailed 3D continuum finite element models are developed, modeling mass concrete and reinforcement steel in a segregated manner. Concrete is represented with a constitutive concrete model CSCM (Murray et al., 2007), that has an inelastic behaviour, with different tension and compression response, hardening, cracking and compression damage and failure. The steel is represented with an elastic-plastic bilinear model with failure. The actual geometry of the concrete is modeled with 3D continuum finite elements and every and each of the reinforcing bars with beam-type finite elements, with their exact position in the concrete mass. The mesh of the model is generated by the superposition of the concrete continuum elements and the beam-type elements of the segregated reinforcement, which are made to follow the deformation of the solid in each point by means of a penalty algorithm, reproducing the behaviour of reinforced concrete. In this work these models will be called continuum FE models as a simplification. With these continuum FE models the response of construction elements (columns, slabs and frames) under explosive actions are analysed. They have also been compared with experimental results of tests on beams and slabs with various explosive charges, verifying an acceptable coincidence and allowing a calibration of the calculation parameters. These detailed models are however not advised for the analysis of complete buildings, as the high number of finite elements necessary raises its computational cost, making them unreliable for the current calculation resources. In addition to that, structural finite elements (beams and shells) models are developed, which, while having a reduced computational cost, are able to reproduce the global behaviour of the structure with a similar accuracy. Mass concrete and reinforcing steel are also modeled segregated. Concrete is represented with the concrete constitutive model EC2 (Hallquist et al., 2013), which also presents an inelastic behaviour, with a different tension and compression response, hardening, compression and cracking damage and failure, and is used in shell-type finite elements. Steel is represented once again with an elastic-plastic bilineal with failure constitutive model, using beam-type finite elements. An equivalent geometry of the concrete and the steel is modeled, considering the relative position of the steel inside the concrete mass. The meshes of both sets of elements are bound with common nodes, therefore producing a joint response. These models will be called structural FE models as a simplification. With these structural FE models the same construction elements as with the continuum FE models are simulated, and by comparing their response under explosive actions a calibration of the former is carried out, resulting in a similar response with a reduced computational cost. It is verified that both the continuum FE models and the structural FE models are also accurate for the analysis of the phenomenon of progressive collapse of a structure, and that they can be employed for the simultaneous study of an explosion damage and the resulting collapse. Both models are validated with an experimental full-scale test in which a six column, two floors module collapses after the removal of one of its columns. The computational cost of the continuum FE model for the simulation of this test is a lot higher than that of the structural FE model, making it non-viable for its application to full buildings, while the structural FE model presents a global response accurate enough with an admissible cost. Finally, structural FE models are used to analyze explosions on several story buildings, and two scenarios are simulated with explosive charges for a full building, with a moderate computational cost.

Computer simulation of form-roll design

Changes in modern structural design have created a demand for products which are light but possess high strength. The objective is a reduction in fuel consumption and weight of materials to satisfy both economic and environmental criteria. Cold roll forming has the potential to fulfil this requirement. The bending process is controlled by the shape of the profile machined on the periphery of the rolls. A CNC lathe can machine complicated profiles to a high standard of precision, but the expertise of a numerical control programmer is required. A computer program was developed during this project, using the expert system concept, to calculate tool paths and consequently to expedite the procurement of the machine control tapes whilst removing the need for a skilled programmer. Codifying the expertise of a human and the encapsulation of knowledge within a computer memory, destroys the dependency on highly trained people whose services can be costly, inconsistent and unreliable. A successful cold roll forming operation, where the product is geometrically correct and free from visual defects, is not easy to attain. The geometry of the sheet after travelling through the rolling mill depends on the residual strains generated by the elastic-plastic deformation. Accurate evaluation of the residual strains can provide the basis for predicting the geometry of the section. A study of geometric and material non-linearity, yield criteria, material hardening and stress-strain relationships was undertaken in this research project. The finite element method was chosen to provide a mathematical model of the bending process and, to ensure an efficient manipulation of the large stiffness matrices, the frontal solution was applied. A series of experimental investigations provided data to compare with corresponding values obtained from the theoretical modelling. A computer simulation, capable of predicting that a design will be satisfactory prior to the manufacture of the rolls, would allow effort to be concentrated into devising an optimum design where costs are minimised.

Finite element analysis of particle impact problems

Particle impacts are of fundamental importance in many areas and there has been a renewed interest in research on particle impact problems. A comprehensive investigation of the particle impact problems, using finite element (FE) methods, is presented in this thesis. The capability of FE procedures for modelling particle impacts is demonstrated by excellent agreements between FE analysis results and previous theoretical, experimental and numerical results. For normal impacts of elastic particles, it is found that the energy loss due to stress wave propagation is negligible if it can reflect more than three times during the impact, for which Hertz theory provides a good prediction of impact behaviour provided that the contact deformation is sufficiently small. For normal impact of plastic particles, the energy loss due to stress wave propagation is also generally negligible so that the energy loss is mainly due to plastic deformation. Finite-deformation plastic impact is addressed in this thesis so that plastic impacts can be categorised into elastic-plastic impact and finite-deformation plastic impact. Criteria for the onset of finite-deformation plastic impacts are proposed in terms of impact velocity and material properties. It is found that the coefficient of restitution depends mainly upon the ratio of impact velocity to yield Vni/Vy0 for elastic-plastic impacts, but it is proportional to [(Vni/Vy0)*(Y/E*)]-1/2, where Y /E* is the representative yield strain for finite-deformation plastic impacts. A theoretical model for elastic-plastic impacts is also developed and compares favourably with FEA and previous experimental results. The effect of work hardening is also investigated.

The fracture toughness of a range of cast steels

A range of plain carbon, carbon-manganese and low alloy cast steels were tested in order to determine their various fracture toughness values under elastic and elastic-plastic conditions. The main fracture toughness parameters which are considered are (1) Linear Elastic Fracture Mechanics (LEFM), (2) the J-Contour Integral, and (3) Crack Opening Displacement (COD). Results are obtained from fracture toughness specimens of various dimensions and the relevance of the validity criteria to cast steels is considered in some detail. In addition, the effect of casting position on specimen toughness values was noted. Valid KIC results according to LEFM, were obtained for three of the eight cast steels tested. Although KIC values from LEFM were not obtained from the remaining five steels, critical COD and J-integral values were determined. It is postulated that these values and particularly the critical J values can be used, with confidence for material selection or in defect tolerance calculations using these steels. Toughness values were found to vary with casting position in several of the steels tested and the possible reasons for such variations are discussed in the Thesis.

Computer simulation for tube-making by the cold roll-forming process

The conventional design of forming rolls depends heavily on the individual skill of roll designers which is based on intuition and knowledge gained from previous work. Roll design is normally a trial an error procedure, however with the progress of computer technology, CAD/CAM systems for the cold roll-forming industry have been developed. Generally, however, these CAD systems can only provide a flower pattern based on the knowledge obtained from previously successful flower patterns. In the production of ERW (Electric Resistance Welded) tube and pipe, the need for a theoretical simulation of the roll-forming process, which can not only predict the occurrence of the edge buckling but also obtain the optimum forming condition, has been recognised. A new simulation system named "CADFORM" has been devised that can carry out the consistent forming simulation for this tube-making process. The CADFORM system applied an elastic-plastic stress-strain analysis and evaluate edge buckling by using a simplified model of the forming process. The results can also be visualised graphically. The calculated longitudinal strain is obtained by considering the deformation of lateral elements and takes into account the reduction in strains due to the fin-pass roll. These calculated strains correspond quite well with the experimental results. Using the calculated strains, the stresses in the strip can be estimated. The addition of the fin-pass roll reduction significantly reduces the longitudinal compressive stress and therefore effectively suppresses edge buckling. If the calculated longitudinal stress is controlled, by altering the forming flower pattern so it does not exceed the buckling stress within the material, then the occurrence of edge buckling can be avoided. CADFORM predicts the occurrence of edge buckling of the strip in tube-making and uses this information to suggest an appropriate flower pattern and forming conditions which will suppress the occurrence of the edge buckling.

Three-dimensional finite element analysis of the effects of anisotropy on bone mechanical properties measured by nanoindentation

A three-dimensional finite element analysis (FEA) model with elastic-plastic anisotropy was built to investigate the effects of anisotropy on nanoindentation measurements for cortical bone. The FEA model has demonstrated a capability to capture the cortical bone material response under the indentation process. By comparison with the contact area obtained from monitoring the contact profile in FEA simulations, the Oliver-Pharr method was found to underpredict or overpredict the contact area due to the effects of anisotropy. The amount of error (less than 10% for cortical bone) depended on the indentation orientation. The indentation modulus results obtained from FEA simulations at different surface orientations showed a trend similar to experimental results and were also similar to moduli calculated from a mathematical model. The Oliver-Pharr method has been shown to be useful for providing first-order approximations in the analysis of anisotropic mechanical properties of cortical bone, although the indentation modulus is influenced by anisotropy.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.