On the importance of thermo-elastic cooling in the fracture of glassy polymers at high rates

In a previous thermo-mechanical analysis [Estevez, R., Basu, S., van der Giessen, E., 2005. Analysis of temperature effects near mode I cracks in glassy polymers. Int. J. Fract. 132, 249–273] in which shear yielding of the bulk and failure by crazing were accounted for, we examined which of these two viscoplastic processes contributed to heat in mode I fracture. The present study completes this work by investigating the conditions for thermo-elastic cooling prior to crack propagation as reported experimentally by Rittel [Rittel, D., 1998. Experimental investigation of transient thermo-elastic effects in dynamic fracture. Int. J. Solids Struct. 35, 2959–2973] and Bougaut and Rittel [Bougaut, O., Rittel, D., 2001. On crack tip cooling during dynamic crack propagation. Int. J. Solids Struct. 38, 2517–2532] on high strain rate loading of PMMA. To this end, coupled thermo-mechanical finite element simulations are carried out by accounting for the thermo-elastic source, in addition to the heat sources related to shear yielding and crazing. The bulk as well as cohesive zone parameters for crazing realistically describe PMMA as they are obtained from detailed calibration experiments. Our results show that if significant thermo-elastic cooling has to be observed in the vicinity of the crack tip of a polymeric material, suppression of shear yielding as well as suppression of crazing is necessary. It seems that at these high strain rates a brittle fracture mechanism activated at very high stresses takes over from crazing, or at least that craze initiation occurs for stress levels very different to those for quasi-static conditions.     

A generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials

In this paper, a generalized anisotropic hardening rule based on the Mroz multi-yield-surface model for pressure insensitive and sensitive materials is derived. The evolution equation for the active yield surface with reference to the memory yield surface is obtained by considering the continuous expansion of the active yield surface during the unloading/reloading process. The incremental constitutive relation based on the associated flow rule is then derived for a general yield function for pressure insensitive and sensitive materials. Detailed incremental constitutive relations for materials based on the Mises yield function, the Hill quadratic anisotropic yield function and the Drucker–Prager yield function are derived as the special cases. The closed-form solutions for one-dimensional stress–plastic strain curves are also derived and plotted for materials under cyclic loading conditions based on the three yield functions. In addition, the closed-form solutions for one-dimensional stress–plastic strain curves for materials based on the isotropic Cazacu–Barlat yield function under cyclic loading conditions are summarized and presented. For materials based on the Mises and the Hill anisotropic yield functions, the stress–plastic strain curves show closed hysteresis loops under uniaxial cyclic loading conditions and the Masing hypothesis is applicable. For materials based on the Drucker–Prager and Cazacu–Barlat yield functions, the stress–plastic strain curves do not close and show the ratcheting effect under uniaxial cyclic loading conditions. The ratcheting effect is due to different strain ranges for a given stress range for the unloading and reloading processes. With these closed-form solutions, the important effects of the yield surface geometry on the cyclic plastic behavior due to the pressure-sensitive yielding or the unsymmetric behavior in tension and compression can be shown unambiguously. The closed form solutions for the Drucker–Prager and Cazacu–Barlat yield functions with the associated flow rule also suggest that a more general anisotropic hardening theory needs to be developed to address the ratcheting effects for a given stress range.     

考虑拉压强度差效应的厚壁圆筒承载能力分析

应用双剪统一强度理论，考虑材料的拉压异性和同性，推导了在内压力和轴力联合作用下的厚壁圆筒的塑性极限载荷表达式．在该表达式中，当反映中间主应力效应的系数取不同的值时，就能得到按Tresca屈服准则、线性逼近的Mises屈服准则和双剪应力屈服准则的计算结果，并且绘制了在相应准则下的极限应力线图．从而可知：在三维应力状态下，应用该理论，可以获得极限载荷分析的精确解；极限载荷线图与三种屈服准则的屈服曲线是相吻合的；计算的结果可以用于拉压异性和同性的材料，为工程应用提供了理论依据．     

Experimental determination of cohesive failure properties of a photodegradable copolymer

