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):Work1,Integer:i,j,k,
2、实型和整型工作变量实型和整型工作变量 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 xyyxxFyFxFyFydxdydxd yxyxcossinsincosFFFF yxyxcossin-s
3、incosdddd 10000CSSC 00sincosTSC xyyxxFyFxFyFydxdydxd T1TT )(TeeeePFTPF eedTd xyyxxFyFxFyFydxdydxd eeeeeeeedTkTdkTPFTPFT )(eeeedkPF TkTkeeT eeeedkPF 1234 2F 1F2 2dP 12F 21F21orj 各各杆杆交交ijiFP1234 314d32223d21122d111dFPFFPFFPFP TT4dT3dT2dT1ddPPPPP TT32T31T22T21T12T11FFFFFFF 2 2dP 12F 21F FIIIIIIP 000000
4、000000000000d 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)(3,4,5)(6,7,8)(6,7,9)(10,11,12)123456
5、(0,0,1)(0,0,2)(3,4,5)(6,7,8)(6,7,9)(10,11,12)123456(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)1234 56(0,0)(0,0)(1,2,3)(4,5)(6,7,8)12345 555453525145444342413534333231252423222115141312113k87654543215432187654(0,0
6、,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)1234587600 4543214 PPPPPP54321 nininnnjninninijiiiinjinjiRRRRKKKKKKKKKKKKKKKKKKKK2121212122222211111211 aKRaaKRaKRKKKKKKKKKKKKniniininnnjnnnjnj2211212122222111121100010000 xy nnnininiiiniKKKKKKKKK111111136101453606672 1)有关有
7、关 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 TLDLA byL yxLD T ULLDLA T nnnnijiinjnjaaaaaaaaaaaaaa
8、2121222221111211 iidjiijll T 1111213231210nnijllllll nnjiiinndldldlddldldldd02223222221113111211111jiiiijldu ijaijl nknkjkkkikkjikijldlula11 iklikiklikik10j,ik )j,i(kikkkjkikijdllamin11)(1111da 1111ad 211112lda 121121adl 222211122dlda 221112222ldad 12a22aiijiikkkjkikijdldlla 11iiikkkjkikijjiddllal 1
9、1121 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 nndddD11122111iiinikkkiidyxlx 1 nikkkiiiiixldyx1nni)n(i)n(ixlyy 1iiinidyy )1(nx)1(nnnyx)1(11 nnnyxnx1 nx 21 1 )(111)1()(,n,njyxj,iixlyyjjjjjijiji PdkF PdTk dTk dkMNQ1F2F3F
