Code_Aster

Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 1/10 Cl : U2.08.04 Rvision : 1145

Notice de calcul au flambage

Rsum : Lobjectif de cette documentation est de prsenter un guide mthodologique pour une analyse de flambage non linaire dune structure. On y aborde principalement deux fonctionnalits de Code_Aster : lanalyse de flambement linaire, dite dEuler, au travers de MODE_ITER_SIMULT, (option TYPE_RESU : MODE_FLAMB), le calcul de lvolution quasi-statique (oprateur STAT_NON_LINE) de la structure qui prsente des non linarits gomtriques et comportementales, dont on cherche un point limite, voire la rponse post-critique. La premire tape est, gnralement, un calcul de flambage dEuler, qui permettra de connatre les modes de flambement et les charges critiques correspondantes. Du point de vue du concepteur, la connaissance du premier mode et de sa charge critique est souvent suffisante, afin de se dfinir une marge de fonctionnement par rapport au chargement impos : le coefficient multiplicateur entre le chargement impos et la charge critique la plus faible donne la marge de scurit. Remarques La connaissance du premier mode de flambement peut aussi servir dindication pour optimiser la gestion du calcul incrmental non linaire men par la suite. En effet, lapproche de la charge critique, on peut alors dcider de modifier le pilotage ou de rduire le pas de temps, voire daugmenter le nombre ditrations de vrification de lquilibre dans la mthode de rsidu, chaque pas de charge. Lallure du mode de flambement dEuler peut aussi servir pour imposer un dfaut gomtrique initial sur la structure, afin de sassurer, entre autre, que le calcul non linaire incrmental bifurquera bien sur ce mode. Lanalyse dEuler tant par dfinition linaire, elle ne permet pas de prendre en compte des relations de comportement inlastiques ou du contact. Il est alors ncessaire de faire un calcul non linaire, qui en quasi-statique sappuiera sur la commande STAT_NON_LINE de Code_Aster. Cest la mthode classique incrmentale par rsidu en quilibre. Les points particuliers de son utilisation seront abords par la suite.

Manuel d'utilisation Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Fascicule u2.08 : Fonctions avances et contrle des calculs

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 2/10 Cl : U2.08.04 Rvision : 1145

Analyse de flambement dEuler

Le calcul des modes de flambement au sens dEuler [bib5] peut se faire par loprateur de rsolution de problmes aux valeurs propres MODE_ITER_SIMULT (ou bien MODE_ITER_INV). Dans le cadre du flambage, on a la syntaxe typique suivante : MODP1 = MODE_ITER_SIMULT ( MATR_A = RAMEP1 , MATR_B = RAGEP1 , TYPE_RESU = 'MODE_FLAMB' , CALC_FREQ = _F( OPTION = 'BANDE' , CHAR_CRIT = ( -2.4 ,-2.2 , ) , DIM_SOUS_ESPACE = 80 , NMAX_ITER_SOREN = 80 , ) , ) Largument du mot cl MATR_A doit tre la matrice de rigidit dite matrielle, alors que le mot cl MATR_B attend la matrice de rigidit gomtrique. Si on avait employ loprateur MODE_ITER_INV, les arguments des mots cls MATR_A et MATR_B seraient les mmes. Pour rappel, les modes de flambement sont les modes propres du problme aux valeurs propres suivant :

K K x=0 K x= K g x
avec :

K : matrice de rigidit K g : matrice de rigidit gomtrique : valeur propre= avec coefficient multiplicateur du chargement

La rigidit matrielle (ou lastique) se calcule avec loption RIGI_MECA de CALC_MATR_ELEM. La rigidit gomtrique se calcule partir du champ de contrainte solution du problme linaire (option RIGI_GEOM de CALC_MATR_ELEM). Il faut donc avoir effectu un calcul linaire statique pralablement lutilisation de MODE_ITER_SIMULT pour le flambement. Si le chargement est compos dune partie fixe (non pilote) et dune partie variable, le coefficient multiplicateur du chargement ne doit, bien sr, porter que sur la partie variable. La contribution de lautre partie du chargement se retrouve dans le premier membre. Notons f c le chargement fixe et f v le chargement pilot (proportionnel ). Le problme aux valeurs propres devient :