In this paper we present a methodology to measure the material traction-separation relation for a poly(ethylene carbon monoxide) copolymer, ECO. This material exhibits a ductile-to-brittle transition when subjected to ultraviolet irradiation and undergoes a change of failure mechanism, from shear yielding to crazing, in the process. Single edge notched tension specimens of ECO irradiated for 50 h were used to generate slow-speed stable crack growth, predominantly from material crazing. Full-field measurements of the in-plane deformation around the growing crack tip were performed using the optical technique of digital image correlation. A multicamera setup was used in which simultaneous measurement of both the far-field displacement and that directly surrounding the craze was performed. The far-field results were used to obtain a value of the energy release rate supplied to the crack tip region. The near-tip results were used to extract a material traction-separation law in the regime of steady-state crack growth. A softening traction-separation relation was measured. The area under the traction-separation curve was then compared with the simultaneous far-field measurements and the agreement was very good (within 6.5%) validating the experimental methodology used.     

饱和超固结黏性土的三剪弹塑性本构模型研究

针对饱和超固结黏性土现有下加载面修正剑桥模型中破坏应力比为定值、土体黏聚力为零,以及不能准确反映不同应力状态下土的强度差异这些问题,基于三剪统一强度准则以及应力坐标平移法得到了扩展破坏应力比,其特点是能更好地反映应力状态变化以及土体黏聚力的影响。在此基础上提出了饱和超固结黏性土的三剪弹塑性本构模型,该模型的特点是能描述土体受力时的中间主应力效应,应力区间效应和拉压差影响,同时也能更好地考虑土体黏聚力的影响。基于该模型对ABAQUS软件进行了二次开发,并利用其模拟了饱和超固结黏性土在排水和不排水条件下的真三轴和常规三轴压缩试验特性。对常规三轴压缩条件下土体力学特性作了模拟和试验结果对比。结果表明所提模型能很好地反映不同超固结比下土体的变形、剪胀、孔隙水压力变化特性。     

Refined Dugdale plastic zones of an external circular crack

The refined Dugdale-type plastic zones ahead of an external circular crack, subjected to a uniform displacement at infinity, are evaluated both analytically and numerically. The analytical method utilizes potential theory in classical linear elasticity with emphasis on the contrast from the internal crack problem. A closed-form solution to the mixed boundary problem is obtained to predict the length of the plastic zone for a Tresca yield condition. The analytical solution is also used to benchmark the results obtained from the numerical method, which show good agreement. Through an iterative scheme, the numerical technique is able to estimate the size of crack tip plasticity zone, which is governed by the non-linear von Mises criterion. The relationships between the applied displacement and the length of the plastic zone are compared for the different yielding conditions. Computational modeling has demonstrated that the plastic constraint effect based on the true yield condition can significantly influence the load-bearing capacity. It is also discovered from the comparative study that the stress components predicted by the three different yield conditions may differ notably; however, the stress triaxiality in the ligament region has only small deviations. The proposed study may find applications in predicting the plastic flow in a circumferentially notched round bars under tension.     

屈服准则及切线模量修正的弹塑性计算模型

为提高不同应力状态下金属弹塑性行为的模拟精度,采用含应力三轴度修正的von Mises屈服准则,材料弹塑性本构关系在等效应力$\!$-$\!$-$\!$等效应变曲线基础上,提出一个切线模量为主应力函数的理论修正项. 将这两个修正项结合,以子程序形式编程嵌入ABAQUS主程序,以此模拟几种不同形状试样的弹塑性行为,并将其他屈服准则在单一曲线假设下编程与之对比. 模拟结果与真实试验结果对比发现,对于屈服而言,含应力三轴度修正的vonMises屈服准则比其他屈服理论准确;对于弹塑性阶段计算而言,提出的切线模量为主应力函数这一假设比单一曲线假设更加接近真实试验.      

Localization in elasto-plastic materials: Influence of the plasticity yield surface in biaxial loading conditions

