1-A袁驷f90程序内容.ppt

2006－7－8 哈尔滨工业大学 王焕定 * * 清华袁驷教授有关的程序内容 定义导出类型 导出类型（自定义类型）——结点 type :: typ_Joint real :: x,y ！坐标 integer :: GDOF(3)！整体位移码 end type typ_Joint 1) 有关 F90 语法 摘自清华袁驷教授教案 导出类型 新特性 清华袁驷教授有关的程序内容 用结点导出类型作为成员导出单元类型: type :: typ_Element integer :: JointNo(2) ！结点编号 type(typ_Joint) :: Node(2) ！结点 integer :: GlbDOF(6) ! 定位向量 real :: EA,EI end type typ_Element type (typ_Element) :: Elem(5) ！定义5个单元类型 … ! 对单元i的端点j的x,y,GDOF(1:3)的赋值 Elem(i)%Node(j) = Joint(Elem(i)%JointNo(j)) … 由导出类型定义新类型 由导出类型定义变量 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 （待续） 有关 F90 语法 数组运算与赋值： matmul(matrix_a,matrix_b) 矩阵相乘 如：locEDisp = matmul(T,glbEDisp) transpose(matrix) 矩阵转置 如：glbEDisp = matmul(transpose(T),locEDisp) size(array,dim) 求数组第dim维的长度 dim为可选变元: size(a,dim=2) 若array为一维时，可不用dim。 sum(array,dim,mask) 数组元素求和 dim,mask为可选变元；mask=条件表达式 sum(a(1:10)) 对a的1到10元素求和 sum(a(1:10),mask=a0) 对a(1:10)中大于0的元素求和 （续） 有关 F90 语法 where结构 新特性 例 where (C0) C = 0 A = B*D end where where (C0) A = B end where 定义 where (数组关系表达式) 数组赋值语句 … elsewhere 数组赋值语句 ... end where 规则: 1)同形数组; 2)不许嵌套; 3)最多二分叉 有关 F90 语法 cycle和exit语句 新特性 用在do循环中 cycle —— 作下一个循环步 exit —— 跳出循环，执行end do后一条语句 等效例 do ... if(.not.cond1)then ... if(cond2) goto 5 ... end if end do 5 ... 用法 do ... if (cond1) cycle ... if (cond2) exit ... end do ... 有关 F90 语法 数组构造函数spread 语法 spread(数组名,dim,ncopies) 将数组沿dim维方向复制ncopies形成新数组 dim,ncopies — 整型、位置变元、关键字变元 （若按位置引用，可略关键字） 例: （仅限一维数组） 1） spread(one,d