An analytical and computational study of crack initiation under transient creep conditions

An analytical method has been developed to predict creep crack initiation (CCI), based on the accumulation of a critical level of damage at a critical distance. The method accounts for the re-distribution of stress from the elastic or elastic–plastic field, experienced on initial loading, to a steady state creep stress distribution, via a transient creep region. The method has been applied to predict CCI times in a fracture specimen of type 316H stainless steel at 550 °C. The failure model has been also been implemented into a finite element (FE) framework. Reasonable and conservative predictions of CCI time can be obtained from the analytical solution relative to FE solutions. Conservative predictions of experimental CCI times are obtained when stress redistribution is taking into account. However, CCI times predicted from a steady state creep model are found to be non-conservative.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Stress analysis of notched anisotropic finite plates under mechanical and hygrothermal loads

Summary Notch-induced stress concentrations in anisotropic composite materials depend on their directional material properties, especially for uniaxially reinforced composites with high-modulus fibres. The design of notched high-performance composites requires therefore a special proof of their notched strength, which includes the structural parameters of the fibre/matrix combination, fibre orientation and layer arrangement. The assessment of the effects of the finite outer boundary is of practical importance when dimensioning critical notched regions. An anisotropic plate with finite dimensions and a hole in its center will be used here to model stress concentrations. The calculation is based on conformal mappings combined with complex-valued stress functions. The outer boundary is described using point-matching and the least-squares method. The solutions are conducive to the assessment of the essential influencing factors of material properties, geometry and loads. Notched finite plates made of fibre/matrix composites, mainly carbon-fibre reinforced polymers, will be presented as illustrations. Received 29 June 1998; accepted for publication 22 October 1998     

Prevention of thermal bending of multilayered beams and plates

A method is described to prevent bending in multilayered beams and plates of different isotropic materials with uniform and nonuniform temperature distribution through the thickness. The method involved the addition of an extra layer to the multilayered beams or plates. With the proper selection of the thermoelastic properties, the added layer would eliminate the curvature produced prior to this addition. A complete analysis for the determination of the various thermoelastic parameters of the extra layer was made. In addition, to ensure that the multilayered beams and plates actually remained straight, a thermoelastic analysis was performed for the solution of thermal stresses and strains in the laminate. The results gave assurance to the straightness of the laminate since the calculated strains have the same value throughout the thickness. The solutions are valid for any given uniform temperature change and for any given nonlinear temperature distribution through the thickness of the multilayered beams and plates. Several numerical examples are presented that illustrate the application of the method for various temperature distributions. A simple experiment was conducted that showed the validity of the analytical method. A brass strip was added to a bimetalic strip made of aluminum and steel at room temperature. The thickness of the brass strip was calculated from the theory to prevent bending. The trimetal strip was placed in a furnace and, as expected, it remained straight for varying temperatures.     

计入横向剪切变形时层合板位移差分解法

在考虑横向剪切变形对层合板弹性解的影响时，本文提出一种数值计算方法。由边界条件给出边界结点位移的表达式，根据薄板的经典理论和一阶横向剪切变形理论导出位移增量所满足的平衡微分方程，引用经典理论计算的横向剪力修正了荷载列阵。致使在较粗的网格划分时、宽广的层合板长厚比范围内，仍能得到与解析解颇为一致的数值解。     

Elastic equilibrium of a half plane containing a finite array of elliptic inclusions

An accurate analytical method has been proposed to solve for stress in a half plane containing a finite array of elliptic inclusions, the last being a model of near-surface zone of the fibrous composite part. The method combines the Muskhelishvili’s method of complex potentials with the Fourier integral transform technique. By accurate satisfaction of all the boundary conditions, a primary boundary-value elastostatics problem for a piece-homogeneous domain has been reduced to an ordinary well-posed set of linear algebraic equations. A properly chosen form of potentials provides a remarkably simple form of equations and thus an efficient computational algorithm. The theory developed is rather general and can be applied to solve a variety of elastostatics problems. Up to several hundred interacting inclusions can be considered in this way in practical simulations which makes the model of composite half plane realistic and flexible enough to account for the microstructure statistics. The stress concentration factors and effective thermoelastic properties of random structure composites with dilute concentration of fibers are estimated in the vicinity of a free edge. The numerical examples are given showing accuracy and numerical efficiency of the developed method and disclosing the way and extent to which the nearby free or loaded boundary influences the local and mean stress concentration in the fibrous composite.     

层状横观各向同性地基上异性矩形薄板的弯曲解析解

选用更具广泛性的层状横观各向同性弹性地基模型,来分析四边自由各向异性矩形地基板的弯曲解析解。先基于直角坐标下横观各向同性体的静力胡海昌通解,借助双重傅里叶变换及矩阵传递法,获得层状横观各向同性地基的静力位移场和应力场;然后将异性薄板的弯曲控制方程,与基于层状横观各向同性弹性地基的位移解建立的板与地基变形协调方程相结合,先按对称性分解,再用三角级数法,得出层状横观各向同性弹性地基上四边自由各向异性矩形薄板的弯曲解析解,包括地基反力、板的挠度及板的内力的解析表达式。克服了数值法的弊端,取消了对地基反力的假设,且避免了矩阵指数函数的计算;同时考虑了地基的层状性及板和地基的各向异性,从而得到板的内力及地基反力更切实际的分布规律。算例结果与文献的有限元结果吻合良好,证明本文方法是切实可行的。     

A layerwise boundary integral equation model for layers and layered media

