- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
END
石家庄铁道学院
单元定位向量回顾
节点从小号到大号编； 结点位移按
编 号 原 则
u v
顺序编。
先编可动结点，后编不动结点
根据位移连续条件，各单元在结构坐标系 下的杆端位移与它的始、末端所对应的结点位 移是相等的，因此当对应于所有结点的结点位 移分量的总体编号已知时，根据各单元的始、 末端的结点号，即可确定该单元6个杆端位移 分量的总体编号。
把每个单元刚度矩阵的4个子块按其下标 的号码送到结构原始刚度矩阵中相应的位置上 去（对号入座）；各单元具有相同下标的子块 被送到总刚中同一位置上叠加起来；而 k0 中 没有子块入座的子块以0子块补入。
单元的子块搬入总刚度矩阵中的位置，完全取 决于结构结点编号。对同一结构，如果改变了 结点的编号，则总刚度矩阵完全不同。
石家庄铁道学院
(1,2,3) 1
3
1
(4,5,6) 2
4
5
2
(13,14,15)
(10,11,12) 4
② k C K 22 ② k 52
② ②
(7,8,9) 3
在R中 4 5 6 13 14 15
② 4 k 5 25 6 ② 13 k 14 55 15
q2(i)
任意值
① ① ② ③
1 2 4 1
L/2 0 L/2 L/2
L
任意值 任意值
-10
任意值 任意值
石家庄铁道学院
单元固端力的计算PAD设计
局部变量定义
efix(i) j=1,6 ff(j)=0
1 Case 2 2 3
Case 3 Case 1
集中荷载
k mf(i),sl al(k)
中 间 变 S1B1*CX*CX+B2*CY*CY，S1(B1－B2)*CX*CY, S3B3*CY， 量
S4B1*CY*CY+B2*CX*CX, S5B3*CX, S6B4
给 上 三 角 赋 值
C(1,1)S1, C(1,2)S2, C(1,3)S3, C(1,4)-S1, C(1,5)-S2, C(1,6)S3, C(2,2)S4,C(2,3)-S5,C(2,4)-S2,C(2,5)-S4, C(2,6)-S5, C(3,3)2*S6,C(3,4)-S3,C(3,5)S5,C(3,6)S6,C(4,4)S1,C(4,5)S2, C(4,6)-S3, C(5,5)S4, C(5,6)S5, C(6,6)2*S6 I=2, 6 J=1, I－1 C(I,J) C(J,I) 给下三角赋值
综 合 结 点 荷 载
形成等效结点荷 载
坐标转换。
反号叠加成[PE]。
形成综合 结点荷载
[P]=[PD]+[PE]
涉及到4个子程序efix、trans、eload、load
石家庄铁道学院
直接结点荷载向量[PD]的形成
与直接结点荷载有关的变量和数组
npj mj（npj）
具有直接结点荷载作用的结点数
具有直接结点荷载作用的结点所 对应的结点整体编号数组
石家庄铁道学院
示例
10kN
20kN.m
1
2 5
25kN 15kN
3 6
注意：可动支座 7、8上结点荷载 的处理。
Y
4 9
X
30kN.m
35kN
7
8
结点荷载值 XD YD MD
有荷载作用 对应的结点总 对应的荷载分量 的结点序号 编号k=mj (i) 号3k-2,3k-1,3k
1 2 3 4
1 3 7 8
3 6 3 6
石家庄铁道学院
单元定位向量理论回顾
如内部结点全部是刚结点，则节点位移编 号和相应的节点位移分量的编号有一种简 单的对应关系： 节点编号 节点位移编号
i
ui 3i 2
vi 3i 1
i 3i
对于任意单元，均可根据其始末端的结点号 确定其杆端位移分量所对应的总体编号，也即可 以确定单刚中的元素在总刚中的位置。
直接给单刚元素赋值， 形成[C]。 根据单元始末结点号 实现。
对号入座的过程
由3个子程序Stiff、Locat、Wstiff实现
石家庄铁道学院
子程序Stiff——形成单 元刚度矩阵[C]
EA EA 0 0 0 0 L L 12 EI 6 EI 12 EI 6 EI 2 0 3 2 0 3 L L L L 6 EI 4 EI 6 EI 2 EI 0 2 0 2 L L L L EA EA 0 0 0 0 L L 12 EI 6 EI 12 EI 6 EI 0 0 L3 2 L L3 L2 6 EI 2 EI 6 EI 4 EI 2 0 0 2 L L L L
石家庄铁道学院
单元定位向量的PAD设计
LOCAT (IE)
数组，变量定义
局部变量的定义
将始末端结点号赋给I、J
I JL(IE)
J JR(IE)
II(1) 3*I-2， II(2) 3*I-1 II(3) 3*I， II(4) 3*J-2
计算II(6)
II(5) 3*J-1，II(6) 3*J
单刚矩阵的元素在总刚中的位置由什么确定？
石家庄铁道学院
组集总刚的程序实现
对于第IE单元（ IE 1,2,3, NE ），根据定位向 量数组 II (6) 把单刚矩阵中处于第 I 行 J 列 （ i, j 1, 2,, 6）的元素 Cij 送到总刚 [ R ] 中的 第 II (i )行II ( j )列的位置上。
qj（npj,1）
qj（npj,3） qj（npj,2）
直接结点荷载XD
直接结点荷载YD
qj（npj,3）
直接结点荷载MD
石家庄铁道学院
直接结点荷载向量[PD]的形成
k=mj(i)
第i个直接 结点荷载作 用的结点
结点整 体编号
结点荷 载分量
qj(i,1)=XD qj(i,2)=YD qj(i,3)=MD
石家庄铁道学院
(1,2,3) 1
3
1
(4,5,6) 2
4
5
2
(13,14,15)
4 (10,11,12) 单元
1 2 3 4
始端 末端
(7,8,9) 3Βιβλιοθήκη vi2 5 2 5
I
J
ui
1 4 1 4
i
结点位移分量编号
uj
4 13 10 7
vj
5 14 11 8
j
6 15 12 9
1 2 1 2
2 5 4 3
END
石家庄铁道学院
组集总刚回顾
1.分块形式的总刚 k0 的行数和列数与结 构的节点数相等。 2.每一子块均为3X3的方阵，结构的阶 数等于3倍的节点数。 3.各单元对结构原始刚度矩阵有影响的 子块的两个下标与结构原始刚度矩阵中 同一子块的两个下标完全相同。
石家庄铁道学院
对号入座组集总刚矩阵
Wstiff
局部变量的定义
i=1,n ie=1,ne
j=1,n
r(i,j) 0
call stiff(ie) call locat(ie) 前后处理结合法PAD
i=2, n
END
j=1, i－1
r(i,j) r(j,i)
石家庄铁道学院
综合结点荷载向量的形成
形成直接结点 荷载
由直接结点荷载信息 形成[PD]。 计算单元固端力。
cos sin
STIFF(IE)
数组、变量定义
IJL(IE), JJR(IR) CX(X(J)－X(I))/AL(IE) CY(Y(J)－Y(I))/AL(IE)
B1EA(IE)/AL(IE), B212*EI(IE)/AL(IE)**3 B36*EI(IE)/AL(IE)**2, B42*EI(IE)/AL(IE)
1．理论回顾
平 面 刚 架 单 元 的 单 刚 矩 阵
k
e
石家庄铁道学院
单刚矩阵理论回顾
K
e
T K T
T
e
结构坐标系下的单元刚度矩阵。
具体形式参见课本P9-P10 平面刚架单元的单刚矩阵的性质？
石家庄铁道学院
子程序Stiff的PAD设计
局部变量定义，全局变量通过调用 相应的数组名或变量名传递。
P(3K-2)
P(3K-1)
P(3K)
荷载列阵中的总分量号
石家庄铁道学院
直接结点荷载向量[PD]的PAD设计
i=1,n npj=0 P(i) 0 i=1,npj k mj(i) P(3k-2) qj(i,1) P(3k-1) qj(i,2) P(3k) qj(i,3) 直接结点荷载向 量[PD]存放在数 组P（n）中。
石家庄铁道学院
前后结合法的程序实现
实质：对固定支 座先处理，对非 固定支座后处理。
ie=1,ne k1=1,6
i ii(k1) i>n
stop
k2=1,6
j ii(k2)
j>n stop r(i,j) r(i,j)+c(k1,k2)
石家庄铁道学院
组集总刚自程序的PAD设计
数组，变量定义
Cij
具 体 实 现
RII (i ), II ( j )
对每一单元，取I,J从1到6循环； 对所有单元进行循环；
石家庄铁道学院
前后处理结合法
前 处 理 法
对所有结点的未知位移进行统一编号，将与 支撑约束相对应的位移分量编为0号，而各单 元刚度矩阵中凡与0位移项对应的元素不参加 组集总刚。这种在建立刚度刚度矩阵过程中 便引入支撑条件的方法称为前处理法。 把全部节点位移分量（包括支座在内）都先 看作是未知量而依次编号，每一单刚的所有 元素均按照定位向量所指示的位置对号入座 形成总刚，然后再处理边界条件。这种在建 立结构刚度矩阵过程中不引入支撑条件的方 法称为后处理法。
