×

Notice

The forum is in read only mode.

periods and forms of natural vibrations

More
13 years 2 weeks ago - 13 years 2 weeks ago #6269 by bolt
beam.7z
beam.zip
Hi! Iam want calculate periods and forms of natural vibrations and form and reserves the loss of stability
#Beam tutorial - Claus Andersen - Feb. 2011 - For CAELinux
#
#Units in mm, MPa, N

DEBUT();

#Read mesh

Mesh=LIRE_MAILLAGE(FORMAT='MED',);

Mesh=DEFI_GROUP(reuse =Mesh,
                MAILLAGE=Mesh,
                CREA_GROUP_MA=_F(NOM='TOUT',
                                 TOUT='OUI',),);

#Assign a beam model - POU_D_E denotes a Euler/Bernoulli

Model=AFFE_MODELE(MAILLAGE=Mesh,
                  AFFE=_F(TOUT='OUI',
                          PHENOMENE='MECANIQUE',
                          MODELISATION='POU_D_E',),);

Material=DEFI_MATERIAU(ELAS=_F(E=66000,
                               NU=0.3,RHO=1200),);

MatField=AFFE_MATERIAU(MAILLAGE=Mesh,
                       AFFE=_F(TOUT='OUI',
                               MATER=Material,),);
A = 27400.0;

I_y = 712000000.0;

I_z = 283000000.0;

J = 6180000.0;

RY = 197;

RZ = 187;

AY = 1.46472;

AZ = 4.64168;











CARA_POU=AFFE_CARA_ELEM(MODELE=Model,
                        POUTRE=_F(GROUP_MA='TOUT',
                                  SECTION='GENERALE',
                                  CARA=('A','IY','IZ','JX','RY','RZ','AY','AZ',),
                                  VALE=(A,I_y,I_z,J,RY,RZ,AY,AZ,),),);

#Boundary conditions
#Point 1 : Fixed

Hold=AFFE_CHAR_MECA(MODELE=Model,
                    DDL_IMPO=_F(GROUP_NO='fixx2',
                                LIAISON='ENCASTRE',),);

#Load
#Point 2 : Point load

Load=AFFE_CHAR_MECA(MODELE=Model,
                    FORCE_NODALE=_F(GROUP_NO='loadd',
                                    FZ=-10000,),);

# Linear elastic solution

RESU1=MECA_STATIQUE(MODELE=Model,
                    CHAM_MATER=MatField,
                    CARA_ELEM=CARA_POU,
                    EXCIT=(_F(CHARGE=Hold,),
                           _F(CHARGE=Load,),),);

# Calculate forces etc.

RESU1=CALC_ELEM(reuse =RESU1,
                RESULTAT=RESU1,
                OPTION=('EFGE_ELNO_DEPL','SIPO_ELNO_DEPL','SIGM_ELNO_DEPL',),);

RESU1=CALC_NO(reuse =RESU1,
              RESULTAT=RESU1,
              OPTION=('REAC_NODA','EFGE_NOEU_DEPL','SIPO_NOEU_DEPL','FORC_NODA',),);
			  



#######################

Me=CALC_MATR_ELEM(OPTION='MASS_MECA',
                  MODELE=Model,
                  CHAM_MATER=MatField,
				  CARA_ELEM=CARA_POU,
                  CHARGE=Load,);

Ke=CALC_MATR_ELEM(OPTION='RIGI_MECA',
                  MODELE=Model,
                  CHAM_MATER=MatField,
				  CARA_ELEM=CARA_POU,
                  CHARGE=Load,);

NU=NUME_DDL(MATR_RIGI=Ke,);

M=ASSE_MATRICE(MATR_ELEM=Me,
               NUME_DDL=NU,);

K=ASSE_MATRICE(MATR_ELEM=Ke,
               NUME_DDL=NU,);

MODE=MODE_ITER_SIMULT(MATR_A=K,
                      MATR_B=M,
					  METHODE='SORENSEN',
#					  PREC_SOREN=0.00005,
                      NMAX_ITER_SOREN=80,
                      CALC_FREQ=_F(OPTION='PLUS_PETITE',
                                   NMAX_FREQ=20,
#                             DIM_SOUS_ESPACE=1,
								   COEF_DIM_ESPACE=20,
								   ),);
								   

								   
