#include
#include
#include
#define M 5
int n,nc,nn,m,e,f;//节点总数,固定节点数,自由度数,杆件数int io,jo;//单根杆对号指示数
int ihl[M],ihr[M];//杆件左右节点号
double a[M];//各杆截面积
double mm[M];//杆件质量
double ea[M];//杆件EA的值
double x[M],y[M];//节点坐标
double dp[M];//总体系下的节点载荷
double t[2];//0,1分别为坐标转换矩阵的cos(),sin()
double c[2][2];//总体系下的单刚
double clxy[3];//0,1,2分别为杆长,正弦,余弦
double h[M];//杆件轴力
double r[M][M];//总刚度阵
double rd;//桁架轴力杆局部系单刚
double u[M];//桁架节点位移
double v[2];//存放节点位移差
double d[M];//LDLT分解时的D矩阵的对角线元素
double l[M][M];////LDLT分解时的D矩阵的对角线元素double fdp[M];//总体系下支座反力
void iojo(int k)//计算对号指示数io,jo
{
int i,j;
i=ihl[k-1];//k号杆左节点号进入i
j=ihr[k-1];//k号杆节点右号进入i
io=2*(i-nc-1);//uxi前未知位移的个数
jo=2*(j-nc-1);//uyi前未知位移的个数
}
void ch(int k)//计算杆长与方向余弦函数
{
int i,j;
i=ihl[k-1];//k号杆左节点进入i
j=ihr[k-1];//k号杆右节点进入j
clxy[1]=x[j-1]-x[i-1];//k号杆x坐标差
clxy[2]=y[j-1]-y[i-1];//k号杆y坐标差
clxy[0]=sqrt(clxy[1]*clxy[1]+clxy[2]*clxy[2]);//k号杆长
clxy[1]=clxy[1]/clxy[0];//k号杆件x轴余弦
clxy[2]=clxy[2]/clxy[0];//k号杆件y轴余弦
void stif(int k)//计算k号杆件总体系下的单元刚度阵
{
int i,j;
ch(k);//调用ch(),计算k号杆件杆长与余弦
t[0]=clxy[1];
t[1]=clxy[2];
rd=ea[k-1]/clxy[0];
for(i=0;i<2;i++)
{
for(j=0;j<2;j++)
{
c[i][j]=t[i]*t[j]*rd;
}
}
}
void dor()//总体系下的总刚度阵的组集
{
int i,j,k;
for(i=0;i { for(j=0;j { r[i][j]=0.0;//总刚度阵清0 } } for(k=1;k<=m;k++) { iojo(k);//调用k的对号指示函数,从而确定组集位置 stif(k);//第k号杆件的总体系下的单刚 if(io>=0) { for(i=io+1;i<=io+2;i++) { for(j=io+1;j<=io+2;j++) { r[i-1][j-1]+=c[i-io-1][j-io-1];//在r中io+1,io+2行以及io+1,io+2列位置累加k的单刚 } for(j=jo+1;j<=jo+2;j++) { r[i-1][j-1]-=c[i-io-1][j-jo-1];//在r中io+1,io+2行以及jo+1,jo+2列位 置累加k的负单刚 r[j-1][i-1]=r[i-1][j-1]; } } for(i=jo+1;i<=jo+2;i++) { for(j=jo+1;j<=jo+2;j++) { r[i-1][j-1]+=c[i-jo-1][j-jo-1];//在r中jo+1,jo+2行以及jo+1,jo+2列位置累加k的单刚 } } } else if(jo>=0)//如果io<0,即左节点为固定节点,jo>=0,右端为可动节点,则只在jo+1,jo+2对角分块位置累加 { for(i=jo+1;i<=jo+2;i++) { for(j=jo+1;j<=jo+2;j++) { r[i-1][j-1]+=c[i-jo-1][j-jo-1];//在r中jo+1,jo+2行以及jo+1,jo+2列位置累加k的单刚 } } } } } void choldlt()//总刚度阵的LDLT分解 { int i,j,k; double m,t[M][M]; for(i=0;i { l[i][i]=1.0; } d[0]=r[0][0];//d[0]=r[0][0] for(i=1;i { for(j=0;j { m=0.0; for(k=0;k { m+=t[i][k]*l[j][k]; } t[i][j]=r[i][j]-m; l[i][j]=t[i][j]/d[j]; //计算l[i][j] m=0.0; for(k=0;k { m+=t[i][k]*l[i][k]; d[i]=r[i][i]-m; //计算d[i] } l[j][i]=l[i][j]; } } } void trildlt()//回代求节点位移 { int i,k; double m,y[M]; y[0]=dp[0]; for(i=1;i { m=0.0; for(k=0;k { m+=l[i][k]*y[k]; y[i]=dp[i]-m; //计算y[i] } } u[nn-1]=y[nn-1]/d[nn-1]; for(i=nn-1;i>=0;i--) { m=0.0; for(k=i+1;k { m+=l[k][i]*u[k]; u[i]=y[i]/d[i]-m; //计算u[i] } } } void doh()//计算杆件的轴力 { int i,k; for(k=1;k<=m;k++) { iojo(k);//调用第k号杆件的左右端点的位移指示数 for(i=0;i<2;i++)//计算每个节点2个自由度循环 { if(io<0)//把右节点的2个位移存入v[0],v[1] { v[i]=u[jo+i]; } else//把右节点的2个位移存入v[0],v[1] { v[i]=u[jo+i]-u[io+i]; } } stif(k);//计算第k号杆件总体系单刚,存入[2] h[k-1]=0.0;//数组h[k-1]清零 for(i=1;i<=2;i++)//对两个位移循环 { h[k-1]=h[k-1]+t[i-1]*v[i-1]*rd;//轴力存入h[k-1] } } } void dowt() //考虑自重,且规定y轴竖直向上 { int k,ko,i; double g; //g为重力 printf("请输入杆件质量:\n"); for(i=0;i { printf("mm[%d]=",i+1); scanf("%lf",&mm[i]); } for(k=1;k<=m;k++) //对桁架杆件循环 { iojo(k); //调用函数 ch(k); g=mm[k-1]*9.80665; //重力计算公式 if(io>=0) //左节点为自由节点 { dp[io+1]=dp[io+1]-(g/2); //左节点的y轴荷载减少二分之一重力 dp[jo+1]=dp[jo+1]-(g/2); //右节点的y轴荷载减少二分之一重力 } else if(jo>=0) //若右节点为自由节点, 则仅有右节点做如下处理 { ko=io+2*nc; //定义反力指示数ko 等于2*(固定节点号-1) dp[jo+1]=dp[jo+1]-(g/2); fdp[ko+1]=fdp[ko+1]-(g/2); //支座反力叠加重力的一半 } } } void dofanli() //计算反力 { int k,ko; for(k=1;k<=m;k++) //对杆件循环 { iojo(k); //引用函数 ch(k); ko=io+2*nc; //记录左节点 if(io<0) //左节点为固定节点 { fdp[ko]=fdp[ko]-h[k-1]*clxy[1]; //为第i个节点x轴向力加杆件轴力反力 fdp[ko+1]=fdp[ko+1]-h[k-1]*clxy[2]; //为第i个节点y轴向力加杆件轴力反力 } } } void verify() //强度校核函数 { int k,i; double sigema,sigema0,sigema1; //定义应力,拉伸许用应力,压缩许用应力 printf("请输入各杆截面积:\n"); //截面面积输入 for(i=0;i { printf("a[%d]=",i+1); scanf("%lf",&a[i]); } printf("请输入杆件拉伸许用应力:\n"); //杆件拉伸许用应力输入scanf("%lf",&sigema0); printf("请输入杆件压缩许用应力(输入正数):\n"); //杆件压缩许用应力输入scanf("%lf",&sigema1); for(k=1;k<=m;k++) //对杆件循环 { sigema=h[k-1]/a[k-1]; //计算公式轴力与面积之商 if(sigema>sigema0||sigema<-1*sigema1) //对应力,许用应力进行比较(注:压应力为负值,所以不小于压缩许用应力) { printf("第%d根杆件超过许用应力,为危险杆件,请增加横截面积或更换其他材料\n",k); } } } void assemble() //装配应力计算 { int k,ko; double l; //定义杆件被拉长l printf("请输入存在装配应力的杆件号:\n"); scanf("%d",&k); printf("请输入杆件装配时的拉长长度:\n"); scanf("%lf",&l); iojo(k); //引用函数 ch(k); h[k-1]=h[k-1]+ea[k-1]*l/clxy[0]; //储存装配杆件的应力值if(io>=0) { dp[io]=dp[io]+ea[k-1]*l*clxy[1]/clxy[0]; //左节点x轴方向附加载荷增加P=△l*EA/l乘其方向余弦 dp[io+1]=dp[io+1]+ea[k-1]*l*clxy[2]/clxy[0]; //y轴方向同上操作 dp[jo]=dp[jo]-ea[k-1]*l*clxy[1]/clxy[0]; //注:右节点与左节点附加装配应力相反 dp[jo+1]=dp[jo+1]-ea[k-1]*l*clxy[2]/clxy[0]; } else { ko=io+2*nc; dp[jo]=dp[jo]-ea[k-1]*l*clxy[1]/clxy[0]; //注:右节点与左节点附加装配应力相反 dp[jo+1]=dp[jo+1]-ea[k-1]*l*clxy[2]/clxy[0]; printf("%f,%f\n",dp[jo],dp[jo+1]); fdp[ko]=fdp[ko]+ea[k-1]*l*clxy[1]/clxy[0]; //固定节点ihl反力叠加装配应力 fdp[ko+1]=fdp[ko+1]+ea[k-1]*l*clxy[1]/clxy[0]; } } void tem() //计算温度应力 { int k,ko; double t0,arf,t1,t2,detal; //定义变量,温差,热膨胀系数,初始温度,最终温度,温变引起的长度变化 printf("请输入初始温度\n"); //变量输入 scanf("%lf",&t1); printf("请输入最终温度\n"); scanf("%lf",&t2); printf("请输入杆件的热膨胀系数\n"); scanf("%lf",&arf); t0=t2-t1; for(k=1;k<=m;k++) { iojo(k); //引用函数 ch(k); detal=-1*arf*t0*clxy[0]; //等效为装配应力杆件受压 h[k-1]=h[k-1]+ea[k-1]*detal/clxy[0]; if(io>=0) { dp[io]=dp[io]+ea[k-1]*detal*clxy[1]/clxy[0]; //左节点x轴方向附加载荷增加P=△l*EA/l乘其方向余弦 dp[io+1]=dp[io+1]+ea[k-1]*detal*clxy[2]/clxy[0]; //y轴方向同上操作 dp[jo]=dp[jo]-ea[k-1]*detal*clxy[1]/clxy[0]; //注:右节点与左节点附加温度应力相反 dp[jo+1]=dp[jo+1]-ea[k-1]*detal*clxy[2]/clxy[0]; } else { ko=io+2*nc; dp[jo]=dp[jo]-ea[k-1]*detal*clxy[1]/clxy[0]; //注:右节点与左节点附加温度应力相反 dp[jo+1]=dp[jo+1]-ea[k-1]*detal*clxy[2]/clxy[0]; fdp[ko]=fdp[ko]+ea[k-1]*detal*clxy[1]/clxy[0]; //固定节点ihl反力叠加温度应力 fdp[ko+1]=fdp[ko+1]+ea[k-1]*detal*clxy[1]/clxy[0]; } } } int main() { int i; printf(" *****************求解平面桁架节点位移与杆端力程序 *****************\n\n"); printf("\n ------------------输入数据时请用空格或回车符间隔------------------\n"); printf("\n\n请输入桁架节点总数n,固定节点数nc,杆件数m:\n"); scanf("%d,%d,%d",&n,&nc,&m); nn=2*(n-nc); printf("请输入节点坐标:\n"); for(i=0;i { printf("x[%d]=",i+1); scanf("%lf",&x[i]); printf("y[%d]=",i+1); scanf("%lf",&y[i]); } printf("输入杆件左右节点号:\n"); for(i=0;i { printf("ihl[%d]=",i+1); scanf("%d",&ihl[i]); printf("ihr[%d]=",i+1); scanf("%d",&ihr[i]); } printf("请输入杆件的EA值[ea]:\n"); for(i=0;i { printf("ea[%d]=",i+1); scanf("%lf",&ea[i]); } printf("请输入节点载荷[dp]:\n"); for(i=0;i { printf("dp[%d]=",i+1); scanf("%lf",&dp[i]); } dor(); choldlt(); printf("是否需要考虑自重,需要请输入‘1’,不需要请输入‘0’。\n"); scanf("%d",&e); if(e==1) { dowt(); e=0; } printf("是否需要考虑装配应力,需要请输入‘1’,不需要请输入‘0’。\n"); scanf("%d",&e); if(e==1) { assemble(); e=0; } printf("是否需要考虑温度应力,需要请输入‘1’,不需要请输入‘0’。\n"); scanf("%d",&e); if(e==1) { tem(); e=0; }//插入结束 trildlt(); printf("\n\n请输出各个节点位移列向量[U]:\n"); for(i=0;i { printf("u[%d]=%f\n",i+1,u[i]); } doh(); printf("输出杆件轴力:\n"); for(i=0;i { printf("h[%d]=%8.4f\n",i+1,h[i]); } printf("是否需要考虑强度校核,需要请输入‘1’,不需要请输入‘0’。\n"); scanf("%d",&e); if(e==1) { verify(); e=0; } printf("是否需要计算支座反力,需要请输入‘1’,不需要请输入‘0’。\n"); scanf("%d",&f); if(f==1) { dofanli(); printf("输出节点反力:\n"); for(i=0;i<=2*nc-1;i++) { printf("fdp[%d]=%8.4f\n",i,fdp[i]); } } return 0; } #include 15-1 多跨静定梁 031=+-=+'=qx qa qx y Q D X a x 3 1 = 2 当l X = α cos 2 l q Q B -= αα0sin sin =--qx y N A X 因在梁上的总载不变:ql l q =11 αcos 11 111q l l q q l l q === ()()()111221122111 1 1 d p l V f H M H H x a p a p l V M b p b p l V A A C B A B A A -?= ===+==+= ∑∑∑ f M H V V V V C A B B A A = = = f=0时,H A =∞,为可弯体系。 简支梁: ① 1 P V Q A - = ()a x P V A- - 1 H=+H A ,(压为正) ②()y H a x p x V M A A - - - = 1 1 即y H M M A - = D截面M、Q、N ()y H a x p x V M A A x ? - - - = 1 1 即y H M M A x - = ? ? ? ? sin sin sin cos H Q N H Q Q x x + = - = 说明:?随截面不同而变化,如果拱轴曲线方程()x f y=已知的话,可利用 dx dy tg= ?确定?的值。 二.三铰拱的合理轴线(拱轴任意截面 = = Q M ) 据:y H M M A ? - = 当0 = M时, A H M y = M是简支梁任意截面的弯矩值,为变值。 说明:合理拱轴材料可得到充分发挥。 f M H c A =(只有轴力,正应力沿截面均匀分布) c M 为简支跨中弯矩。 空间桁架静力分析程序及算例1、变量及数组说明 2、空间桁架结构有限元分析程序源代码 !主程序(读入文件,调用总计算程序,输出结果) CHARACTER IDFUT*20,OUTFUT*20 WRITE(*,*) 'Input Data File name:' READ (*,*)IDFUT OPEN (11,FILE=IDFUT,STATUS='OLD') WRITE(*,*) 'Output File name:' READ (*,*)OUTFUT OPEN(12,FILE=OUTFUT,STATUS='UNKNOWN') WRITE(12,*)'*****************************************' WRITE(12,*)'* Program for Analysis of Space Trusses *' WRITE(12,*)'* School of Civil Engineering CSU *' WRITE(12,*)'* 2012.6.25 Designed By MuZhaoxiang *' WRITE(12,*)'*****************************************' WRITE(12,*)' ' WRITE(12,*)'*****************************************' WRITE(12,*)'*************The Input Data****************' WRITE(12,*)'*****************************************' WRITE(12,100) READ(11,*)NF,NP,NE,NM,NR,NCF,ND WRITE(12,110)NF,NP,NE,NM,NR,NCF,ND 100 FORMAT(6X,'The General Information'/2X,'NF',5X,'NP',5X,'NE',5X,'NM',5X,'NR',& 5X,'NCF',5X,'ND') 110 FORMAT(2X,I2,6I7) NPF=NF*NP NDF=ND*NF CALL ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF) END !******************************************************************** !总计算程序 SUBROUTINE ANALYSE(NF,NP,NE,NM,NR,NCF,ND,NPF,NDF) 平面桁架 程序: #include "stdio.h" #include "math.h" #define unitmax 30 #define pointmax 30 #define matermax 10 #define Pmaxnum 20 #define bindmax 20 main() { int unit[unitmax][4],bind[bindmax][3],number[5]; float point[pointmax][2],material[matermax][2],P[Pmaxnum][3]; float allarray[pointmax*2][pointmax*2]; void readdata(int unit[unitmax][4],float point[pointmax][2],int number[5], int bind[bindmax][3],float material[matermax][2],float P[Pmaxnum][3]); void all(float point[pointmax][2],int unit[unitmax][4],float material[matermax][2], int number[5],float allarray[][pointmax*2]); void equa(float allarray[][pointmax*2],float P[Pmaxnum][3],int bind[][3], int number[]); void result(int unit[unitmax][4],float point[pointmax][2],int number[5], float material[matermax][2],float allarray[][pointmax*2]); readdata(unit,point,number,bind,material,P); all(point,unit,material,number,allarray); equa(allarray,P,bind,number); result(unit,point,number,material,allarray); } /******************************************************************/ void readdata(int unit[unitmax][4],float point[pointmax][2],int number[5], int bind[bindmax][3],float material[matermax][2],float P[Pmaxnum][3]) 常见问题解答七 第七章平面桁架 1.屋盖结构体系分为哪两大类?各有何适用? 答:屋盖结构体系分为无檩屋盖与有檩屋盖。 (1)无檩屋盖:无檩屋盖一般用于预应力混凝土大型屋面板等重型屋面,将屋面板直接放在屋架或天窗架上。 预应力混凝土大型屋面板的跨度通常采用6m,有条件时也可采用12m。当柱距大于所采用的屋面板跨度时,可采用托架支承中间屋架。 (2)有檩屋盖:有檩屋盖常用于轻型屋面材料的情况。如压型钢板、压型铝合金板、石棉瓦、瓦楞铁皮等。石棉瓦和瓦楞铁皮屋面,屋架间距通常为6m;当柱距大于或等于12m时,则用托架支承中间屋架。钢板和压型铝合金板屋面,屋架间距常大于或等于12m。 2.常用的屋架形式有哪几种?各有何特点和适用? 答:屋架外形常用的有三角形、梯形、平行弦和人字形等。 三角形桁架适用于陡坡屋面(i>1/3)的有檩屋盖体系,通常与柱子铰接,房屋的整体横向刚度较低,简支屋架弦杆支座处内力大,跨中内力小,截面不能充分发挥作用。支座处上下弦杆交角过小内力又较大,使支座节点构造复杂。 三角形屋架的腹杆布置有芬克式和人字式。芬克式的腹杆较多,但压杆短、拉杆长,受力相对合理,且可分为两个小桁架制作与运输,较为方便。人字式腹杆的节点较少,但受压腹杆较长,适用于跨度较小(L≤18m)的情况。人字式屋架的抗震性能好,所以在强地震烈度地区,,常用人字式腹杆的屋架。单斜式腹杆的屋架,腹杆和节点数目较多,只适用于下弦需要设置天棚的屋架,一般情况较少采用。 梯形屋架适用于屋面坡度较为平缓的无檩屋盖体系,与简支受弯构件的弯矩图形比较接近,弦杆受力较为均匀。屋架与柱的连接可以铰接也可以刚接。刚接可提高建筑物的横向刚度。腹杆体系可采用单斜式、人字式和再分式。屋架支承分为下承式和上承式。与柱刚接的屋架宜采用下承式。下承式使柱计算高度减小又便于在下弦设置屋盖纵向水平支撑,以往多采用之,但上承式使屋架重心降低,支座斜腹杆受拉,且安装方便,近年来逐渐推广使用。 图7-2 梯形屋架 人字形屋架的上、下弦可以是平行的,坡度为1/20~1/10,节点构造较为统一;也可以上、下弦具有不同坡度或者下弦有一部分水平段。多用于较大跨度。宜采用上承式。 人字形和梯形屋架的中部高度主要取决于经济要求,一般为(1/10~1/8)L,与柱刚接的梯形屋架,端部高度一般为(1/16~1/12)L,通常取为2.0~2.5m。与柱铰接的梯形屋架,端部高度可按跨中经济高度和上弦坡度来决定。人字形屋架因中高度一般为2.0~2.5m,跨度大于36m时可取较大高度但不宜超过3m;端部高度一般为跨度的1/18~1/12,人字形屋架可适应不同的屋面坡度,但与柱刚接时,屋架轴线坡度大于1/7,就应视为折线横梁进行框架分析;与柱铰接时,即使采用了上承式也应考虑竖向荷载作用下折线拱的推力对柱的不利影响,设计时要求在屋面板及檩条等安装完毕后再将屋架支座焊接固定。 图7-3 人字形屋架和平行弦屋架 平行弦桁架在构造方面有突出的优点,弦杆及腹杆分别等长、节点形式相同、能保证桁架的杆件重复率最大,且可使节点构造形式统一,便于制作工业化。 3.屋盖支撑有何作用? 答:屋盖支撑的作用: ①保证结构的空间整体作用。 ②避免压杆侧向失稳,防止拉杆产生过大的振动。 ③承担和传递水平荷载。 ④保证结构安装时的稳定与方便 3.4 静定平面桁架 教学要求 掌握静定平面桁架结构的受力特点和结构特点,熟练掌握桁架结构的内力计算方法——结点法、截面法、联合法 3.4.1 桁架的特点和组成 3.4.1.1 静定平面桁架 桁架结构是指若干直杆在两端铰接组成的静定结构。这种结构形式在桥梁和房屋建筑中应用较为广泛,如南京长江大桥、钢木屋架等。 实际的桁架结构形式和各杆件之间的联结以及所用的材料是多种多样的,实际受力情况复杂,要对它们进行精确的分析是困难的。但根据对桁架的实际工作情况和对桁架进行结构实验的结果表明,由于大多数的常用桁架是由比较细长的杆件所组成,而且承受的荷载大多数都是通过其它杆件传到结点上,这就使得桁架结点的刚性对杆件内力的影响可以大大的减小,接近于铰的作用,结构中所有的杆件在荷载作用下,主要承受轴向力,而弯矩和剪力很小,可以忽略不计。因此,为了简化计算,在取桁架的计算简图时,作如下三个方面的假定: (1)桁架的结点都是光滑的铰结点。 (2)各杆的轴线都是直线并通过铰的中心。 (3)荷载和支座反力都作用在铰结点上。 通常把符合上述假定条件的桁架称为理想桁架。 3.4.1.2 桁架的受力特点 桁架的杆件只在两端受力。因此,桁架中的所有杆件均为二力杆。在杆的截面上只有轴力。 3.4.1.3 桁架的分类 (1)简单桁架:由基础或一个基本铰接三角形开始,逐次增加二元体所组成的几何不变体。(图3-14a) (2)联合桁架:由几个简单桁架联合组成的几何不变的铰接体系。(图3-14b) (3)复杂桁架:不属于前两类的桁架。(图3-14c) 3.4.2 桁架内力计算的方法 桁架结构的内力计算方法主要为:结点法、截面法、联合法 结点法――适用于计算简单桁架。 截面法――适用于计算联合桁架、简单桁架中少数杆件的计算。 联合法――在解决一些复杂的桁架时,单独应用结点法或截面法往往不能够求解结构的内力,这时需要将这两种方法进行联合应用,从而进行解题。 解题的关键是从几何构造分析着手,利用结点单杆、截面单杆的特点,使问题可解。 在具体计算时,规定内力符号以杆件受拉为正,受压为负。结点隔离体上拉力的指向是离开结点,压力指向是指向结点。对于方向已知的内力应该按照实际方向画出,对于方向未知的内力,通常假设为拉力,如果计算结果为负值,则说明此内力为压力。 ! PTA说明: 1,数据输入文件为DATAIN.TXT中,数据依次为(NE,NJ,NR,NP,NP); (X,Y)*NJ ; ! (IJ1,IJ2,A,EI)*NE;(JR1,JR2,JR3,JR4)*NR;(PJ1,PJ2,PJ3)*NP;(PF1,PF2,PF3,PF4)*NF ! 2,NE单元总数;NJ节点总数;NR约束总数;NP节点荷载总数;NF非节点荷载数;X,Y节点坐标; ! IJ(NE,1),IJ(NE,2) 单元期终点;A 单元面积;ZI 截面惯性矩;JR(NR,1)约束结点号; ! JR(NE,2:4)横向竖向转动约束(1); PJ(NP,1:3)分别为结点号、荷载类型、荷载值; ! PF(NF,1:4)分别为单元号、荷载类型,荷载值、荷载据起点距离 !主程序 PROGRAM Plane_Truss_Analysis implicit none integer NE,NJ,NP,NF,NR,N real E integer,allocatable,DIMENSION(:,:)::IJ,JR real,allocatable,DIMENSION(:)::A,P,X,Y !动态数组定义 real,allocatable,DIMENSION(:,:)::PJ,PF,TK OPEN(1,FILE='datain.TXT',STATUS='OLD') !打开文件大datain.txt并存到标号1中 OPEN(2,FILE='dataout.TXT',STATUS='NEW') !建立文件dataout存储计算结果 READ(1,*)NE,NJ,NR,NP,NF,E N=NJ*2 !位移总数 allocate(X(1:NJ),Y(1:NJ),IJ(1:NE,2),JR(1:NR,4),A(1:NE),PJ(1:NP,3),PF(1:NF,4),TK(1:N,1:N),P(1:N) ) WRITE(2,10)NE,NJ,NR,NP,NF,E !打印表头 10 FORMAT(/1X,' ***********平面桁架力计算PTA***********'//4X,'单元数NE=',I2,12X,'结点数NJ=',I2,14X,'支座数NR=',I2,/4X,'结点荷载数NP=',I2,8X,'非节点荷载数 NF=',I2,8X,'弹性模量E=',E12.4) CALL INPUT(NE,NJ,NR,NP,NF,X,Y,IJ,A,JR,PJ,PF) !数据输入 CALL TSM(NE,NJ,E,X,Y,IJ,A,TK,N) !形成结构原始刚度矩阵 CALL JLP(NE,NJ,NP,NF,X,Y,IJ,PJ,PF,P,N) !形成结构综合节点荷载阵列 CALL ISC(NR,JR,TK,P,N) !引入结构约束条件 CALL GA USS(TK,P,N) !高斯消去法计算结构坐标下的节点位移 CALL MVN(NE,NJ,NF,E,X,Y,IJ,A,PF,P,N) !计算单元杆端内力 CLOSE(1) CLOSE(2) deallocate(X,Y,IJ,JR,A,P,PJ,PF,TK) END PROGRAM Plane_Truss_Analysis !原始数据输入 SUBROUTINE INPUT(NE,NJ,NR,NP,NF,X,Y,IJ,A,JR,PJ,PF) DIMENSION X(NJ),Y(NJ),IJ(NE,2),A(NE),JR(NR,4),PJ(NP,3),PF(NF,4) 简单平面桁架受力分析 仪22 廖嵩松 020854 内容摘要:桁架使工程中常用的一种结构,求解桁架受力有三种方法:虚位移法、节点法、截面法,他们个又有缺点,适用于不同情况,更多情况是将它们综和运用。 关键词:桁架——虚位移法——节点法——截面法——载荷 参考书目: 李俊峰,张雄,任革学,高云峰编.理论力学.北京:清华大学出版社,2001 西北工业大学理论力学教研室编.和兴锁主编.理论力学.西安:西北工业大学出版社,2001 牛学仁主编.理论力学.北京:机械工业出版社,2000 一、 桁架简化模型 桁架是工程中常用的一种结构,各构件在同一平面内的桁架称为平面桁架。简单平面桁架是指在一个基本三角框架上每增家两各杆件的同时增加一个节点而形成的桁架。它始终保持其坚固性,且在这种桁架中除去任何一个杆件都会使桁架失去稳固性。 在简单平面桁架中,杆件的数目m 与节点数目n 之间有确定关系。基本三角框家的杆件书和节点数都等于3。此后增加的杆件数(m-3)节点数(n-3)之间的比例是2:1,故有 323-=-n m 即 n m 23=+ 在计算载荷作用下平面桁架各杆件的所受力时,为简化计算,工程上一般作如下规定: (1) 各杆件都是直杆,并用光滑铰链连接; (2) 杆件所受的外载荷都作用在各节点上,各力作用线都在桁架平面内; (3) 各杆件的自重忽略不计; 在以上假设下,每一杆件都是二力构件,故所受力都沿其轴线,或为拉力,或为压力。为便于分析,在受力图中总是假设杆件承受拉力,若计算结果为负值,则表示杆件承受压力。 二、 计算桁架受力的三种方法 1、 虚位移法 接触所求杆的约束,用约束反力 代替,系统仍处于平衡状态,但有一 个自由度。假设系统沿此自由度的方 向有一微小的运动,可得出各主动力 作用点及所加约束力的一组虚位移, 根据虚位移原理可列出一个方程,解 出约束反力的大小。 例1.求解图1所示平面桁架中1 杆的约束力。 解:去掉1杆,用N1和N1’代替。假 施工平台支撑验算 支架搭设高度为7.4米, 搭设尺寸为:立杆的纵距b=1.20m,立杆的横距l=1.20m,立杆的步距h=1.50m,顶托下部采用2根50*100的方通。方通下方为桁架。 1.立杆计算: (1)荷载计算: 取1个计算单元:(1.2m*1.2m) 立杆自重:7.4m*3.5kg/m=0.26kN; 施工荷载取100kg/m2; 堆放荷载取100kg/m2; 水平杆作用在单根立杆上的重量为(5道双向): 2.4*5* 3.5kg/m=0.42kN; 单根立杆荷载总和为: N=2*1.44+0.26+0.42=3.6kN; (2)立杆稳定性验算: A=4.24cm2,i=1.6cm 计算长度l0=uh=1.75*1.5=2.6m λ= l0/i=260/1.6=162.5, φ=0.294 f=N/ΦA=3.6/(424*0.294)=28.9N/mm2<[f]=215N/mm2 满足要求。 2.方通验算: 按三跨连续梁计算: (1)变形验算: 用SAP 2000进行计算,结果如下: 最大挠度位于1.6m处,(双方通) 挠度为14mm/2=7mm<3600mm/250=14.4mm 满足要求。 (2)刚度验算: 弯矩图如下(kN.m): M max=3.54kN.m,W=15.52cm3; f=M/W=3.54/(2*15.52)=114N/mm2<[f]=215N/mm2 满足要求。 (3)支座反力: 支座反力如下: 3.桁架验算: 计算模型: a.Y-Z平面: 内力计算结果为: 上部横杆计算结果为: 下部横杆计算结果为: 25您的位置:在线学习—>在线教程—>教学内容 上一页返回目录下一页 3.4 静定平面桁架 教学要求 掌握静定平面桁架结构的受力特点和结构特点,熟练掌握桁架结构的内力计算方法——结点法、截面法、联合法 3.4.1 桁架的特点和组成 3.4.1.1 静定平面桁架 桁架结构是指若干直杆在两端铰接组成的静定结构。这种结构形式在桥梁和房屋建筑中应用较为广泛,如南京长江大桥、钢木屋架等。 实际的桁架结构形式和各杆件之间的联结以及所用的材料是多种多样的,实际受力情况复杂,要对它们进行精确的分析是困难的。但根据对桁架的实际工作情况和对桁架进行结构实验的结果表明,由于大多数的常用桁架是由比较细长的杆件所组成,而且承受的荷载大多数都是通过其它杆件传到结点上,这就使得桁架结点的刚性对杆件内力的影响可以大大的减小,接近于铰的作用,结构中所有的杆件在荷载作用下,主要承受轴向力,而弯矩和剪力很小,可以忽略不计。因此,为了简化计算,在取桁架的计算简图时,作如下三个方面的假定:(1)桁架的结点都是光滑的铰结点。 (2)各杆的轴线都是直线并通过铰的中心。 (3)荷载和支座反力都作用在铰结点上。 通常把符合上述假定条件的桁架称为理想桁架。 3.4.1.2 桁架的受力特点 桁架的杆件只在两端受力。因此,桁架中的所有杆件均为二力杆。在杆的截面上只有轴力。 3.4.1.3 桁架的分类 (1)简单桁架:由基础或一个基本铰接三角形开始,逐次增加二元体所组成的几何不变体。(图3-14a) (2)联合桁架:由几个简单桁架联合组成的几何不变的铰接体系。(图3-14b) (3)复杂桁架:不属于前两类的桁架。(图3-14c) 展示设计 https://www.doczj.com/doc/0a11465993.html, 模型云 https://www.doczj.com/doc/0a11465993.html, 桁架计算方法 房屋建筑用的桁架,一般仅进行静力计算;对于风力、地震力、运行的车辆和运转的机械等动荷载,则化为乘以动力系数的等效静荷载进行计算;特殊重大的承受动荷载的桁架,如大跨度桥梁和飞机机翼等,则需按动荷载进步履力分析(见荷载)。 支撑系统有上弦支撑、下弦支撑、垂直支撑和桁架租赁共同组成空间稳定体系。桁架的高度与跨度之比,通常采用1/6~1/12,在设计手册和规范中均有具体规定。计算次应力需考虑杆件轴向变形,可用超静定结构的方法或有限元法求解。 平面桁架一般按理想的铰接桁架进行计算,即假设荷载施加在桁架节点上(如果荷载施加在节间时,可按简支梁换算为节点荷载),并和桁架的全部杆件均在同一平面内,杆件的重心轴在一直线上,节点为可自由动弹的铰接点。 工程用的桁架节点,一般是具有一定刚性的节点而不是理想的铰接节点,由于节点刚性的影响而出现的杆件弯曲应力和轴向应力称为次应力。 从力学方面分析,桁架租赁外形与简支梁的弯矩图相似时,上下弦杆的轴力分布均匀,腹杆轴力小,用料最省;从材料与制造方面分析,木桁架做成三角形,钢桁架采用梯形或平行弦形,钢筋混凝土与预应力混凝土桁架为多边形或梯形为宜。 根据桁架杆件所用的材料和计算所得出的内力,选择合适的截面应能保证桁架租赁的整体刚度和稳定性以及各杆件的强度和局部稳定,以满意使用要求。桁架的使用范围很广,在选择桁架形式时应综合考虑桁架的用途、材料和支承方式、施工条件,其最佳形式的选择原则是在满意使用要求前提下,力求制造和安装所用的材料和劳动量为最小。桁架的整体刚度以控制桁架的最大竖向挠度不超过容许挠度来保证;平面桁架的平面外刚度较差,必须依靠支撑体系保证。 空间桁架由若干个平面桁架所组成,可将荷载分解成与桁架租赁同一平面的分力按平面桁架进行计算,或按空间铰接杆系用有限元法计算。理想状态下的静定桁架,可以将杆件轴力作为未知量,按静力学的数解法或图解法求出已知荷载下杆件的轴向拉力或压力(见杆系结构的静力分析)。 作业一 平面桁架ansys 分析 用ansys 分析图1。设250.1,100.2cm A MPa E =?=。 图1 1 设置计算类型 Preferences →select Structural →OK 2 选择单元类型 Preprocessor →Element Type →Add/Edit/Delete →Add →Link 3D finit stn 180 →OK 3 定义实常数 ANSYS Main Menu: Preprocessor →Real Constants… →Add… →select Type 1→OK →input AREA:1 →OK →Close (the Real Constants Window) 4 定义材料属性 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.0e5, PRXY:0.3 →Material →Exit 5 生成几何模型生成关键点,如图2. 图 2 ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS → 1(3,0),2(0,0),3(0,30) →OK 生成桁架 ANSYS Main Menu: Preprocessor →Modeling →Create →Lines →Lines →Straight Line →依次连接点2→1→3→1→OK如图3. 图3 6 网格划分,如图4. ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →Pick All :OK→input NDIV: 1 →OK →(back to the mesh tool window)Mesh: lines →Mesh→Pick All (in Picking Menu) →Close( the Mesh Tool window) 图 4 桁架承重架设计计算书 The latest revision on November 22, 2020 桁架承重架设计计算书 桁架承重架示意图(类型一) 二、计算公式 荷载计算:1.静荷载包括模板自重、钢筋混凝土自重、桁架自重(×1.2); 2.活荷载包括倾倒混凝土荷载标准值和施工均布荷载(×1.4)。 弯矩计算: 按简支梁受均布荷载情况计算 剪力计算: 挠度计算: 轴心受力杆件强度验算: 轴心受压构件整体稳定性计算: 三、桁架梁的计算 桁架简支梁的强度和挠度计算 1.桁架荷载值的计算. 静荷载的计算值为 q1 = 62.18kN/m. 活荷载的计算值为 q2 = 16.80kN/m. 桁架节点等效荷载 Fn = -39.49kN/m. 桁架结构及其杆件编号示意图如下: 桁架横梁计算简图 2.桁架杆件轴力的计算. 经过桁架内力计算得各杆件轴力大小如下: 桁架杆件轴力图 桁架杆件轴力最大拉力为 Fa = 105.31kN. 桁架杆件轴力最大压力为 Fb = -139.62kN. 3.桁架受弯杆件弯矩的计算. 桁架横梁受弯杆件弯矩图 桁架受弯杆件最大弯矩为M = 2.468kN.m 桁架受弯构件计算强度验算= 18.095N/mm 钢架横梁的计算强度小于215N/mm2,满足要求! 4.挠度的计算. 最大挠度考虑为简支梁均布荷载作用下的挠度 桁架横梁位移图 简支梁均布荷载作用下的最大挠度为 V = 0.425mm. 钢架横梁的最大挠度不大于10mm,而且不大于L/400 = 1.25mm,满足要求! 5.轴心受力杆件强度的计算. 式中 N ——轴心拉力或轴心压力大小; A ——轴心受力杆件的净截面面积。 桁架杆件最大轴向力为139.622kN, 截面面积为14.126cm2 . 轴心受力杆件计算强度 = 98.841N/mm2. 计算强度小于强度设计值215N/mm2,满足要求! 6.轴心受力杆件稳定性的验算. 式中 N ——杆件轴心压力大小; A ——杆件的净截面面积; ——受压杆件的稳定性系数。 轴心受力杆件稳定性验算结果列 表 ------------------------------------------------------------- ---------------- 杆件单元长细比稳定系数轴向压力kN 计算强度N/mm2 ------------------------------------------------------------- ---------------- 1 37.948 0.914 0.000 -------- 2 37.948 0.914 105.310 -------- 3 37.948 0.91 4 -52.65 5 40.770 4 40.046 0.907 -139.622 109.010 5 37.948 0.914 0.000 -------- 6 40.046 0.90 7 83.774 -------- 7 37.948 0.914 -26.327 20.385 8 37.948 0.914 -26.327 20.385 9 37.948 0.914 -39.491 30.577 10 37.948 0.914 -52.655 40.770 第二节平面静定桁架的内力计算 桁架是工程中常见的一种杆系结构,它是由若干直杆在其两端用铰链连接而成的几何形状不变的结构。桁架中各杆件的连接处称为节点。由于桁架结构受力合理,使用材料比较经济,因而在工程实际中被广泛采用。房屋的屋架(见图3-10)、桥梁的拱架、高压输电塔、电视塔、修建高层建筑用的塔吊等便是例子。 图3-10房屋屋架 杆件轴线都在同一平面内的桁架称为平面桁架(如一些屋架、桥梁桁架等),否则称为空间桁架(如输电铁塔、电视发射塔等)。本节只讨论平面桁架的基本概念和初步计算,有关桁架的详细理论可参考“结构力学”课本。在平面桁架计算中,通常引用如下假定: 1)组成桁架的各杆均为直杆; 2)所有外力(载荷和支座反力)都作用在桁架所处的平面内,且都作用于节点处; 3)组成桁架的各杆件彼此都用光滑铰链连接,杆件自重不计,桁架的每根杆件都是二力杆。 满足上述假定的桁架称为理想桁架,实际的桁架与上述假定是有差别的,如钢桁架结构的节点为铆接(见图3-11)或焊接,钢筋混凝土桁架结构的节点是有一定刚性的整体节点, 图3-11 钢桁架结构的节点 它们都有一定的弹性变形,杆件的中心线也不可能是绝对直的,但上述三点假定已反映了实际桁架的主要受力特征,其计算结果可满足工程实际的需要。 分析静定平面桁架内力的基本方法有节点法和截面法,下面分别予以介绍。 一、节点法 因为桁架中各杆都是二力杆,所以每个节点都受到平面汇交力系的作用,为计算各杆内力,可以逐个地取节点为研究对象,分别列出平衡方程,即可由已知力求出全部杆件的内力,这就是节点法。由于平面汇交力系只能列出两个独立平衡方程,所以应用节点法往往从只含两个未知力的节点开始计算。 例3-8 平面桁架的受力及尺寸如图3-12a所示,试求桁架各杆的内力。 有限单元法课程设计 有限单元法是基于连续介质力学基础上发展起来的,目前使用最广泛的数值计算方法。有限单元法解决问题的前提是各单元相邻边界的位移协调。有限单元法将连续的求解域离散为一组有限个单元组成的组合体,由细分单元去逼近求解域,由于单元的不同连接方式和形式各异的单元形状,因此可以适应几何形状复杂的求解区域;第二,利用每一个单元内的近似函数(形函数)来表示全求解域上待求的未知场函数,把一个连续的无限自由度问题变成离散的有限自由度问题,只要求出单元结点的物理量,就可以确定单元组合体上的其他未知场函数,如果选择合适的形函数,随着网格密度的减小,近似解将逐步趋向精确解;第三,有限单元法计算得到的总体刚度矩阵为稀疏带状矩阵,这样借助于电子计算机存储和计算的效率大大提高,便于处理大规模问题。 从上述有限单元法的特性可知,其计算原理简单,但由于单元连接方式和单元形状的多元化,以及近似函数的选择合适与否,使得有限元法在针对具体问题求解时比较烦琐,正是基于这样的应用背景,本论文提出了一种更简单实用的单元模型—平面等效桁架单元模型。最后,编制有限元分析程序,将这种桁架单元模型运用于钢筋混凝土结构中,模型中混凝土采用等效桁架单元,钢筋采用一维杆单元,利用混凝土等效的应力应变关系对各种构件进行弹塑性分析,并试探性的提出了单元破坏准则。用本文方法和商用有限元分析软件Ansys9.0的计算结果进行比较,经验证用本文模型在保证同等工程精度的条件下,是一种简单可行的方法。 关键词:有限单元法;平面桁架;形函数;刚度矩阵;有限元分析软件 一、桁架有限元程序流程 (1) 1、子程序说明: (1) 2、平面桁架内力计算的标识符 (2) BH 二维数组,用于存放单元截面尺寸 (2) NRES 二维数组,用于存放约束的位移值 (2) JP 二维数组,用于存放节点的荷载值 (2) ESTIF 四维数组,用于存放整体坐标系下的单元刚度矩阵 (2) 二、数据准备 (3) 三、平面桁架内力计算程序 (4) 参考文献 (11) 桁架内力的计算 3.4 静定平面桁架 教学要求 掌握静定平面桁架结构的受力特点和结构特点,熟练掌握桁架结构的内力计算方法——结点法、截面法、联合法 3.4.1 桁架的特点和组成 3.4.1.1 静定平面桁架 桁架结构是指若干直杆在两端铰接组成的静定结构。这种结构形式在桥梁和房屋建筑中应用较为广泛,如南京长江大桥、钢木屋架等。 实际的桁架结构形式和各杆件之间的联结以及所用的材料是多种多样的,实际受力情况复杂,要对它们进行精确的分析是困难的。但根据对桁架的实际工作情况和对桁架进行结构实验的结果表明,由于大多数的常用桁架是由比较细长的杆件所组成,而且承受的荷载大多数都是通过其它杆件传到结点上,这就使得桁架结点的刚性对杆件内力的影响可以大大的减小,接近于铰的作用,结构中所有的杆件在荷载作用下,主要承受轴向力,而弯矩和剪力很小,可以忽略不计。因此,为了简化计算,在取桁架的计算简图时,作如下三个方面的假定:(1)桁架的结点都是光滑的铰结点。 (2)各杆的轴线都是直线并通过铰的中心。 (3)荷载和支座反力都作用在铰结点上。 通常把符合上述假定条件的桁架称为理想桁架。 3.4.1.2 桁架的受力特点 桁架的杆件只在两端受力。因此,桁架中的所有杆件均为二力杆。在杆的截面上只有轴力。 3.4.1.3 桁架的分类 (1)简单桁架:由基础或一个基本铰接三角形开始,逐次增加二元体所组成的几何不变体。(图3-14a) (2)联合桁架:由几个简单桁架联合组成的几何不变的铰接体系。(图3-14b) (3)复杂桁架:不属于前两类的桁架。(图3-14c) 3.4.2 桁架内力计算的方法 桁架结构的内力计算方法主要为:结点法、截面法、联合法 结点法――适用于计算简单桁架。 截面法――适用于计算联合桁架、简单桁架中少数杆件的计算。 联合法――在解决一些复杂的桁架时,单独应用结点法或截面法往往不能够求解结构的内力,这时需要将这两种方法进行联合应用,从而进行解题。 解题的关键是从几何构造分析着手,利用结点单杆、截面单杆的特点,使问题可解。 在具体计算时,规定内力符号以杆件受拉为正,受压为负。结点隔离体上拉力的指向是离开结点,压力指向是指向结点。对于方向已知的内力应该按照实际方向画出,对于方向未知的内力,通常假设为拉力,如果计算结果为负值,则说明此内力为压力。 管桁架施工方案 编制人: 初审人: 审核人: 审批人: 1?编制依据 2.工程概况 3.施工总平面布置及说明 4.施工准备 6.主要工程项目的施工程序和施工方法7管桁架相贯线焊接措施 8.工程施工计划及其说明 9.工程质量目标及保证措施 10.施工测量控制方案 11.施工图节点构造深化设计措施及方案 12.服务及保修 1.编制依 据目录 3. 4. 5. . 10 28 30 31 36 37 3.9 1.1编制原则 1.1.1严格执行国家基本建设程序、法律、标准、规范,发挥我公司的优势, 遵守我公司的质量信誉方针,用我们的承诺和智慧雕塑时代的艺术精品。 1.1.2科学的安排施工顺序,保证施工的连续性和均衡性,使各施工阶段之间互相搭接、衔接紧凑,力求达到快速优质、安全高效的管理。 1.1.3全面贯彻执行公司ISO9001: 2000质量管理体系(《质量手册》《程序文件》坚持质量第一、预防为主的原则。 1.1.4加强工程项目优化配置及动态管理,确保现场双文明工程。(文明施工、文明现场) 1.2编制依据 1.2.1投标文件、相关设计文件 1.2.2工程质量验收标准 《建筑项目工程施工质量验收统一标准》(GB50300—2001) 《钢结构工程施工质量验收标准》(GB50205—2001) 1.1.3技术规范 《建筑钢结构焊接技术规程》(JGJ81- 2002) 《钢结构制作安装施工规程》(YB9254 —1995) 《优质碳素结构钢》GB/ T699—1999 《普通碳素结构钢》GB/T700—1998 《低合金高强度结构钢》GB/T1591 —1994 《结构用无缝管》GB / T8162—1999 《直缝焊管》GB/T13793—1992 管桁架施工方案 编制人: 初审人: 审核人: 审批人: 年月日 目录 1.编制依据..................................................... 2. 工程概况.................................................... 3.施工总平面布置及说明......................................... 4.施工准备..................................................... 6.主要工程项目的施工程序和施工方法............................. 7管桁架相贯线焊接措施.......................................... 8.工程施工计划及其说明......................................... 9.工程质量目标及保证措施....................................... 10.施工测量控制方案............................................ 11.施工图节点构造深化设计措施及方案............................ 12.服务及保修 (41) 1.编制依据 1.1编制原则 1.1.1严格执行国家基本建设程序、法律、标准、规范,发挥我公司的优势,遵守我公司的质量信誉方针,用我们的承诺和智慧雕塑时代的艺术精品。 1.1.2科学的安排施工顺序,保证施工的连续性和均衡性,使各施工阶段之间互相搭接、衔接紧凑,力求达到快速优质、安全高效的管理。 1.1.3全面贯彻执行公司ISO9001:2000质量管理体系(《质量手册》《程序文件》坚持质量第一、预防为主的原则。 1.1.4加强工程项目优化配置及动态管理,确保现场双文明工程。(文明施工、文明现场) 1.2编制依据 1.2.1投标文件、相关设计文件 1.2.2工程质量验收标准 《建筑项目工程施工质量验收统一标准》(GB50300-2001) 实验二:平面桁架的静力分析 一、平面桁架静力分析程序框图 平面桁架静力分析程序名为PTSAP(Plane Truss Structural Analysis Program)。其主要表示符说明如下: TL(20)——算例标题。实型数组,输入参数。 NJ——结点总数。整型变量,输入参数。 N——结构的自由度,即整体刚度矩阵阶数。整型变量,输入参数。 NNE——单元总数。整型变量,输入参数。 NMT——单元类型总数。同类型单元E、A 相同。整型变量,输入参数。 NPJ——结点荷载总数。整型变量,输入参数。 JE——(2,100)——单元两端结点号数组。整型变量,输入参数。 JN(2,100)——结点位移号数组。整型变量,输入参数。 X(100)、Y(100)——结点坐标数组,X(I)、Y(I)分别为I 号结点的x 坐标、y 坐标。 JEA(100)——单元类型信息数组。JEA(I)为第 I 单元的类型号,同类型的单元弹性模量横截面积相同。整型变量,输入参数。 EA(2,25)——各类型单元的物理、几何性质数组。EA(1,I)、EA(2,I)分别为第 I 类型单元的弹性模量、截面面积。输入参数,实型数组。 JPJ(50)——结点荷载的位移号数组。JPJ(I)为与第 I 个结点荷载相应位移分量的位移号。整数数组,输入参数。 PJ(50)——结点荷载的位移号数组。PJ(I)为与第I 个结点荷载的数值。实型数组,输入参数。 M(4)——单元定位数组,及单元两端的位移号数组。整型变量,输入参数。 AK(200,200)——存结构整体刚度矩阵的上半带元素。 AKE(4,4)——存整体坐标系下单元刚度矩阵。 T(4,4)——存坐标变换矩阵或其转值矩阵。 P(200)——结点荷载。解方程后,存结点位移。 FE(4)——存整体坐标系下单元刚度矩阵与单元杆端位移的乘积。 F(4)——先存整体坐标系下的单元杆端位移,后存局部坐标下单元杆端力列阵。c语言计算平面桁架内力计算程序
桁架内力计算
空间桁架结构程序设计(Fortran)
计算结构力学自编平面桁架
平面桁架
简单桁架内力计算
平面桁架计算程序
简单平面桁架受力分析
桁架支撑计算
简单桁架内力的计算方法
桁架计算方法
平面桁架ansys分析
桁架承重架设计计算书
第二节 平面静定桁架的内力计算
桁架有限元程序流程(有限元课程设计)
桁架内力的计算3.4静定平面桁架
桁架施工方案
桁架施工方案
平面桁架静力分析