A novel non-bonded interface technique (NBIT) is used to analyze internal residual strain by combining a pre-split sample of AISI 4340 steel with the circular grid residual strain analysis technique. NBIT is compared with an implicit non-linear finite element (FE) model using LS-DYNA. A split FE model was compared with a quarter FE model to determine the split interface that causes an average difference of 9.0% on the residual von Mises strain field from a 588.6 N indentation. The homogeneous FE quarter model was then compared with the experimental split model using 588.6, 981.0, and 1471.5 N indentation forces. An average displacement difference of 3.92 µm was found when comparing the experimental split and FE homogeneous samples from a 588.6 N indentation. The internal residual major and minor principal strains from the split experimental sample and homogeneous FE model were compared for each indentation force. The minor principal strain results show the 588.6, 981.0, and 1471.5 N indentation forces resulted in a difference between the experimental split and homogeneous FE model of 28.5%, 34.8%, and 26.0%, respectively. The difference between the comparisons was explained by the inability of the FE model to simulate local non-homogeneous material properties such as grain composition and orientation whereas NBIT does. NBIT can be used for micro- or macro-scale residual strain analysis as the spatial resolution is highly adjustable.

The study presents the anisotropy in the indentation creep response of zinc. Indentation creep tests were conducted on grains with three different orientations, namely, 〈 0001 〉, ⟨ 11 2 ¯ 0 ⟩ , and ⟨ 10 1 ¯ 0 ⟩ . Indentation creep along the 〈c〉 axis showed a pronounced creep exhibiting a 25% higher creep displacement as compared with indentation perpendicular to the 〈c〉 axis. However, the recovery observed was the highest for basal-oriented grains wherein a recovery of 50% to its indented depth was observed, whereas for ⟨ 11 2 ¯ 0 ⟩ the recovery observed was 30%. The stress exponent, n, was obtained by employing two different approaches for each of the orientations and a difference in the stress exponent value was also observed, highlighting the anisotropic creep response. Basal and non-basal oriented grains showed a lower (n ≅ 1.6 ) and higher stress exponent (n ≅ 5 ), respectively, highlighting the different operable creep mechanism. Surface topography using atomic force microscope (AFM) revealed twinning and sink-in for ⟨ 11 2 ¯ 0 ⟩ and ⟨ 10 1 ¯ 0 ⟩ , whereas uniform pile up was observed for basal grain.

Johnson–Cook (JC) strength and failure models have been widely used in finite element analysis (FEA) to solve a variety of thermo-mechanical problems. There are many techniques to determine the required JC parameters; however, a best practice to obtain the most reliable JC parameters has not yet been proposed. In this paper, a genetic-algorithm-based optimization strategy is proposed to calibrate the JC strength and failure model parameters of AISI/SAE 1018 steel. Experimental data were obtained from tensile tests performed for different specimen geometries at varying strain rates and temperatures. FEA was performed for each tensile test. A genetic algorithm was used to determine the optimum JC parameters that best fit the experimental force-displacement data. Calibrated JC parameters were implemented in FEA to simulate the impact tests of standard V-notch Charpy bars to verify the damage mechanism in the material. Considering good agreement of the experimental and FEA results, the current strategy is suggested for calibration proposes in other kind of materials in which plastic behavior could be represented by the JC strength and failure models.

In addition to processing a troubling agricultural by-product and reducing demands on our landfills, prepared agro-waste composites are suitable for a variety of practical applications. However, enhancing value-added options for these agricultural by-products can necessitate ability to assess their mechanical integrity. This paper accordingly describes the preparation of a cellulosic-manure composite and demonstrates ability to determine stresses in a perforated structure of the material from measured displacement data. Processing digital image correlation (DIC) recorded displacement information with an Airy stress function gives reliable results full-field as well as at the edge of geometric discontinuities without having to differentiate the recorded data. Required constitutive properties are evaluated in situ and results are substantiated independently.

