အကြောင်းအရာသို့ ကရန်

အီလီမင့် နံပါတ်ပေးစနစ်နှင့် ပုံသဏ္ဌာန်ဖန်ရှင် လိုင်ဘရီ

ပုံသဏ္ဌာန်ဖန်ရှင်များနှင့် ဖိုင်နိုက်အီလီမင့် ခန့်မှန်းမှု တွင် မိတ်ဆက်ထားသော ပုံသဏ္ဌာန်ဖန်ရှင် \(N_\alpha^e(\boldsymbol{r})\) သည် အီလီမင့်အမျိုးအစားတစ်ခုစီအတွက် သတ်သတ်မှတ်မှတ် ပုံစံရှိသည်။ ဤအခန်းတွင် FrontISTR က ပံ့ပိုးသော အီလီမင့်အုပ်စုများကို သုံးလုံးကိန်း အီလီမင့်အမျိုးအစားနံပါတ် ဖြင့် အမျိုးအစားခွဲပြီး၊ အီလီမင့်အမျိုးအစားတိုင်းကို တူညီသောနည်းဖြင့် ကိုင်တွယ်နိုင်ရန် ပုံသဏ္ဌာန်ဖန်ရှင် လိုင်ဘရီ၏ အင်တာဖေ့စ် ကို ဖော်ပြသည်။ သုံးဖက်မြင် အစိုင်အခဲအီလီမင့် များ၏ သဘာဝကိုဩဒိနိတ်စနစ် နှင့် ပုံသဏ္ဌာန်ဖန်ရှင် များအပြင် ဒုတိယအဆင့် အီလီမင့် များ၏ ဘေးအလယ်-ဆုံမှတ် အစဉ်ပြောင်းလဲမှုကို သုံးဖက်မြင် အစိုင်အခဲအီလီမင့် များ၏ ပုံသဏ္ဌာန်ဖန်ရှင် များ တွင် ဖော်ပြထားပြီး၊ အီလီမင့်အမျိုးအစား ရွေးချယ်မှုဆိုင်ရာ လမ်းညွှန်ကို အီလီမင့် လိုင်ဘရီ (ဖန်ရှင်များ) တွင် ဖော်ပြထားသည်။

အီလီမင့် အမျိုးအစားခွဲခြားမှုနှင့် အမည်ပေးစည်းမျဉ်းများ

FrontISTR အီလီမင့်အမျိုးအစားများကို သုံးလုံးကိန်း နံပါတ်များဖြင့် ခွဲခြားသတ်မှတ်သည်။ ကိန်းလုံးတစ်လုံးစီ၏ အဓိပ္ပာယ်မှာ အောက်ပါအတိုင်းဖြစ်သည်။

  • ပထမကိန်းလုံး (အီလီမင့်အုပ်စု): 1 = မျဉ်းအီလီမင့်/ထရပ်စ်၊ 2 = ပြင်ညီအီလီမင့်၊ 3 = အစိုင်အခဲအီလီမင့်၊ 5 = အင်တာဖေ့စ်အီလီမင့်၊ 6 = ဘီးမ် အီလီမင့်၊ 7 = ရှဲလ်အီလီမင့်။
  • ဒုတိယကိန်းလုံး (အခြေခံပုံသဏ္ဌာန်): အီလီမင့်အုပ်စုအတွင်းရှိ ဂျီဩမေတြီပုံသဏ္ဌာန် (တြိဂံ၊ စတုဂံ၊ တက်ထရာဟီဒရွန်၊ ပရစ်ဇမ် သို့မဟုတ် ဟက်ဇာဟီဒရွန်) ကို ခွဲခြားသည်။
  • တတိယကိန်းလုံး (ကြားဖြည့်တွက်ချက်မှု အဆင့်): 1 = ပထမအဆင့် (ထောင့်နိုဒ်များသာ)၊ 2 = ဒုတိယအဆင့် (ဘေးအလယ်ဆုံမှတ် များပါဝင်)။ ခြွင်းချက်အဖြစ် 2-ဆုံမှတ် ထရပ်စ် အီလီမင့် သည် နံပါတ် 301 ဖြစ်သည်။

ပံ့ပိုးထားသော အီလီမင့်များကို အောက်ပါဇယားတွင် အီလီမင့်အုပ်စုအလိုက် ဖော်ပြထားသည်။

