平面桁架程序计算原理及程序编制.ppt

结构有限元程序设计基本原理 ;5-01 矩阵位移法 5-02 矩阵位移法算例 5-03 最小总势能原理的应用 5-04 矩阵位移法的求解步骤 5-05 结构计算简图的数据结构 5-06 位移未知数的确定 5-07 单根杆件的分析 5-08 结构总刚度矩阵的形成 5-09 杆件内力的计算 5-10 能量原理和矩阵位移法 ;5-1 矩阵位移法 ;求解方法;5-2 矩阵位移法算例 ;根据虎克定律，对于任意一根杆，有： ;由单根杆件的变形几何关系可得 ;对于每根杆件（以两端节点编号A和B定出 角）应用上述公式，有 ;节点平衡方程为 ;将各杆的内力用位移表示的方程代入上式，有;写成矩阵形式 ;采用高斯消元法求得位移为 ;5-3 最小总势能原理的应用 ;对于整个桁架应变能是所有杆件应变能的叠加，即 ;把上述值代入应变能表达式，得到 ;现在计算外力势能。外力 产生势能的原因是当节点发生位移时，外力要作功。所作功的负值便是它们具有势能的改变量，如果取未变形位置外力的势能为零，有 ;根据最小总势能原理，在所有可能的位移状态中，真正发生的位移状态使总势能为最小。即函数对自变量的偏微商为零，即;现对各位移变量分别取偏微商后，得 ;5-4 矩阵位移法的求解步骤 ;5-5 结构计算简图的数据结构 ;5-6 位移未知数的确定 ;5-7 单根杆件的分析 ;注意⊿L= uB-uA，并且由平衡关系得到 ;设杆端分别得到平行于xoy坐标轴的位移uA，vA，uB，vB；则杆件的伸长量为 ;把⊿L代入上面的公式，得到 ;由于最终得到的平衡方程都是相对于全局坐标系建立的，上面所计算出的内力与杆轴方向一致。假设在杆端作用平行于杆轴的力FA、FB，则由平衡方程 ;把 、 、 、 排列成列向量 ，则有 ;5-8 结构总刚度矩阵的形成 ;5-9 杆件内力的计算 ;平面桁架计算程序的编制;矩阵位移法求解一般步骤：;程序流程;根据平面桁架的特性进一步展开为：;根据平面桁架的特性进一步展开为：;{总体数据结构的设计} 平面桁架程序的总体数据结构设计，包括全局量（标识符）的存储安排与说明。;一个平面桁架结构的数据有三个方面： 结构本体描述数据（NW, IESG, NU, X, Y, HL, HR） 性质数据（EF） 荷载数据（PX, PY）;MODULE PDATA REAL, PUBLIC :: BAR_HL(100),BAR_HR(100), BAR_TENSION(100),BAR_EF(100) REAL, PUBLIC :: NODP_X(80),NODP_Y(80) END MODULE REAL B(4),DP(120),R(1200),C(4,4) INTEGER II(4),V(120),NW,IESG,NU;4 ‘ NW 4 ‘ IESG 2 ‘ NU 0 1 ‘ X;Y 1 1 1 0 0 0 1 4 100 ‘ HI;HL;EF 1 2 100 2 3 100 1 3 100 1 0 ‘ PX;PY 1 0 ;{输入描述平面桁架本体结构和性质的数据} 结构本体描述数据（NW, IESG, NU, X, Y, HL, HR） 性质数据（EF）;READ(1,*) NW READ(1,*) IESG READ(1,*) NU NDISP=NU+NU WRITE(2,*) NW=,NW, IESG=,IESG, NU=,NU WRITE(2,*) NODE X Y‘ DO 10 I=1,NW READ(1,*) NODP_X(I),NODP_Y(I) WRITE(2,*) I,NODP_X(I),NODP_Y(I) CONTINUE WRITE(2,*) ELEMENT HL HR EF‘ DO 20 NT=1,IESG READ(1,*) BAR_HL(NT),BAR_HR(NT),BAR_EF(NT) WRITE(2,*) BAR_HL(NT),BAR_HR(NT),BAR_EF(NT) 20 CONTINUE ! LIST DATA, NOW TO ORGANIZE THE GLOBAL MATRIX;4 ‘ NW 4 ‘ IESG 2 ‘ NU 0 1 ‘ X;Y 1 1 1 0 0 0 1 4 100