A hybrid method is presented for the analysis of layers, plates, and multilayered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multilayered system using a total potential energy formulation. The layerwise laminate theory of Reddy is employed to develop a layerwise, two-dimensional, displacement-based, hybrid boundary element model that assumes piecewise continuous distribution of the displacement components through the system's thickness. A one-dimensional finite element model is used for the analysis of the multilayered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of a typical infinite layer (element) assuming linear displacement distribution through its thickness. This fundamental solution is given in a closed form in the cartesian space, and it can be applied in the two-dimensional boundary integral equation model to analyze layered structures with finite dimensions. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems.     

多层复合壳体三维振动分析的谱——微分求积混合法

对于较厚的多层复合壳体,其振动位移沿厚度方向呈锯齿形变化且层间剪切和拉、压应力呈三维耦合状态,采用传统的等效单层理论分析已不能满足精度要求. 建立不受结构厚度、铺层材料性质和铺层方式限制的三维分析方法具有重要的研究价值. 本文以独立铺层为建模对象,结合广义谱方法与微分求积技术建立了一种适用一般边界条件和铺层方式的多层复合壳体三维分析新方法——谱--微分求积混合法. 该方法应用三维弹性理论对独立铺层进行精确建模,有效克服了二维简化理论对横向变形以及层间应力估计不确切的缺点;引入微分求积技术对铺层进行数值离散,将三维偏微分问题转化为二维偏微分问题,降低了求解维度和难度;应用广义谱方法近似地表述离散计算面上的场变量,将获取的二维偏微分方程转化为以场变量谱展开系数为未知量的线性代数方程组,避免了对超越方程的求解. 数值验证结果表明该方法收敛性好,计算精度高.      

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

改进巴西试验:从平台巴西圆盘到切口巴西圆盘

巴西试验是测试混凝土、陶瓷、岩石等脆性或准脆性材料拉伸强度的标准方法,在材料科学和土木工程中有广泛的应用,研究该方法的改进在科学和工程界都受到日益增长的关注.我们对巴西试验的原始试样即巴西圆盘(Brazilian disc,BD)提出一种新的构型:切口巴西圆盘(grooved Brazilian disc,GBD),利用GBD加载直径两端的窄而浅的切口,消除完整巴西圆盘对径压缩加载时可能非中心起裂的缺点,这一点与我们过去提出的平台巴西圆盘(flattened Brazilian disc,FBD)相同.GBD的构型更普适,因为它涵盖了巴西圆盘和平台巴西圆盘.对于求解GBD切口顶点及其正前方近邻的应力,提出两个近似解析模型,分别求出对应的压应力和拉应力,用叠加原理最后推出的公式能够定性地预测应力的变化趋势,以及试样几何参数的影响.近似解析公式,有限元数值计算,对比实验的结果都证实,对径压缩圆盘的几何与表面边界条件的一个微小的改变,对巴西试验试样的几何稳定性和断裂过程产生显著和有利的影响.利用推出的切口应力公式进行试样几何参数的优化选择,可以使GBD比BD和FBD更具优越性.     

含孔曲板弹性波散射与动应力分析

基于敞口浅柱壳弹性波动方程及摄动方法，对无限大含孔曲板弹性波散射及动应力问题进行了分析研究，将经典薄板弯曲波动问题的分析解作为本问题的主项，给出了在稳态波下孔洞附近散射波的零阶渐近解。建立了求解含孔曲板弹性波散射与动应力问题的边界积分方程法，利用积分方程法可获得问题的近似分析解。并给出了无限大曲板圆孔附近动应力集中系数的数值结果，且对计算结果进行了分析与讨论。     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.     

Three-dimensional closed-form analysis of the stress field at rectangular corners of layered plates

Summary A refined and expanded closed-form analytic tool for the calculation of the stress field at free corners of thermally loaded cross-ply and angle-ply laminate structures is developed. Based on adequate approximations for the inplane stresses, the integrated equilibrium conditions yield a full-scale 3D stress field in the vicinity of a free corner. After adjusting the stress field to stress-free conditions at the free edges and facings and obeying conditions of continuity at the interfaces, the remaining free parameters in the stress representation are calculated by minimization of the complementary potential of the laminate. Comparison of the closed-form results with finite element computations shows a good agreement.     

Properties degradation induced by transverse cracks in general symmetric laminates

This paper presents the details of a methodology for predicting the thermoelastic properties degradation in general symmetric laminates with uniform ply cracks in some or all of the 90° layers. First, a stress transfer method is derived by using the concept of state space equation. The laminate can be subjected to any combination of in-plane biaxial and shear loading, and the uniform thermal loading is also taken into account. The method takes into account all independent material constants and guarantees continuous fields of all interlaminar stresses across interfaces between material layers. By this method, a laminate may be composed of an arbitrary number of monoclinic layers and each layer may have different material property and thickness. Second, the concept of the effective thermoelastic properties of a cracked laminate is introduced. Based on the numerical solutions of specially designed loading cases, the effective thermoelastic constants of a cracked laminate can be obtained. Finally, the applications of the methodology are shown by numerical examples and compared with numerical results from other models and experiment data in the literature. It is found that the theory provides good predictions of the thermoelastic properties degradation in general symmetric laminates.     

圆外Stokes问题的边界积分公式

针对圆外区域的Stokes问题,利用Fourier展开法,通过自然边界归化得到了一个只与边界速度有关的Stokes问题的边界积分公式.根据此公式及边界速度值,求得区域内速度及压力分布的解析表达式,并通过数值积分的方法进行求解,计算量小,所得速度及压力分布图曲线光滑.最后借助流体软件进行数值计算,结果验证了边界积分公式的正确性、可行性.