MODE_ITER_SIMULT
- MGolbs
-
Topic Author
- Offline
- Platinum Member
-
ich möchte gern mal etwas an Grundverständnis für die Modalanalyse und die möglichen Optionen in Code Aster gewinnen.
Basis sind ja die zwei Matritzen der Steigikeit und der Masse. Aber schon bei den Methoden fehlt mir der Überblick. Derzeit verwenden wir in Nastran Real Solution Method Lanczos. In Code Aster gibt es Lanczos so erst einmal nicht (nur im Unterpunkt Tri_Diag), aber eben Tri_Diag, Jacobi, Sorensen und QZ. MODE_ITER_SIMULT
Im Kern möchte ich die ersten 20 Moden (ohne RBM), diese bis maximal 200Hz rechnen. Da sollte FREQ=(1,200,), helfen. NMAX_FREQ=10 scheint nur mit BANDE nicht zu gehen.
Wann sollte man welches Verfahren einsetzen - Tri_Diag, Jacobi, Sorensen und QZ? Ist Sorensen die Wahl für große Stahlbauteile?
Habe mit Sorensen die besten Erfahrungen gemacht. Unser Code-Aster Dienstleister arbeitet auch mit Tri_Diag für 1D+2D Elementmodelle. Tri_DIAG hat ja intern auch was von LANCZOS, würde Nastranergebnisse wohl am besten entsprechen?
Als Basis habe ich vorliegen:
[python]
MODE_ITER_SIMULT(MATR_A=stiff,
MATR_B=mass,
METHODE='TRI_DIAG',
CALC_FREQ=_F(OPTION='BANDE',
FREQ=(2,20,),
COEF_DIM_ESPACE=2,),
VERI_MODE=_F(SEUIL=1e-03,),);
[/python]der Solver arbeitet dann mit:
[python]
MODE_ITER_SIMULT(OPTION='SANS',
PREC_LANCZOS=1.E-08,
NMAX_ITER_QR=30,
STOP_FREQ_VIDE='OUI',
NMAX_ITER_ORTHO=5,
PREC_ORTHO=9.9999999999999998E-13,
TYPE_RESU='DYNAMIQUE',
INFO=1,
CALC_FREQ=_F(OPTION='BANDE',
PREC_SHIFT=0.050000000000000003,
APPROCHE='REEL',
NMAX_ITER_SHIFT=5,
COEF_DIM_ESPACE=2,
FREQ=(2,20,),
SEUIL_FREQ=1.E-2),
VERI_MODE=_F(PREC_SHIFT=5.0000000000000001E-3,
STOP_ERREUR='OUI',
STURM='OUI',
SEUIL=1.E-3),
SOLVEUR=_F(RENUM='METIS',
METHODE='MULT_FRONT',
NPREC=8),
METHODE='TRI_DIAG',
MATR_B=mass,
MATR_A=stiff,
);
[/python]Habe mit den Einstellungen folgenden Error:
JDC.py : ERREUR A L'INTERPRETATION DANS ACCAS - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR phase d'initialisation
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! -> La borne minimale de la bande de fréquences est une valeur propre ! !
! Malgré la stratégie de décalage, la matrice de raideur est numériquement !
! singulière (modes de corps rigide). !
! -> Risque & Conseil : !
! Augmenter (ou diminuer) la fréquence (ou la charge critique dans le cas du !
! calcul de !
! flambement) qui définit la borne minimale de la bande de fréquence. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR phase d'initialisation
LE NOMBRE DE DDL
TOTAL EST: 318462
DE LAGRANGE EST: 48
LE NOMBRE DE DDL ACTIFS EST: 318390
DE LAGRANGE EST: 48, habe ich mit 48 RBM's zu rechnen? Eigentlich sollten diese "Freiheitsgrade???" von der Gesamtstruktur (Netzaufbau) nicht als RMB's erscheinen.
Über Tipps und Infos würde ich mich freuen.
Gruß Markus
Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège
Please Log in or Create an account to join the conversation.
- RichardS
- Offline
- Elite Member
-
Grüße,
Richard
SimScale - Engineering Simulation in your browser!
Please Log in or Create an account to join the conversation.
- MGolbs
-
Topic Author
- Offline
- Platinum Member
-
ich habe nun ein Ergebnis:
LE NOMBRE DE DDL
TOTAL EST: 318564
DE LAGRANGE EST: 150
LE NOMBRE DE DDL ACTIFS EST: 318339
L'OPTION CHOISIE EST: BANDE
VERIFICATION DU SPECTRE DE FREQUENCES (SUITE DE STURM)
LE NOMBRE DE FREQUENCES DANS LA BANDE ( 2.00000E+00, 2.00000E+01) EST 9
VALEUR_MIN EN FREQUENCE EST : 2.00000E+00
VALEUR_MAX EN FREQUENCE EST : 2.00000E+01
LA VALEUR DE DECALAGE EN FREQUENCE EST : 1.42127E+01
INFORMATIONS SUR LE CALCUL DEMANDE:
NOMBRE DE MODES RECHERCHES : 9
LA DIMENSION DE L'ESPACE REDUIT EST : 0
ELLE EST INFERIEURE AU NOMBRE DE MODES, ON LA PREND EGALE A
18
LES FREQUENCES CALCULEES INF. ET SUP. SONT:
FREQ_INF : 7.58019E+00
FREQ_SUP : 1.99295E+01
LA PREMIERE FREQUENCE INFERIEURE NON RETENUE EST: -1.27084E-03
LA PREMIERE FREQUENCE SUPERIEURE NON RETENUE EST: 2.06864E+01
CALCUL MODAL: METHODE D'ITERATION SIMULTANEE
METHODE DE LANCZOS
NUMERO FREQUENCE (HZ) NORME D'ERREUR ITER_QR
2 7.58019E+00 3.20465E-04 6
3 8.33254E+00 3.56920E-03 6
4 9.52931E+00 8.28950E-05 6
5 1.33804E+01 3.73212E-06 6
6 1.47234E+01 2.86830E-06 6
7 1.51906E+01 8.23823E-07 6
8 1.54184E+01 9.81097E-07 6
9 1.68361E+01 6.67754E-06 6
10 1.99295E+01 2.59162E-04 6
NORME D'ERREUR MOYENNE: 0.47187E-03
Meine Probleme waren etwas vielfältig.
# Ich habe in den Gruppenzuordnungen des IDEAS Files Modifikationen vornehmen müssen, wie die Gruppen auf ein notwendiges Minimum reduziert.
# In den Gelenkdefinitionen war ich noch unvollständig "...LIAISON_DDL=(_F(GROUP_NO=(..."
auch habe ich "VERI_MODE=_F(SEUIL=1e-01,),); von 1e-03 auf 1e-01 erhöht, ohne zu wissen was ich mir dabei einhandeln kann. Bei 1e-03 wurden die Frequenzen zwar schon gerechnet, aber kein rmed angelegt und mit Error der Solverlauf beendet.
Gruß und Dank Markus
Dem Überflüssigen nachlaufen, heißt das Wesentliche verpassen.
Jules Saliège
Please Log in or Create an account to join the conversation.