Comment on “Stress boundary conditions for plate bending” by F.Y.M. Wan [Int. J. Solids Struct. 40 (2003) 4107–4123]

By using the reciprocal theorem of elasticity, the author obtained the appropriate stress boundary conditions for the Levy solution for plate bending accurate to all order for plates of general edge geometry and loading. Two special cases of k = 0 (axisymmetric deformation of a circular plate) and k ⩾ 2 (unsymmetric deformation of a circular plate) were discussed in detail in the paper.     

Effect of projectile nose shape,impact velocity and target thickness on deformation behavior of aluminum plates

The present study deals with the experimental and numerical investigations of aluminum target plates impacted by blunt, ogive and hemispherical nosed steel projectiles. The projectiles were normally impacted on the target plates of 0.5, 0.71, 1, 1.5, 2, 2.5 and 3 mm thicknesses at different velocities with the help of a pneumatic gun. Effect of projectile nose shape, impact velocity and plate thickness on the deformation of the target plates was studied. Hemispherical nosed projectile caused highest global deformation (dishing) of the target plates. Ogive nosed projectiles were found to be the most efficient penetrator for the case of plates of thicknesses 0.5, 0.71, 1.0 and 1.5 mm. For the case of plates of thicknesses 2.0, 2.5 and 3.0 mm however, blunt nosed projectiles required least energy to perforate the target plates. The ballistic limit velocity of hemispherical nosed projectiles was found to be highest as compared to the other two projectiles. Finite element analysis of the problem was carried out using ABAQUS finite element code. Results of the numerical analysis were compared with the experiments and good correlation between the two was found.     

Analysis of deformation-induced crack tip toughening in ductile single crystals by nano-indentation

This paper reports on the experimental examination of the deformation characteristics near a crack tip in a cyclically work-hardened copper single crystal using a 2D surface scans with nano-indentation. The experimental methodology enables the characterization of the primary deformation field near a crack tip via the modulation of the imposed secondary deformation field by a nano-indentation. In a heavily deformed 4-point bend specimen, the measurements showed an existence of an asymptotic field around the crack tip at a distance of R ⩾ 2.5J/σ₀. The measurements also showed the qualitative details of toughness evolution within the high-gradient deformation field around the crack tip. The nature of dislocation distribution (i.e. statistically distributed vs. distributions necessitated by geometry) around the crack tip is quantified. The measurements indicate the dominance of the geometrically necessary dislocation within the finite deformation zone ahead of the tip up to a distance of R ≈ 3J/σ₀. Thereafter, it is confined in radial rays coinciding with the sector boundaries around the crack tip. These measurements elucidate the origin of the inhomogeneous hardening and the size dependent macroscopic response close to crack tip.     

A 1-D theory for extensional deformation of a viscoelastic filament under exponential stretching

We consider a viscoelastic filament placed between two coaxial discs, with the bottom plate fixed and the top plate pulled at an exponential rate. Using a slender rod approximation, we derive a one-dimensional (1-D) model which describes the deformation of a viscoelastic filament governed by the Oldroyd-B constitutive model. It is assumed that the flow is axisymmetric and that inertia and gravity are negligible. One solution of the model equations corresponds to ideal uniaxial elongation. A linear stability analysis shows that this solution is unstable for a Newtonian fluid and for viscoelastic filaments with small Deborah number (De ≤ 0.5). For Deborah number greater than 0.5, ideal uniaxial elongation is linearly stable. Numerical solution of the nonlinear equations confirms the result of the linear stability analysis. For initial conditions close to ideal uniaxial flow, our results show that if De > 0.5, the central portion of the filament undergoes considerable strain hardening. As a result, the sample remains almost cylindrical and the deformation approaches pure uniaxial extension as the Hencky strain increases. For De ≤ 0.5, the Trouton ratio based on the effective extension rate at the mid-plane radius gives a much better approximation to the true extensional viscosity than that based on the imposed stretch rate.     

弹性曲梁在机械和热载荷共同作用下的几何非线性模型及其数值解

基于Euler-Bernoulli梁的几何非线性理论,建立了弹性曲梁在任意分布机械载荷和热载荷共同作用下的几何非线性静平衡控制方程。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的相互耦合效应。应用打靶法数值求解了半圆形曲梁在横向均匀升温作用下的非线性弯曲问题,数值比较了轴向伸长对曲梁变形的影响。     

