平面桁架有限元分析及程序设计..pptVIP

§2.6 程序设计 2、先处理法处理边界条件 单元定位向量： x y F ① ② ③ 1 2 3 (1,2) (3,4) (5,6) 后处理： x y F ① ② ③ 1 2 3 (1,2) (0,0) (3,0) 划行划列法 先处理法： 单元定位向量： void force() {int i,j,ie,m; float dx,dy,dz,l,cx,cy,cz,ea,w[7]; for(ie=1; ie=ne; ie++) { i=jm[ie][1]; j=jm[ie][2]; m=jm[ie][0]; w[1]=f[2*i-2]; w[2]=f[2*i-1]; w[3]=f[2*j-2]; w[4]=f[2*j-1]; dx=xy[j][1]-xy[i][1]; dy=xy[j][2]-xy[i][2]; 例题 3、程序说明 下面为一个平面桁架计算程序段，试在左端有编号的程序右面写出其注释。（10分） l=sqrt(dx*dx+dy*dy); cx=dx/l; cy=dy/l; ea=EA[m]/l; dx=w[3]-w[1]; dy=w[4]-w[2]; l=ea*(cx*dx+cy*dy); F[ie]=l; } } ? 1 2 4 5 3 1、运用有限单元法，计算图示桁架：杆件截面积为A，弹性模量为E，结点2作用集中力，结点3给定水平位移b，要求写出整体平衡方程及边界条件处理方法。 作业 1 2 3 4 Y2 X2 l l b 2、运用有限单元法程序,计算程序书第4页作业题，要求打印输入文件、输出文件，并画出各个杆件的轴力图。 3、计算分析题（15分） 平面桁架结构网格如图所示，已知EA=6000 kN，采用填0置1法（0、1置换法）引入支撑条件，试求后处理法引入支撑条件后的满阵存贮的整体刚度矩阵 。结构坐标系下单元刚度矩阵计算公式如下： 式中： 1 3 2 ① ② ③ 3m y x 4m 作业 4、程序说明（10分） 下面为一个平面桁架计算程序段，左端有编号的程序部分有两处错误，请改正，并在没有错误的程序右面写出其注释。（10分） int ekzk(int ie ) { int i1,j1,i,j,i2,j2,ii,jj,ji; for(i1=1;i1=2;i1++) { for(i2=1;i2=2;i2++) { 1 i=2*(i1-1)+i2; 2 ii=2*jm[ie][i1]+i2; for(j1=1; j1=2; j1++) { for(j2=1; j2=2; j2++) { j=2*(j1-1)+j2; 3 jj=2*(jm[ie][j1]-1)+j2; 4 ji=bw+jj-ii+1; 5 if(ji=bw) K[ii-1][ji-1]=K[ii-1][ji-1]+ek[i][j]; } } } } } 作业 此例说明了桁架有限元分析的基本思路和步骤。 * 建立杆端力与杆端位移的关系。 注意杆端力与应力正负号规定的差别，就可以明白杆端力公式中的正负。 * 结点力-此处提杆端力更合适，公式中正负号来自轴力的正负号与杆端力正负号的规定。 * * * * * * * 注意讲解整体刚度矩阵每个元素的特点。特点1：脚标；特点2：与哪几个杆件连接；特点3：对称；特点4：行列式的值为零。 * 讲清楚为什么可以划去行和列（有几个约束就多几个方程）。 * * * 说明半带宽公式的来历（所有单元的最大结点编号差，就意味着这个单元的元素离对角线位置的距离，2倍表示每个结点有两个自由度） 下半带存储后，元素距对角线多远，决定了这个元素在新矩阵中列的位置。 原来讲过的方法均属后处理方法，后处理法相当于划行划列法。 * * * * 单元①整体坐标系下的刚