အီလီမင့်အုပ်စု အီလီမင့်အမျိုးအစား နိုဒ်အရေအတွက် ဖော်ပြချက်
မျဉ်း အီလီမင့် 111 2 2-ဆုံမှတ် လင့်ခ် အီလီမင့် (စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှု တွင် အသုံးမပြုနိုင်)
112 3 3-ဆုံမှတ် လင့်ခ် အီလီမင့် (စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှု တွင် အသုံးမပြုနိုင်)
ပြင်ညီ အီလီမင့် 231 3 3-ဆုံမှတ် ပထမအဆင့် တြိဂံ အီလီမင့်
232 6 6-ဆုံမှတ် ဒုတိယအဆင့် တြိဂံ အီလီမင့်
241 4 4-ဆုံမှတ် ပထမအဆင့် စတုဂံ အီလီမင့်
242 8 8-ဆုံမှတ် ဒုတိယအဆင့် စတုဂံ အီလီမင့် (Serendipity မျိုးစု)
ဆော်လစ် အီလီမင့် 301 2 2-ဆုံမှတ် ထရပ်စ် အီလီမင့်
341 4 4-ဆုံမှတ် ပထမအဆင့် လေးမျက်နှာ အီလီမင့်
342 10 10-ဆုံမှတ် ဒုတိယအဆင့် လေးမျက်နှာ အီလီမင့်
351 6 6-ဆုံမှတ် ပထမအဆင့် ပရစ်ဇမ် အီလီမင့်
352 15 15-ဆုံမှတ် ဒုတိယအဆင့် ပရစ်ဇမ် အီလီမင့်
361 8 8-ဆုံမှတ် ပထမအဆင့် ခြောက်မျက်နှာ အီလီမင့်
362 20 20-ဆုံမှတ် ဒုတိယအဆင့် ခြောက်မျက်နှာ အီလီမင့် (Serendipity မျိုးစု)
အင်တာဖေ့စ် အီလီမင့် 541 4×2 4-ဆုံမှတ် ပထမအဆင့် စတုဂံ-မျက်နှာ အီလီမင့် (စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှုတွင် အသုံးမပြုနိုင်)
542 8×2 8-ဆုံမှတ် ဒုတိယအဆင့် စတုဂံ-မျက်နှာ အီလီမင့် (စထရက်စ် ခွဲခြမ်းစိတ်ဖြာမှုတွင် အသုံးမပြုနိုင်)
ဘီးမ် အီလီမင့် 611 2 2-ဆုံမှတ် Bernoulli-Euler ဘီးမ် အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 6 ခု)
641 2×2 2-ဆုံမှတ် Bernoulli-Euler ဘီးမ် အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 3 ခု၊ ရောစပ် လွတ်လပ်မှုဒီဂရီ အတွက်)
ရှဲလ် အီလီမင့် 731 3 MITC3 တြိဂံ ရှဲလ်အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 6 ခု)
741 4 MITC4 စတုဂံ ရှဲလ်အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 6 ခု)
743 9 MITC9 စတုဂံ ရှဲလ်အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 6 ခု)
761 3×2 MITC3 တြိဂံ ရှဲလ်အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 3 ခု၊ ရောစပ် လွတ်လပ်မှုဒီဂရီ အတွက်)
781 4×2 MITC4 စတုဂံ ရှဲလ်အီလီမင့် (နိုဒ်တစ်ခုလျှင် လွတ်လပ်မှုဒီဂရီ 3 ခု၊ ရောစပ် လွတ်လပ်မှုဒီဂရီ အတွက်)

ပထမအဆင့် ခြောက်မျက်နှာ အီလီမင့် 361 တွင် ဖော်မြူလေးရှင်း အမျိုးမျိုး (အပြည့်အဝ အင်တီဂရေးရှင်း၊ မကိုက်ညီမုဒ်၊ B-bar နှင့် F-bar) ရှိပြီး ခွဲခြမ်းစိတ်ဖြာမှု ထိန်းချုပ်ချက်မှ ရွေးချယ်သည်။ အသေးစိတ်ကို စွမ်းဆောင်ရည်မြင့် အီလီမင့် ဖော်မြူလာရေးဆွဲမှု တွင် ကြည့်ပါ။

ဘီးမ် အီလီမင့် 641 နှင့် ရှဲလ်အီလီမင့် 761/781 တို့သည် တူညီသောနိုဒ်များတွင် လွတ်လပ်မှုဒီဂရီ 6 ခုရှိ ဖွဲ့စည်းပုံဆိုင်ရာ အီလီမင့် များကို လွတ်လပ်မှုဒီဂရီ 3 ခုရှိ အစိုင်အခဲအီလီမင့် များနှင့် ချိတ်ဆက်ရန် ဖော်မြူလေးရှင်း များဖြစ်သည်။ ဖွဲ့စည်းပုံဆိုင်ရာ ဆုံမှတ် တစ်ခုကို 3-လွတ်လပ်မှုဒီဂရီ ဆုံမှတ် နှစ်ခုဖြင့် ကိုယ်စားပြုသည်။ ဖော်မြူလေးရှင်း ကိုယ်တိုင်သည် သက်ဆိုင်ရာ 6-လွတ်လပ်မှုဒီဂရီ ဗားရှင်းများ (611၊ 731၊ 741) နှင့် တူညီသည်။

