一种新的缺陷压电体奇异准协调单元及其断裂分析

在缺陷压电体断裂理论分析的基础上构造了一种新的缺陷压电体广义位移模式并建立了有限元断裂分析模型;提出了求解压电体含大曲率缺口广义强度因子等断裂参量的数值分析计算方法。考查了不同结构尺寸、电位移与电场对解的影响,本方法所得的结果与理论解一致,说明了本文方法的有效性。     

考虑应力梯度的缺口疲劳寿命预测方法

基于试样缺口根部区域应力分布规律,综合考虑了平均应力、应力梯度及尺寸效应的影响,发展了一种考虑梯度影响的缺口疲劳寿命预测方法.以TC4合金材料含两种形式缺口的试样为例,计算得到了缺口应力梯度影响因子和尺寸效应影响因子,并进行了疲劳寿命预测.寿命预测结果表明:引入缺口局部平均应力、应力梯度影响因子和尺寸效应影响因子可以较为全面地考虑缺口对试样疲劳寿命的影响,基于所发展的缺口疲劳寿命预测方法给出的TC4合金缺口试样疲劳寿命预测结果在2倍分散带以内.      

基于三维裂纹尖端应力场的应力强度因子计算方法

提出一种基于无限大体裂纹尖端弹性应力场理论解的前几项多项式函数,对实际裂纹体弹性应力场有限元解进行拟合来计算应力强度因子的方法.该方法在计算应力强度因子时不需要预先假设裂纹尖端的应力应变状态,应力强度因子计算结果更符合三维裂纹体裂纹尖端实际的应力应变状态.首先基于二维无限大板中心穿透裂纹应力场理论解验证了方法的有效性,探讨了拟合确定应力强度因子需要的多项式应力函数的项数.然后分别以二维大板中心穿透裂纹、三维大体内埋圆裂纹和三维有限厚板中心穿透裂纹的应力强度因子计算为例,通过与无限大板和无限大体应力强度因子理论解以及基于位移外推法和1/4节点张开位移法的应力强度因子有限元解的对比分析,验证了该方法的有效性和合理性.研究表明该方法能够合理反映三维裂纹体裂纹尖端的实际应力应变状态,计算得到的应力强度因子数值更合理.      

运用损伤力学求解金属构件疲劳裂纹形成寿命的附加位移法

从损伤力学出发,建立了对金属构件疲劳裂纹形成过程进行损伤耦合的弹塑性应力分析及寿命估算的附加位移法。该方法较充分地利用了缺口构件疲劳损伤的局部性特点,其计算量远低于普通的全场数值解法。算例表明寿命的理论解与试验值基本吻合。该方法对于分析其他局部耦合问题也具有参考价值     

湍流模型对三维翼梢涡流场数值模拟的影响

 梢涡流场数值模拟对减少梢涡的不良影响意义深远。应用雷诺平均Navier-Stokes(RANS)方法建立了三维翼梢涡流场数值模型,选用的湍流模型包括两种经典两方程模型和一种代数雷诺应力模型。计算内容包括翼梢涡流场速度分布、压力分布及涡核位置,为了准确模拟梢涡涡核内系统旋转和流线曲率影响,将旋转和曲率修正方法应用于选用的湍流模型中。计算时对翼梢涡涡核处网格进行了加密处理。将计算结果与实验结果对比后发现,修正后的湍流模型较标准湍流模型明显提高了预报精度,与实验结果吻合较好,具有工程应用价值。     

反映流动曲率影响的非线性湍流模式

本文从一般曲线坐标系下的雷诺应力模式出发，在雷诺应力模式向代数应力模式简化过程中，保留工率项对雷诺应力输运的影响，得到了曲率修正的代数应力模式和非线性κ-ε模式。并且该模式对一典型曲率流动－－U型管道内的湍流流动进行了数值模拟与分析。计算结果显示,率修正的代数应力模式较好地模拟出了凸曲率地雷诺剪切应力的阻滞作用和凹曲率对剪切应力的增强作用，捕捉到了曲率流动的主要特征，是对雷诺应力模式的一个合理而有效的简化。     

考虑应力松弛的缺口疲劳寿命预测方法

