3杆系静力分析2ppt课件

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 有关有关 F90 语法语法where构造构造 新特性新特性有关有关 F90 语法语法cycle和和exit语句语句 新特性新特性 用在用在do循环中循环中 cycle 作下一个循环步作下一个循环步 exit 跳出循环，执行跳出循环，执行end do后一后一条语句条语句 . . end if end if end do end do5 .5 . . .end doend do.有关

8、有关 F90 语法语法数组构造函数数组构造函数spread 054301054301054301 spread(A(2,2:),dim=1,ncopies=2) spread(A(2,2:),dim=1,ncopies=2) 332313322212312111A. 32223222.有关有关 F90 语法语法指针指针pointer！即：！即： 一切一切 = 3.14 = 3.14 deallocate (b)deallocate (b)有关有关 F90 语法语法用指针建立动态数组用指针建立动态数组.allocate (Kcol(5)%row(3:5) ! allocate (Kcol(5)%

9、row(3:5) ! 第第5 5列只列只用用3 3至至5 5行行Indx 类型类型 pos1 均布均布 长度长度2 集中集中 位置位置3 .(2)(1)(3)24135i = 6i = 4i = 310 kN4 kN/m4m4mEA= 109 N!=!=subroutine SetMatBand (Kcol, Elem) ! subroutine SetMatBand (Kcol, Elem) ! 接口简单接口简单!=!= type(typ_Kcol),intent(in out) : Kcol(:) type(typ_Kcol),intent(in out) : Kcol(:) ！总刚列总刚

10、列 type(typ_Element),intent(in) : Elem(:) type(typ_Element),intent(in) : Elem(:) ！单元单元 integer(ikind) : minDOF,ELocVec(6) integer(ikind) : minDOF,ELocVec(6) integer(ikind) : Row1(size(Kcol,dim=1) integer(ikind) : Row1(size(Kcol,dim=1) ！Row1Row1为自动数组，子程序终了后自动释放。为自动数组，子程序终了后自动释放。 ！这样做可使接口简单，减少了数组的控制变量。

11、！这样做可使接口简单，减少了数组的控制变量。 integer(ikind) : ie,j,NGlbDOF,NElem integer(ikind) : ie,j,NGlbDOF,NElem NGlbDOF = size(Kcol,dim=1) NGlbDOF = size(Kcol,dim=1) ！使接口简单！使接口简单621,.,j , i IJ)ie(ijKk 0)(0)()ie(j)ie(iJI do ie=1,NElem ! 计算单刚 EK(6,6),ELocVec(6) do j=1,6 ! 对单元逐列集成 JGDOF = ELocVec(j) ! 取出位移码 if (JGDOF =

12、 0) cycle ! 作下一循环步 where (ELocVec 0 . And . ELocVec = JGDOF) ! 位移码非零同时小于第j个部分码对应的位移码 Kcol(JGDOF)%row(ELocVec) = & Kcol(JGDOF)%row(ELocVec) + EK(:,j) ! 集成 end where end do end do4) 变带宽矩阵的分解求解变带宽矩阵的分解求解4-1) LDLT分解法求解分解法求解 Ax = b nn A TLDLA =单位上三角阵对角阵 byLxLDL T TLDLA byL yxLD T ULLDLA T nnnnijiinjn

13、jaaaaaaaaaaaaaa2121222221111211iidjiijll T 1111213231210nnijllllll nnjiiinndldldlddldldldd02223222221113111211111jiiiijldu ijaijl nknkjkkkikkjikijldlula11 iklikiklikik10j , ik )j , i(kikkkjkikijdllamin11)(1111da 1111ad 211112lda 121121adl 222211122dlda 221112222ldad 12a22aiijiikkkjkikijdldlla 11iiik

14、kkjkikijjiddllal 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 nnnnnndydydyx