!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!PRE-PROCESSAMENTO
/PREP7
/title,Viga - Geracao Automatica em Volumes
!
! Material
! 
mp,ex,1,210e9       * E = 210GPa
mp,nuxy,1,0.3      * Coeficiente de Poisson
!
! Tipo Elemento
!
et,1,SOLID45         * Elemento Solid45: 3D Isopara. Solid  
!
!Parametros definindo a secao em metros
a=1e-2
b=4.5e-2
c=1.0e-2
d=8.0e-2
e=1.0e-2
f=2.5e-2

!Comprimento da viga em m
len=4

!Valor da carga distribuida em N/m^2
q=100.0

!Posicoes inicial e final da carga distribuida
xi=0
xf=len

!Secao do engaste (x=0)
cc=0

! Keypoints
!
k, 1,   0,  0         * base
k, 2,   0,  a
k, 3,  2*b+c,  0
k, 4,  2*b+c ,  a
k, 5, b,  a
k, 6, b+c,  a
k, 7, b,  0
k, 8, b+c,  0
k, 9, b,  d+a         * alma
k,10, b+c,  d+a
k,11, b, d+a+e         * topo
k,12, b+c, d+a+e
k,13,   b-f,  d+a
k,14,   b+c+f,  d+a
k,15,   b-f, d+a+e
k,16,   b+c+f, d+a+e
kplot
! 
! Linha para o Vdrag (numero 1)
!
k,17,0,0,len  * kpoint
l,1,17       * linha (no. 1)
/view,1,1,1,1         * Observador em (1,1,1)
kplot
! 
! Areas
!
a,1,2,5,7      * base
a,5,6,8,7
a,6,4,3,8
a,5,9,10,6     * alma
a,13,15,11,9   * topo
a,11,12,10,9
a,12,16,14,10
aplot
!
! Vdrag (gera volume por extrusao de areas)
!
vdrag,all,,,,,,1   * areas 1-6 ao longo da linha 1
!
! Ajusta o modelo gerado
!
numm,all              * executa um merge das entidades
!
! Controle/Geracao Automatica
!
lesize,1,,,20       *numero de divisoes ao longo do comprimento
esize,0.02          *tamanho medio do elemento igual a len/10
eshape,2            *geracao com apenas cubos
vmesh,all           *gera malha
eplot               *grafico da malha gerada
!
save         *Salva dados
fini         *Abandona prep7

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!SOLVER
/solu
!
! Restricoes
!
nsel,s,loc,z,cc      *seleciona todos os noz com z=0
d,all,all,0         *viga engastada numa extremidade
nall
! 
! Carregamento
!
! (Carga distribuida de q kgf/cm^2 no topo)
! (areas: 26, 30 e 33)
nsel,s,loc,y,a+d+e
nsel,r,loc,z,xi,xf
sf,all,press,q         * pressao nas areas selecionadas
nall              * seleciona todas as areas
!
! Transporta B.C. para o mesh
!
sbctran
!
! Visualiza
!
/pbc,all,1
/psf,press,2
eplot
!
!Resolve modelo
solve
fini                **Abandona SOLVER

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Pos-Processamento
!
/post1
set                **le resultados
pldisp,1           **deformada
plnsol,s,eqv       **tensao equivalente de von Mises
plnsol,s,x        **tensao normal em x
prnsol,s,comp      **imprime componentes do tensor de tensoes para os nos
prnsol,s,prin      **imprime tensoes principais e de von Mises para os nos
fini               **abandona pos-processador