אנליזת מגע (חלק 1)¶
לניתוח זה משתמשים בנתונים מתוך tutorial/09_contact_hertz.
מושא האנליזה¶
הניתוח הוא בעיית מגע הרצית. הגאומטריה של מושא האנליזה מוצגת באיור 4.9.1 ונתוני הרשת באיור 4.9.2.
| פריט | תיאור | הערות | הפניה |
|---|---|---|---|
| סוג האנליזה | אנליזה סטטית לא־ליניארית (אלסטו־פלסטיות, מגע) | !SOLUTION,TYPE=NLSTATIC !CONTACT | |
| מספר צמתים | 408 | ||
| מספר אלמנטים | 168 | ||
| סוג אלמנט | אלמנט הקסהדרלי ליניארי בעל 8 צמתים | !ELEMENT,TYPE=361 | |
| שם חומר | MAT1 | !MATERIAL,NAME=MAT1 | |
| תכונות חומר | ELASTIC | !ELASTIC | |
| תנאי שפה | אילוץ, הזזה מוכתבת | ||
| פותר מטריצה | CG/SSOR | !SOLVER,METHOD=CG,PRECOND=1 |

איור 4.9.1: גאומטריית מושא האנליזה

איור 4.9.2: נתוני הרשת של מושא האנליזה
תוכן האנליזה¶
בשיטת מכפילי לגראנז' המורחבת מתבצעת אנליזת מגע שבה מוחלת הזזה מוכתבת בכיוון לחיצה על המשטח העליון של מודל רבע־דיסקה. נתוני בקרת האנליזה מוצגים להלן.
נתוני בקרת אנליזה cgs3.cnt¶
# קובץ בקרה עבור FISTR
## בקרת אנליזה
!VERSION
3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## בקרת פותר
### תנאי שפה
!BOUNDARY, GRPID=1
ALL, 3, 3, 0.0
BOTTOM, 2, 2, 0.0
CENTER, 1, 1, 0.0
UPPER, 2, 2, -0.306
!CONTACT_ALGO, TYPE=ALAGRANGE
!CONTACT, GRPID=1
CP1, 0.0
### צעד
!STEP, SUBSTEPS=5, CONVERG=1.0e-5
BOUNDARY, 1
CONTACT, 1
### חומר
!MATERIAL, NAME=MAT1
!ELASTIC
1100.0, 0.0
### הגדרות פותר
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES
1000, 1
1.0e-10, 1.0, 0.0
## בקרת עיבוד לאחר
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
נוהל האנליזה¶
מריצים את פקודת ההפעלה של FrontISTR, fistr1.
(הרצה עם 4 תהליכונים)
תוצאות האנליזה¶
עבור תוצאות תת־השלב החמישי, REVOCAP_PrePost משמש ליצירת תצוגת עיוות עם קווי מתאר של ההזזה בכיוון y, המוצגת באיור 4.9.3. בנוסף, כחומר מספרי של תוצאות הניתוח, מוצג להלן חלק מקובץ יומן תוצאות הניתוח.