An experimental and numerical study of the localization behavior of tantalum and stainless steel

In general, the shear localization process involves initiation and growth. Initiation is expected to be a stochastic process in material space where anisotropy in the elastic–plastic behavior of single crystals and inter-crystalline interactions serve to form natural perturbations to the material’s local stability. A hat-shaped sample geometry was used to study shear localization growth. It is an axi-symmetric sample with an upper “hat” portion and a lower “brim” portion with the shear zone located between the hat and brim. The shear zone length is 870–890 μm with deformation imposed through a split-Hopkinson pressure bar system at maximum top-to-bottom velocity in the range of 8–25 m/s. We present experimental results of the deformation response of tantalum and 316L stainless steel samples. The tantalum samples did not form shear bands but the stainless steel sample formed a late stage shear band. We have also modeled these experiments using both conductive and adiabatic continuum models. An anisotropic elasto-viscoplastic constitutive model with damage evolution was used within the finite element code EPIC. A Mie-Gruneisen equation of state and the rate and temperature sensitive MTS flow stress model together with a Gurson flow surface were employed. The models performed well in predicting the experimental data. The numerical results for tantalum suggested a maximum equivalent strain rate on the order of 7 × 10⁴ s⁻¹ in the gage section for an imposed top surface displacement rate of 17.5 m/s. The models also suggested that for an initial temperature of 298 K a temperature in the neighborhood of 900 K was reached within the shear section. The numerical results for stainless steel suggest that melting temperature was reached throughout the shear band shortly after peak load. Due to sample geometry, the stress state in the shear zone was not pure shear; a significant normal stress relative to the shear zone basis line was developed.     

Relating inhomogeneous deformation to local texture in zirconium through grain-scale digital image correlation strain mapping experiments

The effect of local texture on inhomogeneous plastic deformation is studied in zirconium subjected to uniaxial compression. Cross-rolled commercially pure Zr 702 plate that had a strong basal (0 0 0 1) texture through the plate thickness, and a non-basal texture in cross-section, was obtained. At a compressive strain rate of 1 s^?1, samples loaded either in the through-thickness or in-plane directions exhibited significant differences in yield strength, hardening response and failure mechanisms. These macroscopic differences are related to microstructural features by combining information from electron backscattered diffraction with real time in situ imaging and subsequent full-field strain measurements obtained using digital image correlation. Experimental results indicate that the through-thickness loaded zirconium samples, which show a strong basal-texture in the loading direction, do not deform homogeneously – implying the lack of a representative volume element. The detailed surface deformation fields provided by digital image correlation allow for a qualitative and quantitative analysis of the relationship between grain orientation and patterns of deformation bands that form as the precursors to development of an adiabatic shear band in the through-thickness loaded sample. For the in-plane loaded samples, inhomogeneities still exist at the microscale, but the collective behavior of several grains leads to a homogeneous response at the macroscale. It is observed that local texture for hcp polycrystals, which are significantly slip restricted, can directly affect both local and global response, even at low to moderate plastic strains.     

Effect of elastic anisotropy on surface pattern evolution of epitaxial thin films

Stress-induced surface instability and evolution of epitaxial thin films leads to formation of a variety of self-assembled surface patterns with feature sizes at micro- and nanoscales. The anisotropy in both the surface and bulk properties of the film and substrate has profound effects on the nonlinear dynamics of surface evolution and pattern formation. Experimentally it has been demonstrated that the effect of anisotropy strongly depends on the crystal orientation of the substrate surface on which the film grows epitaxially. In this paper we develop a nonlinear model for surface evolution of epitaxial thin films on generally anisotropic crystal substrates. Specifically, the effect of bulk elastic anisotropy of the substrate on epitaxial surface pattern evolution is investigated for cubic germanium (Ge) and SiGe films on silicon (Si) substrates with four different surface orientations: Si(0 0 1), Si(1 1 1), Si(1 1 0), and Si(1 1 3). Both linear analysis and nonlinear numerical simulations suggest that, with surface anisotropy neglected, ordered surface patterns form under the influence of the elastic anisotropy, and these surface patterns clearly reflect the symmetry of the underlying crystal structures of the substrate. It is concluded that consideration of anisotropic elasticity reveals a much richer dynamics of surface pattern evolution as compared to isotropic models.     

