散斑光弹性干涉法

全息干涉法和普通光弹性法相结合产生了全息光弹性法。与此类似,本文建议将散斑干涉法和光弹性相结合形成具有一定特点的散斑光弹干涉法。近十几年来,散斑干涉法得到了很快的发展,可以获得面内或离面位移,错位散斑干涉法则可获得位移梯度场。本文建议将光弹性透明模型置于双光束散斑干涉光路中,通过两次曝光法可以获得模型的应力条纹图案。但和全息光弹性法相比,只要选择合适的漫反射器,用散斑光弹干涉法可以获得不受等差线干扰的、单纯的等和线条纹图。同样的理由,也可以获得三维应力模型冻结切片的等和线条纹图。更有意义的是采用双光片散光穿透模型的散斑光弹性干涉法将有可能获得整体三维应力模型内部的绝对光程应力条纹图,这将为三维应力分析提供一种新的实验方法。     

基于自然单元法的功能梯度板固有频率优化

对功能梯度板的固有频率优化问题进行研究.从无网格方法自然单元法出发,建立了功能梯度板的一阶剪切变形理论的自由振动分析格式.自然单元法是一种基于自然邻近插值的无网格数值方法.自然邻近插值不涉及矩阵求逆运算,也不需要任何人为的参数.较之于传统复合材料,功能梯度板材料属性沿厚度方向呈梯度连续变化,因此材料加工时产生的残余应力小很多,此外,功能梯度板具有可设计的优点.以第一阶固有频率为优化目标建立优化模型,通过单纯形法搜索较优的功能梯度板的成分分布.最后通过矩形功能梯度板在不同边界条件下的基频优化的实例,验证了整个算法的可行性和有效性.     

板壳畸变单元的无网格和有限元自动耦合方法

板壳大变形时单元的严重畸变会使计算精度降低。无网格局部Petrov-Galerkin法是一种真正的无网格方法,能够消除网格畸变,但比有限元法计算效率低。根据板壳网格畸变的局部性特点,利用过渡单元法,基于板壳网格质量,建立了板壳的网格严重畸变区域由有限元分析切换为无网格分析的自动耦合算法,实现了有限元法和无网格局部彼得罗夫．迦辽金法的耦合。应用实例表明：通过自适应耦合,既能发挥有限元法计算效率高的特点,又能发挥无网格法适合大变形分析、没有网格畸变造成计算困难的特点。     

基于弹性理论的残余应力反演和变形计算

毛坯在制造过程中,材料力学性能的非均匀性导致其内部产生残余应力.残余应力会造成结构破坏,在工件的切削去除过程中,残余应力会逐渐释放并引起变形.采用有限元的“生死单元”技术模拟材料的切削去除过程,并转化为残余应力的释放,基于板壳理论、小变形理论、弹性理论和叠加原理,将径向基函数插值法与几何方程、物理方程相结合,开发出一种新的解析方法,反演残余应力场并计算变形.结果表明解析方法的理论解与有限元解高度吻合,能够用于判断残余应力分布并计算残余应力弹性释放引起的变形.     

12.7 mm弹侵彻不同强度钢靶的数值模拟

针对12.7 mm弹侵彻不同强度钢靶时可能出现子弹保持完整或发生破碎的情况,过去的数值模拟仅限于模拟单一模式的子弹侵彻行为。为了克服这种数值模拟的局限性,开展了模型算法、网格尺寸对模拟结果影响的研究,并将模拟结果与实验结果进行了对比,提出了一种能够用于模拟子弹保持完整或破碎的弹靶模型。研究结果表明,为模拟子弹保持完整状态,子弹和靶板应分别采用基于Lagrange算法的有限元法和光滑粒子算法,而且子弹网格尺寸和靶板粒子间距之比应至少保持在5.3左右,否则弹头会产生与实验结果不符合的异常变形。但是,在模拟子弹发生破碎侵蚀时,该比例的网格/粒子尺寸比会引起计算中止。为了克服该问题,进一步建立了一种弹体表面采用大尺寸网格、内部采用细化小尺寸网格的有限元/光滑粒子法耦合弹靶模型。计算结果表明,改进的弹靶模型可模拟子弹保持完整或者发生破碎的情况。     