This work documents the development of a tool to perform automated parameter fitting of constitutive material models. Specific to this work is the fitting of a Swift hardening rule and isotropic linear plasticity model to aluminum 2024-T351, C36000 brass, and C10100 copper. Material characterization was conducted through the use of compressive, cold upsetting tests. A noncontact, optical displacement measurement system was applied to measure the axial and radial deformation of the test specimens. Nonlinear optimization techniques were then applied to tune a finite element model to match experimental results through the optimization of material model parameters as well as frictional coefficient. The result is a system, which can determine constitutive model parameters rapidly and without user interaction. While this tool provided material parameters for each material and model tested, the quality of the fit varied depending on how appropriate the constitutive model was to the material's actual plastic behavior. Aluminum's behavior proved to be an excellent match to the Swift hardening rule while the behavior of brass and copper was described better by the linear plasticity model.

Important from exploitation point of view mechanical properties of single-layer, double-layer, and mixed alumina and ceria films and their stainless steel (SS) substrate were investigated by means of nanoindentation experiments. As a result, we obtained the experimental load–displacement curves and calculated the indentation hardness (H IT ) and indentation modulus (E IT ), by means of Oliver and Pharr approximation method. Numerical simulations of the process of nanoindentation by means of finite element method were performed as well, in order to obtain more information about the plastic properties of the investigated films. The obtained results show that the mixed Al 2 O 3 +Ce 2 O 3 film, obtained at dominant concentration of cerium ions in the working electrolyte, has the highest indentation hardness and modulus, followed by the single Ce 2 O 3 -CeO 2 film, the mixed Al 2 O 3 +Ce 2 O 3 film, obtained at dominant concentration of aluminum ions in the working electrolyte, the double Ce 2 O 3 -CeO 2 /Al 2 O 3 layer, and single Al 2 O 3 layer.

Micro-computed tomography (CT) was used as a tool to investigate the deformation behavior of particulate-filled composite materials. Three different shapes of glass fillers (spherical, flake, and fiber) and filler mass fractions (5%, 10%, and 15%) were introduced to the epoxy resin. Rockwell hardness H scale indentation test was used to deform the composite material. The composite materials were scanned before and after the indentation test by using micro-CT. Displacement field for each filler type and mass fraction were measured through correlation of before and after scan data. The effects of filler type and mass fraction on the internal displacement field were investigated. It was also demonstrated that micro-CT can be used as a tool to create realistic representative volume elements (RVEs) for particulate-filled composite materials instead of randomly distributed particles within the matrix material.

This study explores the application of a proper orthogonal decomposition (POD) and radial basis function (RBF)-based surrogate model to identify the parameters of a nonlinear viscoelastic material model using nanoindentation data. The inverse problem is solved by reducing the difference between finite element simulation-trained surrogate model approximation and experimental data through genetic algorithm (GA)-based optimization. The surrogate model, created using POD–RBF, is trained using finite element (FE) data obtained by varying model parameters within a parametric space. Sensitivity of the model parameters toward the load–displacement output is utilized to reduce the number of training points required for surrogate model training. The effect of friction on simulated load–displacement data is also analyzed. For the obtained model parameter set, the simulated output matches well with experimental data for various experimental conditions.

Displacement and stress fields in a functionally graded (FG) fiber-reinforced rotating disk of nonuniform thickness subjected to angular deceleration are obtained. The disk has a central hole, which is assumed to be mounted on a rotating shaft. Unidirectional fibers are considered to be circumferentially distributed within the disk with a variable volume fraction along the radius. The governing equations for displacement and stress fields are derived and solved using finite difference method. The results show that for disks with fiber rich at the outer radius, the displacement field is lower in radial direction but higher in circumferential direction compared to the disk with the fiber rich at the inner radius. The circumferential stress value at the outer radius is substantially higher for disk with fiber rich at the outer radius compared to the disk with the fiber rich at the inner radius. It is also observed a considerable amount of compressive stress developed in the radial direction in a region close to the outer radius. These compressive stresses may prevent any crack growth in the circumferential direction of such disks. For disks with fiber rich at the inner radius, the presence of fibers results in minimal changes in the displacement and stress fields when compared to a homogenous disk made from the matrix material. In addition, we concluded that disk deceleration has no effect on the radial and hoop stresses. However, deceleration will affect the shear stress. Tsai–Wu failure criterion is evaluated for decelerating disks. For disks with fiber rich at the inner radius, the failure is initiated between inner and outer radii. However, for disks with fiber rich at the outer radius, the failure location depends on the fiber distribution.

