入口整型参数

1、Subroutine Elem_Stiff() 说明说明 Stiff=0.0 ! 单元刚度清零单元刚度清零 Select Case (Type) Case (1) 平面杆系结构单元平面杆系结构单元 Case (2) 空间杆系结构单元空间杆系结构单元 Case Default 出错信息出错信息 End SelectEnd Subroutine Elem_Stiff Integer,Intent (in) : 入口整型参数入口整型参数 Real(8),Intent(in) : 入口实型参数入口实型参数 Real(8),Intent(out) : 出口实型参数出口实型参数 Real(8) : Wor

2、k1, Integer : i,j,k, 实型和整型工作变量实型和整型工作变量 Select Case (Plane) Case (1) 平面桁架元素赋值平面桁架元素赋值 Case (2) 平面梁柱元素赋值平面梁柱元素赋值 Case (3) Case Default 出错信息出错信息 End Select Select Case (Space) Case (1) 空间桁架元素赋值空间桁架元素赋值 Case (2) 空间梁柱元素赋值空间梁柱元素赋值 Case (3) 交叉梁元素赋值交叉梁元素赋值 Case Default 出错信息出错信息 End Select xyyxxFyFxFyFydxdy

3、dxd yxyxcossinsincosFFFF yxyxcossin-sincosdddd 10000CSSC 00sincosTSC xyyxxFyFxFyFydxdydxd T1TT )(TeeeePFTPF eedTd xyyxxFyFxFyFydxdydxd eeeeeeeedTkTdkTPFTPFT )( eeeedkPF TkTkeeT eeeedkPF 1234 2F 1F2 2dP 12F 21F21orj 各各杆杆交交ijiFP1234 314d32223d21122d111dFPFFPFFPFP TT4dT3dT2dT1ddPPPPP TT32T31T22T21T12T1

4、1FFFFFFF 2 2dP 12F 21F FIIIIIIP 000000000000000000d FAP d123 000000000000000000IIIIIIA1234 TT4T1 TT32T31T22T21T12T11dddddd TA 1234 31kkdiagK eqdPKF;FAP TT3T2T1eqPPPP )(eqdPKAP eqeqPAP 记记 TeqdAKAPP KAKA;RPP Teqd记记 KR 则则 31TTiiiiAkAAKA3211234 000000000000000000IIIIIIA iikkkkk 22211211rssr(0,0,1)(0,0,2

5、)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0,1)(0,0,2)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0,1)(0,0,2)(3,4,5)(6,7,8)(6,7,9)(10,11,12)1234 56(0,0)(0,0)(1,2,3)(4,5)(6,7,8)12345(0,0,1)(0,0,2)(3,4,5)(6,7,8) (6,7,9)(10,11,12)123456(0,0)(0,0)(1,2,3)(4,5)(6,7,8)12345 555453525145444342413534333231252423

6、222115141312113k87654543215432187654(0,0,1)(0,0,2)(3,4,5)(6,7,8) (6,7,9)(10,11,12)123456(0,0)(0,0)(1,2,3)(4,5)(6,7,8)1234587600 4543214 PPPPPP54321 nininnnjninninijiiiinjinjiRRRRKKKKKKKKKKKKKKKKKKKK2121212122222211111211 aKRaaKRaKRKKKKKKKKKKKKniniininnnjnnnjnj2211212122222111121100010000 xy nnninini

7、iiniKKKKKKKKK111111136101453606672 1) 有关有关 F90 语法语法导出类型 新特性新特性 有关有关 F90 语法语法 eedTd 有关有关 有关有关 有关有关 054301054301054301 332313322212312111A. 32223222.有关有关 有关有关 Indx (2)(1)(3)24135i = 6i = 4i = 310 kN4 kN/m4m4mEA= 109 N!=!=!=!= 621,.,j , i IJ)ie(ijKk 0)(0)()ie(j)ie(iJI 4) nn A TLDLA =单位上三角阵对角阵 byLxLDL T

8、 TLDLA byL yxLD T ULLDLA T nnnnijiinjnjaaaaaaaaaaaaaa2121222221111211 iidjiijll T 1111213231210nnijllllll nnjiiinndldldlddldldldd02223222221113111211111jiiiijldu ijaijl nknkjkkkikkjikijldlula11 iklikiklikik10j , ik )j , i(kikkkjkikijdllamin11)(1111da 1111ad 211112lda 121121adl 222211122dlda 22111222

9、2ldad 12a22aiijiikkkjkikijdldlla 11iiikkkjkikijjiddllal 11121 j,ijjjkkkjkjjddla 112n.,j32 112jkkkjkjjjjdladiiiikkkjkikijjiddllKl 11121 j,ii11 ik 121jikkkjkjjjjdlKdii1ji1)max(111jii ,ii n.,j32 776657565546544436534333275242322241311100000KKKKdKldKlldKllldlld nnnnijbbbyyyllllll2121213231211111022121byyl iikikikbyyl 11kikikiiylby 11n,i21 12122ylby 11by 2bkiikikiiylby 111i iyiyibib nnnnnndydydyxxxlllll222111212321312111110nnnndyx 1nx yDxL1T nndddD11122111iiinikkkiidy