איור 4.9.3: תוצאות ניתוח של עיוות והזזה בכיוון y
יומן תוצאות האנליזה 0.log¶
fstr_setup: OK
#### Result step= 0
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 0.0000E+00 1 0.0000E+00 1
//U2 0.0000E+00 1 0.0000E+00 1
//U3 0.0000E+00 1 0.0000E+00 1
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Global Summary @Node :Max/IdMax/Min/IdMin####
//U1 0.0000E+00 1 0.0000E+00 1
//U2 0.0000E+00 1 0.0000E+00 1
//U3 0.0000E+00 1 0.0000E+00 1
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 0.0000E+00 1 0.0000E+00 1
//E22 0.0000E+00 1 0.0000E+00 1
//E33 0.0000E+00 1 0.0000E+00 1
//E12 0.0000E+00 1 0.0000E+00 1
//E23 0.0000E+00 1 0.0000E+00 1
//E31 0.0000E+00 1 0.0000E+00 1
//S11 0.0000E+00 1 0.0000E+00 1
//S22 0.0000E+00 1 0.0000E+00 1
//S33 0.0000E+00 1 0.0000E+00 1
//S12 0.0000E+00 1 0.0000E+00 1
//S23 0.0000E+00 1 0.0000E+00 1
//S31 0.0000E+00 1 0.0000E+00 1
//SMS 0.0000E+00 1 0.0000E+00 1
...
#### Result step= 5
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 1.3079E-02 31 -4.1442E-02 54
//U2 4.9597E-03 1008 -3.0604E-01 69
//U3 0.0000E+00 1 0.0000E+00 1
//E11 7.1171E-03 1049 -4.7627E-02 1
//E22 2.0644E-03 1012 -5.3428E-02 50
//E33 9.3301E-03 50 -4.0103E-03 54
//E12 4.3231E-02 1000 -2.7666E-02 54
//E23 1.3028E-12 1046 -1.2515E-12 2
//E31 3.5660E-13 48 -2.1028E-12 51
//S11 7.9455E+00 1049 -4.9987E+01 1
//S22 2.2819E+00 1012 -5.3282E+01 50
//S33 1.1359E+01 50 -4.4937E+00 54
//S12 2.0802E+01 1000 -1.4668E+01 54
//S23 6.7682E-10 1046 -6.3666E-10 2
//S31 2.0235E-10 48 -1.1497E-09 51
//SMS 6.0944E+01 1000 1.0198E-01 69
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.6311E-03 1024 -3.5443E-02 8
//E22 1.7091E-03 1031 -5.4262E-02 8
//E33 1.5903E-04 7 -1.6631E-05 101
//E12 1.8590E-02 1009 -2.6073E-02 14
//E23 1.5204E-13 1000 -4.1521E-13 16
//E31 1.3902E-13 5 -8.2391E-13 7
//S11 5.2516E+00 1024 -3.8099E+01 8
//S22 1.8872E+00 1031 -5.5262E+01 8
//S33 2.2824E-01 7 -1.8263E-02 101
//S12 9.8516E+00 1009 -1.3352E+01 14
//S23 7.8767E-11 1000 -2.1075E-10 16
//S31 7.7553E-11 5 -4.4780E-10 7
//SMS 4.9470E+01 8 2.4786E-01 1007
##### Global Summary @Node :Max/IdMax/Min/IdMin####
//U1 1.3079E-02 31 -4.1442E-02 54
//U2 4.9597E-03 1008 -3.0604E-01 69
//U3 0.0000E+00 1 0.0000E+00 1
//E11 7.1171E-03 1049 -4.7627E-02 1
//E22 2.0644E-03 1012 -5.3428E-02 50
//E33 9.3301E-03 50 -4.0103E-03 54
//E12 4.3231E-02 1000 -2.7666E-02 54
//E23 1.3028E-12 1046 -1.2515E-12 2
//E31 3.5660E-13 48 -2.1028E-12 51
//S11 7.9455E+00 1049 -4.9987E+01 1
//S22 2.2819E+00 1012 -5.3282E+01 50
//S33 1.1359E+01 50 -4.4937E+00 54
//S12 2.0802E+01 1000 -1.4668E+01 54
//S23 6.7682E-10 1046 -6.3666E-10 2
//S31 2.0235E-10 48 -1.1497E-09 51
//SMS 6.0944E+01 1000 1.0198E-01 69
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.6311E-03 1024 -3.5443E-02 8
//E22 1.7091E-03 1031 -5.4262E-02 8
//E33 1.5903E-04 7 -1.6631E-05 101
//E12 1.8590E-02 1009 -2.6073E-02 14
//E23 1.5204E-13 1000 -4.1521E-13 16
//E31 1.3902E-13 5 -8.2391E-13 7
//S11 5.2516E+00 1024 -3.8099E+01 8
//S22 1.8872E+00 1031 -5.5262E+01 8
//S33 2.2824E-01 7 -1.8263E-02 101
//S12 9.8516E+00 1009 -1.3352E+01 14
//S23 7.8767E-11 1000 -2.1075E-10 16
//S31 7.7553E-11 5 -4.4780E-10 7
//SMS 4.9470E+01 8 2.4786E-01 1007