Télécharger fluendo3d_helmholtz.dgibi
* fichier : fluendo3d_helmholtz.dgibi ************************************************************************ ************************************************************************ * * test de la formulation non-locale differentielle avec fluendo3d * ---------------------------------------------------------------- * * Alain Sellier, Alain Millard mercredi 4 janvier 2023 * *----------------------------------------------------------------------- * Exemple de calcul non local avec une formulation d'Helmholtz * a plusieurs variables, en 3D sur un tirant de section variable. * 6 variables d'Helmholtz sont utilisees dans un contexte a coefficients * constant (sans sous iteration de HelmHoltz). Les variables * correspondent aux increments du tenseur de plasticite en traction * du modele FLUENDO3D. Dans ce modele des variables internes nommees * DNL1 a DNL6 permettent de diffuser les increments avec une longueur * caracteristique lc0, les increments sont sommes dans ENL1 a ENL6, * ces dernieres variables sont alors utilisees pour calculer * un endommagement orthotrope *----------------------------------------------------------------------- * test formulation Helmholtz avec Fluendo3d * graphique en fin de calcul si VRAI complet=FAUX; arret=VRAI; test=VRAI; L1=0.05; H1=0.05; H2=0.04; *nombre d element sur 1/2 longueur n1=10; DH12=(H1-H2)/2; P1=0. 0. 0.; P2=L1 Dh12 0.; P3=L1 (DH12+H2) 0.; P4=0. H1 0.; P5=L1 DH12 (DH12+H2); surf1=daller l12 l23 l34 l41; vol3=vol1 et vol2; *trac vol3; *blocage blo4 = blo0 et blo1 et blo2 et blo3 ; *deplcement impose *modele mod1=MODELE vol3 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'VISCOPLASTIQUE' 'FLUENDO3D' 'NON_LOCAL' 'HELM' 'V_MOYENNE' *materiau rt0=3.8; rc0 = 38.; *ep0=1.5e-4; E0=32000.; ept0=1.0*rt0/e0; gft0 = 1.10e-4 ; *Beton * endommagement HYDS 0.2 RT rt0 EPT ept0 GFT gft0 RC rc0 EPC 2.0e-3 DELT 0.75 BETA 0.15 EKDC 5.0e-3 REF 3. GFR 2.0e-4 ALTC 0. ; * effets de l eau CSHR 0. TTKW 40. TTRW 20. KWRT 0.25 KWRC 0.1; * fluage XFLU 1. NRJM 20000. DFMX 0. ; *Weibull scale effect : note cf below to prepar Vmax field *matw1=MATE mod1 VREF 2.1e-4 VMAX 2.1e-4 CVRT 0.15 DHI1 0.5; *---------NON LOCAL ------------ * 1ere variable non locale ***** *longueur de diffusion lc0=0.02; *coeff de diffusion dh0=(lc0**2)/2.d0; *diffusion des composantes du tenseur des increment de deformations plastiques * DNL1=dexx 'DH11' dh0 'V111' 1. 'V112' 0. 'V113' 0. 'DH12' 0. 'V121' 0. 'V122' 1. 'V123' 0. 'DH13' 0. 'V131' 0. 'V132' 0. 'V133' 1.; * DNL2=deyy 'DH21' 0. 'V211' 1. 'V212' 0. 'V213' 0. 'DH22' dh0 'V221' 0. 'V222' 1. 'V223' 1. 'DH23' 0. 'V231' 0. 'V232' 0. 'V233' 1.; * DNL3=dezz 'DH31' 0. 'V311' 1. 'V312' 0. 'V313' 0. 'DH32' 0. 'V321' 0. 'V322' 1. 'V323' 0. 'DH33' dh0 'V331' 0. 'V332' 0. 'V333' 1.; * DNL4=dexy 'DH41' (dh0/2.) 'V411' 1. 'V412' 0. 'V413' 0. 'DH42' (dh0/2.) 'V421' 0. 'V422' 1. 'V423' 0. 'DH43' 0. 'V431' 0. 'V432' 0. 'V433' 1.; * DNL5=dexz 'DH51' (dh0/2.) 'V511' 1. 'V512' 0. 'V513' 0. 'DH52' 0.0 'V521' 0. 'V522' 1. 'V523' 0. 'DH53' (dh0/2.) 'V531' 0. 'V532' 0. 'V533' 1.; * DNL6=deyz 'DH61' 0.0 'V611' 1. 'V612' 0. 'V613' 0. 'DH62' (dh0/2.) 'V621' 0. 'V622' 1. 'V623' 0. 'DH63' (dh0/2.) 'V631' 0. 'V632' 0. 'V633' 1.; *----------------------------------------------------------------------- *material chracteristics assemblage mat1=matb1 ET matd1 ET matf1 ET matgw1 ET matnl1 ET matnl2 ET matnl3 ET matnl4 ET matnl5 ET matnl6 ; tab1=table; tab1.modele=mod1; tab1.caracteristiques=mat1; tab1.blocages_mecaniques=blo4; tab1.chargement=char1; *liste des temps calcules tab1.temps_calcules=list3 ; tab1.precision=1.0E-4; TAB1.processeurs= CHAINE 'AUTOMATIQUE'; *controleurs pour les sous iteration de Helmholtz TAB1.MAXSOUSITERATION = 500; TAB1.PRECSOUSITERATION = 1.0E-6; pasapas tab1; ************************************************************************ * exploitation des resultats ************************************************************************ *force deplacement si (test); n1=20; fref1=7.14727E-03; n2=40; fref2=4.51279E-03; err1 = abs ((fref1 - fcal1)/fref1) ; err2 = abs ((fref2 - fcal2)/fref2) ; si ((err1 < 1.e-3) et (err2 < 5.e-3)) ; sinon ; finsi ; finsi; *exemple de traces si (complet); dess evolfx1; trac cmod1 mod1 def1 ; finsi; si(arret); fin; finsi;
© Cast3M 2003 - Tous droits réservés.
Mentions légales