abaqus随机振动米塞斯应力理论.docxVIP

一种在随机振动环境中计算米塞斯均方根 应力 RMS 的高效方法 1.介绍 提出了一种严谨高效的方法来计算稳态随机载荷下线性结构的米塞斯应力均方根值。 均 方根应力 RMS 值通过载荷的零时间间隔的协方差矩阵来表示，在大多数结构分析中，我们 在频域中通过应力分量转换矩阵和所加激励载荷的互功率谱密度矩阵来计算激励载荷的协 方差矩阵。在以前的文献中提出了关键的相关性，但是在之前这种相关性并没有被应用过。 在一般的设计中所用的米塞斯均方根应力的准确确定方法可能是保守或不保守的。 最终， 由 于准确计算米塞斯均方根应力 RMS 的高效率，分析工程师现在可以通过研究米塞斯应力来 对工程设计的可靠性进行深入的探索研究。 我们越来越多的采用数值计算方法，尤其是有限元方法，来预测结构和构建的可靠性。 数值计算工具先进性。 多样性的数值计算能力及误差的降低推动了这种趋势的发展。 在本专 论中研究了经受一类随机激励载荷的结构可靠性的问题， 失效总是与接近屈服强度的米塞斯 应力值有关。 这个问题本质比静态和动态问题的确定要困难的多， 因为输入载荷仅具有统计 学意义， 并且响应本质上也是统计学意义上的。 另外， 由于冯米塞斯应力是应力分量的非线 性函数， 因此用于计算加速度、 位移、 应力分量随机响应的计算方法及理论不能直接用于计 算冯米塞斯应力。 线性系统非确定性响应的深入讨论可以参考文献 【1 】，关于这些问题的更 多讨论可以参考文献【 2】【 3】。 从频率数据中计算冯米塞斯应力的最直接方法要求对线性应力分量的很长的时间序列 进行计算。 可以在每个时间步计算应力不变量并且通过时间积分来确定均方根。 这种计算流 程的代价使得这种方法在冯米塞斯应力的大量研究中显得不切实际。 涉及近似值计算的简化 计算方法，例如 Miles 关系【 4】，可能不够保守【 5】。 最终的方法就是找到冯米塞斯应力在施加的激励载荷的统计学性质上的概率分布。 本专 论将指出， 怎样从输入的激励载荷的统计转换到冯米塞斯应力的概率分布。 特别将给出频域 内直接计算的冯米塞斯应力均方根值的显式表达式。 这种新方法可以使分析工程师对冯米塞 斯应力的研究常规化，通过考虑结构全频域的响应实现对工程设计可靠性的深入研究。 2 问题 一个典型的随机振动试验， 对一个结构施加一个单独的载荷源， 例如一个振动器， 受到 了一个用输入的加速度功率谱密度来表征的振动载荷。 将加速度计和应变计布置在测试试件 上可以用来对有限元模态分析的准确度或精度进行校核。 结构分析工程师所面临的问题就是 去评估在结构上产生的应力场是否达到失效。 对于延性材料， 冯米塞斯应力提供了这种度量， 分析工程师需要一种工具将输入的激励载荷的统计学内容转换正所导致的冯米塞斯应力的 某些统计学参量。 这个模型在为了说明这个问题，我们采用壳单元创建了一个铝合金圆筒的有限元模型, 底部受到横向的随机振动。 图1和图2分别显示了圆筒模型和施加在其底部的加速度功率谱 密度（PSD）。 这个模型在 当前的标准分析流程就是： 选择一个单一的模态频率 （这个模态频率是典型的一个， 它 与输入带宽内的最高阶模态等效质量相联系【 6】并且按照Miles相关性在这个频率上用模 态阻尼和加速度功率谱去计算一个 等效静态重力场”。执行一个静态应力分析，就好像是结 构受到了一 g水平的重力。由其他结构频率所贡献的结构响应被忽略了并且没有用来形成任 何模态形状。如果我们选择的模态在所有模态中起主导作用并且如果这个模态与静态变形相 类似，这种方法就可以得出一个较合理的应力分布估计。 对于确定整体的随机应力响应， 上 述方法通常都是不准确的。 本文中提出了一种方法可以准确的、高效的获得贯穿于整个结构的所有模态和所有关注 频率上的冯米塞斯应力均方根。 3冯米塞斯应力均方值 冯米塞斯应力 p t是被定义成一个标量，因此它的平方式应力发二次函数： TOC \o 1-5 \h \z 2 T , p t y [t A；「t （ 1） 这里坊 ^^xx^yy^zz^ ° yz，① zx ▽xy]是应力张量不重复的分量， 并且, - 1 -1/2 -1/ 2 1 -1/2 1 -1/ 2 A = -1/2 -1/ 2 1 3 3 （2） 3 一 TOC \o 1-5 \h \z 按照模态坐标来表示应力， 二 t 八 qk t ?k「 （3） k 在这里T气丁是所关注位置的k阶模态下的应力分量的列向量， qk t是那一时刻的第k阶广 义坐标（可能是主坐标或正则坐标 ）。公式中求和是对所有模态进行求和。将公式（ 3）带入 公式（1），模态冯米塞斯应力的平方可以用模态坐标的方式表达： 2 t P t q t qj t 讥-A?： （4） i j 按照模态坐标写出的冯米塞斯应力