Commit 8332fa34 authored by mozul's avatar mozul
Browse files

Squashed 'examples/' changes from 104cc2e..791c9d3


git-subtree-dir: examples
git-subtree-split: 791c9d3e7fc6911753deb056894cc1238c7a9a4c
parent 8e599252
This diff is collapsed.
This diff is collapsed.
! File BEHAVIOUR
!
! The symbol '$' preceeds a keyword used in scanning files.
!
! The symbol 'behav' stands for the nickname of a bulk or
! contact behaviour law, character(LEN=5).
!
! The symbol 'lawty' stands for the name of a bulk or
! contact behaviour law, character(LEN=30).
!
! The symbol 'seety' stands for description of a candidate
! 'cdbdy' type of body, 'cdtac' type of contactor, 'cdcol' color
! ready to meet with the contact behaviour law 'behav' an antagonist
! 'anbdy' type of body, 'antac' type of contactor, 'ancol' color.
!
! Candidate antagonist objects are considered only within some distance
! 'alert'.
!
! STANDARD PACKAGE of bulk or contact behaviour laws
!
! 123456789012345678901234567890:
! :
! bulk behaviour :
! :
!
$gravy
grv1= 0.0000000e+00 grv2= 0.0000000e+00 grv3=-9.8100000e+00
$behav lawty
Rigid RIGID Umas= 2.8000000e+03
$steps 0 time= 0.0000000D+00
$bdyty
RBDY3 1
$nodty
NO6xx 1 X(1)= 0.0000000e+00 X(2)= 0.0000000e+00 X(3)= 0.0000000e+00
X(4)= 0.0000000e+00 X(5)= 0.0000000e+00 X(6)= 0.0000000e+00
V(1)= 0.0000000e+00 V(2)= 0.0000000e+00 V(3)= 0.0000000e+00
V(4)= 0.0000000e+00 V(5)= 0.0000000e+00 V(6)= 0.0000000e+00
a(1)=-1.3555327e-01 a(2)= 9.6245403e-02 a(3)= 9.8608424e-01
b(1)=-2.2765982e-02 b(2)=-9.9531032e-01 b(3)= 9.4016348e-02
c(1)= 9.9050847e-01 c(2)=-9.7049534e-03 c(3)= 1.3710869e-01
$$$$$$
$bdyty
RBDY3 2
$nodty
NO6xx 1 X(1)= 0.0000000e+00 X(2)= 0.0000000e+00 X(3)= 0.8000000e+00
X(4)= 0.0000000e+00 X(5)= 0.0000000e+00 X(6)= 0.0000000e+00
V(1)= 0.0000000e+00 V(2)= 0.0000000e+00 V(3)= 0.0000000e+00
V(4)= 0.0000000e+00 V(5)= 0.0000000e+00 V(6)= 0.0000000e+00
a(1)= 0.1000000e+01 a(2)= 0.0000000e+00 a(3)= 0.0000000e-00
b(1)= 0.0000000e-00 b(2)=-0.1000000e+01 b(3)= 0.0000000e-00
c(1)= 0.0000000e-00 c(2)= 0.0000000e-00 c(3)=-0.1000000e+01
$$$$$$
$steps 0 time= 0.0000000D+00
$bdyty
RBDY3 1
$nodty
NO6xx 1 X(1)= 0.0000000e+00 X(2)= 0.0000000e+00 X(3)= 0.0000000e+00
X(4)= 0.0000000e+00 X(5)= 0.0000000e+00 X(6)= 0.0000000e+00
V(1)= 0.0000000e+00 V(2)= 0.0000000e+00 V(3)= 0.0000000e+00
V(4)= 0.0000000e+00 V(5)= 0.0000000e+00 V(6)= 0.0000000e+00
a(1)=-1.3555327e-01 a(2)= 9.6245403e-02 a(3)= 9.8608424e-01
b(1)=-2.2765982e-02 b(2)=-9.9531032e-01 b(3)= 9.4016348e-02
c(1)= 9.9050847e-01 c(2)=-9.7049534e-03 c(3)= 1.3710869e-01
$$$$$$
$bdyty
RBDY3 2
$nodty
NO6xx 1 X(1)= 0.0000000e+00 X(2)= 0.0000000e+00 X(3)= 0.75000000e+00
X(4)= 0.0000000e+00 X(5)= 0.0000000e+00 X(6)= 0.0000000e+00
V(1)= 0.0000000e+00 V(2)= 0.0000000e+00 V(3)= 0.0000000e+00
V(4)= 0.0000000e+00 V(5)= 0.0000000e+00 V(6)= 0.0000000e+00
a(1)= 0.1000000e+01 a(2)= 0.0000000e+00 a(3)= 0.0000000e-00
b(1)= 0.0000000e-00 b(2)=-0.1000000e+01 b(3)= 0.0000000e-00
c(1)= 0.0000000e-00 c(2)= 0.0000000e-00 c(3)=-0.1000000e+01
$$$$$$
! DOF
$bdyty
RBDY3 1
$nodty
NO6xx 1
$dofty numbr [CT......+......AMP..*..cos.(..OMEGA.*.time.+.PHI..)]...*...[RAMPI.....+.....RAMP.*.time]
vlocy 1 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 2 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 3 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 4 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 5 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 6 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
$$$$$$
! DOF
$bdyty
RBDY3 2
$nodty
NO6xx 1
$dofty numbr [CT......+......AMP..*..cos.(..OMEGA.*.time.+.PHI..)]...*...[RAMPI.....+.....RAMP.*.time]
vlocy 1 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 2 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 3 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 4 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 5 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
vlocy 6 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 1.0000000e+00 0.0000000e+00
$$$$$$
! File MODELS
!
! The symbol '$' preceeds a keyword used in scanning files.
! \
! The symbol 'model' stands for the nickname of a model (LEN=5).
!
! The symbol 'mdlty' stands for the nickname of a physical model (LEN=5),
! meca, therm
!
!
! The symbol 'finel' stands for the nickname of a finite element used for
! the discretisation of the body for the done behaviour character(LEN=5).
!
! The symbol 'lawty' stands for the name of a bulk or character(LEN=30)
!
! The symbol 'eleop' stands for the name of the elementary options
!
! STANDARD PACKAGE of bulk model options
!
!
!
! 123456789012345678901234567890
$model mdlty finel eleop value
Rigid
MECAx Rxx3D
$behav
gapG0 GAP_SGR_CLB fric= 5.0000000e-02
$seety
cdbdy cdtac cdcol behav anbdy antac ancol alert
RBDY3 POLYR REDxx gapG0 RBDY3 POLYR REDxx 2.0000000e-03
$behav
gapG0 GAP_SGR_CLB fric= 5.0000000e-02
$seety
cdbdy cdtac cdcol behav anbdy antac ancol alert
RBDY3 POLYR REDxx gapG0 RBDY3 POLYR BLUEx 2.0000000e-03
# importing chipy module
from pylmgc90.chipy import *
import numpy as N
import pylab as P
# Initializing
Initialize()
# checking/creating mandatory subfolders
checkDirectories()
#utilities_DisableLogMes()
# ------------------------------------------------------
# modeling hypothesis ( 1 = plain strain, 2 = plain stress, 3 = axi-symmetry)
mhyp = 1
dim= 3
# Set space dimension
SetDimension(dim,mhyp)
POLYR_SkipAutomaticReorientation()
#POLYR_FlatnessAngle(5.)
#POLYR_TopologyAngle(80.)
PRPRx_UseNcDetection(0.01)
PRPRx_LowSizeArrayPolyr(1000)
# theta integrator parameter
theta = 0.5
# time evolution parameters
dt = 1e-2
# Defined the time discretisation and integrator
t_final = 20*dt
utilities_logMes('INIT TIME STEPPING')
TimeEvolution_SetTimeStep(dt)
# Initialize theta integrator for all physical model
Integrator_InitTheta(theta)
# write parameter
freq_write = 10
# display parameters
freq_display = 1
ref_radius = 2.5e-2
# interaction parameters
freq_detect = 1
Rloc_tol = 5.e-2
# nlgs parameters
tol = 1e-4
relax = 1.0
norm = 'Quad '
gs_it1 = 50
gs_it2 = 10
type = 'Stored_Delassus_Loops '
nlgs_3D_DiagonalResolution()
#
utilities_logMes('READ BEHAVIOURS')
ReadBehaviours()
ReadModels()
#
utilities_logMes('READ BODIES')
ReadBodies()
#
utilities_logMes('LOAD BEHAVIOURS')
LoadBehaviours()
LoadModels()
#
utilities_logMes('READ INI DOF')
ReadIniDof()
#
utilities_logMes('READ INI GPV')
ReadIniGPV()
#
utilities_logMes('READ DRIVEN DOF')
ReadDrivenDof()
#
utilities_logMes('LOAD TACTORS')
LoadTactors()
#
utilities_logMes('READ INI Vloc Rloc')
ReadIniVlocRloc()
#
# open display & postpro
#
utilities_logMes('DISPLAY & WRITE')
OpenDisplayFiles()
#OpenPostproFiles()
# since constant compute elementary mass matrices once
utilities_logMes('COMPUTE MASS')
ComputeMass()
WriteDisplayFiles(freq_display,ref_radius)
##################### PARTIE THERMIQUE ####################
while TimeEvolution_GetTime() < t_final :
#
utilities_logMes('INCREMENT STEP')
IncrementStep()
#
utilities_logMes('COMPUTE Fext')
ComputeFext()
#
utilities_logMes('COMPUTE Bulk')
ComputeBulk()
#
utilities_logMes('COMPUTE Free Vlocy')
ComputeFreeVelocity()
#
utilities_logMes('SELECT PROX TACTORS')
SelectProxTactors(freq_detect)
#
utilities_logMes('RESOLUTION' )
RecupRloc()
#
ExSolver(type, norm, tol, relax, gs_it1, gs_it2)
UpdateTactBehav()
#
StockRloc()
#
utilities_logMes('COMPUTE DOF, FIELDS, etc.')
ComputeDof()
#
utilities_logMes('UPDATE DOF, FIELDS')
UpdateStep()
#
utilities_logMes('WRITE OUT DOF')
WriteOutDof(freq_write)
#
utilities_logMes('WRITE OUT Rloc')
WriteOutVlocRloc(freq_write)
#
utilities_logMes('VISU & POSTPRO')
WriteDisplayFiles(freq_display,ref_radius)
#~ WritePostproFiles()
# close display & postpro
#
CloseDisplayFiles()
#ClosePostproFiles()
# this is the end
#~ Finalize()
! File BODIES
!
! The symbol '$' preceeds a keyword used in scanning files.
!
! the symbol 'bdyty' stands for 'body type data'.
! These data are distributed according to some species.
!
! the specy 'blmty' stands for 'bulk element type data',
! i.e. part or total bulk geometric description,
! and bulk behaviour laws;
!
! the specy 'nodty' stands for 'node type data',
! i.e. degrees of freedom data;
!
! the specy 'tacty' stands for 'contactor type data';
!
! the keyword '$$$$$$' ends a body record.
!
$bdyty
RBDY3 1
$blmty
PLAIN 1 behav PLAN1 avrd= 0.4211376D+01
I1 = 0.1382191D+06 I2 = 0.1337255D+06 I3 = 0.2719432D+06
$nodty
NO6xx 1 coo1= 0.6079758D+01 coo2= 0.6363961D+01 coo3=-0.1000000D-01
coo4= 0.0000000D+00 coo5= 0.0000000D+00 coo6= 0.0000000D+00
$tacty
PLANx 1 color BLEUx axe1= 0.2532044D+02 axe2= 0.2574235D+02 axe3= 0.6000000D-01
localframe
alp1= 0.1000000D+01 alp2= 0.0000000D+00 alp3= 0.0000000D+00
bet1= 0.0000000D+00 bet2= 0.1000000D+01 bet3= 0.0000000D+00
gam1= 0.0000000D+00 gam2= 0.0000000D+00 gam3= 0.1000000D+01
coo1= 0.0000000D+00 coo2= 0.0000000D+00 coo3= 0.0000000D+00
$$$$$$
$bdyty
RBDY3 2
$blmty
PLAIN 1 behav PLEX1 avrd= 0.2751069D+01
I1 = 0.2491754D+03 I2 = 0.2905847D+03 I3 = 0.3279970D+03
$nodty
NO6xx 1 coo1= 0.4242641D+01 coo2= 0.4242641D+01 coo3= 0.4242641D+01
coo4= 0.0000000D+00 coo5= 0.0000000D+00 coo6= 0.0000000D+00
$tacty
POLYR 1 color BLEUx nb_vertex= 8 nb_faces= 12
coo1= 0.3000000D+01 coo2=-0.3000000D+01 coo3= 0.0000000D+00
coo1= 0.3000000D+01 coo2= 0.3000000D+01 coo3= 0.0000000D+00
coo1=-0.3000000D+01 coo2= 0.3000000D+01 coo3= 0.0000000D+00
coo1=-0.3000000D+01 coo2=-0.3000000D+01 coo3= 0.0000000D+00
coo1= 0.1267949D+01 coo2= 0.0000000D+00 coo3= 0.3000000D+01
coo1=-0.1267949D+01 coo2= 0.0000000D+00 coo3= 0.3000000D+01
coo1= 0.1267949D+01 coo2= 0.0000000D+00 coo3=-0.3000000D+01
coo1=-0.1267949D+01 coo2= 0.0000000D+00 coo3=-0.3000000D+01
ver1= 1 ver2= 2 ver3= 5
ver1= 3 ver2= 4 ver3= 6
ver1= 2 ver2= 6 ver3= 5
ver1= 2 ver2= 3 ver3= 6
ver1= 5 ver2= 4 ver3= 1
ver1= 5 ver2= 6 ver3= 4
ver1= 1 ver2= 7 ver3= 2
ver1= 8 ver2= 4 ver3= 3
ver1= 2 ver2= 7 ver3= 8
ver1= 8 ver2= 3 ver3= 2
ver1= 7 ver2= 1 ver3= 8
ver1= 8 ver2= 1 ver3= 4
$$$$$$
$bdyty
RBDY3 3
$blmty
PLAIN 1 behav PLEX1 avrd= 0.2751069D+01
I1 = 0.2491754D+03 I2 = 0.2905847D+03 I3 = 0.3279970D+03
$nodty
NO6xx 1 coo1= 0.4242641D+01 coo2= 0.4242641D+01 coo3= 0.1024264D+02
coo4= 0.0000000D+00 coo5= 0.0000000D+00 coo6= 0.0000000D+00
$tacty
POLYR 1 color BLEUx nb_vertex= 8 nb_faces= 12
coo1= 0.3000000D+01 coo2=-0.3000000D+01 coo3= 0.0000000D+00
coo1= 0.3000000D+01 coo2= 0.3000000D+01 coo3= 0.0000000D+00
coo1=-0.3000000D+01 coo2= 0.3000000D+01 coo3= 0.0000000D+00
coo1=-0.3000000D+01 coo2=-0.3000000D+01 coo3= 0.0000000D+00
coo1= 0.1267949D+01 coo2= 0.0000000D+00 coo3= 0.3000000D+01
coo1=-0.1267949D+01 coo2= 0.0000000D+00 coo3= 0.3000000D+01
coo1= 0.1267949D+01 coo2= 0.0000000D+00 coo3=-0.3000000D+01
coo1=-0.1267949D+01 coo2= 0.0000000D+00 coo3=-0.3000000D+01
ver1= 1 ver2= 2 ver3= 5
ver1= 3 ver2= 4 ver3= 6
ver1= 2 ver2= 6 ver3= 5
ver1= 2 ver2= 3 ver3= 6
ver1= 5 ver2= 4 ver3= 1
ver1= 5 ver2= 6 ver3= 4
ver1= 1 ver2= 7 ver3= 2
ver1= 8 ver2= 4 ver3= 3
ver1= 2 ver2= 7 ver3= 8
ver1= 8 ver2= 3 ver3= 2
ver1= 7 ver2= 1 ver3= 8
ver1= 8 ver2= 1 ver3= 4
$$$$$$
! File BEHAVIOUR
!
! The symbol '$' preceeds a keyword used in scanning files.
!
! The symbol 'behav' stands for the nickname of a bulk or
! contact behaviour law, character(len=5).
!
! The symbol 'lawty' stands for the name of a bulk or
! contact behaviour law, character(len=30).
!
! The symbol 'seety' stands for description of a candidate
! 'cdbdy' type of body, 'cdtac' type of contactor, 'cdcol' color
! ready to meet with the contact behaviour law 'behav' an antagonist
! 'anbdy' type of body, 'antac' type of contactor, 'ancol' color.
!
! Candidate antagonist objects are considered only within some distance
! 'alert'.
!
! STANDARD PACKAGE of bulk or contact behaviour laws
!
! 123456789012345678901234567890:
! :
! bulk behaviour :
! :
! RIGID : Rigid body
! :
$gravy
grv1= 0.0000000D+01 grv2= 0.0000000D+01 grv3=-0.9810000D+01
$behav lawty
BLEUx RIGID Umas= 0.2800000D+04
$behav lawty
PLAN1 RIGID Umas= 0.2800000D+04
$behav lawty
PLEX1 RIGID Umas= 0.5800000D+04
$behav lawty
TETEy RIGID Umas= 0.5800000D+04
$behav lawty
TETEz RIGID Umas= 0.5800000D+04
#COMMAND.DAT :
#23456789012345678901234567890:
#NO CHECK CHIC COMMAND :
##############################
ECHO OFF :
DIME :
3D :
TIME STEP :
2.D-03 :
THETA :
0.5D0 :
NEW ROTATION SCHEME :
##############################
READ BODIES :
READ BEHAVIOURS :
READ INI DOF :
READ INI Vloc Rloc :
READ DRIVEN DOF :
##############################
WRITE BODIES :
WRITE BEHAVIOURS :
WRITE DRIVEN DOF :
##############################
FD INIT SOURCE POINT :
2 :
0.3000000D+01 :
12. 12. 4.0 :
#
XPERIODIC CONDITION :
25.00
YPERIODIC CONDITION :
25.00
#
COMPUTE BOX :
#
##############################
DISPLAY POSITION :
DISPLAY AVERAGE VELOCITY :
DISPLAY STRESS :
DISPLAY FORCES :
#DISPLAY HEAT :
INIT POST 3D :
START POSTPRO :
WRITE OUTPUT GMV :
##############################
CUNDALL ITERATION :
300
LOW SIZE ARRAY POLYR :
70
#DIAGONAL RESOLUTION :
##############################
COMPUTE MASS :
##############################
loading step :
TOTALISATEURS A ZERO :
INCREMENT STEP :
DISPLAY TIMES :
##############################
COMPUTE Fext :
COMPUTE BULK :
COMPUTE FREE VELOCITY :
##############################
WCP SELECT PROX TACTORS :
##############################
.doing Signorini Coulomb : Some label to draw attention
RECUP Rloc : from verlt into this
###ECHO OFF :
NLGS CHECK TYPE :
Maxm 0.1666D-03 : Maxm,Quad,QM/16
1.D0 : RELAX
EX PREP NLGS :
Exchange_Local_Global :
#
>more :
quickSCRAMBLE :
>iteration :
EX ITER NLGS :
REPETE 10 FOIS :
DEPUIS >iteration :
AFTER ITER CHECK :
# DISPLAY AFTER ITER CHECK :
IOK = FLAG 1 :
INTO AFTER ITER CHECK :
STOP REPETE SI IOK = 0 :
REPETE 100 FOIS :
DEPUIS >more :
EX POST NLGS :
###ECHO ON :
STOCK Rloc : from this into verlt
##############################
##############################
COMPUTE DOF :
UPDATE DOF :
##############################
#WRITE LAST DOF :
#WRITE LAST Vloc Rloc :
#WRITE LAST GMV :
WRITE OUT DOF STEP 100 :
WRITE OUT Vloc Rloc STEP 100 :
UPDATE POST 3D :
POSTPRO DURING COMPUTATION :
WRITE OUTPUT GMV STEP 100 :
##############################
CHECK SOURCE POINT :
##############################
###AFFICHE TOTALISATEURS :
###AFFICHE COMPTEURS :
REPETE 100000 FOIS :
DEPUIS loading step :
##############################
CLOSE POSTPRO FILES :
##############################
WRITE OUT DOF STEP 1 :
WRITE OUT Vloc Rloc STEP 1 :
##############################
##############################
FIN DU PROGRAMME :
FIN DU FICHIER :
! DOF
$steps 0 time= 0.0000000D+00
!-----------------------------------------------------------------------
$bdyty