张系斌结构力学简明教程讲义3杆系静力分析.pptxVIP

3.4 杆系结构单元分析子程序;3.4 杆系结构单元分析子程序;4.1 杆系结构整体分析;4.1 杆系结构整体分析;4.1 杆系结构整体分析;4.1 杆系结构整体分析;4.1 杆系结构整体分析;;4.1 杆系结构整体分析;;;4.1 杆系结构整体分析;①;定位向量; 根据单元局部位移码和定位向量的对应关系用定 位向量位移码送元素，定位向量元素为零时不送。;4.1 杆系结构整体分析; 5) 综合荷载由两部分组成，因此首先要将直接作 用结点的荷载按结点位移码送入，如果还有单 元等效荷载，再按定位向量集装、累加。;4.1 杆系结构整体分析;;4.1 杆系结构整体分析;最大 半带宽;2.5 杆系结构整体分析;;real :: A(5),B(5,10),C(5) B=0.0 !对B清零 A=1.0 !对A赋1: A(i)=1.0, i=1,5 C=A+2 !数组与标量运算: A(1:5)+(/ 2,2,2,2,2 /) A=C+A !数组与数组运算（同形） C=sqrt(A) !数组的函数运算: C(i)=sqrt(A(i), i=1,5 数组内部函数： dot_product(vector_a,vector_b) !点积 如： dot_product((/1,2,3/),(/2,3,4/)) 的值为20 （待续）;;;;;;;;;;指针与allocatable数组的区别 具备allocatable数组的所有功能 还可以用在导出类型中，例如整体刚度矩阵的变带宽存储： type :: typ_Kcol !整体刚度矩阵K的列 real(8),pointer :: row(:) !该列的行元素 end type ... type (typ_Kcol),allocatable :: Kcol(:) … allocate (Kcol(NGlbDOF)) ! 分配了NGlbDOF列 ... allocate (Kcol(5)%row(3:5)) ! 第5列只用3至5行;;;read(5,*) NElem,NJoint,NGlbDOF,NJLoad,NELoad allocate (Joint(NJoint)) allocate (Elem(NElem)) allocate (JLoad(NJLoad)) allocate (ELoad(NELoad)) ... read(5,*) (Joint(i),i=1,NJoint) read(5,*) (Elem(ie)%JointNo,Elem(ie)%EA, ELem(ie)%EI,ie=1,NElem) if (NJLoad0) read(5,*) (JLoad(i),i=1,NJLoad) if (NELoad0) read(5,*) (ELoad(i),i=1,NELoad);;;;;;;;;;;;;;;;;;