结构力学 矩阵位移法课件.ppt

土木工程专业系列教材—结构力学 出版社 科技分社 9.9.2 子程序流程图 1）子程序tsm的流程图 总刚数组tk置零 m=0,ne-1 调用lsc 调用esm 调用elv l=0,5 i=lv(l) i==0 k=0,5 j=lv(k) j==0 tk[i][j]=tk[i][j]+ek[l][k] 返回 对单元循环 求单元常数 求ke 形成单元定位向量 对ke的行循环 ke的第l行在总刚中的行号 是 否 是 否 对ke的列循环 ke的第k列在总刚中的列号 单刚元素向总刚中对号入座 土木工程专业系列教材—结构力学 出版社 科技分社 2）子程序jlp的流程图 数组p[n]置零 j=int(pj[i][0]) k=int(pj[i][1]) l=jn[j][k] np0 i=0,np-1 l!=0 p[l]=pj[i][2] nf0 i=0,nf-1 m=int(pf[i][0]) 调用lsc 调用eff 调用elv 求 j=0,5 l=lv[j] l!=0 p[l]= p[l]+pe[j] 返回 是 否（无结点荷载） 是 否 是 否（无非结点荷载） 1 1 第一步：处理结点荷载 第二步：处理非结点荷载 是 否 对非结点荷载循环 将 叠加到F 排除定位向量的零分量 的第j个分量在F中的地址 对 的分量循环 存入数组pe[6] 荷载作用单元号 求单元m的常数 求 形成单元定位向量 对结点荷载循环 荷载作用结点号 荷载方向代码 荷载定位编码 结点荷载送入F 土木工程专业系列教材—结构力学 出版社 科技分社 2）子程序jlp的流程图 数组p[n]置零 j=int(pj[i][0]) k=int(pj[i][1]) l=jn[j][k] np0 i=0,np-1 l!=0 p[l]=pj[i][2] nf0 i=0,nf-1 m=int(pf[i][0]) 调用lsc 调用eff 调用elv 求 j=0,5 l=lv[j] l!=0 p[l]= p[l]+pe[j] 返回 是 否（无结点荷载） 是 否 是 否（无非结点荷载） 1 1 第一步：处理结点荷载 第二步：处理非结点荷载 是 否 对非结点荷载循环 将 叠加到F 排除定位向量的零分量 的第j个分量在F中的地址 对 的分量循环 存入数组pe[6] 荷载作用单元号 求单元m的常数 求 形成单元定位向量 对结点荷载循环 荷载作用结点号 荷载方向代码 荷载定位编码 结点荷载送入F 土木工程专业系列教材—结构力学 出版社 科技分社 9.9.3 算例 1）算例 【例9.6】已知各杆E=3×107kN/m2、A=0.16m2、I=0.002m4。为图示刚架准备PFF程序使用的原始数据文件。 土木工程专业系列教材—结构力学 出版社 科技分社 【解】进行离散化。对单元、结点和结点位移未知量进行编号，如图(b)所示。接下来建立一个文本文件，按上述原始数据文件的格式，输入原始数据如下（右侧蓝色文字为注释，不输入） 。 5,7,13,2,2,3.0e7 ne,nj,n,np,nf,e 0,0 (x[i],y[i])(i=0;inj;i++) 6,0 0,6 0,6 6,6 0,12 6,12 土木工程专业系列教材—结构力学 出版社 科技分社 1,2,0.16,0.002 (ij[i][0],ij[i][1],a[i],zi[i])(i=