Télécharger snap_non_associe.dgibi
* fichier : snap.dgibi ************************************************************************ ************************************************************************ ************************************************************************* * TEST SNAP * * EXEMPLE D UTILISATION DE LA PROCEDURE PASAPAS ( ET INCREME) * * PROBLEME DE GRANDS DEPLACEMENTS * * PROBLEME DU SNAP une seule barre est maillee et calculee * * Utilisation d'une relation non associee entre le deplacement du point de controle * et la force de chargement * * Utilisation de K_SIGMA pour eviter la raideur tangente nulle au declanchement du snap * * * || * || * \/ F * /\ * / \ * | / \ | * |/ \| * | | ************************************************************************* * * TITRE ' SNAP ' ; * *TEMPS ; * *------------- geometrie ligne ST formee d' 1 SEG2 ------------------- * PAR SYMETRIE, ON N'ETUDIE QUE LA MOITIE * P1 = 0. 1. ; P2 = 10. 0. ; ST = P1 D 1 P2 ; * *------------ calcul mecanique . ------------------ * * * MATERIAU ET CARACTERISTIQUES * * MACA= MA1 ET CAR1; * *----------- calcul de la rigidite -------------------------------------- * * *----------- definition des conditions aux limites ---------------------- * CL = CL1 ET CL2 ET CL3 ; * *----------- definition du chargement ----------------------------------- * FP11 = FORCE ( 0 -12.5e5 ) p1; * *----------- resolution par la procedure NONLIN ------------------------- * CALCUL EN GRANDS DEPLACEMENTS * TAB2 = TABLE ; **CHA1 = CHAR MECA FP11 EV1 ; TAB2.'GRANDS_DEPLACEMENTS'=VRAI; * * relation non associee pour avoir un chargement en deplacement impose TAB2.'BLOCAGES_MECANIQUES' = CL et cls; TAB2.'MODELE' = MO; TAB2.'CARACTERISTIQUES' = MACA; TAB2.'TEMPS_CALCULES' = LIS11; TAB2.'PRECISION' = 1d-8; TAB2.'K_SIGMA' = vrai; PASAPAS TAB2 ; * *----------- resultats -------------------------------------------------- * courbe de snap through : montee descente montee REPETER TBOU2 NDIM ; LEDEP = TAB2 . DEPLACEMENTS. (&TBOU2); FIN TBOU2 ; * pgx = pgx * -1.; sref = -2.e+6; smp= sref * 1.05;smb = sref * 0.95; mess ' borne min ' smp 'valeur trouvée' strou 'borne max' smb; * si (( strou > smb ) ou ( strou < smp) ) ; erreur (5); finsi; * DESS EV5 ; fin;
© Cast3M 2003 - Tous droits réservés.
Mentions légales