文档介绍：Subroutine Elem_Stiff(······)
说明
Stiff= ! 单元刚度清零
Select Case (Type)
Case (1)
平面杆系结构单元
Case (2)
空间杆系结构单元
Case Default
出错信息
End Select
End Subroutine Elem_Stiff
杆系结构单元分析子程序
单元刚度总体设计
说明部分设计
Integer,Intent (in) :: · · ·
入口整型参数
Real(8),Intent(in) :: · · ·
入口实型参数
Real(8),Intent(out) :: · · ·
出口实型参数
Real(8) :: Work1, · · ·
Integer :: i,j,k, · · ·
实型和整型工作变量
杆系结构单元分析子程序
平面杆系结构设计
Select Case (Plane)
Case (1)
平面桁架元素赋值
Case (2)
平面梁柱元素赋值
Case (3)
· · · · · ·
Case Default
出错信息
End Select
空间杆系结构设计
Select Case (Space)
Case (1)
空间桁架元素赋值
Case (2)
空间梁柱元素赋值
Case (3)
交叉梁元素赋值
Case Default
出错信息
End Select
杆系结构单元分析子程序
有关单元等效结点荷载设计和进一步的考虑
1) 单元等效结点荷载设计同仿单元刚度。
2) 从各类单元刚度元素的计算,可看到要用到长度、
单元弹性特性、单元截面特性等数据。因此,要
确定存放它们的数据结构。要将它们作为出口。
3) 为计算单元等效结点荷载元素,首先要建立各种
荷载情况等效荷载表达式,它们可由积分或载常
数表得到。然后要解决荷载信息的存放结构,也
要将它们作为出口量。
4) 单元刚度矩阵、等效结点荷载矩阵都应先清零。
杆系结构整体分析
首先就全刚结点平面刚架进行讨论,然后推广。
总的思路
在单元特性搞清后,将单元拼装回去。在结点处
位移自动协调基础上,如果全部结点平衡,则求得
的结点位移将是实际结构的解。因此,整体分析就
是设法建立结点平衡方程。
坐标转换
组成结构的杆件可以各个
方向,单元分析对局部坐标,
因此,必须将物理量转为统
一坐标——整体坐标。
1) 力的转换关系
杆系结构整体分析
2) 位移转换关系
3) 转换矩阵
转换矩阵是正交矩阵。
杆系结构整体分析
4) 杆端力转换
5) 杆端位移转换
6) 刚度方程的转换
如果记称为整体单元刚度矩阵
则
这就是整体坐标下的单元刚度方程。
本节以后的讨论认为
都是对整体坐标的
杆系结构整体分析
结点平衡方程的建立
1) 一简单例子(如图)
图中有两套编号,红的
是单元杆端编号,黑的是
结构整体编号。
1-1) 结点示意
1
2
1
2
2
1
①
②
③
图中蓝色的表示结点荷载(已知),红色的表示
杆端力(未知的), 、分别1、2单元杆端力
子矩阵。对1、4结点“荷载”含有未知反力。
2
1-2) 结点平衡
由示意图可见,结构
结点的平衡方程为
杆系结构整体分析
从例图可见,其全部结点
平衡方程为
1
2
1
2
2
1
①
②
③
若记
2
杆系结构整体分析
式中[I]、[0]分别为单位和零矩阵。
若引入矩阵记号
,则结点平衡方程可改写作
这一结论虽然是由一个例子得到的,但是显然
对一切结构都是成立的。问题在于不同结构,[A]
矩阵是不同的。
杆系结构整体分析
杆端位移用结点位移来表示
1
2
1
2
2
1
①
②
③
仍以简单例子来说明
若记
由结点、杆端位移的协调条件,可得[]、[]
的对应关系为
式中[A]T是前面力关系[A]的转置,因此[A]T称
为位移转换矩阵。