تحليل التلامس (الجزء 1)¶
يستخدم هذا التحليل البيانات الموجودة في tutorial/09_contact_hertz.
نموذج التحليل¶
هذه مسألة تلامس Hertz. يوضح الشكل 4.9.1 هندسة النموذج المراد تحليله، ويوضح الشكل 4.9.2 بيانات الشبكة.
| البند | الوصف | ملاحظات | المرجع |
|---|---|---|---|
| نوع التحليل | تحليل ساكن غير خطي (مرونة لدنة، تلامس) | !SOLUTION,TYPE=NLSTATIC !CONTACT | |
| عدد العقد | 408 | ||
| عدد العناصر | 168 | ||
| نوع العنصر | عنصر سداسي السطوح من الرتبة الأولى بثماني عقد | !ELEMENT,TYPE=361 | |
| اسم المادة | MAT1 | !MATERIAL,NAME=MAT1 | |
| خاصية المادة | ELASTIC | !ELASTIC | |
| الشروط الحدّية | قيود، إزاحة مفروضة | ||
| المحلل الخطي | CG/SSOR | !SOLVER,METHOD=CG,PRECOND=1 |
الشكل 4.9.1 هندسة نموذج التحليل
الشكل 4.9.2 بيانات شبكة نموذج التحليل
وصف التحليل¶
يُجرى تحليل تلامس باستخدام طريقة مضاعف Lagrange المعزَّز، وذلك بتطبيق إزاحة مفروضة في اتجاه الضغط على السطح العلوي لنموذج ربع قرص. وتظهر بيانات التحكم في التحليل أدناه.
بيانات التحكم في التحليل cgs3.cnt¶
# Control File for FISTR
## Analysis Control
!VERSION
3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## Solver Control
### Boundary Condition
!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
!STEP, SUBSTEPS=5, CONVERG=1.0e-5
BOUNDARY, 1
CONTACT, 1
### Material
!MATERIAL, NAME=MAT1
!ELASTIC
1100.0, 0.0
### Solver Setting
!SOLVER,METHOD=CG,PRECOND=1,ITERLOG=YES,TIMELOG=YES
1000, 1
1.0e-10, 1.0, 0.0
## Post Control
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
إجراء التحليل¶
شغّل أمر FrontISTR التنفيذي fistr1.
(تشغيل باستخدام 4 خيوط)
نتائج التحليل¶
بالنسبة إلى الخطوة الفرعية الخامسة، يُنشأ مخطط للشكل المتشوه مع خطوط كنتور لإزاحة اتجاه y باستخدام REVOCAP_PrePost، ويُعرض في الشكل 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