K K g f c f v x=0 K K g f c x= K g f v x
avec :

K: K g f c : K g f v : :

matrice de rigidit matrielle matrice de rigidit gomtrique pour le chargement non pilot matrice de rigidit gomtrique pour le chargement variable valeur propre=

Dans ce cas, il faut donc rsoudre deux problmes lastiques linaires pralables, pour pouvoir calculer les deux matrices de rigidit gomtriques diffrentes. Afin dtre exhaustif, la prsentation portera sur une structure soumise des dplacements imposs ainsi que des efforts, qui seront la combinaison dun chargement fixe et dun chargement variable que lon pilotera avec un coefficient croissant pouvant conduire au flambage.

Manuel d'utilisation

Fascicule u2.08 : Fonctions avances et contrle des calculs

Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 3/10 Cl : U2.08.04 Rvision : 1145

1.1

Etape 1 : Calcul(s) linaire(s) pralable(s)

On va se servir de MECA_STATIQUE. La structure, maille en lments de type coque (lments de type coques volumiques [bib3]), est soumise des conditions aux limites de Dirichlet (CONDLIM) et de Neumann. Ces dernires se dcomposent en : PESA : champ de pesanteur, PRESPH : champ de pression impos non pilot, PRESPS1 : champ de pression impos variable. Pour lanalyse de flambage, il faut sparer les efforts constants de ceux qui sont variables (pilots par un coefficient). On va donc faire deux calculs statiques linaires. Le premier sera le cas de la structure soumise aux dplacements imposs et aux efforts constants, le second verra la structure soumise aux dplacements imposs et aux efforts variables. Chargement pilot : RESC11P1 = MECA_STATIQUE ( MODELE = MODELE , CHAM_MATER = CHMAT , CARA_ELEM = CARAELEM , EXCIT = ( _F( CHARGE = CONDLIM , ) , _F( CHARGE = PRESPS1 , ) , ) , OPTION = 'SIEF_ELGA_DEPL' , PLAN = 'MOY' , ) Chargement non pilot : RESC12P1 = MECA_STATIQUE ( MODELE = MODELE , CHAM_MATER = CHMAT , CARA_ELEM = CARAELEM, EXCIT = ( _F( CHARGE = CONDLIM , ) , _F( CHARGE = PESA , ) , _F( CHARGE = PRESPH , ) , ) , OPTION = 'SIEF_ELGA_DEPL' , PLAN = 'MOY' , ) On va utiliser le champ de contrainte pour calculer les matrices de rigidit gomtrique associes, pour les deux chargements : SIGC11P1 = CREA_CHAMP ( TYPE_CHAM = 'ELGA_SIEF_R' , OPERATION = 'EXTR' , RESULTAT = RESC11P1 , NOM_CHAM = 'SIEF_ELGA_DEPL' , TYPE_MAXI = 'MINI' , TYPE_RESU = 'VALE' , ) # REGC11P1 = CALC_MATR_ELEM ( OPTION = 'RIGI_GEOM' , MODELE = MODELE , CARA_ELEM = CARAELEM , SIEF_ELGA = SIGC11P1 , ) REGC11P1 est donc la matrice de raideur gomtrique associe au cas de chargement variable (PRESPS1).

Manuel d'utilisation Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Fascicule u2.08 : Fonctions avances et contrle des calculs

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 4/10 Cl : U2.08.04 Rvision : 1145