ပုံသဏ္ဌာန်ဖန်ရှင် လိုင်ဘရီ

အီလီမင့်တွက်ချက်ရေး လုပ်ငန်းစဉ် များ (တောင့်တင်းမှုမက်ထရစ် နှင့် အတွင်းအားဗက်တာ တည်ဆောက်ခြင်း) တွင် အီလီမင့်အမျိုးအစား တစ်ခုစီအတွက် သီးခြား လုပ်ငန်းစဉ် မရှိပါ။ ထိုအစား အီလီမင့်အမျိုးအစားနံပါတ် fetype ကို အာဂျူမင့် အဖြစ် လက်ခံသော တစ်ပြေးညီ အင်တာဖေ့စ် မှတစ်ဆင့် ပုံသဏ္ဌာန်ဖန်ရှင် လိုင်ဘရီကို ခေါ်ယူသည်။ elementInfo မော်ဂျူး (fistr1/src/lib/element/element.f90) သည် အောက်ပါ ဖန်ရှင် များကို ပံ့ပိုးသည်။

  • NumOfQuadPoints(fetype): အီလီမင့်အမျိုးအစားအတွက် အင်တီဂရေးရှင်းအမှတ် အရေအတွက် \(n_q\) ကို ပြန်ပေးသည်။
  • getQuadPoint(fetype, np, pos): အင်တီဂရေးရှင်းအမှတ် \(i\) ၏ သဘာဝကိုဩဒိနိတ် \(\boldsymbol{r}_i\) ကို ပြန်ပေးသည်။
  • getWeight(fetype, np): အင်တီဂရေးရှင်းအမှတ် \(i\) ၏ အလေးချိန် \(w_i\) ကို ပြန်ပေးသည်။
  • getShapeFunc(fetype, r, N): သဘာဝကိုဩဒိနိတ် \(\boldsymbol{r}\) ရှိ ပုံသဏ္ဌာန်ဖန်ရှင် တန်ဖိုး \(N_\alpha^e(\boldsymbol{r})\) ကို ပြန်ပေးသည်။
  • getShapeDeriv(fetype, r, dN): ပုံသဏ္ဌာန်ဖန်ရှင် များ၏ သဘာဝကိုဩဒိနိတ် ဆင်းသက်ချက် \(\partial N_\alpha^e/\partial \boldsymbol{r}\) ကို ပြန်ပေးသည်။
  • getGlobalDeriv(fetype, n_e, r, X^e, J, dN_x): အီလီမင့်ဆုံမှတ် ကိုဩဒိနိတ် \(\boldsymbol{X}^e\) ကို ပေးထားလျှင် ရုပ်ပိုင်းဆိုင်ရာ ကိုဩဒိနိတ် ဆင်းသက်ချက် \(\partial N_\alpha^e/\partial \boldsymbol{x}\) နှင့် ယာကိုဘီယန် ဒီတာမီနန့် \(J\) ကို ပြန်ပေးသည်။
  • getNumberOfSubface(fetype), getSubFace(fetype, k, ...): အီလီမင့် နယ်နိမိတ်မျက်နှာပြင် အရေအတွက်နှင့် မျက်နှာပြင်တစ်ခုစီကို ဖွဲ့စည်းသော ဆုံမှတ်နံပါတ် များကို ပြန်ပေးသည်။ မျက်နှာပြင်ဝန် နှင့် ထိတွေ့မျက်နှာပြင် များတည်ဆောက်ရာတွင် အသုံးပြုသည်။

လုပ်ဆောင်ချက် တစ်ခုစီသည် အီလီမင့်အမျိုးအစားနံပါတ် အလိုက် အတွင်းပိုင်းတွင် ခွဲခွာမှု ခွဲပြီး အီလီမင့်အမျိုးအစား တစ်ခုစီအတွက် ပံ့ပိုးထားသော အောက်အဆင့် မော်ဂျူး (tet4n.f90, tet10n.f90, prism6n.f90, prism15n.f90, hex8n.f90, hex20n.f90, tri3n.f90, tri6n.f90, quad4n.f90, quad8n.f90, quad9n.f90, line2n.f90, line3n.f90) ရှိ ပုံသဏ္ဌာန်ဖန်ရှင် နှင့် ပုံသဏ္ဌာန်ဖန်ရှင် ဆင်းသက်ချက် ဆပ်ရူတင်း များကို ခေါ်သည်။ အင်တီဂရေးရှင်းအမှတ် ကိုဩဒိနိတ် နှင့် အလေးချိန် များ၏ ကိန်းဂဏန်းဇယားများကို quadrature.f90 တွင် စုစည်းထားသည်။ အီလီမင့်အမျိုးအစား အသစ် ထည့်ရန် သက်ဆိုင်ရာ အောက်အဆင့် မော်ဂျူး ကို အကောင်အထည်ဖော်ပြီး elementInfo ရှိ ဖန်ရှင် တစ်ခုစီတွင် ခွဲခွာမှု ထည့်ရမည်။

