Fracture toughness of open-cell Kelvin foam

Brittle fracture behavior of a perfect open-cell Kelvin foam is considered. The foam is modeled as a spatial lattice consisting of brittle elastic struts rigidly connected to each other at the nodal points. The fracture toughness is determined from the analysis of a quasi-plane problem for a slice of the foam with an embedded finite length crack generated by broken struts. The crack plane is chosen on the basis of a previous study of crack nucleation phenomenon, and the crack length, which assures the self-similar K-field in the tip vicinity, is established by numerical experiments. For the considered densities range the crack includes several hundreds of broken struts and, consequently, the portion of the foam to be considered in the analysis has a very large number of nodal degrees of freedom. The computational cost is reduced significantly by using for the analysis the representative cell method based on the discrete Fourier transform. As a result, the initial problem for the foam slice is reduced to the problem for the repetitive cell which includes 12 struts.     

通孔泡沫铝的动态压缩行为

在SHPB装置上对渗流法制备的通孔泡沫铝进行了动态压缩实验,研究了相对密度为0.341~0.419的通孔泡沫铝在10-3~2000 s-1应变率范围内的压缩响应特征和应变率相关性,并用扫描电镜（scanning electron microscope,SEM)分析了泡沫铝的压缩变形特征。实验结果表明,通孔泡沫铝有明显应变率效应,随应变率上升,泡沫铝流动应力提高。SEM观察结果揭示,在动态压缩下,通孔泡沫铝宏观上均匀变形,微观变形机制以泡孔横向伸展坍塌为主。     

Constitutive modeling of aluminum foam including fracture and statistical variation of density

An existing constitutive model applicable for aluminum foam was evaluated. The model was implemented in LS-DYNA, and several test cases were analyzed and compared to experimental data. The numerical analyses consisted of foam subjected to both simple and more complex loading conditions where fracture was of varying importance. Therefore, simple fracture criteria were added to the model. Additionally, the inhomogeneities in the foam were modeled by incorporating the possibility of statistical variation of the foam density. The implemented model is efficient and robust, and gives satisfactory results compared with the experimental results.     

Multi-scale modeling of tensile behavior of carbon nanotube-reinforced composites

A multi-scale representative volume element (RVE) for modeling the tensile behavior of carbon nanotube-reinforced composites is proposed. The RVE integrates nanomechanics and continuum mechanics, thus bridging the length scales from the nano- through the mesoscale. A progressive fracture model based on the modified Morse interatomic potential is used for simulating the behavior of the isolated carbon nanotubes and the FE method for modeling the matrix and building the RVE. Between the nanotube and the matrix a perfect bonding is assumed until the interfacial shear stress exceeds the corresponding strength. Then, nanotube/matrix debonding is simulated by prohibiting load transfer in the debonded region. Using the RVE, a unidirectional nanotube/polymer composite was modeled and the results were compared with corresponding rule-of-mixtures predictions. A significant enhancement in the stiffness of the polymer owing to the adding of the nanotubes is predicted. The effect of interfacial shear strength on the tensile behavior of the nanocomposite was also studied. Stiffness is found to be unaffected while tensile strength to significantly decrease with decreasing the interfacial shear strength.     

A constitutive model of aluminum foam for crash simulations

A new constitutive model for metallic foams is developed to overcome the deficiency of existing models in commercial finite element codes. The proposed constitutive model accounts for volume changes under hydrostatic compression and combines the hydrostatic pressure and von Mises stress into one yield function. The change of the compressibility of the metallic foam is handled in the constitutive model by allowing for shape changes of the yield surface in the hydrostatic pressure-von Mises stress space. The backward Euler method is adopted to integrate the constitutive equations to achieve numerical accuracy and stability. The model is implemented into LS-DYNA as a user-defined subroutine, verified with existing solutions, and validated with foam testing data. The verified and validated model is then utilized in the crushing simulations of foam-filled columns with square and hexagonal cross-sections. Two constitutive models are studied: the first using an exponential function to describe the relationship of plastic Poisson's ratio with respect to true strain and the second using linear interpolation function as an alternative approximation. The new foam model provides satisfactory prediction of crushing forces and deformed shapes compared to experimental results. Additionally, the new foam model was shown to have better numerical stability and accuracy than existing LS-DYNA built-in material models.     

Analysis and simulation for tensile behavior of anisotropic open-cell elastic foams

Based on the elongated Kelvin obtained to investigate the tensile behavior Kelvin model＇s periodicity and symmetry in model, a simplified periodic structural cell is of anisotropic open-cell elastic foams due to the whole space. The half-strut element and elastic deflection theory are used to analyze the tensile response as done in the previous studies. This study produces theoretical expressions for the tensile stress-strain curve in the rise and transverse directions. In addition, the theoretical results are examined with finite element simulation using an existing formula. The results indicate that the theoretical analysis agrees with the finite element simulation when the strain is not too high, and the present model is better. At the same time, the anisotropy ratio has a significant effect on the mechanical properties of foams. As the anisotropy ratio increases, the tensile stress is improved in the rising direction but drops in the transverse direction under the same strain.     