This study concerns the development of peridynamic (PD) strain energy density functions for a Neo-Hookean type membrane under equibiaxial, planar, and uniaxial loading conditions. The material parameters for each loading case are determined by equating the PD strain energy density to that of the classical continuum mechanics. The PD equations of motion are derived based on the Neo-Hookean model under the assumption of incompressibility. Numerical results concern the deformation of a membrane with a defect in the form of a hole, a crack, and a rigid inclusion under equibiaxial, planar, and uniaxial loading conditions. The PD predictions are verified by comparison with those of finite element analysis.

We investigate a domain decomposition method (DDM) of finite element method (FEM) using Intel's many integrated core (MIC) architecture in order to determine the most effective MIC usage. For this, recently introduced high-scalable parallel method of DDM is first introduced with a detailed procedure. Then, the Intel's Xeon Phi MIC architecture is presented to understand how to apply the parallel algorithm into a multicore architecture. The parallel simulation using the Xeon Phi MIC has an advantage that traditional parallel libraries such as the message passing interface (MPI) and the open multiprocessing (OpenMP) can be used without any additional libraries. We demonstrate the DDM using popular libraries for solving linear algebra such as the linear algebra package (LAPACK) or the basic linear algebra subprograms (BLAS). Moreover, both MPI and OpenMP are used for parallel resolutions of the DDM. Finally, numerical parallel efficiencies are validated by a two-dimensional numerical example.

The paper presents the derived equations for calculations of the initial wall thickness g 0 of a tube bent to elbow. The expressions for calculating g 0 are presented in a suitable measure of the “great active actual radius R j ” in the bending zone for an exact-generalized solution (continuous fields) and for three formal simplifications (discontinuous fields) of the first-, second-, and third-orders. The expressions to calculate the components of deformation for a generalized solution (continuous fields) are obtained on the basis of kinematically admissible fields of plastic deformations. In any case, a value of initial tube thickness depends on the radius and on the angle of bending α b on the external diameter of the tube, on the displacement of the neutral axis, and on the allowable (required) elbow thickness according to European, American, or other national technical standard or regulations. The initial thickness also depends on the coordinates of the point where the allowable thickness was determined and on the technological–material coefficient of the bending zone range k (defined during the tests). The obtained calculation results are presented in the form of graphs and in table.

In real physical experiments, three typical deformation stages including elastic deformation stage, symmetric deformation stage, and asymmetric deformation stage appear step by step when the stainless steel hemispherical shell structure is under axial compression loading. During the asymmetric deformation stage, the rolling-plastic-hinge-radius which characterizes the size of the deformation area evolves along the circumferential direction with the compressive displacement. For the hemispherical shell structures with apparent asymmetric deformation stage, the double-buckling phenomenon of the structures in experiments can be clearly detected. The traditional theoretical analysis based on the assumption with circumferentially constant rolling-plastic-hinge-radius is not suitable to predict this phenomenon. For these hemispherical shell structures, load capacity and absorbed energy predicted by the traditional analysis are usually higher than experimental results in the asymmetric deformation stage. In this paper, a new description based on experimental observation for the evolution of rolling-plastic-hinge-radius has been proposed. Minimum energy principle was employed to obtain the postbuckling behavior. The energy evolution of different buckling stages during compression loading is investigated to evaluate the structure load capacity. Stainless steel hemispherical specimens with different sizes are tested under axial compression between two rigid plates to verify the theoretical modification. Good agreement is achieved between proposed model and experimental results. The theoretical model proposed in this paper can be used in prediction of postbuckling behavior for different deformation patterns in the asymmetric deformation stage. It also provides higher flexibility and efficiency for the postbuckling behavior prediction of hemispherical shell structures.