The influence of the plasticity yield surface on the development of instabilities in plane plates in biaxial loading is analyzed in order to understand and simulate the localization pattern observed in an expanding hemisphere experiment. First, a criterion for the activation of slip bands is formulated in the form of a critical hardening coefficient: it is particularized to the Von Mises and Tresca surfaces. In the Von Mises case, the criterion gives a strongly negative hardening coefficient in biaxial loading conditions different from the ones of plane strain. In the Tresca case, the criterion is fulfilled for a perfectly plastic material in uniaxial and biaxial loading; besides, in equi-biaxial loading, two possible orientations for slip bands are exhibited; this can be understood, with a few approximations, by the existence of a vertex point on the Tresca yield surface which give additive degrees of freedom for the direction of the plastic strain rate. Second, the development of localization in the loading conditions met in an expanding hemisphere experiment is simulated using both plasticity yield surfaces; whereas the Von Mises simulation does not localize, the Tresca simulation exhibits a pattern composed of a network of shear bands of different orientations; this pattern is not far from the pattern observed experimentally.     

基于D-P准则的压力相关材料结构拓扑优化

基于描述材料力学行为的Drucker-Prager(D-P)屈服准则, 研究了压力相关材料连续体结构拓扑优化设计问题的数学模型和数值算法. 以单元材料人工密度为设计变量, 结合SIMP惩罚模型和多孔微结构局部应力插值模型, 建立了以材料体积最小化为目标、考虑材料D-P屈服条件约束的优化问题数学模型. 利用\varepsilon-松弛方法消除奇异解现象, 采用伴随法有效推导约束函数灵敏度计算公式, 运用基于梯度的连续变量优化算法迭代求解优化问题. 数值算例验证了优化模型的正确性及数值算法的有效性, 并通过与von Mises应力约束优化结果的比较, 说明了材料的压力相关特性会对结构最优拓扑产生重要影响. 该方法设计出的最优拓扑由于充分利用了压力相关材料的抗压能力, 因而更为合理和实际.      

The yield surface of viscoelastic and plastic fluids in a vane viscometer

Vane viscometers are often used to investigate the low shear rate properties of plastic fluids. The shear stress is determined by assuming that the material is held in the space between the vane blades so that it behaves like a rigid cylinder. Experimental evidence supports this assumption and the aim of the present study is to model numerically the yield process in a vane rheometer using viscoelastic and plastic fluids. The finite element method has been used to model the behavior of Herschel-Bulkley (Bingham), Casson and viscoelastic (Maxwell type) fluids. The penalty function approach for the pressure approximation and a rotating reference frame are used together with fine meshes containing more than 1300 elements. The results show that for Herschel-Bulkley (Bingham), and Casson fluids a rotating rigid cylinder of fluid is trapped inside the periphery of the vane, the shear stress is uniformly distributed over the surface of the cylinder. Finally a modified second order fluid is used to simulate the viscoelastic behaviour, anticipated to be an intermediate between the elastic deformation and the plastic flow, to provide a more realistic simulation of the yield process about a vane. In this case, contrast with the concentration of the elastic strain rate at the blade tips, a nearly uniform distribution of the plastic shear rate is still found. This implies that the plastic shear always distributes uniformly during the entire yielding process. Evidently the assumption of uniform shear on a rotating cylinder of material occluded in the blades of a vane is a valid and useful model for many types of fluid possessing a yield stress.     

Multiaxial yield surface of transversely isotropic foams: Part I—Modeling

A new yield criterion is proposed for transversely isotropic solid foams. Its derivation is based on the hypothesis that the yielding in foams is driven by the total strain energy density, rather than a completely phenomenological approach. This allows defining the yield surface with minimal number of parameters and does not require complex experiments. The general framework used leads to the introduction of new scalar measures of stress and strain (characteristic stress and strain) for transversely isotropic foams. Furthermore, the central hypothesis that the total strain energy density drives yielding in foams ascribes to the characteristic stress an analogous role of von Mises stress in metal plasticity. Unlike the overwhelming majority of yield models in literature the proposed model recognizes the tension–compression difference in yield behavior of foams through a linear mean stress term. Predictions of the proposed yield model are in excellent agreement with the results of uniaxial, biaxial and triaxial FE analyses implemented on both isotropic and transversely isotropic Kelvin foam models.     

基于SPH-FEM耦合方法的柔性导爆索分离装置爆炸分离过程数值模拟

