xuefeng03 发表于 2008-11-1 15:44

问一个有限元程序里的问题

是别人的程序,一个求悬臂板响应的有限元程序,坐标系建立如图所示

在组装总矩阵的时候,有这样一段我不太明白
       forI=1:4
          if(YE>=0)
      Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))=Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))+Kuv((I-1)*5+1:I*5);
         else
      Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))=Ksu(5*(ELE(M1,I)-1)+1:5*ELE(M1,I))-Kuv((I-1)*5+1:I*5);
          end
       end
其中YE为节点的坐标,Ksu为整体的节点载荷,Kuv为单元节点载荷,我不明白的就是为什么YE<0时组装单元节点载荷时用减号

xuefeng03 发表于 2008-11-1 20:59

自己给顶起来
页: [1]
查看完整版本: 问一个有限元程序里的问题