Resul=EXTR_MODE(FILTRE_MODE=_F(MODE=MODE,
                               NUME_MODE=(1,2,3,4,5,6,7,8,),),);
							   
IMPR_RESU(MODELE=Model,
          FORMAT='MED',
          RESU=_F(MAILLAGE=Mesh,
                  RESULTAT=Resul,
                  NOM_CHAM=('DEPL','EFGE_NOEU_DEPL','REAC_NODA','SIPO_NOEU_DEPL','FORC_NODA',),),);
				  							   
					   
											   
#############		  

# Write to .med file

IMPR_RESU(MODELE=Model,
          FORMAT='MED',
          RESU=_F(MAILLAGE=Mesh,
                  RESULTAT=RESU1,
                  NOM_CHAM=('DEPL','EFGE_NOEU_DEPL','REAC_NODA','SIPO_NOEU_DEPL','FORC_NODA',),),);
				  
				  
				  

#Write to .resu file

IMPR_RESU(MODELE=Model,
          FORMAT='RESULTAT',
          RESU=_F(RESULTAT=RESU1,
                  NOM_CHAM=('DEPL','EFGE_NOEU_DEPL','SIPO_NOEU_DEPL','REAC_NODA','SIGM_ELNO_DEPL',),
                  VALE_MAX='OUI',
                  VALE_MIN='OUI',),);

FIN();


I problems.
   !-------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_23>                      !
   !                                                 !
   ! pour le concept  MODE                           !
   !                                                 !
   !   dans l'intervalle  ( -0.010025  ,  0.010025 ) !
   !   il y a thйoriquement  176 frequence(s)        !
   !   et on en a calculй  10                        !
   !-------------------------------------------------!
   
   
   !----------------------------!
   ! <EXCEPTION> <ALGELINE2_74> !
   !                            !
   ! erreur de vйrification     !
   !----------------------------!



If I increase NMAX_FREQ=176,

I have many problems
!----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  98                                        !
   !                                                                                  !
   !   de frequence  0.002285                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000001  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  108                                       !
   !                                                                                  !
   !   de frequence  0.002793                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000003  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  109                                       !
   !                                                                                  !
   !   de frequence  0.002797                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000002  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  110                                       !
   !                                                                                  !
   !   de frequence  0.002797                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000003  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  111                                       !
   !                                                                                  !
   !   de frequence  0.002800                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000003  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  115                                       !
   !                                                                                  !
   !   de frequence  0.003259                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000002  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  116                                       !
   !                                                                                  !
   !   de frequence  0.003260                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000004  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  117                                       !
   !                                                                                  !
   !   de frequence  0.003262                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000001  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  118                                       !
   !                                                                                  !
   !   de frequence  0.003262                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000007  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  121                                       !
   !                                                                                  !
   !   de frequence  0.003556                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000011  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).                                                                  !
   !----------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------!
   ! <EXCEPTION> <ALGELINE5_15>                                                       !
   !                                                                                  !
   ! pour le concept  MODE  le mode numйro  122                                       !
   !                                                                                  !
   !   de frequence  0.003556                                                         !
   !                                                                                  !
   !   a une norme d'erreur de  0.000008  supйrieure au seuil admis  0.000001         !
   !                                                                                  !
   !   Conseil: Si vous utiliser la mйthode 'TRI_DIAG' ou                             !
   ! 'SORENSEN' vous                                                                  !
   !   pouvez amйliorer cette norme en:                                               !
   !     - augmentant la valeur de COEF_DIM_ESPACE (par dйfaut valeur 4 pour TRI_DIAG !
   !                                                                                  !
   !      et 2 pour SORENSEN),                                                        !
   !     - rйduire le nombre de valeurs propres recherchйes (NMAX_FREQ ou taille de   !
   !      la BANDE).
Last edit: 13 years 2 weeks ago by bolt.
Moderators: catux
Time to create page: 0.146 seconds
Powered by Kunena Forum