为了能够快速准确获得循环载荷条件下缺口局部应力应变并用于缺口疲劳寿命预测,基于粘塑性本构方程及光滑试样循环应力应变试验结果,发展了一种应力应变简化计算方法,采用该方法计算得到的光滑试样循环应力应变结果与试验结果的误差在5%以内,并将该方法结合寿命预测方程采用迭代的方式用于粉末高温合金FGH97单边圆形缺口试样疲劳寿命预测。结果表明:所发展的方法充分考虑了FGH97合金单边圆形缺口试样缺口根部区域的应力松弛,取得了较好的寿命预测结果,其分散带基本在2倍以内。     

海绵/橡胶适配器应力和变形模式的轴对称平面应变解析

海绵/橡胶适配器的应力和变形模式的解析方法是其力学理论研究的重要前提。基于轴对称平面应变假设,推导了海绵圆筒和橡胶圆筒的径向应力、切向应力和径向位移的解析公式;建立了位移边界条件下适配器受压问题的非线性常微分解析方程组,并通过数值模型仿真验证了解析公式推导的正确性;调整海绵/橡胶粘合面半径,研究了海绵层厚度比对适配器受压力学特性的影响。结果表明:数值模型和解析公式计算结果基本一致;适配器海绵层径向应力值和切向应力值相差不大,橡胶层切向应力高于径向应力约1个数量级,径向应力最大值位于适配器外表面,切向应力和径向位移最大值均位于海绵/橡胶粘合面;随海绵层厚度比增加,适配器应力和位移减小且减幅越来越小。     

基于转角模态曲率的损伤识别研究

利用测量得到的模态参数,应用中心差分法求取振型曲率,用损伤前后同一位置的振型曲率变化作为结构损伤的识别指标。采用有限元方法,通过数值仿真,对不同位置发生损伤前后的悬臂梁进行仿真计算,得到其模态参数。利用构造的损伤识别指标对结构损伤情况进行识别研究。在前人研究的基础上,特别提出了利用角位移模态曲率变化来识别损伤位置。仿真算例表明：利用角位移模态曲率比平动位移模态曲率识别损伤的效果更好,精度更高。     

复合材料叠层结构剖面翘曲修正理论

 首先假设复合材料叠层结构的面内位移沿厚度方向为直线分布,然后根据剪切变形进行连续修正,获得逐步升高的各阶翘曲位移分布函数。文中以叠层梁为例,应用最小势能原理导出其静力方程和力的边界条件。在长度方向则采用升阶谱位移函数,进行了数值计算。所得叠层梁的位移和应力分布的数值结果与精确解相比,表明本文方法有很好的收敛性。     

三维钝裂纹的边界元模拟和剩余腐蚀疲劳寿命估算方法

提出了一个三维钝裂纹的计算模型,并建立了适用于这类模型的边界元法程序。还利用现有的疲劳裂纹扩展模型,对含三维钝裂纹构件的剩余疲劳寿命估算方法进行了讨论,并对三维钝裂纹在空气和海水中的扩展进行了计算。研究表明,本方法和边界元程序,既适用于三维尖裂纹的剩余寿命估算,也适用于三维钝裂纹的剩余疲劳寿命估算。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

涡轮盘销钉孔损伤容限分析新方法及其应用

本文提出温度及离心载荷作用下三维构件损伤容限分析方法。该方法采用新型双重边界元法分析温度及离心力载荷作用下的涡轮盘三维裂纹应力强度因子 (SIF) ,并结合 Paris公式及 Euler法获得涡轮盘销钉孔边裂纹扩展形状及裂纹扩展寿命 ,进而采用二次估计方法获得涡轮盘裂纹扩展寿命修正值。利用新型双重边界元法分析了销钉载荷、温度场、裂纹形状对轮盘孔边三维裂纹应力强度因子的影响 ,并对以往涡轮盘销钉孔边裂纹扩展寿命分析工程方法的近似性进行了讨论。文中算例表明采用本文方法分析复杂载荷三维裂纹扩展寿命具有效率高、建模方便的优点。     

