user_post.log

fortran编的有限元平面问题分析程序！

!-----------------------------user_post.log----------begin----------------------
!以下为2D问题平面三节点(四节点)单元上用户定义物理量的信息传递过程
!用户提供的信息数据见FOR_POST.DAT文件
!FOR_POST.DAT文件的格式如下
!PART1(共1行两个数，格式2f9.4):节点总数INFO(1)、单元总数INFO(2)
! PART2(基于节点的信息，共INFO(1)行，每行7个数，格式7f9.4):
!节点x坐标、节点y坐标、节点物理量1(x方向位移)、节点物理量2(y方向位移)、节点物理量3(x方向应力)、
!节点物理量4(y方向应力)、节点物理量5(剪应力)
! PART3(基于单元的信息，共INFO(2)行，每行7个数，格式7f9.4):
!单元节点1、单元节点2、单元节点3、单元节点4（对于三角形单元，第4个节点号与第3个相同）、
!单元物理量1（x应力）、单元物理量2（y应力）、单元物理量3（z应力）
!--------------------------------------
/PREP7
ET,1,PLANE42 
DOF,rotx               !define new DOF at node which should be ploted           
DOF,roty
DOF,rotz
*dim,INFO,,2    !define an array INFO which refers to:  INFO(1): number of node, INFO(2): number of element    
*vread,INFO(1),FOR_POST,dat
(2f9.4)
*dim,ND_INFO,,INFO(1),7         !define an array ND_INFO which refers to the information of nodes
*vread,ND_INFO(1,1),FOR_POST,dat,,JIK,7,INFO(1),,1
(7f9.4)
*dim,ELE_INFO,,INFO(2),7        !define an array ele_info which refers to the information of elements
*vread,ELE_INFO(1,1),FOR_POST,dat,,JIK,7,INFO(2),,INFO(1)+1
(7f9.4)
*do,i,1,INFO(1)
N,i,ND_INFO(i,1),ND_INFO(i,2)      !creat the nodes by first two columns of nd_info
*enddo
*do,i,1,INFO(2)
E,ELE_INFO(i,1),ELE_INFO(i,2),ELE_INFO(i,3),ELE_INFO(i,4)   !creat elements by first four columns of ele_info
*enddo
/post1
*do,i,1,INFO(1)
dnsol,i,u,x,ND_INFO(i,3)
dnsol,i,u,y,ND_INFO(i,4)
dnsol,i,rot,x,ND_INFO(i,5)
dnsol,i,rot,y,ND_INFO(i,6)
dnsol,i,rot,z,ND_INFO(i,7)
*enddo
*do,i,1,INFO(2)
*do,j,1,3
num=ELE_INFO(i,j)
desol,i,num,s,x,ELE_INFO(i,5)
desol,i,num,s,y,ELE_INFO(i,6)
desol,i,num,s,z,ELE_INFO(i,7)
*enddo
*enddo!-----------------------------user_post.log----------end----------------------