计算力学平面桁架程序说明选编.docx

题目：编制一个能够计算任意形状的平面桁架（线弹性）的有限元程序一、基本思想1、先建立一个平面桁架的力学模型，再把结构整体拆开，分解成若干个单元（在杆件结构中就把每个杆件取作一个单元，基础数据有：结点坐标、结点编号、单元编号、单元的节点号、单元的属性、节点荷载、节点自由度）。2、单元分析：采用刚度影响系数建立单元节点的平衡方程（即单元刚度方程）。3、整体分析：把所有单元的刚度方程组合成整体的刚度方程，这是一组以节点物理量为未知量的线性方程组，引入边界条件求解该方程组。4、计算位移和内力二、程序实现整个程序的流程：1、主程序%整个执行过程（总控）process.mclear;InFileName=input(请输入基础数据文件名(默认为：model data.txt）：,s);ifisempty(InFileName)InFileName=model data.txt;end[ Joint,Elem,fidout ] = Struss_ReadData(InFileName );fori=1:length(Elem.Def(:,1)) [Len(i),Orient(i,:)]=LCS(Elem.Def(i,:),Joint.Coord);endStructStiff=strus_Stiff( Joint,Elem,Len,Orient);%计算刚度Node_Disp=Node_D