ການວິເຄາະການສຳຜັດ (ສ່ວນ 2)¶
ການວິເຄາະນີ້ໃຊ້ຂໍ້ມູນໃນ tutorial/10_contact_2tubes.
ແບບຈຳລອງວິເຄາະ¶
ການວິເຄາະແມ່ນບັນຫາການກົດກະບອກເຂົ້າຫາກັນ; ຮູບຮ່າງຂອງວັດຖຸວິເຄາະສະແດງໃນຮູບ 4.10.1 ແລະ ຂໍ້ມູນເມຊໃນຮູບ 4.10.2.
| ລາຍການ | ເນື້ອຫາ | ໝາຍເຫດ | ອ້າງອີງ |
|---|---|---|---|
| ປະເພດການວິເຄາະ | ການວິເຄາະສະຖິດແບບບໍ່ເປັນເສັ້ນ(ອີລາສຕິກ,ການສຳຜັດ) | !SOLUTION,TYPE=NLSTATIC !CONTACT | |
| ຈຳນວນໂຫນດ | 4,000 | ||
| ຈຳນວນເອເລເມັນ | 2,888 | ||
| ປະເພດເອເລເມັນ | ເອເລເມັນຫົກໜ້າອັນດັບໜຶ່ງ 8 ໂຫນດ | !ELEMENT,TYPE=361 | |
| ຊື່ຄຸນສົມບັດວັດສະດຸ | M1 | !MATERIAL,NAME=M1 | |
| ຄຸນສົມບັດວັດສະດຸ | ELASTIC | !ELASTIC | |
| ເງື່ອນໄຂຂອບເຂດ | ການຍຶດ, ການກະຈັດທີ່ກຳນົດ | ||
| ວິທີແກ້ເມທຣິກ | ວິທີໂດຍກົງ | !SOLVER,METHOD=MUMPS |
ຮູບ 4.10.1 ຮູບຮ່າງຂອງວັດຖຸວິເຄາະ
ຮູບ 4.10.2 ຂໍ້ມູນເມຊຂອງວັດຖຸວິເຄາະ
ເນື້ອຫາການວິເຄາະ¶
ດຳເນີນການວິເຄາະການສຳຜັດດ້ວຍວິທີຕົວຄູນ Lagrange ໂດຍໃຫ້ການກະຈັດທີ່ກຳນົດໃນທິດກົດແກ່ພື້ນຜິວບັງຄັບທີ່ສະແດງໃນຮູບ 4.10.1. ຂໍ້ມູນຄວບຄຸມການວິເຄາະສະແດງດ້ານລຸ່ມ.
ຂໍ້ມູນຄວບຄຸມການວິເຄາະ 2tubes.cnt¶
# ໄຟລ໌ຄວບຄຸມສຳລັບ FISTR
## ການຄວບຄຸມການວິເຄາະ
!VERSION
3
!SOLUTION, TYPE=NLSTATIC
!WRITE,RESULT
!WRITE,VISUAL
## ການຄວບຄຸມຕົວແກ້ສົມຜົນ
### ເງື່ອນໄຂຂອບເຂດ
!BOUNDARY, GRPID=1
X0, 1, 3, 0.0
Y0, 2, 2, 0.0
Z0, 3, 3, 0.0
!BOUNDARY, GRPID=2
X1, 1, 1, 0.0
!BOUNDARY, GRPID=3
X1, 1, 1, -1.0
!CONTACT_ALGO, TYPE=SLAGRANGE
!CONTACT, GRPID=1, INTERACTION=FSLID, NPENALTY=1.0e+2
CP1, 0.0, 1.0e+5
### STEP
!STEP, SUBSTEPS=4, CONVERG=1.0e-5
BOUNDARY, 1
BOUNDARY, 3
CONTACT, 1
### ວັດສະດຸ
!MATERIAL, NAME=M1
!ELASTIC
2.1e+5, 0.3
### ການຕັ້ງຄ່າຕົວແກ້ສົມຜົນ
!SOLVER,METHOD=MUMPS
## ການຄວບຄຸມຫຼັງການປະມວນຜົນ
!VISUAL,method=PSR
!surface_num=1
!surface 1
!output_type=VTK
!END
ຂັ້ນຕອນການວິເຄາະ¶
ເອີ້ນໃຊ້ຄຳສັ່ງ fistr1 ຂອງ FrontISTR.
(ເອີ້ນໃຊ້ດ້ວຍ 4 ເທຣດ)
ຜົນການວິເຄາະ¶
ສຳລັບຜົນຂອງຂັ້ນຍ່ອຍທີ 4, ຮູບປ່ຽນຮູບທີ່ຊ້ອນຄອນທົວຄວາມເຄັ້ນ von Mises ສ້າງດ້ວຍ REVOCAP_PrePost ແລະ ສະແດງໃນຮູບ 4.10.3. ນອກຈາກນັ້ນ ບາງສ່ວນຂອງໄຟລ໌ບັນທຶກຜົນການວິເຄາະສະແດງດ້ານລຸ່ມເປັນຂໍ້ມູນຕົວເລກ.
ຮູບ 4.10.3 ຜົນການວິເຄາະການປ່ຽນຮູບ ແລະ ຄວາມເຄັ້ນ von Mises
ບັນທຶກຜົນການວິເຄາະ 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= 4
##### Local Summary @Node :Max/IdMax/Min/IdMin####
//U1 8.7136E-04 32 -1.0021E+00 2006
//U2 8.7633E-03 104 -7.0546E-03 2006
//U3 8.7634E-03 4 -7.0546E-03 1901
//E11 7.5985E-04 1901 -4.1250E-04 105
//E22 9.8385E-04 2 -9.3128E-04 2058
//E33 9.8386E-04 102 -9.3120E-04 3843
//E12 5.3544E-04 133 -2.8361E-04 278
//E23 1.2509E-03 1901 -1.4195E-03 4
//E31 5.3554E-04 33 -2.8366E-04 1678
//S11 7.7172E+01 103 -8.9880E+01 101
//S22 2.0158E+02 2 -2.2918E+02 1905
//S33 2.0159E+02 102 -2.2920E+02 2010
//S12 4.3210E+01 133 -2.2913E+01 278
//S23 1.0057E+02 1901 -1.1446E+02 4
//S31 4.3219E+01 33 -2.2918E+01 1678
//SMS 3.0051E+02 1901 3.1670E+00 2454
##### Local Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.3715E-04 2809 -3.1965E-04 1293
//E22 8.8407E-04 73 -7.9864E-04 1369
//E33 8.8408E-04 1365 -7.9867E-04 2809
//E12 5.3226E-04 1324 -1.9554E-04 188
//E23 1.0836E-03 2149 -1.0974E-03 684
//E31 5.3223E-04 32 -1.9559E-04 1176
//S11 6.6688E+01 1296 -7.9313E+01 1293
//S22 1.8533E+02 73 -2.1547E+02 1369
//S33 1.8533E+02 1365 -2.1548E+02 2809
//S12 4.2951E+01 1324 -1.5800E+01 188
//S23 8.7338E+01 2149 -8.8751E+01 684
//S31 4.2949E+01 32 -1.5804E+01 1176
//SMS 1.9064E+02 1369 3.2782E+00 2126
##### Global Summary @Node :Max/IdMax/Min/IdMin####
//U1 8.7136E-04 32 -1.0021E+00 2006
//U2 8.7633E-03 104 -7.0546E-03 2006
//U3 8.7634E-03 4 -7.0546E-03 1901
//E11 7.5985E-04 1901 -4.1250E-04 105
//E22 9.8385E-04 2 -9.3128E-04 2058
//E33 9.8386E-04 102 -9.3120E-04 3843
//E12 5.3544E-04 133 -2.8361E-04 278
//E23 1.2509E-03 1901 -1.4195E-03 4
//E31 5.3554E-04 33 -2.8366E-04 1678
//S11 7.7172E+01 103 -8.9880E+01 101
//S22 2.0158E+02 2 -2.2918E+02 1905
//S33 2.0159E+02 102 -2.2920E+02 2010
//S12 4.3210E+01 133 -2.2913E+01 278
//S23 1.0057E+02 1901 -1.1446E+02 4
//S31 4.3219E+01 33 -2.2918E+01 1678
//SMS 3.0051E+02 1901 3.1670E+00 2454
##### Global Summary @Element :Max/IdMax/Min/IdMin####
//E11 4.3715E-04 2809 -3.1965E-04 1293
//E22 8.8407E-04 73 -7.9864E-04 1369
//E33 8.8408E-04 1365 -7.9867E-04 2809
//E12 5.3226E-04 1324 -1.9554E-04 188
//E23 1.0836E-03 2149 -1.0974E-03 684
//E31 5.3223E-04 32 -1.9559E-04 1176
//S11 6.6688E+01 1296 -7.9313E+01 1293
//S22 1.8533E+02 73 -2.1547E+02 1369
//S33 1.8533E+02 1365 -2.1548E+02 2809
//S12 4.2951E+01 1324 -1.5800E+01 188
//S23 8.7338E+01 2149 -8.8751E+01 684
//S31 4.2949E+01 32 -1.5804E+01 1176
//SMS 1.9064E+02 1369 3.2782E+00 2126