Cellular materials' two important properties—structure and mechanism—can be selectively used for materials design; in particular, they are used to determine the modulus and yield strain. The objective of this study is to gain a better understanding of these two properties and to explore the synthesis of compliant cellular materials (CCMs) with compliant porous structures (CPSs) generated from modified hexagonal honeycombs. An in-plane constitutive CCM model with CPSs of elliptical holes is constructed using the strain energy method, which uses the deformation of hinges around holes and the rotation of links. A finite element (FE) based simulation is conducted to validate the analytical model. The moduli and yield strains of the CCMs with an aluminum alloy are about 4.42 GPa and 0.57% in one direction and about 2.14 MPa and 20.9% in the other direction. CCMs have extremely high positive and negative Poisson's ratios (NPRs) (νxy* ∼ ±40) due to the large rotation of the link member in the transverse direction caused by an input displacement in the longitudinal direction. A parametric study of CCMs with varying flexure hinge geometries using different porous shapes shows that the hinge shape can control the yield strength and strain but does not affect Poisson's ratio which is mainly influenced by rotation of the link members. The synthesized CPSs can also be used to design a new CCM with a Poisson's ratio of zero using a puzzle-piece CPS assembly. This paper demonstrates that compliant mesostructures can be used for next generation materials design in tailoring mechanical properties such as moduli, strength, strain, and Poisson's ratios.

Instrumented indentation is commonly used for determining mechanical properties of a range of materials, including viscoelastic materials. However, most—if not all—studies are limited to a flat substrate being indented by various shaped indenters (e.g., conical or spherical). This work investigates the possibility of extending instrumented indentation to nonflat viscoelastic substrates. In particular, conical indentation of a sphere is investigated where a semi-analytical approach based on “the method of functional equations” has been developed to obtain the force–displacement relationship. To verify the accuracy of the proposed methodology selected numerical experiments have been performed and good agreement was obtained. Since it takes significantly less time to obtain force–displacement relationships using the proposed method compared to conducting full finite element simulations, the proposed method is an efficient substitute of the finite element method in determining material properties of viscoelatic spherical particles using indentation testing.

In this paper, the effect of surface damage induced by focused ion beam (FIB) fabrication on the mechanical properties of silicon (Si) nanowires (NWs) was investigated. Uniaxial tensile testing of the NWs was performed using a reusable on-chip tensile test device with 1000 pairs of comb structures working as an electrostatic force actuator, a capacitive displacement sensor, and a force sensor. Si NWs were made from silicon-on-nothing (SON) membranes that were produced by deep reactive ion etching hole fabrication and ultrahigh vacuum annealing. Micro probe manipulation and film deposition functions in a FIB system were used to bond SON membranes to the device's sample stage and then to directly fabricate Si NWs on the device. All the NWs showed brittle fracture in ambient air. The Young's modulus of 57 nm-wide NW was 107.4 GPa, which was increased to 144.2 GPa with increasing the width to 221 nm. The fracture strength ranged from 3.9 GPa to 7.3 GPa. By assuming the thickness of FIB-induced damage layer, the Young's modulus of the layer was estimated to be 96.2 GPa, which was in good agreement with the literature value for amorphous Si.