断裂测试系列讲座 第二讲 云纹法及其在断裂力学中的应用

1.云纹法(Moire Method) 由粘牢于试件(模型或构件)表面并可随其一起变形的栅片和栅板重迭(图1),栅线间因几何干涉而产生条纹(图2),这种干涉条纹称为云纹。云纹法就是根据这类云纹测定和分析试件的位移场或应变场的实验应力分析方法。     

矩形加肋板线性弯曲分析的移动最小二乘无网格法

提出了一种求解矩形加肋板线性弯曲问题的移动最小二乘无网格法。将矩形加肋板模拟成平板与肋条组成的复合结构。先基于一阶剪切变形理论,由移动最小二乘近似建立板和肋条的位移场,再利用板与肋条叠合处的位移协调条件,推导肋条与板的节点参数转换方程,最后利用此方程将板与肋条的应变能叠加,推导出整个加肋板的刚度方程。由于本文提出的加肋板无网格模型中不涉及到网格,肋条不必像有限元那样必须沿网格线布置,肋条位置的改变也不会导致板网格的重新划分。文末算例表明由本文方法得到的解与采用实体单元得到的ANSYS有限元解吻合良好,证明了本文方法的准确性。     

高分子材料平面应变拉伸变形局部化

将基于应变软化玻璃状高分子材料微观特征建立的ＢＰＡ８－链分子网络模型引入ＵｐｄａｔｉｎｇＬａｇｒａｎｇｅ有限元方法，建立了适于变形局部化分析的大变形弹塑性有限元驱动应力法．在此基础上，数值模拟了初始各向同性高分子材料平面应变拉伸变形局部化的传播过程．探讨了ＢＰＡ模型对具有加工硬化特性的结晶性高分子材料变形分析的适应性；分析了局部化传播过程中颈缩截面的非均匀应力三轴效应；最后，讨论了网格尺寸以及初始几何不均匀性对颈缩扩散以及应力三轴效应的影响     

有限元分析中的信息处理——单元延拓里兹法

有限元分析给出连续场在若干离散点(网格结点)的信息(例如平面应力板的位移)。如何利用这些离散点信息获得全场信息的连续分布(变形曲面)和进一步推算其它性态信息(应力),一般采用分片插值和分片拟合方法。本文提出了单元延拓里兹法,并将其应用于有限元分析中的信息处理。因利用了单元邻近区域的信息来处理单元内部的信息,从而提高了精度;利用较粗的网格作有限元分析,又节约了计算时间。     

基于径向基点插值法的旋转Mindlin板高次刚柔耦合动力学模型

将无网格径向基点插值法(radial point interpolation method,RPIM)用于中心刚体?旋转柔性板的动力学分析.基于浮动坐标系方法和一阶剪切变形理论即Mindlin板理论,考虑剪切变形的影响,并计入板面内变形的非线性耦合变形项,采用径向基点插值法描述板的变形场,保留动能中有关非线性耦合变形项...     

光弹性应力分析的条纹灰度法

条纹灰度法在光弹性中能自动采集数据、分析和显示等色线条纹图。这种方法组合了光弹性和数字图象分析两种技术,这篇论文描述了一个Apple—Ⅱ微型计算机图象分析系统,它把等色线图分成256×256个图象元素(象素),每一象素将亮度转换成8—bit分辨率的视频讯号,在黑到亮之间产生256级灰度,灰度Z与条纹级数N在暗场时,它能写成N=(1/π)arcsin(AZ~(1/2r))其中A、r是常数。这种方法对相对迟后小于半波长的双折射特别有效,这篇论文中用两片玻璃模型研究弹性应力场,具有很好的分辨率,其结果表明对应力分析主应力差(σ_1-σ_2)有较高精度。     

聚碳酸酯弹丸穿甲实验与数值模拟

利用空气炮为加载设备,开展了不同头部形状聚碳酸酯弹丸的穿甲实验,获得了弹体穿甲时靶板 的变形时间历程曲线及弹靶相互作用时间,观察到不同弹体在穿甲过程中产生的断裂、内部损伤及头部塑性 变形等失效特征。采用光塑性方法,分析了截锥型聚碳酸酯弹丸穿甲后的变形特征,通过用户子程序将DSGZ 模型引入商用有限元程序,对截锥型聚碳酸酯弹丸的穿甲过程进行了数值模拟,获得了不同时刻聚碳酸 酯弹丸应力分布及靶板的变形特征,模拟结果与实验结果吻合较好。     

