02-07 边界约束的处理_共6页
- 格式：pdf
- 大小：142.67 KB
- 文档页数：6
§2-7 边界约束的处理
一、边界约束
●由于总体刚度矩阵是一个奇异矩阵，在求得总刚矩阵和总体载荷列阵之后，还不能立即
求解整体节点平衡方程组。
◎从数学上讲，此时的总刚矩阵无逆矩阵，方程组没有确定的解。
◎从其物理意义来说，是由于整个结构未引入边界约束，为一自由结构，对于一个定
常力系的作用，没有定常的位移。
◎因此，为进一步解得结构位移，必须引入足够的几何边界约束，以消除结构的刚体
位移。
对于同一结构，在受相同载荷的条件下，由于不同的边界约束，求得的结构位移、应力
等会大不相同。因此，引入正确的边界条件是获得较高精度解的前提。
●根据结构的实际情况，离散出现的边界约束大致可分为如下三种：
1．基础刚性支承
◎大多数结构要支承在基础上。当基础的刚性很大时，根据不同的支承类型，可以认
为结构和基础相连的节点的一个或几个方向的自由度受到了限制，即位移分量为零。
☆如一简支梁，可以认为其支承点处的一个或二个方向的位移分量为零。
2．对称结构的对称部分支承
◎当结构和外载荷均对称于某些轴线时，为减少工作量或提高计算精度，可只计算结
构的1/2或1/4。此时，为保持原有结构特性，要在对称剖分面的节点上施加垂直
于剖分面的刚性约束，以限制该方向的位移。
☆如轧机机架。
3．允许产生给定位移的支承
◎由于结构本身或安装的需要，在支承和结构之间存在给定的间隙，在结构受到实际
约束之前，此节点处允许产生该距离的位移。
☆如高炉下降管的多余支承。
★ 从数学意义上来讲，上述三种支承(几何约束)可以归纳为零位移约束和给定
位移约束二种，而前者则又是后者的一个特例。
二、边界约束的处理
根据边界约束的类型及后续处理方法和要求的不同，边界约束处理大致采用如下方法：
1.
划行划列法●这种方法适用于预定边界位移为零的约束条件。
●具体做法：在用矩阵表示的线性方程组中，划去相应于己知为零的节点位移分量的行和
列，以消除刚度位移。
◎如图2-13所示的单元组合体，其边界条件为，足0
654421vvvuuu
以消除结构的刚体位移。
◎处理时，则是将以上各为零位移分量相应的行与列划掉，这样，原来12阶的线性
方程组及其12×12阶的总体刚度矩阵，就变成了6阶的线性方程组及其6×6阶的总体刚度矩阵，即
1111
21222
3
313233
3132334
33
5
52535553
636
6365660
0
0
00
000yKvR
KKv
KKKu
KKKKv
KKKKu
KKKKu
对
称
●这样约束处理是必要的。
（1）因为总体刚度矩阵在约束处理前是一个奇异矩阵，而经过约束处理划掉某几行
和几列后变为非奇异矩阵，即约束处理后的总体刚度矩阵的行列式不等于零。
（2）另外，如果不进行约束处理，那么包括在总体节点载荷列矩阵中的约束反力必
须事先求出，作为已知节点载荷。
☆然而，对于形状较复杂一点的单元组合体，在高次超静定情况下，约束反
力很难求出。
☆经过约束处理后，在划去总体节点位移列矩阵与总体刚度矩阵中相应于已
知节点位移分量为零行与列的同时，总体节点载荷列矩阵中未知的约束反
力的行也都被划掉。这样一来，无论次数多高的超静定问题，约束反力都
不必事先求出。
●这种约束处理也是可行的。
（1）因为线性方程组是由各节点平衡方程建立起来的，而方程组的未知量就是节点
位移分量，那么受约束的节点有一个或两个位移分量已知为零，就不必再去求
它，因此该节点的一个或两个平衡方程就可不要，即可以把它们所在的行划去；
（2）同时，在其它方程中，与已知零位移分量和相应的载荷分量，即相应刚度矩阵
元素和此位移的乘积也为零，所以该列也可划去。
◎由此可见，划行划列的约束处理方法是完全可行的，并不影响计算结果。
●划行划列约束处理使总体刚度矩阵发生了两个变化：
（1）总体刚度矩阵的阶数下降。若单元组合体有n个节点和r个约束，则总体刚度
矩阵在约束处理前为2n×2n阶，约束处理后变为(2n-r)(2n-r)
阶。（2）总体刚度矩阵的奇异性发生变化。约束处理前是奇异矩阵；约束处理后变为非
奇异性矩阵。而对总体刚度矩阵的对称性，稀疏性和带形分布等特性并无影响。
由于约束处理时在划去某行的同时划去同序号的列，所以总体刚度矩阵仍保持
其对称性；另外一般单元组合体的r/2n比值是很小的，所以约束处理后总体刚
度矩阵仍保持稀疏性和带形分布的特点。
●经过约束处理后，所建立起来的线性方程组的个数与要求解的未知节点位移分量的个数
都是2n-r个。
●特点：
◎这种处理方法，由于舍弃了相应于已知位移分量为零的行与列各元素，这样就改变
了各方程及元素的编排序号；
◎另外，若是求出各节点位移｛δ｝之后，需计算约束反力，则需重新计算相应行中
各刚度矩阵元素。以上二点是利用此法在编写程序时要注意的。
2. 划0置1法
●适用：这种方法适用于边界节点位移分量为已知(含为0)的各种约束。
●做法：
（1）将总刚矩阵〔K〕中相应于已知位移行主对角线元素置1，其他元素改为零；同
时将载荷列阵｛R｝中相应元素用已知位移置换。
◎这样，由该方程求得的此位移值一定等于已知量。
（2）将〔K〕中已知位移相应的列的非主对角成元素也置0，以保持〔K〕的对称性。
◎当然，在已知位移分量不为零的情况下，这样做就改变了方程左端的数值，
为保证方程成立，须在方程右端减去已知位移对该方程的贡献——已知位移
和相应总刚元素的乘积。
◎若约束为零位移约束时，此步则可省去。
●举例：
为具体说明，现举一具有四个方程(二个节点)的简例。其节点平衡方程为
1
111111212
1
111211112
2
221212222
2
212222122x
y
x
yR
KKKKu
R
KKKKv
R
KKKKu
R
KKKKv
设结构在1点受到约束=,=，则上式中、为未知的约束反力。利用划0置
1u
1
1v
2
xR
1yR
1
1
的约束处理方法，上式变为
41
1
2
1
212
212122222
21212
22212221000
0100
00
00x
yu
v
RKK
KKu
RKK
KKv
●特点：
（1）经以上处理同样可以消除刚性位移(约束足够的前提下)，去掉未知约束反力。
（2）但这种方法不改变方程阶数，利于存贮。
（3）不过，若是要求出约束反力，仍要重新计算各个划去的总刚元素。
3. 乘大数法
●适用：这种方法同样适用于边界节点位移分量为已知(含为0)的各种约束。
●做法：
（1）将整体刚度矩阵中与给定节点位移相应的主对角线元素乘上一个大数，如；20
10
（2）再将方程右端载荷列阵中的相应元素用己知位移和该大数及主对角线元素的乘
积来置换。其余各项均保持不变。
●举例：如上例用此法进行约束处理后，节点平衡方程组变成
20
20
11
1211
111112
20
20
2111
1112
1112
22
21212222
22
2122
212210
10
10
10
x
yKu
KKKK
Kv
KKKK
Ru
KKKK
Rv
KKKK
●特点：
（1）使用此一方法，只要大数选得足够大，就可保证求得的位移有足够的精度。
（2）由于在处理过程中，不失去总刚矩阵的任一行(列)及各个元素，便于进行程序
处理及约束反力计算。
★小结：
经过约束处理，最终建立了系数矩阵正定的2n-r阶(划行划列法)或是2n阶(划零置1法
和乘大数法)方程组。
三、后续工作
●下一步即求解此方程组，最终获得2n-r个未知的位移分量。
●线性方程
的解法有直接法和和迭代法两大类：