On calcule, de mme, la matrice de raideur gomtrique pour le chargement constant (PESA et PRESPH), partir de RESC12P1 : SIGC12P1 = CREA_CHAMP ( TYPE_CHAM = 'ELGA_SIEF_R' , OPERATION = 'EXTR' , RESULTAT = RESC12P1 , NOM_CHAM = 'SIEF_ELGA_DEPL' , TYPE_MAXI = 'MINI' , TYPE_RESU = 'VALE' , ) # REGC12P1 = CALC_MATR_ELEM ( OPTION = 'RIGI_GEOM' , MODELE = MODELE , CARA_ELEM = CARAELEM , SIEF_ELGA = SIGC12P1 , ) Il reste calculer la matrice de rigidit matrielle pour le chargement total : REMEP1 = CALC_MATR_ELEM ( OPTION = 'RIGI_MECA' , MODELE = MODELE , CHAM_MATER = CHMAT , CARA_ELEM = CARAELEM , CHARGE = ( CONDLIM , PESA , PRESPH , PRESPS1 , ) , ) Toutes les matrices lmentaires sont calcules, ltape suivante est donc leur assemblage : NUP1 = NUME_DDL ( MATR_RIGI = REMEP1 , ) # RAMC1P1 = ASSE_MATRICE ( MATR_ELEM = REMEP1 , NUME_DDL = NUP1 , ) # RAGEP1 = ASSE_MATRICE ( MATR_ELEM = REGC11P1 , NUME_DDL = NUP1 , ) # RAGC12P1 = ASSE_MATRICE ( MATR_ELEM = REGC12P1 , NUME_DDL = NUP1 , ) On somme ensuite les matrices de rigidit matrielle (RAMC1P1) et gomtrique (RAGC12P1) correspondant au cas de chargement constant : RAMEP1 = COMB_MATR_ASSE ( COMB_R = ( _F( MATR_ASSE = RAMC1P1 , COEF_R = 1.0 , ) , _F( MATR_ASSE = RAGC12P1 , COEF_R = 1.0 , ) , ) , ) Les deux matrices ncessaires au calcul des modes de flambement sont donc construites.

Manuel d'utilisation

Fascicule u2.08 : Fonctions avances et contrle des calculs

Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 5/10 Cl : U2.08.04 Rvision : 1145

1.2

Etape 2 : Calcul des modes dEuler

Il peut tre utile de faire des tests de STURM (oprateur IMPR_STURM) sur lintervalle de recherche sur lequel on veut trouver les cas de flambement. Ainsi, cela permettra doptimiser la taille de lintervalle et de contrler le bon droulement du calcul modal ultrieur puisquon connatra davance le nombre de modes existants. La syntaxe est : IMPR_STURM ( MATR_A = RAMEP1 , MATR_B = RAGEP1 , TYPE_RESU = 'MODE_FLAMB' , CHAR_CRIT_MIN = -2.4 , CHAR_CRIT_MAX = -2.2 , )

Une fois lintervalle de recherche de charge critique de flambage choisi, on peut alors mettre en uvre MODE_ITER_SIMULT comme suit : MODP1 = MODE_ITER_SIMULT ( MATR_A = RAMEP1 , MATR_B = RAGEP1 , TYPE_RESU = 'MODE_FLAMB' , CALC_FREQ = _F( OPTION = 'BANDE' , CHAR_CRIT = ( -2.4 ,-2.2 , ) , DIM_SOUS_ESPACE = 80 , NMAX_ITER_SOREN = 80 , ) , ) Remarque Si lalgorithme ne converge pas ou si le nombre de modes nest pas celui prdit par IMPR_STURM , il peut tre utile daugmenter les valeurs de DIM_SOUS_ESPACE et NMAX_ITER_SOREN . On norme les modes [bib6], uniquement en se servant des degrs de libert de translation : MODP1 = NORM_MODE ( reuse = MODP1 MODE = MODP1 , NORME = 'TRAN' , ) Les modes peuvent ensuite tre post-traits. Remarques Il est indispensable de vrifier que la raideur gomtrique du modle choisi est bien une option disponible dans Code_Aster (par exemple, ce nest pas le cas des DKT). Une discrtisation plus fine conduit normalement une baisse des charges critiques. La discrtisation doit tre apte capter les modes de flambement, sachant que ces modes peuvent engendrer des dformations localises (plis). Le calcul pralable des modes dynamiques peut constituer une premire indication sur la qualit du maillage, bien que ces modes puissent tre trs diffrents des modes de flambement. Les charges critiques des diffrents modes sont proportionnelles au module dYoung E.

