《跨连续梁内力计算程序程序.docx》由会员分享，可在线阅读，更多相关《跨连续梁内力计算程序程序.docx（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、跨连续梁内力计算程序程序六跨连续梁内力计算程序讲明文档一程序适用范围本程序用来解决六跨连续梁在荷载作用下的弯矩计算。荷载能够是集中力Fp作用于跨中、分布荷载q分布全垮、集中力偶m作用于结点的任意组合情况。端部支承可为铰支或固支。二程序编辑方法使用TurboC按矩阵位移法的思路进行编辑，用TurboC中的数组来完成矩阵的实现，关键的求解K=P的步骤用高斯消元法。三程序使用方法运行程序后，根据提示，依次输入结点编号，单元编号，单元长度，抗弯刚度EI的倍数，集中力，均部荷载，集中力偶，各个数据间用空格隔开，每一项输入完毕后按回车键，所有数据输入完毕后按任意键输出结果。输出结果中包括输入的数据以便校核
2、，角位移的值以1/EI为单位以及每个单元的左右两端弯矩值。四程序试算1算例1算力图示：输入数据：结点：1234560；单元：123456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(11.383738,-1.434142,-8.980504,14.053733,-10.192107,10.048027,0)EI单元编号123456左端弯矩右端弯矩2算例2算例图示：6EI8kN/m4m3m2m8mkN/m12365474kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEIEIEI1.5
3、EI1.52EI输入数据：结点：0123456；单元：123456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(0,1.686721,-10.080218,14.871010,-12.183221,17.195206,-26.597603)EI单元编号123456左端弯矩右端弯矩弯矩图如下：3算例3算例图示：6EI8kN/m4m3m2m8mkN/m1236544kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEI1.5EIEIEI2EI1.5输入数据：结点：0123450；单元：12
4、3456；长度：466846；EI：1121；Fp：0128060；q：800406；m：00-801000运行程序如下：结果为：角位移为：1(0,1.653925,-9.949034,14.264283,-10.248253,10.062063,0)EI单元编号123456左端弯矩右端弯矩弯矩图如下：4算例4算例图示：6EI8kN/m4m3m2m8mkN/m12365474kN/m3m3m3m2m6m12kN8kN8kN.m6kN10kN.mEIEIEI1.5EI1.52EI输入数据：结点：1234567；单元：123456；长度：466846；EI：1121；Fp：0128060；q：80
5、0406；m：00-801000运行程序如下：结果为：角位移为：1(11.364772,-1.396211,-9.113262,14.660626,-12.126579,17.179023,-26.589511)EI单元编号123456左端弯矩右端弯矩弯矩图如下：五源程序#include#include定义变量intjoint7;结点编号intunit6;单元编号floatlength6;长度floatEI6;抗弯刚度floatP6;集中力floatq6;均部荷载floatm7;集中力偶doubleI6;线刚度intnumber=7,i,j;doubleK77=0;整体刚度矩阵doublek6
6、22;单元刚度矩阵doubleMP62,Mq62,Mlast62,M7=0,tatleM7,mm7,Mqueue7=0,antiM7=0;依次为：集中力、均布荷载引起的固端弯矩，最终杆端弯矩，全部荷载引起的弯矩，总的结点荷载，按结点编号排列的结点集中力偶，按结点编号排列的固端弯矩，等效结点弯矩doubleangle7=0;角位移voidinput();输入函数voidsolve();求解函数voidoutput();输出函数voidGauss();高斯消元法voidmain()input();solve();output();voidinput()clrscr();printf(Pleaseinputdata.nnThejointnumber:);for(i=0;i=1;i+)Mqueuejointi-1=Mi;mmjointi-1=mi;for(i=0;i=1&jointi+1=1&i=j;m-)Klm+=-Kjm*Klj;tatleMnumber-1/=Knumber-1number-1;Knumber-1number-1=1;for(i=0;i=0;i-)for(j=number-1;ji;j-)anglei=anglei-Kij*anglej;