Four different thermoviscoplastic relations, namely, the Litonski-Batra, the Johnson-Cook, the Bodner-Partom and the power law are used to model the thermoviscoplastic response of a material. Each one of these relations accounts for strain hardening, strain-rate hardening and thermal softening of the material. The material parameters in these relations are found by solving an initial-boundary-value problem corresponding to simple shearing deformations so that the computed effective stress vs. the effective plastic strain curves match closely with the experimental data of Marchand and Duffy who tested thin-walled HY-100 steel tubes in torsion. These four viscoplastic relations are used to analyze dynamic thermomechanical deformations of a prenotched plate impacted on the notched side by a cylindrical projectile made of the same material as the plate. The impact loading on the contact surface is simulated by prescribing the time history of the normal component of velocity and null tangential tractions. A plane strain state of deformation is assumed to prevail in the plate and its deformations are studied for different values of the impact speed. The in-house developed finite element code employs constant strain triangular elements, one point integration rule, and a lumped mass matrix. The Lagrangian description of motion is used to describe deformations of the plate. The coupled nonlinear partial differential equations are first reduced to coupled nonlinear ordinary differential equations (ODEs) by using the Galerkin approximation. The ODEs are integrated by using the stiff solver, LSODE, which adaptively adjusts the time step size and computes the solution within the prescribed accuracy. Results computed with the four constitutive relations are found to be qualitatively similar to each other and the general trends agree with the experimental observations in the sense that at low speed of impact, a brittle failure ensues at a point on the upper surface of the notch tip. However, at high impact speeds, a ductile failure in the form of a shear band initiates first from a point on the lower surface of the notch tip. The predicted speed at which the failure mode transitions from brittle to ductile is different for the four viscoplastic relations. Results have been computed using the Bodner-Partom law to study the effects of the notch tip radius and the presence of a circular hole ahead of the notch-tip. For sharp elliptic notch tips, it is found that there is no failure transition speed and the ductile failure always preceeded the brittle failure for the range of the impact speeds studied. For the hole located on the axis of the circular notch tip, the brittle failure always preceeded the ductile failure and it initiated at a point on the lower surface of the circular hole. / Ph. D.

Plane strain deformations of an elastoplastic material are studied using numerical methods. In the first chapter, a meshless formulation of the static small strain elastic-plastic problem is formulated using the meshless local Petrov-Galerkin method. The code is validated against the small strain plasticity routines in the commercial finite element code ABAQUS for two basic configurations with loading, unloading, and reloading. The results are found to agree within 5%. The validated code is then used to analyze the stress intensity factor (SIF) in a double edge-cracked plate. Deformations of the plate are studied both with and without exploiting the symmetry conditions. The penalty method is used to enforce the essential boundary condition in the former case. When analyzing the deformations of the entire plate, the diffraction method is employed in order to introduce the discontinuity in the displacement field across the crack faces. The log-log and a higher order extrapolation technique due to Dally and Berger (1996) are used to calculate the SIF. It is found that the penalty method was inadequate to enforce the essential boundary conditions in the vicinity of the crack tip and that in this region the deformations were oscillatory. Consequently, the SIF calculation using the higher order technique was not accurate. It is also found that for a small plastic zone (3% of the cracked length) the SIFs do not differ significantly from their values for the corresponding linear elastic problem. In the second chapter, a finite element formulation of the dynamic deformations of a micro-porous thermoviscoplastic solid is formulated. The heat conduction in a material is assumed to be governed by a hyperbolic heat equation; thus thermal and mechanical waves propagate with finite speeds. The formation and propagation of an adiabatic shear band (ASB) inplane strain tensile deformations is studied for eleven materials. The ASB is assumed to form when the maximum shear stress has been reduced to 80% of its peak value at a point and it is deforming plastically. The materials are ranked according their susceptibility to the formation of an ASB. A parametric study of the effect of the initial defect strength where the defect is assumed through an initially inhomogeneous distribution of porosity, the thermal conductivity, the thermal wave speed, and the applied strain-rate upon the ASB initiation and propagation is conducted. It is found that the susceptibility ranking for this configuration differs somewhat from that previously found for simple shear and torsion of thin-walled tubes. It is also found that thermal conductivity influences ASB initiation and propagation only for materials with large values of · and that for such materials an adiabatic model may not be adequate. The effects of initial defect strength and the nominal strain-rates are both found to be consistent with simple shearing studies except that the ASB propagation speed was found to decrease with increasing nominal strain-rate. It is found that the criterion employed for ASB initiation accurately predicts the onset of the collapse of the total axial load applied to the body. In the final chapter, the formulation from the previous chapter is modified to permit the formation and propagation of brittle and ductile fracture. Deformations of the impact loaded double edge-crack specimen of Kalthoff and Winkler (1987) are studied. The brittle to ductile failure mode transition with increasing impact speed was found. Previous studies have focused on identifying the transition speed and did not allow for crack propagation. In this study, crack propagation is achieved through a nodal release algorithm and interpenetration of the crack surfaces is prevented using stiff-spring contact elements. Brittle fracture is assumed to occur when the maximum tensile principal stress achieves a critical value and the ductile fracture is assumed to occur when the effective plastic strain reaches a critical value. It is found that the transition speed for 4340 steel is approximately 54 m/s. For the brittle failure, the stress field is found to be significantly modified by the propagating crack and in the vicinity of the propagating crack the field is mode-I dominant. The crack formed through brittle fracture is found to completely propagate through the plate. For the ductile failure, the distribution of effective plastic strain about the crack tip is not significantly altered by the formation of the crack. The temperature rise in the vicinity of the ductile crack is found to be approximately 45% of the melting temperature of the material. / Ph. D.

The investigation of failure in brittle materials, subjected to dynamic transient loading conditions, represents one of the ongoing challenges in the mechanics community. Progresses on this front are required to support the design of engineering components which are employed in applications involving extreme operational regimes. To this purpose, this thesis is devoted to the development of a framework which provides the capabilities to model how crack patterns form and evolve in brittle materials and how they affect the quantitative description of failure. The proposed model is developed within the context of diffusive interfaces which are at the basis of a new class of theories named phase field models. In this work, a set of additional features is proposed to expand their domain of applicability to the modelling of (i) rate and (ii) pressure dependent effects. The path towards the achievement of the first goal has been traced on the desire to account for micro-inertia effects associated with high rates of loading. Pressure dependency has been addressed by postulating a mode-of-failure transition law whose scaling depends upon the local material triaxiality. The governing equations have been derived within a thermodynamically-consistent framework supplemented by the employment of a micro-forces balance approach. The numerical implementation has been carried out within an updated lagrangian finite element scheme with explicit time integration. A series of benchmarks will be provided to appraise the model capabilities in predicting rate-pressure-dependent crack initiation and propagation. Results will be compared against experimental evidences which closely resemble the boundary value problems examined in this work. Concurrently, the design and optimization of a complimentary, improved, experimental characterization platform, based on the split Hopkinson pressure bar, will be presented as a mean for further validation and calibration.