Manuel d'utilisation Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Fascicule u2.08 : Fonctions avances et contrle des calculs

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 6/10 Cl : U2.08.04 Rvision : 1145

Etude non linaire quasistatique de la structure

Cette tape se justifie si la structure prsente de fortes non linarits, dont lanalyse dEuler ne peut tenir compte. Loprateur de rsolution des problmes non linaires en quasi-statique se nomme STAT_NON_LINE [bib7]. Ces non linarits peuvent tre lies au matriau qui peut avoir un comportement lastoplastique [bib8], comme dans lexemple qui va suivre. La prise en compte du contact, voire du frottement, est une autre source de non linarits. On peut aussi citer le cas des chargements suiveurs, comme la pression ([bib1] et [bib2] pour les lments de type coques volumiques), qui ncessitent une approche non linaire. Pour ltude dune structure potentiellement instable ou susceptible de connatre un point limite, qui risque donc de rencontrer une bifurcation en solution au cours de lvolution du chargement, il est souvent utile de pouvoir choisir une branche de solution particulire (souvent la solution physique quand elle est dfinie a priori sans ambiguts). Pour cela, lutilisateur peut avoir introduire un dfaut initial qui va forcer la structure bifurquer sur la branche de solution particulire. Plusieurs mthodes existent pour dfinir ce dfaut. Lune des plus adapte est de prdformer lgrement la structure suivant lallure du mode dEuler de flambement correspondant la branche que lon veut suivre. Lamplitude de cette prdformation doit tre faible, par exemple moins de 1/10me de lpaisseur pour une structure mince. Lidal tant de trouver le dfaut minimal qui est compatible avec une performance satisfaisante de lalgorithme de rsidu en quilibre. En effet, un dfaut trop faible peut entraner une difficult de convergence du rsidu, principalement dans le cas dun pilotage en effort. Le dfaut gomtrique peut aussi tre dfini par mesures exprimentales de la pice relle dont la gomtrie ne saurait tre parfaite. Le dfaut peut aussi prendre la forme dune perturbation du chargement (dfaut dalignement, rajout dun chargement localis, ...) ou des caractristiques mcaniques du matriau (affaiblissement local du module dYoung, par exemple). Il peut nanmoins tre alors plus difficile dadapter le dfaut au mode de flambage dsir, surtout si la structure prsente des modes relativement voisins. Remarque Dans certains cas, mme sur le problme non perturb, le chargement est tel quil provoque la bifurcation dsire. Un des autres points particuliers, lis linstabilit, est le choix de la technique de pilotage de lalgorithme STAT_NON_LINE. En effet, le pilotage classique en effort nest plus adapt car il ne peut capter une branche instable de solution. De mme, lapproche dun point limite, la convergence avec le pilotage en effort deviendra de plus en plus difficile, la matrice de rigidit tangente devenant singulire. Il est alors ncessaire de rduire lincrment de charge et daugmenter le nombre maximal ditration pour continuer le calcul. Il existe des techniques de pilotage [bib9] permettant de contourner ces difficults numriques. Parmi les mthodes proposes par Code_Aster, celle dite par longueur darc [bib12] (option TYPE=LONG_ARC du mot cl PILOTAGE dans STAT_NON_LINE), qui est la plus adapte pour les instabilits de type flambage, dans le cas de snap-backs ventuels doux [bib13]. Pour les cas de snap-backs plus brutaux, Crisfield propose une variante [bib13], non disponible dans Code_Aster. Dautres mthodes existent, comme celle de Riks [bib14] (non disponible non plus), qui traite aussi le cas dynamique. Si lon ne veut quobtenir le point limite, y compris avec une bonne prcision, un pilotage en chargement peut suffire, condition de bien grer les paramtres de pas dincrment de charge (SUBD_PAS et SUBD_PAS_MINI du mot cl INCREMENT) et de nombre ditrations maximal autoris (ITER_GLOB_MAXI de CONVERGENCE). Il peut aussi tre utile, lapproche du point limite, de ne plus employer la matrice tangente ractualise pour le solveur, puisquelle est quasi-singulire. On peut alors se contenter de ne pas ractualiser cette matrice chaque calcul (paramtres REAC_INCR et REAC_ITER) ou, dans le pire des cas, adopter la matrice lastique de base (PREDICTION=ELASTIQUE et MATRICE=ELASTIQUE du mot cl NEWTON).