金属材料延性断裂韧度测试中的柔度方法研究

 对紧凑拉伸(CT)试样的柔度转动修正方法进行了研究,理清了现行柔度转动修正方法存在的问题,给出了新的转动修正方法,并提出了考虑转动效应的CT试样的裂纹嘴张开位移(CMOD,V₀)与加载线张开位移(V_LL)的转换公式。弹塑性有限元分析表明,CT试样产生刚性转动的中心并不在试样的剩余韧带中心,而是在偏靠裂纹尖端的位置;转动半径 R 基本不受材料本构关系的影响,仅与CT试样的裂纹长度a与宽度W之比(a/W)有关,从而提出了R与a/W之间的单调多项式;当裂纹尖端产生较大程度的塑性变形时,需要考虑转动效应对CT试样裂纹张开位移(COD)测量的影响, J积分的塑性功计算须采用经转动修正后的COD。采用两种延性材料Cr2Ni2MoV和16MnR对COD转换公式进行了实验验证,结果表明,基于转动分析的COD转换公式由于考虑了裂纹尖端附近区域的弹塑性变形行为,用于断裂韧性测试中的塑性功计算更加符合实际,且同实验结果符合良好;直通型紧凑拉伸(FFCT)试样的刚性转动对较大裂纹扩展情形下的J_R阻力曲线影响较大。     

机械连接件中螺栓孔角裂纹应力强度因子数值计算(英文)

The three-dimensional finite element method is used to solve the problem of the quarter-elliptical comer crack of the bolt-hole in mechanical joints being subjected to remote tension. The square-root stress singularity around the corner crack front is simulated using the collapsed 20-node quarter point singular elements. The contact interaction between the bolt and the hole boundary is considered in the finite element analysis. The stress intensity factors （SIFs） along the crack front are evaluated by using the displacement correlation technique. The effects of the amount of clearance between the hole and the bolt on the SIFs are investigated. The numerical results indicate that the SIF for mode I decrease with the decreases in clearance, and in the cases of clearance being present, the corner crack is in a mix-mode, even if mode I loading is dominant.     

幂硬化材料平面应力裂纹稳定扩展的有限元模拟

对幂硬化材料中平面应力Ⅰ型裂纹的稳定扩展进行了有限元模拟。计算中采用了逐步释放裂尖节点力法和子结构方法并配合裂尖附近高密度的网格划分,获得了有关裂尖应力变形场、塑性区尺寸、裂纹开口位移、裂纹扩展阻力等大量资料。计算结果还证实了部分已有的裂尖场渐近解。     

Research on the Plane Multiple Cracks Stress Intensity Factors Based on Stochastic Finite Element Method

Taylor stochastic finite element method (SFEM) is applied to analyze the uncertainty of plane multiple cracks stress intensity factors (SIFs) considering the uncertainties of material properties, crack length, and load. The stochastic finite element model of plane multiple cracks are presented. In this model, crack tips are meshed with six-node triangular quarter-point elements; and other area is meshed with six-node triangular elements. The partial derivatives of displacement and stiffness matrix with respect to all the random variables obtained by Taylor SFEM are derived. Meanwhile, the mean and variance expressions of SIF under uncertain factors are also derived. Parallel double-crack illustrative example of using the proposed method is given. The calculation results indicate that the uncertainty of SIF is influenced by the distance of the cracks, the smaller the distance of cracks is, the greater the SIF uncertainty is, and the uncertainties of material elastic modulus, load and crack length markedly affect the uncertainty of SIF; and the Poisson ratio of material has little influence. Among the variables, the elastic modulus has the greatest effect on SIF uncertainty. The next is external load. The crack length has lower effect on SIF uncertainty than both the elastic modulus and external load.     

航空管路补偿器耐久性计算方法对比分析

本文研究了管路补偿器补偿位移疲劳寿命和耐久振动寿命的有限元分析计算方法。通过对某型管路补偿器有限元仿真计算结果和传统工程计算法计算的结果进行对比分析,研究管路补偿器补偿位移疲劳寿命的有限元计算方法,探索耐久振动寿命的有限元计算流程,并将结果同试验数据对比,以验证其设计方法的计算精度。结果表明,补偿位移的疲劳寿命有限元法的计算精度优于工程设计法;耐久振动寿命有限元法能够给出确定的振动应力分布情况,并预测出耐久振动寿命时间,工程设计法仅能计算出管路补偿器的自振频率,不能明确振动应力和寿命。