Modeling of the role of defects in sintered FeCrAlY foams

The metal sintering approach offers a cost- effective means for the mass-production of open-cell foams from a range of materials, including high-temperature steel alloys, which offer novel mechanical and acoustic properties. In a separate experimental study, the mechanical properties of open-celled steel alloy (FeCrAlY) foams have been characterized under uniaxial compression and shear loading. Compared to predictions from established models, a significant knockdown in material properties was observed. This knockdown was attributed to the presence of defects throughout the microstructure that result from the unique fabrication process. In the present paper, the microstructure of sintered FeCrAlY foams was modeled by using a finite element (FE) model. In particular, microstructural variations were introduced to a base lattice, and the effects on the strength and stiffness calculated. A range of defects identified under scanning electronic microscope (SEM) imaging were considered including broken ligaments, thickness variations, and pore blockages, which are the three primary imperfections observed in sintered foams. The corresponding levels of defect present in the material were subsequently input into the FE model, with the resulting predictions correlating well with experimental data. The project supported by the National Basic Research Program of China (2006CB601202), the National Natural Science Foundation of China (10328203, 10572111, 10632060), the National 111 Project of China (B06024), and the US Office of Naval Research (N000140210117).     

Influence of skin/core debonding on free vibration behavior of foam and honeycomb cored sandwich plates

The dynamic behavior of partially delaminated at the skin/core interface sandwich plates with flexible cores is studied. The commercial finite element code ABAQUS is used to calculate natural frequencies and mode shapes of the sandwich plates containing a debonding zone. The influence of the debonding size, debonding location and types of debonding on the modal parameters of damaged sandwich plates with various boundary conditions is investigated. The results of dynamic analysis illustrated that they can be useful for analyzing practical problems related to the non-destructive damage detection of partially debonded sandwich plates.     

Growth limit of carbon onions – A continuum mechanical study

The growth of carbon onions is simulated using continuum mechanical shell models. With this models it is shown that, if a carbon onion has grown to a critical size, the formation of an additional layer leads to the occurrence of a structural instability. This instability inhibits further growth of carbon onions and, thus, can be a reason for the limited size of such particles. The loss of stability is mainly evoked by van der Waals interactions between misfitting neighboring layers leading to self-equilibrating stress states in the layers due to mutual accommodation. The influence of the curvature induced surface energy and its consequential stress state is investigated and found to be rather negligible. Furthermore, it is shown that the nonlinear character of the van der Waals interactions has to be considered to obtain maximum layer numbers comparable to experimental observations. The proposed model gives insight into mechanisms which are assumed to limit the size of carbon onions and can serve as basis for further investigations, e.g., of the formation of nanodiamonds in the center of carbon onions.     

Modeling of anisotropic wound healing

Biological soft tissues exhibit non-linear complex properties, the quantification of which presents a challenge. Nevertheless, these properties, such as skin anisotropy, highly influence different processes that occur in soft tissues, for instance wound healing, and thus its correct identification and quantification is crucial to understand them. Experimental and computational works are required in order to find the most precise model to replicate the tissues' properties. In this work, we present a wound healing model focused on the proliferative stage that includes angiogenesis and wound contraction in three dimensions and which relies on the accurate representation of the mechanical behavior of the skin. Thus, an anisotropic hyperelastic model has been considered to analyze the effect of collagen fibers on the healing evolution of an ellipsoidal wound. The implemented model accounts for the contribution of the ground matrix and two mechanically equivalent families of fibers. Simulation results show the evolution of the cellular and chemical species in the wound and the wound volume evolution. Moreover, the local strain directions depend on the relative wound orientation with respect to the fibers.     

Three-dimensional modeling of the mechanical property of linearly elastic open cell foams

Three-dimensional Voronoi models are developed to investigate the mechanical behavior of linearly elastic open cell foams. Dependence of the Young’s modulus, Poisson’s ratio and bulk modulus of the foams on the relative density is evaluated through finite element analysis. Obtained results show that in the low density regime the Young’s modulus and bulk modulus of random Voronoi foams can be well represented by those of Kelvin foams, and are sensitive to the geometric imperfections inherent in the microstructure of foams. In contrast, the compressive plateau stress of the foams is less sensitive to the imperfections. Failure surface of the foams subject to multi-axial compression is determined and is found to comply with the maximum compressive principal stress criterion, consistent with available experimental observations on polymer foams. Numerical results also show that elastic buckling of cell edges at microscopic level is the dominant mechanism responsible for the compressive failure of elastic open cell foams.     