Manuel d'utilisation

Fascicule u2.08 : Fonctions avances et contrle des calculs

Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 7/10 Cl : U2.08.04 Rvision : 1145

Voici un exemple dutilisation de STAT_NON_LINE pour un calcul lastoplastique en grands dplacements ([bib4] pour les lments employs, qui sont de type coques volumiques), avec pilotage en efforts : RESU = STAT_NON_LINE ( MODELE = MODELE , CHAM_MATER = CHMAT , CARA_ELEM = CARAELEM , EXCIT = ( _F( CHARGE = CONDLIM , TYPE_CHARGE = 'FIXE_CSTE' , ) , _F( CHARGE = PESA , TYPE_CHARGE = 'FIXE_CSTE' , ) , _F( CHARGE = PRESPH , FONC_MULT = FONCMUL2 , TYPE_CHARGE = 'SUIV' , ) , _F( CHARGE = PRESPS1 , FONC_MULT = FONCMUL , TYPE_CHARGE = 'SUIV' , ) , ) , COMP_INCR = ( _F( RELATION = 'VMIS_ISOT_TRAC' , DEFORMATION = 'GREEN_GR' , GROUP_MA = ( 'VIROLE' , 'TOIT' , 'ANNEAUX' , 'SGOU') , ) , ) , COMP_ELAS = _F( RELATION = 'ELAS' , DEFORMATION = 'GREEN_GR' , GROUP_MA = 'LTIGE' , ) , INCREMENT = _F( LIST_INST = L_INST1 , NUME_INST_FIN = 14 , SUBD_PAS = 4 , SUBD_PAS_MINI = 1.E-9 , ) , NEWTON = _F( REAC_INCR = 1 , PREDICTION = 'TANGENTE' , MATRICE = 'TANGENTE' , REAC_ITER = 1 , ) , CONVERGENCE = _F( RESI_GLOB_RELA = 1.E-06 , ITER_GLOB_MAXI = 40 , ARRET = 'OUI' , ) , SOLVEUR = _F( METHODE = 'MULT_FRONT' , RENUM = 'METIS' , ) , ) Remarques On utilise la matrice tangente ractualise chaque calcul, en autorisant le sousdcoupage du pas de charge. Les pressions imposes sont des efforts suiveurs ( TYPE_CHARGE='SUIV' ). Dans le cas dune modlisation en lments massifs, le tenseur de dformation recommand en grands dplacements est SIMO_MIEHE .

Manuel d'utilisation Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Fascicule u2.08 : Fonctions avances et contrle des calculs

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 8/10 Cl : U2.08.04 Rvision : 1145

