杆系静力分析阅读.pptx

会计学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(:)!该列的行元素
endtype
...
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(NJLoad>0)read(5,*)(JLoad(i),i=1,NJLoad)
if(NELoad>0)read(5,*)(ELoad(i),i=1,NELoad)