基于近场动力学的薄板弯曲大变形及断裂分析

薄板结构仅在较小的荷载下就能产生较大的位移、旋转,甚至引发结构产生裂纹并扩展,进而发生结构整体断裂,因此,建立薄板结构在大变形过程中的裂纹扩展及断裂仿真模型,具有重要的工程实际意义.文章建立了用于薄板结构几何大变形和断裂分析的近场动力学(PD)和连续介质力学(CCM)耦合模型.首先基于冯·卡门假设,采用更新的拉格朗日法得到薄板在几何大变形增量步下的虚应变能密度增量公式,并利用虚功原理和均质化假设求出几何大变形微梁键的本构模型参数;接着分别建立几何大变形薄板PD模型与CCM模型的虚应变能密度增量,并建立了薄板几何大变形PD-CCM耦合模型;最后模拟了薄板结构在横向变形作用下的渐进断裂过程,得到与实验结果高度一致的仿真结果,验证了所提出的几何非线性PD-CCM耦合模型的精度.结果表明:本文所提出的薄板PD-CCM耦合模型具有简单高效,无需考虑材料参数限制和边界效应的特点,可以很好地用于预测薄板结构在几何大变形过程中的局部损伤和结构断裂,有利于薄板结构的断裂安全评价和理论发展.     

Thermo-viscoelastic and viscoplastic behavior of high-density polyethylene

Observations are reported on high-density polyethylene in uniaxial tensile tests with constant strain rates and relaxation tests at various temperatures ranging from 25 to 90 °C. A constitutive model is derived for the nonlinear viscoelastic and viscoplastic behavior of semi-crystalline polymers at three-dimensional deformations. Adjustable parameters in the stress–strain relations are found by fitting the experimental data. It is demonstrated that (i) the model correctly approximates the observations and (ii) material parameters are independent of strain rate and change consistently with temperature.     

Dual-element directional shock wave attenuators

This study explores the interaction between shock waves and dual-element porous plates used to ameliorate the hazardous effects of these waves. Tests were performed in a shock tube to determine the effects that a pair of porous plates with directional resistance properties had on the initial peak pressure and impulse amelioration experienced by an end wall. Mild steel test specimens, ranging in porosity values from 6.6% to 41.1%, were mounted two at a time at different spacings in the test section. Each plate had directional properties, i.e. resistance to flow was different for flow coming from either side. Four plates were used, and 48 plate configurations were tested. Side wall and end wall pressure measurements and schlieren photographs were taken of the interactions. Tests were run at Mach numbers of 1.23, 1.35 and 1.42. The separation distances between the plate specimens were varied between 30 mm and 60 mm; however, the distance between the downstream plate and the end wall was kept constant at 140 mm for all tests.Both the initial peak pressure and impulse amelioration values were found to be dependant on the plate combination porosity. As the porosity of the combination increased, the amelioration values decreased. Complementary plate combinations produced differing results as different wave interactions occur when plate positions were interchanged. The porosity of the combined plates was found to have an overriding influence on the end wall peak pressure and impulse amelioration values when compared to the effect that plate arrangement (i.e. geometrical influences) had. Impulse amelioration values were found to increase as the separation distance between plates were increased. The amplitude of the end wall pressure trace was found to increase as the incident Mach number was increased. The significant attenuation of the incident shock wave obtained during this study is attributed to the system of multiple reflected and transmitted waves that are produced by the presence of the plate specimens in series. This increases the frequency of shock wave and barrier interactions, when compared to just using a single barrier.     

Coupled annealing temperature and layer thickness effect on strengthening mechanisms of Ti/Ni multilayer thin films

A systematic study was performed on mechanical and microstructural properties of Ti/Ni multilayers with layer thickness from 200 nm to 6 nm and annealing temperature from room temperature to 500 °C. Based on the observed hardness evolution, a coupled layer-thickness and annealing-temperature dependent strengthening mechanism map is proposed. For as-deposited films, the deformation behavior follows the traditional trend of dislocation mediated strengthening to grain boundary mediated softening with decreasing layer thickness. For annealed films, grain boundary relaxation is considered to be the initial strengthening mechanism with higher activation temperature required for thicker layers. Under further annealing, solid solution hardening, intermetallic precipitation hardening, and fully intermixed alloy structure continue to strengthen the thin layered films, while recrystallization and grain-growth lead to the eventual softening of thick layered films. For the films with intermediate layer thickness, a strong orientation dependent hardness behavior is exhibited under high temperature annealing due to mechanism switch from grain growth softening to intermetallic precipitation hardening when changing the loading orientation from perpendicular to parallel to the layer interfaces.     

