FINISH
/CLEAR
/COM,  Structural
!
/FILNAME,geom_simple,1
!
! Ce programme pour tester parametres visco
! sur une eprouvette de geometrie simple.
!
FINISH
/Prep7
ET,1,shell281
MPTEMP,1,0
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Proprietes du materiau 																						 !
! Partie elastique																								 !
!																												 !
MPDATA,EX,1,,??????????, !!!! Entrer la valeur de E0															 !
MPDATA,PRXY,1,,???????????, !!!! Entrer la valeur de nu0													     !
!   																											 !
!                                                                                                                !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!													     !
! Module de relaxation en compressibilite : alpha_k_ansys														 !
!																												 !
TB,PRONY,1,1,??,BULK     !!!! Ne pas oublier d'entrer le nombre de paires				                		 !
TBTEMP,0																										 !
TBDATA,,alpha_1_k_ans,tau_1_k_ans,alpha_2_k_ans,tau_2_k_ans,alpha_3_k_ans,tau_1_3_ans,  						 !
! Rajouter autant de paires que necessaire, 3 paires par lignes maximimum                                        !
TBDATA,,alpha_4_k_ans,tau_4_k_ans,alpha_5_k_ans,tau_5_k_ans,              			            		         !
!                                                                                                                !
!                                                                                                                !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!														 !
! Module de relaxation en cisaillement	: alpha_mu_ansys														 !
!																												 !
TB,PRONY,1,1,??,SHEAR    !!!! Ne pas oublier d'entrer le nombre de paires					                	 !
TBTEMP,0																										 !																										 !
TBDATA,,alpha_1_mu_ans,tau_1_mu_ans,alpha_2_mu_ans,tau_2_mu_ans,alpha_3_mu_ans,tau_1_3_ans,  					 !
! Rajouter autant de paires que necessaire, 3 paires par lignes maximimum                                        !
TBDATA,,alpha_4_mu_ans,tau_4_mu_ans,alpha_5_mu_ans,tau_5_mu_ans,                                                 !
!                                                                                                                !
!																												 !
!fin des donnees pour les proprietes																			 !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Charge 283 N <=> 5 MPa
LT=165-25*2
LG=57
RD=76
WT=19
WG=12.6

SECTYPE,1,shell
secdata, 4.49,1,0,5
BLC4, -WG/2, -LG/2, WG, LG,
BLC4, -WT/2, LG/2, WT, (LT-LG)/2,
BLC4, -WT/2, -LT/2, WT, (LT-LG)/2,
CYL4,(WG/2+RD),LG/2,76
CYL4,-(WG/2+RD),LG/2,76
CYL4,-(WG/2+RD),-LG/2,76
CYL4,(WG/2+RD),-LG/2,76
FLST,2,3,5,ORDE,2
FITEM,2,1
FITEM,2,-3
FLST,3,4,5,ORDE,2
FITEM,3,4
FITEM,3,-7
ASBA,P51X,P51X
AGLUE,ALL
NUMCMP,ALL

smrtsize,1
AATT,1,,1,,1
AMESH,1
AMESH,2,3,1

nsel,S,loc,y,-LT/2
D,all,ALL,0
nsel,all
nsel,S,loc,y,LT/2
CM,charge,NODE
*get,nnode,node,,count,
nsel,all

FINISH
/SOLU
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
!NSUBST,20,30,
NSUBST,50,2000,
OUTRES,ERASE
OUTRES,ALL,1
KBC,0 !0 rampe 1 step
TIME,5.5 !temps fin du LS
F, charge, FY, 283/nnode,      ! Vérifier force pour obtenir 5MPA max.
lswrite,1

!NSUBST,20,700,
NSUBST,50,2000,
OUTRES,ERASE
OUTRES,ALL,1
KBC,1
F, charge, FY, 283/nnode,
TIME,905 !temps fin du LS
lswrite,2

!NSUBST,20,30,
NSUBST,50,2000,
OUTRES,ERASE
OUTRES,ALL,1
KBC,0
F, charge, FY, 0/nnode,
TIME,910.5 !temps fin du LS
lswrite,3

!NSUBST,20,700,
NSUBST,50,2000,
OUTRES,ERASE
OUTRES,ALL,1
KBC,1
F, charge, FY, 0/nnode,
TIME,4000 !temps fin du LS
lswrite,4

lssolve,1,4,1

FINISH
/Post26

ANSOL,2,895,EPEL,X,strainx
ANSOL,3,895,EPEL,Y,strainy
ANSOL,4,895,EPEL,XY,strainxy

PLVAR,3,

LINES,100
PRVAR, 2, 3, 4