Determining a material constitutive law that is representative of the extreme conditions found in the cutting zone during machining operations is a very challenging problem. In this study, dynamic shear tests, which reproduce, as faithfully as possible, these conditions in terms of strain, strain rate, and temperature, have been developed using hat-shaped specimens. The objective was to identify the parameters of a Johnson–Cook material behavior model by an inverse method for two titanium alloys: Ti6Al4V and Ti555-3. In order to be as representative as possible of the experimental results, the parameters of the Johnson–Cook model were not considered to be constant over the total range of the strain rate and temperature investigated. This reflects a change in the mechanisms governing the deformation. The shear zones observed in hat-shaped specimens were analyzed and compared to those produced in chips during conventional machining for both materials. It is concluded that the observed shear bands can be classified as white-etching bands only for the Ti555-3 alloy. These white bands are assumed to form more easily in the Ti555-3 alloy due to its predominately β phase microstructure compared to the Ti6Al4V alloy with a α + β microstructure.

The indentation properties of pure aluminum (99.9%, 3N aluminum) and high purity aluminum (99.999%, 5N aluminum) with respect to the strain rate dependence of strength are experimentally investigated in order to clarify the effect of strain rate on the micro-indentation test. A micro-indentation test using a Berkovich indenter was performed at loading rates of 0.7, 7, and 70 mN/s. In all of the specimens, the indenter was loaded to a maximum value of 1200 mN, and then was maintained for 30 s. In the 3N specimen, the dependence of the loading rate on the load was slight at loading rates of 0.7 and 7 mN/s, whereas the load at the loading rate of 70 mN/s was higher than the loads at loading rates of 0.7 and 7 mN/s. On the other hand, the load for the 5N specimen increased with the increasing loading rate. Thus, the effect of the loading rate on the load-displacement curve for the 3N and 5N specimens was similar to the strain rate dependence of strength for theses metals. In addition, the micro-indentation test was demonstrated to be strongly affected by high strain rate at a loading rate of 70 mN/s.

Finite element analysis techniques were used to study the release step in a nanofabrication process. These calculations employed a novel adhesion/atomistic friction surface interaction model to define how the glassy polymer interacts with the hard mold. This model is applicable to solids that interact via relatively weak, van der Waals forces and is applicable to intentionally weakened interfaces (e.g., when a mold release is used). The goal of this effort is to simulate the entire separation process. The release step was studied by performing unit cell calculations for a pattern composed of identical, parallel channels. The interface between the mold and the glassy polymer did not unzip in a continuous, quasi-static manner in these simulations. Instead, there was a complex failure sequence that included multiple dynamic separations and arrest events as well as adhesive reattachment. The sensitivity of the release process to interface and bulk material properties, polymer shrinkage, and feature geometry was then quantified by examining variations from a baseline configuration. Finally, the feasibility of a hierarchical analysis that represents the nanometer-scale pattern by a pattern traction–separation (T–U) relationship, which is defined by a unit cell analysis, was demonstrated.

In finite element analysis of sheet metal forming the use of combined isotropic-kinematic hardening models is advisable to improve stamping simulation and springback prediction. This choice becomes compulsory to model recent materials such as high strength steels. Cyclic tests are strictly required to evaluate the parameters of these constitutive models. However, for sheet metal specimens, in case of simple axial tension-compression tests, buckling occurrence during compression represents a serious drawback. This is the reason why alternative set-ups have been devised. In this paper, two experimental arrangements (a cyclic laterally constrained tension-compression test and a three-point fully reversed bending test) are compared so as to point out the advantages and the disadvantages of their application in tuning the well-known Chaboche’s hardening model. In particular, for tension-compression tests, a new clamping device was specifically designed to inhibit compressive instability. Four high strength steel grades were tested: two dual phases (DP), one transformation induced plasticity (TRIP) and one high strength low alloy material (HSLA). Then, the Chaboche’s model was calibrated through inverse identification methods or by means of analytical expressions when possible. The proposed testing procedure proved to be successful in all investigated materials. The achieved constitutive parameters, obtained independently from the two experimental techniques, were found to be consistent. Their accuracy was also been assessed by applying the parameter set obtained from one test to simulate the other one, and vice versa. Clues on what method provides the better transferability are given.