abaqus稳态动力学多载荷施加方法.pdf

1 Abaqus 稳态动力学分析中多个载荷的定义 张智渊 摘要：本文通过一个例子研究使用 Abaqus 进行稳态动力学分析时，考虑不同载荷相位的定 义方式。通过 Matlab 程序对于 Abaqus 计算结果进行验证,得出了 Abaqus 稳态动力学计算中 多个相位不同载荷施加方式和结果获取的途径。 1．模型建立和模态分析 如图 1 所示的 5 自由度质点弹簧体系。五个质点，每个质点只有 1 个水 平自由度。各质点质量均为 1，质点之间用无质量的弹簧单元连接，弹簧刚度为 100。弹簧和质点均无阻尼。容易列出该体系的质量矩阵与刚度矩阵，求解特征 值可得到该体系的各阶振型与自振频率，如表 1 和图 2 所示。 在 Abaqus 中通过弹簧单元和质点建立计算模型(*.inp 文件见附件 1)。得到 模型自振频率见表一。 表 1: 体系的自振频率（角频率） 模态编号 1 2 3 4 5 自振频率（MATLAB） 2.8463 8.3038 13.0972 16.8251 19.1899 自振频率（ABAQUS） 2.8463 8.3083 13.0973 16.8251 19.1901 模型受两个载荷 F1，F2，分别作用在 1，2 两点处。其中：F1=50cos10t， F2=40cos10t+30sin10t。可见 F2 与 F1 相比幅值相同，相位相差-36.8698 度。 2 2．公式推导分析 该系统的振动方程为： ? ?Mx Kx F?? 其中， 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? M ， 200 100 0 0 0 100 200 100 0 0 0 100 200 100 0 0 0 100 200 100 0 0 0 100 100 ?? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? ? ??? ? K 激励力 50cos10 40cos10 30sin10 0 0 0 t t i t ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? F ，可以表示成另一种形： 50 40 30 0 0 0 i ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? F 。 对于第一种激励力可以把激励力分解为正弦部分和余弦部分，并分别进行求 解。对于正弦部分设 1 sin10t?X X ，对于余弦部分设 2 cos10t?X X ，分别求解然后 线性叠加得到最后结果为： 0.4 0.3 0.5cos(10 -143.1) 0.9 0.3 0.9487cos(10t-161.6) cos10 sin100.9 0.3 0.9487cos(10t-161.6) 0 0 0 0.9 0.3 0.9487cos(10t+18.43) t t t ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? X 对于第二种激励，假设 F 为三角函数激励，角频率为 10，假设响应的也为 三角函数形式，角频率为 10，这样带入振动方程中求解可得 0.4 0.3 0.9 0.3 0.9 0.3 0 0.9 0.3 i i i i ? ?? ? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ? ??? ? X 当例子不考虑相位差进行幅值计算时， 激励近似为为 50cos10 50cos10 0 0 0 t t ? ? ? ? ? ? ? ?? ? ? ? ? ? ?? ? F ，响应为 0.5 1.0 cos101.0 0 1.0 t ?? ? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? X 。 3 该例题中相位差较小，但是结果相比较还是有一定的差距。当激励相位差较 大时，不考虑相位差的计算结果将会产生很大的误差。因此很有必要寻找在 Abaqus 中正确添加存在相位差的多载荷的方法。 3．Abaqus 分析及与公式推导结果对比 将三角函数激励转换为复数进行计算有其方便之处，在此并没有完全体现出 来，当遇到积分或导数等情况时，复数形式则带来很大方便。Abaqus 中也是采 用复数形式进行计算，输出为幅值和相角。 在 Abaqus 中采用复数形式计算，input 文件载荷语句如下。Cload 为集中力， load case=1 说明为实部，load case=2 说明为虚部，_PickedSet29