Modeling segmentation cracking of a brittle coating due to underneath periodic eigenstrains

The formation mechanism for a segmented Cr coating by the hybrid technique of laser pre-quenching steel substrate plus post-electroplating was investigated. The discrete laser quenched zones (LQZs) were modeled as multiple inclusions with prescribed eigenstrains. The stress field was determined to account for the onset of segmentation cracking. Also addressed was segmentation crack growth through the evaluation of the stress intensity factor (SIF). The computations were implemented by using the finite element (FE) method. The dependence of a wide range of dimensionless variables of interest on both the stress field and SIF was assessed through a detailed parametric study.     

Thermal vibration and apparent thermal contraction of single-walled carbon nanotubes

It is of fundamental value to understand the thermo-mechanical properties of carbon nanotubes. In this paper, by using molecular dynamics simulation, a systematic numerical investigation is carried out to explore the natural thermal vibration behaviors of single-walled carbon nanotubes and their quantitative contributions to the apparent thermal contraction behaviors. It is found that the thermo-mechanical behavior of single-walled carbon nanotubes is exhibited through the competition between quasi-static thermal expansion and dynamic thermal vibration, while the vibration effect is more prominent and induces apparent contraction in both radial and axial directions. With increasing temperature, the anharmonic interatomic potential helps to increase the bond length, which leads to thermally induced expansion. On the other hand, the higher structural entropy and vibrational entropy of the system cause the carbon nanotube to vibrate, and the apparent length of nanotube decreases due to various vibration modes. Parallel analytical and finite element analyses are used to validate the vibration frequencies and provide helpful insights. The unified multi-scale study has successfully decoupled and systematically analyzed both thermal expansion and contraction behaviors of single-walled carbon nanotube from 100 to 800 K, and obtained detailed information on various vibration modes as well as their quantitative contributions to the coefficient of thermal expansion in axial and radial directions. The results of this paper may provide useful information on the thermo-mechanical integrity of single-walled carbon nanotubes, and become important in practical applications involving finite temperature.     

Finite element computation of trapped and leaky elastic waves in open stratified waveguides

Elastic guided waves are of interest for inspecting structures due to their ability to propagate over long distances. In numerous applications, the guiding structure is surrounded by a solid matrix that can be considered as unbounded in the transverse directions. The physics of waves in such an open waveguide significantly differs from a closed waveguide, i.e. for a bounded cross-section. Except for trapped modes, part of the energy is radiated in the surrounding medium, yielding attenuated modes along the axis called leaky modes. These leaky modes have often been considered in non destructive testing applications, which require waves of low attenuation in order to maximize the inspection distance. The main difficulty with numerical modeling of open waveguides lies in the unbounded nature of the geometry in the transverse direction. This difficulty is particularly severe due to the unusual behavior of leaky modes: while attenuating along the axis, such modes exponentially grow along the transverse direction. A simple numerical procedure consists in using absorbing layers of artificially growing viscoelasticity, but large layers may be required. The goal of this paper is to explore another approach for the computation of trapped and leaky modes in open waveguides. The approach combines the so-called semi-analytical finite element method and a perfectly matched layer technique. Such an approach has already been successfully applied in scalar acoustics and electromagnetism. It is extended here to open elastic waveguides, which raises specific difficulties. In this paper, two-dimensional stratified waveguides are considered. As it reveals a rich structure, the numerical eigenvalue spectrum is analyzed in a first step. This allows to clarify the spectral objects calculated with the method, including radiation modes, and their dependency on the perfectly matched layer parameters. In a second step, numerical dispersion curves of trapped and leaky modes are compared to analytical results.     

Continuum Mechanics Modeling and Simulation of Carbon Nanotubes

The understanding of the mechanics of atomistic systems greatly benefits from continuum mechanics. One appealing approach aims at deductively constructing continuum theories starting from models of the interatomic interactions. This viewpoint has become extremely popular with the quasicontinuum method. The application of these ideas to carbon nanotubes presents a peculiarity with respect to usual crystalline materials: their structure relies on a two-dimensional curved lattice. This renders the cornerstone of crystal elasticity, the Cauchy–Born rule, insufficient to describe the effect of curvature. We discuss the application of a theory which corrects this deficiency to the mechanics of carbon nanotubes (CNTs). We review recent developments of this theory, which include the study of the convergence characteristics of the proposed continuum models to the parent atomistic models, as well as large scale simulations based on this theory. The latter have unveiled the complex nonlinear elastic response of thick multiwalled carbon nanotubes (MWCNTs), with an anomalous elastic regime following an almost absent harmonic range.     

Experiments and modeling of carbon nanotube-based NEMS devices

