jingli.txt

生成网架结构的程序

finish
/clear
M=15
N=15
S1=3
S2=3
H=3
N1=(M+1)*(N+1)
/title,Zheng jiao zheng fang
/prep7
et,1,link8
mp,ex,1,2.06e11
r,1,1.717e-03,0
r,2,1.308e-03,0
r,3,7.900e-04,0
r,4,7.900e-04,0
!形成节点
*DO,J,1,N+1
  *DO,I,1,M+1
    K=(J-1)*(M+1)+I
    N,K,(I-1)*S1,(J-1)*S2,0
  *ENDDO
*ENDDO
*DO,J,1,N+1
  *DO,I,1,M+1
    K=N1+(J-1)*(M+1)+I
    N,K,(I-1)*S1,(J-1)*S2,-H
  *ENDDO
*ENDDO
!形成上弦杆
REAL,1
*DO,J,1,N+1
  *DO,I,1,M
    E1=(J-1)*(M+1)+I
    E2=E1+1
    E,E1,E2
  *ENDDO
*ENDDO
*DO,J,1,N
  *DO,I,1,M+1
    E1=(J-1)*(M+1)+I
    E2=J*(M+1)+I
    E,E1,E2
  *ENDDO
*ENDDO
!形成下弦杆
REAL,2
*DO,J,1,N+1
  *DO,I,1,M
    E1=N1+(J-1)*(M+1)+I
    E2=E1+1
    E,E1,E2
  *ENDDO
*ENDDO
*DO,J,1,N
  *DO,I,1,M+1
    E1=N1+(J-1)*(M+1)+I
    E2=N1+J*(M+1)+I
    E,E1,E2
  *ENDDO
*ENDDO
!形成竖杆
REAL,3
*DO,J,1,N+1
  *DO,I,1,M+1
    E1=(J-1)*(M+1)+I
    E2=N1+(J-1)*(M+1)+I
    E,E1,E2
  *ENDDO
*ENDDO
!形成斜杆
REAL,4
*DO,J,1,N+1
  *DO,I,1,NINT(M/2)
    E1=(J-1)*(M+1)+I
    E2=N1+(J-1)*(M+1)+I+1
    E,E1,E2
  *ENDDO
*ENDDO
*DO,J,1,N+1
  *DO,I,NINT((M+1)/2)+1,M+1
    E1=(J-1)*(M+1)+I
    E2=N1+(J-1)*(M+1)+I-1
    E,E1,E2
  *ENDDO
*ENDDO
*DO,I,1,M+1
  *DO,J,1,NINT(N/2)
    E1=(J-1)*(M+1)+I
    E2=N1+J*(M+1)+I
    E,E1,E2
  *ENDDO
*ENDDO
*DO,I,1,M+1 
  *DO,J,NINT((N+1)/2)+1,N+1 
    E1=(J-1)*(M+1)+I
    E2=N1+(J-2)*(M+1)+I
    E,E1,E2
  *ENDDO
*ENDDO
FINISH
/SOLU
*DO,J,1,N+1,N
  *DO,I,1,M+1
    K=(J-1)*(M+1)+I
    D,K,UX,0
    D,K,UZ,0
  *ENDDO
*ENDDO
*DO,J,1,N+1
  *DO,I,1,M+1,M
    K=(J-1)*(M+1)+I
    D,K,UY,0
    D,K,UZ,0
  *ENDDO
*ENDDO
*DO,I,1,N1
  F,I,FZ,-1
*ENDDO