C语言编写的计算平面桁架的代码

C语言编写的计算平面桁架的代码，注释详细

//Space truss structural analysis program
//
//注:本程序由FORTRAN程序直接改编过来.
//（要注意FORTRAN主子程序是按列传数据而C是按行传数据）
//
/*	（1）整型变量
	NJ---------结点总数；
	NE---------单元总数；
    NS---------支座结点总数；
	NW---------最大半带宽；
    NPJ--------载荷作用的结点总数；
    N----------结点位移分量总数，N=3*NJ；
    M----------单元序号。

    （2）实型变量
    BL----------单元长度l;
    SI----------单元的sin（a）值；
    CO----------单元的cos（a）值；
    RX----------x方向的支座反力；
    RY----------y方向的支座反力；
    SG----------单元横截面上的正应力；
    EAL---------单元常数EA/l；
    S1S2S3----单刚ke中的三个常数。

    （3）整型数组
    JE（2NE）-----单元杆端结点编号数组；
    JS（4NS）-----支座结点约束信息数组；

    （4）实型数组
    X（NJ）Y（NJ） Z（NJ）--分别为结点的x坐标y坐标Z坐标数组；
    E（NE）--------弹性模量数组；
    PJ（4NPJ）----结点载荷数组；
    KE（66）------单元刚度矩阵ke；
    KB（NNW）-----带形存储的整体刚度矩阵K；
    P（N）---------结点载荷向量，后存结点

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....      8569  2004-11-01 11:07  stsap三维桁架有限元分析源程序.cpp

----------- ---------  ---------- -----  ----

                 8569                    1