In this paper, carbon nanotube-based nanoelectromechanical systems (NEMS) are nanofabricated and tested. In-situ scanning electron microscopy measurements of the deflection of the cantilever under electrostatic actuation are reported. In particular, a cantilever nanotube suspended over an electrode (nanoswitch), or two symmetric cantilever nanotubes (nanotweezers), from which a differential in electrical potential is imposed, are studied. The finite deformation regime investigated here is the first of its kind. An analytical model based on the energy method in both small deformation and finite kinematics (large deformation) regimes is used to interpret the measurements. The theory overcomes limitations of prior analysis reported in the literature towards the prediction of the structural behavior of NEMS. Some of the simplifying hypotheses have been removed. Furthermore, the theory takes into account the cylindrical shape of the deflected nanotube in the evaluation of its electrical capacitance, the influence of the van der Waals forces as well as finite kinematics. In addition, tip charge concentration and a quantum correction of the electrical capacitance are also considered. The energy-based method is used to predict the structural behavior and instability of the nanotube, corresponding to the on/off states of the nanoswitch, or to the open/close states of the nanotweezers—at the so-called pull-in voltage. Accuracy of the derived formulas is assessed by comparison of the theoretical prediction and experimental data in both small deformation and finite kinematics regimes. The results reported in this work are particularly useful in the characterization of the electromechanical properties of nanotubes as well as in the optimal design of nanotube-based NEMS devices.     

泡沫铝率相关性能的有限元模拟

构建了三维随机分布球形泡孔模型,模拟开、闭孔混合结构泡沫铝材料的微细观结构,并通过有限元方法计算了10~104 s-1应变率范围内、孔隙率35%~65%泡沫铝材料的率相关性以及应变率和相对密度变化对泡沫铝动态压缩力学性能的影响。研究表明:中、低应变率下,泡沫铝材料率相关性能主要取决于基体材料的应变率敏感性;高应变率下,泡沫铝材料率相关性能受基体材料的应变率敏感性以及微结构惯性联合作用,且相对密度较低泡沫铝材料的微结构惯性效应更显著。     

空心微珠/Al复合材料的动态压缩力学性能和吸能特性

利用分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）系统对空心微珠体积分数为0.4的空心微珠/1199Al复合泡沫在1 700~2 900s^-1应变率范围内的动态压缩力学性能、吸能性能进行了研究,还利用SEM扫描电镜对压缩试件断口进行微观组织分析,与准静态条件下材料的压缩力学性能及压缩变形机制进行了对比。结果表明,空心微珠/1199Al复合泡沫是一种应变率敏感材料,与准静态结果相比,在高应变率下复合材料的流动应力和塑性应变有明显的增大,应变率硬化效应对复合材料的流动应力的影响明显大于应变硬化的影响。复合材料的准静态和动态压缩变形机制存在一定差异,动态载荷作用下,空心微珠/1199Al复合泡沫内部空心微珠的压缩和基体材料的充填同时发生,组分之间具有良好的协调变形能力。

     

The influence of mechanical deformation on the electrical properties of single wall carbon nanotubes

Recent experimental studies and atomistic simulations have shown that carbon nanotubes (CNTs) display strong interplay between the mechanical deformation and electrical properties. We have developed a simple and accurate method to determine atom positions in a uniformly deformed CNT via a continuum analysis based on the interatomic potential. A shift vector is introduced to ensure the equilibrium of atoms. Such an approach, involving only three variables for the entire CNT, agrees very well with the molecular mechanics calculations. We then study the effect of mechanical deformation on the band gap change of single wall CNTs under tension, torsion, and combined tension/torsion via the k-space tight-binding method. Prior studies without this shift vector lead to significant overestimation of the band gap change. It is established that the conducting CNTs may easily become semi-conducting ones subject to mechanical deformation, but the semi-conducting CNTs never become conducting ones upon deformation.     

The role of plastic deformation on the impact behaviour of high aspect ratio aluminium foam-filled sections

The main objective of this work is to investigate the role of the plastic deformation of metal foams on the dynamic behaviour of aluminium foam-filled columns with respect to their energy absorbing capabilities. The influence of the cross-section shape as well as other parameters is thoroughly studied. A comparison with correspondent hollow-sections is performed concerning the dissipation of kinetic energy and the obtained deformed profiles. For this particular purpose, three-dimensional finite element modelling dynamic analyses are carried out using ABAQUS/Explicit in order to achieve an in-depth study of the structural crash behaviour, during which energy needs to be absorbed in a controlled manner. A comprehensive numerical study of the crush behaviour of aluminium foam-filled sections undergoing axial compressive loading is performed. The results obtained are also analysed with respect to the reduction in the length of the structural element and impact time, the effect of friction between the foam and the outer skin, the energy decomposition, the role of plastic deformation, the influence of the skin material and impact velocity, and the influence of the shape of the cross-section on the impact behaviour. A comparison with existing analytical expressions is made in order to corroborate the numerical results.