Si on veut remplacer le pilotage en effort par une mthode par longueur darc, il suffit dcrire : RESU = STAT_NON_LINE ( MODELE = MODELE , CHAM_MATER = CHMAT , CARA_ELEM = CARAELEM , EXCIT = ( _F( CHARGE = CONDLIM , TYPE_CHARGE = 'FIXE_CSTE' , ) , _F( CHARGE = PESA , TYPE_CHARGE = 'FIXE_CSTE' , ) , _F( CHARGE = PRESPH , FONC_MULT = FONCMUL2 , TYPE_CHARGE = 'SUIV' , ) , _F( CHARGE = PRESPS1 , TYPE_CHARGE = 'FIXE_PILO' , ) , ) , COMP_INCR = ( _F( RELATION = 'VMIS_ISOT_TRAC' , DEFORMATION = 'GREEN_GR' , GROUP_MA = ( 'VIROLE' , 'TOIT' , 'ANNEAUX' , 'SGOU') , ) , ) , COMP_ELAS = _F( RELATION = 'ELAS' , DEFORMATION = 'GREEN_GR' , GROUP_MA = 'LTIGE' , ) , INCREMENT = _F( LIST_INST = L_INST1 , NUME_INST_FIN = 14 , SUBD_PAS = 4 , SUBD_PAS_MINI = 1.E-9 , ) , NEWTON = _F( REAC_INCR = 1 , PREDICTION = 'TANGENTE' , MATRICE = 'TANGENTE' , REAC_ITER = 1 , ) , CONVERGENCE = _F( RESI_GLOB_RELA = 1.E-06 , ITER_GLOB_MAXI = 40 , ARRET = 'OUI' , ) , PILOTAGE = _F( GROUP_NO = 'G' , TYPE = 'LONG_ARC' , NOM_CMP = ( 'DY' , ) , COEF_MULT = 7. ) , ) Remarques Dans Code_Aster, on ne peut pas piloter de forces suiveuses. Pour le pilotage par longueur darc, il est, en gnral, recommand que GROUP_NO contienne toute la structure.

Manuel d'utilisation

Fascicule u2.08 : Fonctions avances et contrle des calculs

Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 9/10 Cl : U2.08.04 Rvision : 1145

Pour finir, citons deux articles de Crisfield qui donnent une bonne vision gnrale des problmes et mthodes lis aux calculs non linaires pouvant prsenter divers types dinstabilits ([bib15] et [bib11]). Quelques cas-tests de Code_Aster traitant du flambage : Modes dEuler : sdls504 sdls505 ssll103 ssll105 ssll403 ssll404 ssls110 Modes dEuler et calcul non linaire : ssnl123 Calcul non linaire : ssnl502 ssnp305 : calcul jusqu' un snap-through

Manuel d'utilisation Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

Fascicule u2.08 : Fonctions avances et contrle des calculs

Code_Aster
Titre : Notice de calcul au flambage Responsable : Nicolas GREFFET

Version default
Date : 29/04/2009 Page : 10/10 Cl : U2.08.04 Rvision : 1145

Bibliographie
1) 2) 3) 4) 5) 6) 7) 8) 9) Efforts extrieurs de pression en grands dplacements [R3.03.04]. Pression suiveuse pour les lments de coques volumiques [R3.03.07]. Elments finis de coques volumiques [R3.07.04]. Elments de coques volumiques en non linaire gomtrique [R3.07.05]. Algorithme de rsolution pour le problme gnralis [R5.01.01]. Paramtres modaux et norme des vecteurs propres [R5.01.03]. Algorithme non linaire quasi statique [R5.03.01]. Intgration des relations de comportement lasto-plastique de Von Mises [R5.03.02]. Mthodes de pilotage du chargement [R5.03.80].

10) Mthode multifrontale [R6.02.02]. 11) M.A. Crisfield, G. Jelenic, Y. Mi, H.-G. ZhongG & Z. Fan : Some aspects of the non-linear finite element method, Finite Elements in Analysis and Design, Vol. 27, 19-40, 1997. 12) M.A. Crisfield : A fast incremental iterative solution procedure that handles snap through, Computers & Structures, Vol. 13, 55-62, 1981. 13) H.-B. Hellweg & M.A. Crisfield : A new arc-length method for handling sharp snap-backs, Computers & Structures, Vol. 66, 705-709, 1998. 14) E. Riks, C.C. Rankin & F.A. Brogan : On the solution of mode jumping phenomena in thinwalled shell structures , Comp. Meth. In Applied Mech. And Engrg., Vol. 1367, 59-92, 1996. 15) J. Shi & M.A. Crisfield : Combining arc-length and line searches in path-following , Comm. Numer. Meth. Engrg, Vol. 11, 793-803, 1995.

Manuel d'utilisation

Fascicule u2.08 : Fonctions avances et contrle des calculs

Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)