弹性力学仿真软件：MSC Nastran：谱分析与随机振动技术教程.pdfVIP

弹性力学仿真软件：MSCNastran：谱分析与随机振动技术

教程

1弹性力学基础

1.11弹性力学基本概念

弹性力学是研究弹性体在外力作用下变形和应力分布的学科。弹性体是指

在外力作用下能够产生变形，当外力去除后，能够恢复原状的物体。在弹性力

学中，我们关注的是物体的内部应力、应变以及位移，这些量可以通过求解弹

性方程来获得。

1.1.1弹性体的分类

一维弹性体：如杆、轴等，主要考虑轴向应力和应变。

二维弹性体：如板、壳等，需要考虑平面内的应力和应变。

三维弹性体：如块体、复杂结构等，需要全面考虑三个方向的应

力和应变。

1.1.2弹性常数

弹性模量：描述材料抵抗弹性变形的能力，分为杨氏模量（E）和

剪切模量（G）。

泊松比：描述材料横向变形与纵向变形的比值，通常用符号ν表

示。

1.22应力与应变分析

1.2.1应力

应力是单位面积上的内力，分为正应力（σ）和剪应力（τ）。在三维空间

中，应力可以用一个3x3的对称矩阵表示，称为应力张量。

1.2.2应变

应变是物体变形的程度，分为线应变（ε）和剪应变（γ）。同样，应变也

可以用一个3x3的对称矩阵表示，称为应变张量。

1

1.2.3应力应变关系

在弹性范围内，应力与应变之间存在线性关系，由胡克定律描述：

=

对于三维情况，胡克定律可以表示为：

=

其中，是应力张量，是应变张量，是弹性常数矩阵。

1.2.4应力应变分析示例

××

假设有一个长方体，其尺寸为111，材料的杨氏模量=

200泊松比=0.3。当在长方体的顶部施加一个均匀的压力=100，

我们可以计算出长方体的应变和位移。

#Python示例代码

importnumpyasnp

#材料属性

E=200e9#杨氏模量，单位：Pa

nu=0.3#泊松比

#应力

P=100e6#压力，单位：Pa

#计算应变

epsilon=P/E

#计算位移

delta=epsilon*1#假设长方体的长度为1m

#输出结果

print(f应变：{epsilon:.6f})

print(f位移：{delta:.6f}m)

1.33弹性方程与边界条件

1.3.1弹性方程

弹性方程是描述弹性体内部应力和应变分布的微分方程，通常包括平衡方

程和本构方程。平衡方程描述了力的平衡条件，而本构方程则描述了材料的物

理性质，即应力与应变之间的关系。

2

1.3.2边界条件

边界条件是弹性方程在边界上的约束，分为位移边界条件和应力边界条件。

位移边界条件规定了边界上的位移，而应力边界条件则规定了边界上的应力。

1.3.3弹性方程与边界条件示例

考虑一个简单的梁，两端固定，中间受到一个集中力的作用。我们可以使

用弹性方程和边界条件来求解梁的变形。

#Python示例代码，使用有限元方法求解梁的变形

importnumpyasnp

fromscipy.sparseimportlil_matrix

fromscipy.sparse.linalgimportspsolve

#定义梁的长度、宽度、高度和材料属性

L=1.0

b=0.1

h=0.1

E=200e9

nu=0.3

#定义网格和节点

n=10

dx=L/n

nodes=np.linspace(0,L,n+1)

#定义单元和刚度矩阵

elements=np.arra