ဒုတိယအဆင့် တြိဂံအီလီမင့် များအတွက် ဘေးအလယ်ဆုံမှတ် အစဉ်ပြောင်းလဲခြင်း

ဒုတိယအဆင့် တြိဂံ အီလီမင့် 232 အတွက် ဘေးအလယ်ဆုံမှတ် အစဉ်သည် HEC-MW မက်ရှ် အဝင်ဒေတာ ဖော်မတ် (အီလီမင့် ပတ်လည် ထောင့်ဆုံမှတ် များကို အရင်စီပြီးနောက် အစွန်း 1-2၊ 2-3၊ … တစ်လျှောက် ဘေးအလယ်ဆုံမှတ် များကို စီသောပုံစံ) နှင့် အီလီမင့် လိုင်ဘရီ ရှိ ပုံသဏ္ဌာန်ဖန်ရှင် များက ယူဆထားသော FrontISTR အတွင်းပိုင်း ဖော်မတ် ကြား ကွာခြားသည်။ မက်ရှ်ကို ဖတ်ပြီးချင်း hecmw2fstr_mesh_conv (fistr1/src/common/hecmw2fstr_mesh_conv.f90 နှင့် hecmw2fstr_connect_conv.c) က ပြောင်းလဲမှုဇယား Table232 အတိုင်း ဆုံမှတ် ခြောက်ခုကို \(\{1,2,3,6,4,5\}\) သို့ ပြန်စီခြင်းဖြင့် ဖော်မတ် နှစ်ခုကြား ကွာခြားမှုကို ဖြေရှင်းသည်။ ပြောင်းလဲမှုကို အီလီမင့် ဆက်သွယ်မှုအစီအစဉ် သိမ်းထားသော အာရေး ပေါ်တွင် တိုက်ရိုက်ပြုလုပ်သဖြင့် နောက်ပိုင်း အီလီမင့် လိုင်ဘရီ ခေါ်ယူမှုအားလုံးသည် FrontISTR အတွင်းပိုင်း ဆုံမှတ် အစဉ် ကို အသုံးပြုသည်။ ရလဒ် အထွက်ဒေတာ တွင် ပြောင်းပြန် ပြောင်းလဲမှု fstr2hecmw_mesh_conv က ရေးသားမီ HEC-MW ဖော်မတ် သို့ ပြန်လည်ထားပေးသဖြင့် အသုံးပြုသူများအနေဖြင့် အဝင်ဒေတာ နှင့် အထွက်ဒေတာ နှစ်မျိုးစလုံးအတွက် HEC-MW ဖော်မတ် ကိုသာ သိရှိထားရန် လိုအပ်သည်။

ဒုတိယအဆင့် သုံးဖက်မြင် အစိုင်အခဲအီလီမင့် 342 (လေးမျက်နှာ) နှင့် 352 (ပရစ်ဇမ်) တွင် အလားတူ အလိုအလျောက် ပြောင်းလဲမှု ရှိသော်လည်း 362 (ခြောက်မျက်နှာ ဒုတိယအဆင့်) တွင် ပြောင်းလဲမှု မရှိပါ။ ယင်းတို့၏ ကိုင်တွယ်ပုံကို သုံးဖက်မြင် အစိုင်အခဲအီလီမင့် များ၏ ပုံသဏ္ဌာန်ဖန်ရှင် များ တွင် ကြည့်ပါ။ ဒုတိယအဆင့် ပြင်ညီစတုဂံ အီလီမင့် 242 နှင့် အင်တာဖေ့စ်အီလီမင့် 542 အတွက် HEC-MW အဝင်ဒေတာ ဖော်မတ် နှင့် အတွင်းပိုင်း ဖော်မတ် ရှိ ဘေးအလယ်ဆုံမှတ် အစဉ် တူညီသောကြောင့် ပြောင်းလဲမှု မပြုလုပ်ပါ။

ဆက်စပ်အကြောင်းအရာများ