Fortran语言编写的有限元结构程序

算例一计算简图及结果输出 用平面刚架静力计算程序下图结构的内力。各杆EA，EI 相同。 已知：642EA=4.0 10 KN,EI=1.610 KNm 计算简图如下： . （1）输入原始数据 控制参数 3,5,8,7,1,2 （NE,NJ,N,NW,NPJ,NPF） 结点坐标集结点未知量编号 0 . 0 , 0 . 0 , 0 ,0 . 0 , 4 . 0 , 1 , 2 , 30 . 0 , 4 . 0 , 1 , 2 , 44 . 0 , 4 . 0 , 5 , 6 , 74 . 0 , 0 . 0 , 0 , 0 , 8 单元杆端结点编号及单元EA、EI 1,2,4.0E+06,1.6E+043,4,4.0E+06,1.6E+045,4,4.0E+06,1.6E+04 结点荷载 7 . 0 , - 1 5 . 非 结点荷 载 1 . 0 , 2 . 0 , 2 . 0 , -2 . 0 , 1 . 0 , 4 . 0 , - 2 5 . 0 （2）输出结果 NE= 3 NJ= 5 N= 8 NW= 7 NPJ= 1 NPF= 2 NODE X Y XX YY ZZ 1 0.0000 0.0000 0 0 0 2 0.0000 4.0000 1 2 3 3 0.0000 4.0000 1 2 4 4 4.0000 4.0000 5 6 7 5 4.0000 0.0000 0 0 8 ELEMENT NODE-I NODE-J EA EI 1 1 2 0.400000E+07 0.160000E+05 2 3 4 0.400000E+07 0.160000E+05 3 5 4 0.400000E+07 0.160000E+05 CODE PX-PY-PM 7. -15.0000 ELEMENT IND A Q 1. 2. 2.0000 -18.0000 2. 1. 4.0000 -25.0000 NODE U V CETA 1 0.000000E+00 0.000000E+00 0.000000E+00 2 -0.221743E-02 -0.464619E-04 -0.139404E-02 3 -0.221743E-02 -0.464619E-04 0.357876E-02 4 -0.222472E-02 -0.535381E-04 -0.298554E-02 5 0.000000E+00 0.000000E+00 0.658499E-03 ELEMENT N Q M 1 N1= 46.4619 Q1= 10.7119 M1= -6.8477 N2= -46.4619 Q2= 7.2881 M2= 0.0000 2 N1= 7.2881 Q1= 46.4619 M1= 0.0000 N2= -7.2881 Q2= 53.5381 M2= 14.1523 3 N1= 53.5381 Q1= 7.2881 M1= 0.0000 N2= -53.5381 Q2= -7.2881 M2= -29.1523 算例二计算简图及结果输出 用平面刚架静力计算程序下图结构的内力。 已知：桁架单元的抗拉刚度为6EA=2.0 10 KN，平面刚架单元的抗拉刚度为已知：6EA=4.0 10 KN，抗弯刚度为42EI=1.8410 KNm。 计算简图如下： （1）输入原始数据 控制参数 5 , 4 , 8 , 7 , 0（NE,NJ,N,NW,NPJ,NPF） 结点坐标集结点未知量编号 0 . 0 , 0 . 0 , 0 , 04 . ...