ANSYS的仿真及教学日记-第4-8章-Current.docx

ANSYS的仿真及教学日记  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc350021389 4. 动力学有限元分析  PAGEREF _Toc350021389 \h 2  HYPERLINK \l _Toc350021390 4.1. 动力学分析与矩阵特征值  PAGEREF _Toc350021390 \h 2  HYPERLINK \l _Toc350021391 4.2. 播放音频文件  PAGEREF _Toc350021391 \h 5  HYPERLINK \l _Toc350021392 4.3. 模态分析  PAGEREF _Toc350021392 \h 6  HYPERLINK \l _Toc350021393 4.4. 几个教学案例的分析要点：  PAGEREF _Toc350021393 \h 8  HYPERLINK \l _Toc350021394 5. 非线性有限元分析  PAGEREF _Toc350021394 \h 12  HYPERLINK \l _Toc350021395 5.1. 非线性分析中的纯粹增量式近似与牛顿-拉普森近似  PAGEREF _Toc350021395 \h 12  HYPERLINK \l _Toc350021396 5.2. ANSYS仿真中应力刚化(stress stiffening)  PAGEREF _Toc350021396 \h 12  HYPERLINK \l _Toc350021397 5.3. ANSYS对单元几何形状特征的限制  PAGEREF _Toc350021397 \h 18  HYPERLINK \l _Toc350021398 5.4. 习题2(2)圆柱壳体的非线性屈曲分析  PAGEREF _Toc350021398 \h 19  HYPERLINK \l _Toc350021399 6. 热场及热力耦合的有限元分析  PAGEREF _Toc350021399 \h 25  HYPERLINK \l _Toc350021400 6.1. 热场与温度场耦合的计算  PAGEREF _Toc350021400 \h 25  HYPERLINK \l _Toc350021401 7. 其他有限元分析  PAGEREF _Toc350021401 \h 30  HYPERLINK \l _Toc350021402 7.1. 可靠性分析  PAGEREF _Toc350021402 \h 30  HYPERLINK \l _Toc350021403 7.2. 拓扑优化设计实例  PAGEREF _Toc350021403 \h 31  HYPERLINK \l _Toc350021404 7.3. 电磁场分析  PAGEREF _Toc350021404 \h 34  HYPERLINK \l _Toc350021405 8. ANSYS的一些使用技巧  PAGEREF _Toc350021405 \h 35   动力学有限元分析 动力学分析与矩阵特征值 也就是对于某一个向量，它通过??阵变换得到的值相当于其通过代数变换得到的，这个向量就是这个矩阵的特征向量，λ就是这个矩阵的特征值。 如果矩阵是某个向量空间到自身的变换，Ax=λx，就是向量X经过矩阵A的变换后相等于向量X的一个代数变换，即方向不变，只是长度变化而已。特征值λ一般对应于关键值，比如振动频率、稳定性参数、原子的量级等等。 求解的依据是，如果方程组有非平凡解，则矩阵的行列式必须等于0。将行列式展开，得到一个关于λ的多项式。 比如两个质量块m1,m2以弹簧连接于两边墙上，对于质量块m1, m2使用牛顿第二定律写出每个物体的受力平衡关系式为： 方程可以表示为： 由振动定理知道，公式的解是： ，Ai为物体i的振幅，ω是振动频率，这样 带入得到 如果m1=m2=40kg, k=200N/m,计算特征值和特征向量 该矩阵的行列式的值为： 求得=15或者5，物体的振动角频率为3.873 s-1，物体的振动角频率为2.236s-1。 我们尝试着自己去对一下： 观测得到： 采用Matlab求解， 得到特征值和特征向量，和预期的一样。 #向量的特征值和特征向量的计算 A=[10 -5;-5 10]; [V,D]=eig(A) V = -0.7071 -0.7071 -0.7071 0.7071 D = 5 0