Elastic deformation of Wind turbine tower under wind pressure. #apdl #ANSYS #ANSYS_CLASSIC

#apdl #ANSYS #ANSYS_CLASSIC
Elastic deformation of wind turbine tower under wind pressure modeling.
apdl used to modeling 3D-model of tower and load. Load static[simplified].

finish
/CLEAR
/BATCH 
/COM,ANSYS RELEASE Release 18.0      BUILD 18.0      UP20161219       08:12:45
/TITLE, Tower
/FILNAME, tower

!/GRA,POWER
!/GST,ON
!/PLO,INFO,3
!/GRO,CURL,ON
!/CPLANE,1  

/PREP7 
*set,W_Press,XXX000
*set,M_Press,XXX000
*set,K,1 !proportional factor
*SET,H,XXXXX/K
*SET,D_0,XXXX
*SET,S_0,XX
*SET,d_X,XXXX
*SET,s_X,XX
*SET,R_0,D_0/2
*SET,R_1,D_1/2
k,,(R_0-S_0),0,0
k,,R_0,0,0
k,,(R_1-S_1),H,0
k,,R_1,H,0
!klist  
!ksel,,,,1,4,,  
A,1,2,4,3,
   
k,,0,0,0,0 
k,,0,H,0
vrotat,1,,,,,,5,6,,,

!* 
ET,1,SOLSH190   !element type
!ET,1,SOLID185    !element type
!* 
UIMP,1,EX, , ,2.1e8, !
UIMP,1,NUXY, , ,0.3 !

!vmesh,1,4
ESIZE,300,0, ! Element size
vmesh,ALL

!*set,up_f,(asel,s,LOC,Z,H)
asel,s,LOC,Y,0
nsel,s,LOC,Y,0,0+1
!*set,down_f,(asel,s,LOC,Z,0)


d, ALL, ux, 0.             ! Displacement at node down_face_ns in x-dir is zero
d, ALL, uy, 0.             ! Displacement at node down_face_ns in y-dir is zero
d, ALL, uz, 0.

asel,s,LOC,Y,H-1,H+1,0
nsla,s,1
f,ALL,FX,500000! W_Press!   !Wind pressure=Forces at node 1 in X direction.
f,ALL,FY,-600000 !M_Press     ! Mass_Pressure
f,ALL,FZ,500000 !W_Press    !Wind pressure on Z

allsel

antype,static    !0=static
solve

antype,1
solve
antype,2
solve
antype,3
solve

save
finish
/post1                   ! Use 'ssum' to compute sums interactively
!etable, Stress, ls, 1    ! Create a table of element stress values
!etable, Volume, volum    ! Create a table of element volumes
PLNSOL, U,SUM, 0,1.0
/show,png
PLDISP,2
/show,png
!* 
/EFACET,1  
PLNSOL, S,EQV, 0,1.0
/show,png
!PLETAB,STRESS,NOAV 

!/GRAPHICS,power
finish