The effect of initial geometric imperfection on the nonlinear resonance of functionally graded carbon nanotube-reinforced composite rectangular plates

The purpose of the present study is to examine the impact of initial geometric imperfection on the nonlinear dynamical characteristics of functionally graded carbon nanotube-reinforced composite (FG-CNTRC) rectangular plates under a harmonic excitation transverse load. The considered plate is assumed to be made of matrix and single-walled carbon nanotubes (SWCNTs). The rule of mixture is employed to calculate the effective material properties of the plate. Within the framework of the parabolic shear deformation plate theory with taking the influence of transverse shear deformation and rotary inertia into account, Hamilton’s principle is utilized to derive the geometrically nonlinear mathematical formulation including the governing equations and corresponding boundary conditions of initially imperfect FG-CNTRC plates. Afterwards, with the aid of an efficient multistep numerical solution methodology, the frequency-amplitude and forcing-amplitude curves of initially imperfect FG-CNTRC rectangular plates with various edge conditions are provided, demonstrating the influence of initial imperfection, geometrical parameters, and edge conditions. It is displayed that an increase in the initial geometric imperfection intensifies the softening-type behavior of system, while no softening behavior can be found in the frequency-amplitude curve of a perfect plate.     

含冷却孔镍基合金次级取向效应的应变梯度晶体塑性有限元研究

单晶镍基合金具有优异的耐高温、高强、高韧等性能, 这些力学性能受制造过程引入的次级取向和冷却孔的影响. 已有研究大多关注单孔薄板的变形机理和力学性能, 而工程中应用的往往是多孔薄板, 当前亟需阐明多孔的塑性滑移带变形机理、次级取向效应以及冷却孔引起的应变梯度效应. 文章采用基于位错机制的非局部晶体塑性本构模型对含冷却孔镍基单晶薄板的单拉变形进行了数值模拟. 此模型基于塑性滑移梯度与几何必需位错的关系引入了位错流动项, 因此可有效刻画非均匀变形过程中的应变梯度效应. 为了全面揭示含孔镍基薄板的次级取向效应, 系统研究了[100]和[110]取向(两种次级取向)下镍基薄板的单拉变形行为, 并重点探究了在两种次级取向下冷却孔数量对薄板塑性行为的影响. 此外, 还分析了镍基合金板变形过程中各个滑移系上分切应力变化、主导滑移系开动以及几何必需位错密度的演化过程, 并讨论了塑性滑移量及其分布特征对不同次级取向镍基合金板强度的影响. 研究表明, 单孔和多孔的[110]薄板抗拉强度均低于[100]薄板, 多孔薄板的塑性变形过程比单孔薄板更为复杂且受次级取向影响更大, 并且发生滑移梯度位置主要位于冷却孔附近以及塑性滑移带区域. 研究结果可为工程中镍基合金的设计和服役提供理论指导.      

Nonlinear vibrations of FGM rectangular plates in thermal environments

Geometrically nonlinear vibrations of FGM rectangular plates in thermal environments are investigated via multi-modal energy approach. Both nonlinear first-order shear deformation theory and von Karman theory are used to model simply supported FGM plates with movable edges. Using Lagrange equations of motion, the energy functional is reduced to a system of infinite nonlinear ordinary differential equations with quadratic and cubic nonlinearities. A pseudo-arclength continuation and collocation scheme is used and it is revealed that, in order to obtain the accurate natural frequency in thermal environments, an analysis based on the full nonlinear model is unavoidable since the plate loses its original flat configuration due to thermal loads. The effect of temperature variations as well as volume fraction exponent is discussed and it is illustrated that thermally deformed FGM plates have stronger hardening behaviour; on the other hand, the effect of volume fraction exponent is not significant, but modal interactions may rise in thermally deformed FGM plates that could not be seen in their undeformed isotropic counterparts. Moreover, a bifurcation analysis is carried out using Gear’s backward differentiation formula (BDF); bifurcation diagrams of Poincaré maps and maximum Lyapunov exponents are obtained in order to detect and classify bifurcations and complex nonlinear dynamics.     