为深入研究柔性导爆索在爆炸分离装置中的作用过程和机理,提出一种改进的光滑粒子流体动力学方法（smoothed particle hydrodynamics, SPH）与有限单元法（ finite element method, FEM）耦合算法。新方法中不仅包含导爆索模拟的SPH方法与分离装置模拟的FEM方法之间的接触算法,同时将完全损伤失效后的单元采用转化算法动态转化成SPH粒子继续参与计算,转化后的粒子与未转化的有限单元之间采用接触算法计算。采用该方法对环型和平板型两种爆炸分离结构的分离过程进行了数值模拟,验证了新方法的准确性与问题适用性;分析了分离板的变形断裂及损伤碎片的飞溅过程,得到了分离装置表面不同时刻的应力分布、损伤因子的变化趋势、von Mises应力的变化趋势;探讨了炸药在不同比内能情况下单元的屈服损伤速度、碎片的飞溅位移速度。     

On the application of FEM to deformation of high-density polyethylene

This paper presents a study that uses finite element method (FEM) to simulate deformation behaviour of high-density polyethylene (HDPE) when subjected to tensile loading, either without or with the presence of pre-cracks. For the former, dog-bone-shaped model of HDPE was deformed under uni-axial tensile (UT) loading beyond the initial yielding, to generate stable necking in the gauge section. The paper proposes a simple correction factor to determine the equivalent von Mises stress that is needed for the input to the FEM model, in order to generate the same loading level as that observed experimentally for neck propagation. The paper points out that such consistency in the loading level could not be generated in the past mainly because of a misconception that axial stress in the neck was regarded as the equivalent of the von Mises stress. The study also explored the consideration of crack growth in double-edge-notched tensile (DENT) specimen, and showed that the conventional von Mises yield function, with the assumption of isotropic work hardening, cannot be directly applied to simulate the deformation behaviour of DENT specimen. Instead, empirical parameters were employed to reflect the polymer orientation during the necking process. The paper shows that FEM models for both UT and DENT tests can reproduce the experimental load–displacement curves quite accurately, and concludes that with a proper yield function to reflect the deformation involved in the mechanical tests, the deformation behaviour observed experimentally can be accurately mimicked using the FEM simulation.     

Collapse of thick-walled cylinders under external pressure

Hydrostatic collapse tests performed on thick-walled capped cylinders are described. Finite-element predictions which incluce the effects of end-cap stiffening, cross-section ovalities and material strain hardening are compared to experimental results. The analyses correctly predict the sequence of events leading to collapse, but experimental failure pressures are significantly below predictions. It is concluded that the von Mises yield criterion used in the analysis did not accurately represent the yield behavior of the 1018 steel tubing material of the test-cylinders for the triaxial-stress conditions of interest.     

Experimental studies on the yield behavior of ductile and brittle aluminum foams

Specimen size effects are a major cause of the unreliability of foam models in finite element codes. Here, the modified Arcan apparatus is used to investigate the biaxial yielding of ductile and brittle Al foams. This apparatus subjects a central section of a “butterfly-shaped” specimen to a uniform state of plane stress. The stresses have local maxima at the central section, thus ensuring that yielding occurs there. A yield envelope, which directly relates to the crushing process, can then be determined. Size effects are introduced when using conventional methods such as tri-axial or plate-shear tests. In such tests, averages of stress and strain are measured. These measures do not represent the actual yield event, because foam's internal structure is inhomogeneous and so is the deformation field. Strain localization and failure can occur at any weak layer of cells in the bulk. In this study, we have performed a series of biaxial tests on isotropic Alporas and anisotropic Hydro closed-cell Al foams of approximately equal densities. Alporas failed locally by a ductile phenomenon of progressive crushing of cells. It also possessed uniaxial strength asymmetry. Hydro specimens parallel and perpendicular to ‘foam rise’ were investigated. The Hydro foam developed a local, characteristic brittle crack at loads in the vicinity of the yield point. Phenomenological yield surfaces, which incorporate these features are obtained for the foams, and show dependence on both the deviatoric and hydrostatic stresses. We also provide expressions for the shear and hydrostatic strengths in terms of the uniaxial strengths. Finally, the size-independence of the yield surface is verified using the uniaxial compression of tapered specimens.     

弹体贯穿混凝土数值模拟的改进材料模型

