!------------------------------------!
!- Modelisation du pylone en poutre -!
!------------------------------------!

/SOL

!--------------------!
!- Cas de chargement
!--------------------!

!- deplacement impose
DK,1, ,0, ,0,ALL, , , , , , 
DK,2, ,0, ,0,ALL, , , , , , 
DK,3, ,0, ,0,ALL, , , , , , 

!- Gravite
ACEL,0,0,10,

!- FORCES vent
rho=1.2
Da=0.1
VREF=40
ZREF=20
ALPHA=0.16

NNODES=453
*DIM,VFVent,ARRAY,NNODES,1

NEL=516
*DO,i,1,NEL,1
  *GET,nodeI,ELEM,i,NODE,1
  *GET,xI,NODE,nodeI,LOC,X
  *GET,yI,NODE,nodeI,LOC,Y
  *GET,zI,NODE,nodeI,LOC,Z
  *GET,nodeJ,ELEM,i,NODE,2
  *GET,xJ,NODE,nodeJ,LOC,X
  *GET,yJ,NODE,nodeJ,LOC,Y
  *GET,zJ,NODE,nodeJ,LOC,Z
  ZELT=(zI+zJ)/2
  V=VREF*(ZELT/ZREF)**ALPHA
  *SET,long,sqrt((zJ-zI)**2+(yJ-yI)**2)
  FVel=0.5*rho*(Da*long)*(V**2)
  VFVent(nodeI)=VFVent(nodeI)+FVel/2
  VFVent(nodeJ)=VFVent(nodeJ)+FVel/2
*ENDDO
*DO,i,1,NNODES,1
  F,i,FX,VFVent(i)
*ENDDO