Thermo-mechanical response of nylon 101 under uniaxial and multi-axial loadings: Part I,Experimental results over wide ranges of temperatures and strain rates

A comprehensive study of the thermo-mechanical response of a thermoplastic polymer, nylon 101 is presented. Quasi-static and dynamic compression uniaxial and multi-axial experiments (stress states) were performed at a wide range of strain rates (10⁻⁵ to 5000 s⁻¹) and temperatures (−60 to 177 °C or −76 to 350 °F). The material is found to be non-linearly dependent on strain rate and temperature. The change in volume after plastic deformation is investigated and is found to be negligibly small. The relaxation and creep responses at room temperature are found to be dependent on strain rate and the stress–strain level at which these phenomena are initiated. Total deformation is decomposed into visco-elastic and visco-plastic components; these components have been determined at different levels of deformation. Results from non-proportional uniaxial to biaxial compression, and torsion experiments, are also reported for three different strain rates at room temperature. It is shown that nylon 101 has a response dependent on the hydrostatic pressure.     

Nonlinear bending of size-dependent circular microplates based on the modified couple stress theory

The present study proposes a nonclassical Kirchhoff plate model for the axisymmetrically nonlinear bending analysis of circular microplates under uniformly distributed transverse loads. The governing differential equations are derived from the principle of minimum total potential energy based on the modified couple stress theory and von Kármán geometrically nonlinear theory in terms of the deflection and radial membrane force, with only one material length scale parameter to capture the size-dependent behavior. The governing equations are firstly discretized to a set of nonlinear algebraic equations by the orthogonal collocation point method, and then solved numerically by the Newton–Raphson iteration method to obtain the size-dependent solutions for deflections and radial membrane forces. The influences of length scale parameter on the bending behaviors of microplates are discussed in detail for immovable clamped and simply supported edge conditions. The numerical results indicate that the microplates modeled by the modified couple stress theory causes more stiffness than modeled by the classical continuum plate theory, such that for plates with small thickness to material length scale ratio, the difference between the results of these two theories is significantly large, but it becomes decreasing or even diminishing with increasing thickness to length scale ratio.     

Influence of temperature and strain rate on the mechanical behavior of three amorphous polymers: Characterization and modeling of the compressive yield stress

Uniaxial compression stress–strain tests were carried out on three commercial amorphous polymers: polycarbonate (PC), polymethylmethacrylate (PMMA), and polyamideimide (PAI). The experiments were conducted under a wide range of temperatures (−40 °C to 180 °C) and strain rates (0.0001 s⁻¹ up to 5000 s⁻¹). A modified split-Hopkinson pressure bar was used for high strain rate tests. Temperature and strain rate greatly influence the mechanical response of the three polymers. In particular, the yield stress is found to increase with decreasing temperature and with increasing strain rate. The experimental data for the compressive yield stress were modeled for a wide range of strain rates and temperatures according to a new formulation of the cooperative model based on a strain rate/temperature superposition principle. The modeling results of the cooperative model provide evidence on the secondary transition by linking the yield behavior to the energy associated to the β mechanical loss peak. The effect of hydrostatic pressure is also addressed from a modeling perspective.     

Experimental research on heat transfer of natural convection in vertical rectangular channels with large aspect ratio

This work presents the experimental research on the steady laminar natural convection heat transfer of air in three vertical thin rectangular channels with different gap clearance. The much higher ratio of width to gap clearance (60–24) and the ratio of length to gap clearance (800–320) make the rectangular channels similar with the coolant flow passage in plate type fuel reactors. The vertical rectangular channels were composed of two stainless steal plates and were heated by electrical heating rods. The wall temperatures were detected with the K-type thermocouples which were inserted into the blind holes drilled in the steal plates. Also the air temperatures at the inlet and outlet of the channel were detected. The wall heat fluxes added to the air flow were calculated by the Fourier heat conduction law. The heat transfer characteristics were analyzed, and the average Nusselt numbers in all the three channels could be well correlated with the Rayleigh number or the modified Rayleigh number in a uniform correlation. Furthermore, the maximum wall temperatures were investigated, which is a key parameter for the fuel’s integrity during some accidents. It was found that even the wall heat flux was up to 1500 W/m², the maximum wall temperature was lower than 350 °C. All this work is valuable for the plate type reactor’s design and safety analysis.