研究混凝土结构在冲击载荷下的力学特性对武器以及防护结构的设计和评估具有重要意义,而合适的材料模型可以更准确地预测混凝土结构的力学行为和破坏模式。因此,本文中提出了一种改进的混凝土塑性损伤材料模型来描述其在冲击载荷下的力学响应。该改进模型考虑了压力-体积应变关系、应变率效应、洛德角效应和塑性损伤累积对混凝土材料力学特性的影响,并引入了一个与损伤相关的硬化/软化函数来描述压缩状态下的应变硬化和软化行为。随后,通过对3个独立的强度面进行线性插值得到了该改进模型的破坏强度面,并采用部分关联流动法则考虑了混凝土材料的体积膨胀特性。最后,开展了单个单元在不同加载条件下和弹体贯穿钢筋混凝土靶的数值模拟,验证了该改进模型的可行性、准确性以及预测性能提升。     

A yield-surface corner lowers the buckling stress of an elastic-plastic plate under compression

The bifurcation stress is calculated for a simply-supported elastic-plastic plate in uniaxial compression, in the case when the stress is at a vertex of the yield surface which is locally similar to that of Tresca. The effect of coupled hardening between yield-surface facets which meet at the vertex is included. It is found that the bifurcation stress may be substantially lower than that associated with the von Mises yield surface. The latter is known to give results too high to be reconciled with experiment. Reductions of order 10–30 per cent are shown to be unexceptional, even with the retention of an elastic value for the shear modulus.     

岩土的一种强度准则及其变换应力法

基于岩土摩擦性,假设岩土破坏是由其物理空间内特征面上的应力比决定,提出了等效应力比的概念,即物理空间特征面上的剪应力合力与正应力合力的比值.在二维条件下,等效应力比可表示为σ-τ坐标系下与摩尔圆相切的直线扣除截距正切值;在三维条件下,假设在XYZ空间内存在一三维物理空间平面,此三维空间特征平面的等效应力比为影响材料强度特性的决定性因素,基于上述三维空间特征面建立了强度准则并称之为a准则.SMP准则以及广义Mises准则都是a准则的特例,当二维坐标中的截距为零时,则强度准则退化为SMP(spatially mobilized plane)强度准则,而当正切角为零时,则强度准则退化为广义Mises准则.而当截距与外切角均不为零时,则强度准则为介于上述两者之间的一种强度准则,在偏平面上为介于SMP曲边三角形与广义Mises圆形之间的曲边三角形.在子午面上,采用考虑岩土压剪耦合的屈服准则,破坏准则采用幂函数表达式.在偏平面上提出了基于a准则的形状函数,并采用真三维应力状态表示的破坏强度准则表示在三轴压缩路径下以p,q二维应力变量表达的准则公式,推导得到了基于a准则的变换应力公式,可简单地将一般以p,q为基本变量的二维模型转变为三维应力模型.通过强度以及多种应力路径的测试对比,验证了a准则及基于该准则的变换应力公式的合理性.     

Development of an analytical framework for viscoelastic corrugated-core sandwich plates and validation against FEM

In this study, an analytical procedure for the bending problem of a viscoelastic sandwich plate with a corrugated core is presented. Reissner–Mindlin plate theory and N-termed Prony series are employed to define the elastic and time-dependent contributions of the governing equations, respectively. Three different corrugation patterns, i.e., rectangular, trapezoidal, and triangular, are examined. Moreover, the structure is analyzed under both simply support and clamp boundary conditions. The calibrated material parameters of polymethyl methacrylate (PMMA) for the Generalized Maxwell rheological model are employed to show the viscoelastic response of the structure. A 3D finite element simulation of the problem is also conducted to confirm the accuracy of the analytical formulation. The two well-known creep and stress relaxation phenomena of the viscoelastic materials are examined for the mentioned corrugation cores and both boundary conditions analytically and numerically. The time-dependent dimensionless deflection and resultant von Mises stress distributions are provided. Besides, the variation of the results with various rise-times and applied load are studied in detail. The von Mises stress contours of the upper surface of the structure at the end of the creep test are also presented. The finite element method outcomes verify the analytical results with excellent compatibility. The proposed analytical procedure can be used as an efficient tool to study the effects of various parameters such as material, geometrical constants, and corrugation pattern on bending of viscoelastic sandwich plates with corrugated core problems for design and optimization, which involves a high number of simulations.