/CLEAR,NOSTART

!menu pre process
/PREP7  

L1=1.08
L2=0.8

H1=15
H2=25

R1=0.05
EP1=0.01

Pa=100 ! poid antenne

! creation des points

K,1,L1,0,0,  
K,2,-0.5*L1,sqrt(3)/2*L1,0,   
K,3,-0.5*L1,-sqrt(3)/2*L1,0, 

K,4,L2,0,H1,  
K,5,-0.5*L2,sqrt(3)/2*L2,H1,   
K,6,-0.5*L2,-sqrt(3)/2*L2,H1,  

K,7,L2,0,H2,  
K,8,-0.5*L2,sqrt(3)/2*L2,H2,   
K,9,-0.5*L2,-sqrt(3)/2*L2,H2,

! creation des lignes
*DO,i,1,3,1
LSTR,       i,       i+3  
LSTR,       i+3,       i+6 
*ENDDO

! decoupage ligne montants
*DO,i,1,5,2
LDIV,i, , ,3,0
*ENDDO
*DO,i,2,6,2
LDIV,i, , ,2,0
*ENDDO

! decoupage lignes transverses
*DO,i,1,15,1
LDIV,i,0.44, ,,0
*ENDDO


LSEL,S,,,16,20,2
LSEL,A,,,22,26,2
LSEL,A,,,23,27,2
LSEL,A,,,17,21,2
LSEL,A,,,28,30,1

LDIV,ALL,0.9, ,2,0 

LSEL,S,,,ALL

LSTR,      23,      19  
LSTR,      19,      21  
LSTR,      21,      23  
 
LSTR,      38,      34  
LSTR,      34,      36  
LSTR,      36,      38  

LSTR,      29,      25  
LSTR,      25,      27  
LSTR,      27,      29  
 
LSTR,      44,      40  
LSTR,      40,      42  
LSTR,      42,      44  

LSTR,      30,      26  
LSTR,      26,      28  
LSTR,      28,      30  

LSTR,      45,      41  
LSTR,      41,      43  
LSTR,      43,      45  

LSTR,      24,      20  
LSTR,      20,      22  
LSTR,      22,      24  

LSTR,      39,      35  
LSTR,      35,      37  
LSTR,      37,      39  

LSTR,      33,      31  
LSTR,      32,      31  
LSTR,      33,      32  

LSTR,      48,      46  
LSTR,      46,      47  
LSTR,      47,      48  

! creation des lignes transverses diagonales
LSTR,       1,      21  
LSTR,      21,      34  
LSTR,      34,      27  
LSTR,      27,      40   
LSTR,      40,      28  
LSTR,      28,      41  
LSTR,      41,      22  
LSTR,      22,      35   
LSTR,      35,      32  
LSTR,      32,      46   
LSTR,       3,      19  
LSTR,      19,      38   
LSTR,      38,      25   
LSTR,      25,      44   
LSTR,      44,      26   
LSTR,      26,      45  
LSTR,      45,      20  
LSTR,      20,      39   
LSTR,      39,      31  
LSTR,      31,      48  
LSTR,       2,      23  
LSTR,      23,      36  
LSTR,      36,      29   
LSTR,      29,      42  
LSTR,      42,      30  
LSTR,      30,      43  
LSTR,      43,      24   
LSTR,      24,      37   
LSTR,      37,      33   
LSTR,      33,      47






! options affichage
GPLOT
/PNUM,KP,0  
/PNUM,LINE,0

!-----------------!
! MODELISATION EF
!-----------------!

! Type de modele
ET,1,BEAM188

! sections
SECTYPE,   1, BEAM, CTUBE, pieds, 0 
SECOFFSET, CENT 
SECDATA,R1-EP1,R1,0,0,0,0,0,0,0,0,0,0   

SECTYPE,   2, BEAM, CTUBE, traverses, 0 
SECOFFSET, CENT 
SECDATA,0.02,0.025,0,0,0,0,0,0,0,0,0,0

! materiau
MPTEMP,,,,,,,,  
MPTEMP,1,0  

MPDATA,EX,1,,2e11   
MPDATA,PRXY,1,,0.3  
MPDATA,DENS,1,,7800 


poidAnt=100
MPDATA,EX,2,,2e11   
MPDATA,PRXY,2,,0.3  
MPDATA,DENS,2,,7800+poidAnt/(2.5*3.14*0.1*0.1/4)




!Attribution de la modelisation
LSEL,S,,,1,45,1
LSEL,U,,,28,30,1
LATT,1, ,1, , , ,1 !section pieds 

LSEL,S,,,28,30,1
LATT,2, ,1, , , ,1 !section pieds + antenne

NLIGNES=105
LSEL,S,,,46,NLIGNES,1
LATT,1, ,1, , , ,2 !section traverse

LSEL,S,,,ALL

! taille des mailles
LESIZE,ALL,0.5, , , ,1, , ,1,   

! ceation du maillage
LMESH,ALL

! affichage des sections
/ESHAPE,1

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

! Gravite
ACEL,0,0,-10,

! FORCES vent
V=10
rho=1.2
Sa=0.5*0.1
Fvent=0.5*rho*Sa*V*V

Sant=0.5*0.3
Fant=0.5*rho*Sant*V*V

!NNODES=408
!*DO,i,1,NNODES,1
!  *GET,Znode,NODE,i,LOC,Z
!  *IF,Znode,GE,20,THEN
!    F,i,FX,Fant
!  *ELSE
!    F,i,FX,Fvent
!  *ENDIF
!*ENDDO

VREF=10
ZREF=10
ALPHA=0.16

NNODES=408
*DIM,VFVent,ARRAY,NNODES,1
NEL=465
*DO,i,1,NEL,1
  *GET,nodeI,ELEM,i,NODE,1
  *GET,xI,NODE,nodeI,LOC,X
  *GET,zI,NODE,nodeI,LOC,Z
  *GET,nodeJ,ELEM,i,NODE,2
  *GET,xJ,NODE,nodeJ,LOC,X
  *GET,zJ,NODE,nodeJ,LOC,Z
  ZELT=(zI+zJ)/2
  V=VREF*(ZELT/ZREF)**ALPHA
  *SET,long,abs(XJ-XI)
  FVel=0.5*rho*Sa*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






!--
! UTILISATION DU SOLVEUR ("atelier" resolution) 
!--
/SOL

! calcul statique
ANTYPE,0

! resolution
SOLVE

!--
! Post traitement ("atelier" post traitement)
!--
/POST1

! affichage deformee avec option
PLDISP,1