关于三维模型光测“冻结”法的若干问题

一、方法的简介 1.光测法 当偏振光通过受载荷作用的透明塑料模型(光测模型)时,模型各点将产生与应力状态相关联的光学干涉条纹图案——应力光图(等色线及等倾线)。用偏光弹性仪,可以定量地量测应力光图,进而对模型作全面的应力分析。若模型及作用于其上的载荷与真实结构(原型)相似,则根据相似原理,可将模型应力转换为原型应力。上述     

高阶数值流形方法的初应力公式

高阶数值流形方法和高阶DDA方法可以显著提高结构变形的计算精度,但目前涉及几何非线性问题的研究成果大都计算精度差甚至不收敛,这是由高阶初应力公式的不准确或不正确引起的。本文介绍数值流形方法的大变形计算格式,基于平面三角形数学网格和多项式覆盖函数,提出高阶流形法的两种初应力处理方法,首次导出了高阶初应力的准确公式。该公式在分步计算的初应力累加中考虑了大变形结构的构形变化,并将初应力表示成多项式函数形式以满足单纯形积分的要求。文中给出的悬臂梁大变形数值算例与理论解的对比结果证明了方法的正确性。本文的方法和公式也适用于三维四面体数学网格,稍加修改后将可应用于高阶DDA方法和常规的有限元方法。     

界面上应力波传播规律的动态光弹性试验研究

采用动态光弹性方法,对应力波在杆-杆和板-板双材料连接界面上的反射和透射传播规律展开研究．通过对双材料粘接模型施加冲击载荷作用,采集和分析等差线条纹,计算出最大剪应力值．对于两种试验模型应力波在界面上都以透射为主,对应力波的阻隔作用不明显,这与接触式界面对应力波有较强的阻隔作用有较大区别,在应用界面减弱应力波的作用时应注意区分．试验结果与理论计算结果基本保持一致,并由此证明采用动态光弹性法研究应力波在界面上传播规律的可行性．     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为。计算采用二维轴对称应变条件下的绝热模型。钨合金的本构方程采用热粘塑性形式的Johnson Cook模型。为了得到不同尺度的变形信息 ,计算中用了两种网格 ;先用粗糙网格分析试件变形局部化的概貌 ;接着 ,用细密网格 (在变形局部化区域 ,网格尺寸达到 10 m)分析绝热剪切带的形成和发展。有限元模拟得到的绝热剪切带位置和方向与实验一致。计算结果表明 ,绝热剪切带的形成和发展与试件的应力状态密切相关。     

冲击载荷下钨合金圆台试件绝热剪切变形局部化的数值模拟

采用有限元计算编码ABAQUS模拟了钨合金圆台试件在冲击载荷下的变形和剪切局部化行为.计算采用二维轴对称应变条件下的绝热模型.钨合金的本构方程采用热粘塑性形式的Johnson-Cook模型.为了得到不同尺度的变形信息,计算中用了两种网格;先用粗糙网格分析试件变形局部化的概貌;接着,用细密网格(在变形局部化区域,网格尺寸达到10μm)分析绝热剪切带的形成和发展.有限元模拟得到的绝热剪切带位置和方向与实验一致.计算结果表明,绝热剪切带的形成和发展与试件的应力状态密切相关.     

海洋平台K型管节点的疲劳裂纹扩展分析 II:数值分析

对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定。本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度。这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型。同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性。     

无网格局部Petrov-Galerkin法求解板壳弹塑性大变形

无网格局部Petrov-Galerkin法构造的高阶光滑的形函数非常适合建立板壳结构场函数的逼近函数,是一种比较理想的研究板壳问题的方法。基于Mindlin板壳理论,采用更新拉格朗日原理和大变形条件下场量的无网格表达形式,实现了率型无网格局部Petrov-Galerkin方法对板壳弹塑性大变形的求解,算例分析表明了方法的有效性和较高的分析精度。     

海洋平台K型管节点的疲劳裂纹扩展分析Ⅱ:数值分析

对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定.本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度.这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型.同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性.