یک مسئله وابسته به زمان را با استفاده از گام های زمانی ضمنی یا صریح حل کنید.
نحو
model.sol(sname).create(fname”Time”)
model.sol(sname).feature(fname).set(pname,value)
شرح
ویژگی عملیات.
فاصله زمانی و مقادیر زمانی میانی ممکن در ویژگی Tlist آورده شده است . زمان های خروجی توسط ویژگی Tout کنترل می شود .
ویژگی Time خواص و مقادیر زیر را می پذیرد:
ویژگی | ارزش های | پیش فرض | شرح |
algebraicsolveeverydt | دو برابر مثبت | 0 | مقدار دوره Δ t، زمانی که روش حل جبری برای روش های آدامز-باشفورث روی dt تنظیم شده است . |
algebraicsolveeverydtRK | دو برابر مثبت | 0 | مقدار دوره Δ t، زمانی که روش حل جبری برای روش های Runge-Kutta روی dt تنظیم شده است . |
جبری حلواورین | عدد صحیح مثبت | 1 | مقدار N ، زمانی که روش حل جبری برای روش های آدامز-باشفورث روی n تنظیم شده است . |
جبريحلveeverynRK | عدد صحیح مثبت | 1 | مقدار N ، زمانی که روش حل جبری روی n برای روش های Runge-Kutta تنظیم شده است . |
روش حل جبری | n | dt | n | استراتژی مورد استفاده برای حل معادلات جبری: هر مرحله N یا به صورت دوره ای برای روش های آدامز-باشفورث. |
روش حل جبری RK | n | dt | n | استراتژی مورد استفاده برای حل معادلات جبری: هر مرحله N یا به صورت دوره ای برای روش های Runge-Kutta. |
جزیره مرجانی | رشته | خالی | تحمل مطلق در هر زمینه. زیر را ببینید. |
روش اتل | رشته | خالی | نحوه تفسیر مقدار atolfields زیر را ببینید. |
اتلگلوبال | اسکالر مثبت | 1e-3 | تحمل مطلق جهانی، اگر روش atolglobalvalue دستی باشد . |
اتلگلوبال فاکتور | اسکالر مثبت | 0.1 | تحمل مطلق جهانی به عنوان عاملی از تحمل نسبی، اگر روش atolglobalvalue فاکتور باشد . |
روش atolglobal | مقیاس شده | بدون مقیاس | مقیاس شده | نحوه تفسیر مقدار اتلگلوبال |
روش اتلگلوبال | عامل | کتابچه راهنمای | عامل | از ضریب تحمل نسبی یا مقدار تعریف شده توسط کاربر برای تحمل مطلق استفاده کنید. |
بوی | رشته | خالی | تحمل مطلق برای مشتقات زمانی در هر میدان. فقط در صورتی قابل استفاده است که آتولودوتاکتیو روشن باشد . زیر را ببینید. |
اتولودوتاکتیو | رشته | خالی | برای فعال کردن مشخصات دستی تحمل مطلق برای مشتقات زمانی استفاده می شود. زیر را ببینید. |
bdforder | 1-5 | 2 | سفارش BDF برای تنظیمات دستی BDF. |
bwinitstepfrac | اسکالر مثبت | 0.001 | کسری از مرحله اولیه، زمانی که سازگار به bweuler تنظیم شده است . |
bwinitfactor | اسکالر مثبت | 20 | ضریب ایمنی مورد استفاده در خاتمه معادله جبری، زمانی که سازگار بر روی bweuler تنظیم شود . |
اعتبار چک | واقعی | نادرست | نادرست | اگر odesolvertype روی صریح تنظیم شده باشد ، اعتبار سیستم جفت شده را بررسی کنید . |
مجتمع | در | خاموش | خاموش | اجازه دادن به اعداد مختلط |
استوار | خاموش | در | قلع و قمع | قلع و قمع | اولیه سازی مداوم سیستم های DAE. |
فهرست | آرایه رشته ای | با استفاده از یک آرایه رشته ای مقادیر ثابت را به عنوان پارامترهای ورودی ارائه دهید. برای نام های ثابت مربوطه، از cname استفاده کنید . مشخصات clist و cname را ببینید . | |
cname | آرایه رشته ای | با استفاده از یک آرایه رشته ای نام ثابت ها را به عنوان پارامترهای ورودی ارائه دهید. برای مقادیر ثابت مربوطه، از clist استفاده کنید . مشخصات clist و cname را ببینید . | |
کنترل | رشته | کاربر | نام مرحله مطالعه یا کاربر کنترل کننده در صورتی که ویژگی به صورت دستی کنترل شود. |
دوپریگرومکس | اسکالر مثبت | 10 | حداکثر نسبت رشد اندازه گام برای Dormand–Prince 5. |
دوپرگرومین | اسکالر مثبت | 0.2 | حداقل نسبت رشد اندازه گام برای Dormand–Prince 5. |
کنترل دوپری | صاف | سریع | معلول | صاف | رفتار کنترل کنترلر متناسب-انتگرال که میرایی را در تغییرات اندازه پله برای Dormand-Prince 5 اضافه می کند. |
doprisafe | اسکالر مثبت | 0.9 | ضریب ایمنی اندازه پله برای Dormand–Prince 5. |
درون یابی پایان زمان | بولی | درست است، واقعی | آخرین بار مشخص شده در لیست زمان های خروجی را درون یابی کنید، اگر درست است . اگر روی false تنظیم شود ، آخرین زمان خروجی حل می شود و درون یابی نمی شود. به ویژه، حل کننده از آخرین بار گام بر نمی دارد. |
سوابق | عدد صحیح، 1-4 | 4 | نظم کلاسیک رانگ–کوتا. |
لایه | شامل | حذف کنند | عبارتند از | استراتژی تخمین خطا |
exprs | رشته | عبارت برای گام برداشتن زمانی که tstepping=elemexprs . | |
eventtol | اسکالر مثبت | 0.01 | تحمل رویداد برای یافتن ریشه شرایط رویداد هنگام استفاده از رویدادهای ضمنی برای BDF استفاده می شود. |
ewtrescale | در | خاموش | بر | به روز رسانی تحمل مطلق مقیاس شده برای BDF. |
هندسی غیرخطی | در | خاموش | خاموش | شامل غیرخطی بودن هندسی اگر مجوز شامل ماژول مکانیک سازه، ماژول آکوستیک، ماژول MEMS یا ماژول دینامیک چند بدنه باشد، در دسترس است. |
تأخیر ناپذیر | عدد صحیح مثبت | 15 | تعداد مراحل زمانی برای به تاخیر انداختن افزایش گام زمانی. |
غیر تاخیری | در | خاموش | خاموش | از تاخیر در افزایش گام زمانی استفاده کنید. |
startepbdf | اسکالر مثبت | 1e-3 | مرحله زمانی اولیه برای BDF. |
initialstepbdfactive | در | خاموش | خاموش | از یک مرحله زمانی اولیه برای BDF استفاده کنید. |
مرحله اولیه 5 | اسکالر مثبت | 1e-3 | مرحله زمانی اولیه برای Cash–Karp 5. |
inicialstepdopri5 | اسکالر مثبت | 1e-3 | مرحله زمانی اولیه برای Dormand–Prince 5. |
initialsteprk34 | اسکالر مثبت | 1e-3 | مرحله زمانی اولیه برای روش RK34. |
initialstepck5active | در | خاموش | خاموش | از یک مرحله زمانی اولیه برای Cash–Karp 5 استفاده کنید. |
initialstepdopri5active | در | خاموش | خاموش | از یک مرحله زمانی اولیه برای Dormand-Prince 5 استفاده کنید. |
initialsteprk34active | در | خاموش | خاموش | از یک مرحله زمانی اولیه برای RK34 استفاده کنید. |
مرحله اولیه | اسکالر مثبت | 1e-3 | مرحله زمانی اولیه برای آلفای تعمیم یافته. |
مرحله اولیه فااکتیو | در | خاموش | خاموش | از یک مرحله زمانی اولیه برای آلفای تعمیم یافته استفاده کنید. |
firststepfractionbdf- i | اسکالر مثبت | کسری از گام زمانی برای مرحله اولیه یک پله زمانی دستی برای BDF. نام و پیش فرض به ترتیب BDF بستگی دارد. به عنوان مثال، initialstepfractionbdf-2 برای BDF order 2. | |
startstepgrowthratebdf- i | اسکالر مثبت | نرخ رشد برای مراحل اولیه یک پله زمانی دستی برای BDF. نام و پیش فرض به ترتیب BDF بستگی دارد. برای مثال، initialstepgrowthratebdf-2 برای BDF order 2. | |
نگهدارنده | در | خاموش | خاموش | هشدارها را در گزارش ذخیره شده نگه دارید. |
lumpedflux | در | خاموش | خاموش | هنگام محاسبه شارها از برآمدگی استفاده کنید. |
توده ای | بله | شاید | شاید | ماتریس جرم منفرد. |
حداکثر سفارش | عدد صحیح بین 1 و 5 | 5 | حداکثر سفارش BDF. |
maxstepbdf | اسکالر مثبت | 1e-1 | حداکثر گام زمانی برای BDF، زمانی که maxstepconstraintbdf روی const تنظیم شده است . |
maxstepconstraintbdf | خودکار | const | expr | خودکار | حداکثر مرحله زمانی برای BDF: خودکار ( خودکار )، ثابت ( const )، یا یک عبارت ( expr ). |
maxstepconstraintdopri5 | خودکار | const | expr | خودکار | حداکثر مرحله زمانی برای Dormand-Prince 5: خودکار ( خودکار )، ثابت ( const ) یا یک عبارت ( expr ). |
maxstepconstraintgenalpha | خودکار | const | expr | خودکار | حداکثر گام زمانی برای آلفای تعمیم یافته: خودکار ( خودکار )، ثابت ( const ) یا یک عبارت ( expr ). |
maxstepdopri5 | اسکالر مثبت | 1e-1 | حداکثر گام زمانی برای Dormand-Prince 5، زمانی که maxstepconstraintdopri5 ثابت است . |
maxstepexpressionbdf | رشته | عبارت برای حداکثر گام زمانی برای BDF، زمانی که maxstepconstraintbdf expr است . | |
maxstepexpressiongendopri5 | رشته | عبارت برای حداکثر گام زمانی برای Dormand-Prince 5، زمانی که maxstepconstraintdopri5 expr است . | |
maxstepexpressiongenalpha | رشته | بیان حداکثر گام زمانی برای آلفای تعمیم یافته، زمانی که maxstepconstraintgenalpha expr است . | |
maxstepgenalpha | اسکالر مثبت | 1e-1 | حداکثر گام زمانی برای آلفای تعمیمیافته، زمانی که maxstepconstraintgenalpha ثابت است . |
پیام | رشته | پیام ورود به سیستم از آخرین فرآیند راه حل. | |
جزئی | 1 | 2 | 1 | حداقل سفارش BDF |
حل کننده mlin | مستقیم | توده ای | مستقیم | حل کننده ماتریس جرم: مستقیم یا توده ای (اگر timemethodexp rk نباشد ). |
nlsolver | اتوماتیک | کتابچه راهنمای | کتابچه راهنمای | تنظیمات حل کننده غیرخطی |
ntimestepsupdate | عدد صحیح مثبت | 100 | تعداد مراحل زمانی بین بهروزرسانیها، زمانی که updtlvl روی دستی تنظیم میشود ، وقتی timemethodexp روی ab3loc تنظیم نشده است . |
ntimestepsupdateab3loc | عدد صحیح مثبت | 100 | تعداد گامهای زمانی بین بهروزرسانیها، زمانی که updtlvl روی دستی تنظیم میشود وقتی timemethodexp روی ab3loc تنظیم شده است . |
odesolvertype | ضمنی | به صراحت | ضمنی | نوع حل کننده ODE: با استفاده از روش گام به گام زمانی ضمنی یا صریح. |
طرح | در | خاموش | خاموش | طرح در حین حل. |
plotfreq | tsteps | همه | همه | زمان به روز رسانی طرح. |
گروه طرح | رشته | نام گروه طرح برای طرح هنگام حل. | |
فرکانس نمونه | tsteps | همه | tsteps | زمان به روز رسانی پروب. |
کاوشگرها | بردار رشته ها | Probed برای استفاده از if probesel=manual . | |
کاوشگر | همه | هیچکدام | کتابچه راهنمای | همه | کاوشگرها برای محاسبه |
پیشگو | خطی | ثابت | خطی | نوع پیشگوی مورد استفاده (خطی یا ثابت)، |
واکنش نشان دادن | در | خاموش | بر | محاسبه نیروهای واکنش |
rescaleafterinitbw | در | خاموش | خاموش | پس از مقداردهی اولیه برای روشهای BDF و آلفای تعمیمیافته، زمانی که سازگاری روی bweuler تنظیم شده است، مقیاس مجدد را تغییر دهید . |
گرد | عددی | 0.75 | ضریب تقویت برای فرکانس های بالا. |
روش rk | rk34 | ck5 | dopri5 | rk34 | روش Runge–Kutta: RK34، Cash–Karp 5، یا Dormand–Prince 5. فقط زمانی در دسترس است که timemethodexp روی rk تنظیم شده باشد . |
rkstiffcheck | در | خاموش | بر | اگر مشکل از نظر عددی برای حل کننده های Runge-Kutta سخت شد، بررسی کنید و متوقف کنید. |
rktimestep | رشته | 1e-3 | گام زمانی برای گام برداشتن زمان دستی با روش کلاسیک Runge-Kutta. |
rtol | عددی | 0.01 | تحمل نسبی |
rtstepab3loc | عددی | 0.1 | تغییر مرحله زمانی نسبی، چه زمانی |
stabcntrl | در | خاموش | خاموش | از یک کنترل کننده غیرخطی برای کنترل گام-زمان کارآمدتر در روش BDF استفاده کنید. |
storeudot | در | خاموش | بر | مشتقات زمانی را ذخیره کنید. |
فاکتور tderglobal | عددی | 1 | فاکتور مشتق زمان جهانی، اگر روش tderglobal روی دستی تنظیم شود . |
روش tderglobal | خودکار | کتابچه راهنمای | خودکار | روش تحمل مشتق جهانی. |
tderfactor | عددی | 1 | فاکتور مشتق زمانی، اگر روش tder روی دستی تنظیم شده باشد . |
روش tder | خودکار | کتابچه راهنمای | خودکار | روش تحمل مشتق. |
روش زمانی | bdf | جنالفا | init | بی دی اف | روش گام زمانی ضمنی |
timemethodexp | rk | erk | ab3 | ab3loc | rk | حل کننده صریح زمان: رانج–کوتا، رانگ–کوتا کلاسیک، آدامز–باشفورث 3 یا آدامز–باشفورث 3 (محلی). |
timestepbdf | اسکالر عددی | بردار عددی | رشته با بیان | 0.01 | گام زمانی هنگام گام برداشتن دستی با استفاده از روش BDF. |
timestepck5 | اسکالر عددی | بردار عددی | رشته با بیان | 0.01 | گام زمانی هنگام گام برداشتن زمان دستی با استفاده از روش Cash–Karp 5 Runge–Kutta. |
timestepdopri5 | اسکالر عددی | بردار عددی | رشته با بیان | 0.01 | گام زمانی زمانی که گامهای زمانی دستی با استفاده از روش خوابیده-پرینس رانج-کوتا استفاده میشود. |
timestepgenalpha | اسکالر عددی | بردار عددی | رشته با بیان | 0.01 | مرحله زمانی زمانی که گامهای زمانی دستی با استفاده از روش آلفای تعمیمیافته انجام میشود. |
timesteprk34 | اسکالر عددی | بردار عددی | رشته با بیان | 0.01 | گام زمانی هنگام گام برداشتن زمان دستی با استفاده از روش RK34 Runge–Kutta. |
tlist | بردار عددی | لیست زمان های خروجی | |
همه | tlist | tsteps | tstepsclosest | tlist | زمان برای ذخیره: زمان خروجی توسط درون یابی. هر مرحله N که توسط حل کننده برداشته می شود. یا مراحل انجام شده توسط حل کننده نزدیک به زمان خروجی. |
قدم زدن | دفترچه راهنما | elemexprs | کتابچه راهنمای | دستی یا از عبارات زمان پله. |
tstepsbdf | رایگان | متوسط | سختگیرانه | کتابچه راهنمای | رایگان | حالت گام به گام زمانی که timemethod روی bdf تنظیم شده است . |
tstepsck5 | رایگان | متوسط | سختگیرانه | کتابچه راهنمای | رایگان | حالت پله زمانی زمانی که rkmethod روی ck5 تنظیم شده است . |
tstepsdopri5 | رایگان | متوسط | سختگیرانه | کتابچه راهنمای | رایگان | حالت پله زمانی زمانی که rkmethod روی dopri5 تنظیم شده است . |
tstepsgenalpha | رایگان | متوسط | سختگیرانه | کتابچه راهنمای | رایگان | حالت پله زمانی زمانی که timemethod روی genalpha تنظیم شده است . |
tstepsrk34 | رایگان | متوسط | سختگیرانه | کتابچه راهنمای | رایگان | حالت پله زمانی زمانی که rkmethod روی rk34 تنظیم شده است . |
tstepsstore | عدد صحیح مثبت | 1 | مقدار N برای هر N گام از حل کننده تا ذخیره زمانی که tout روی tsteps تنظیم شده است . |
updtlvl | نادرست | دفترچه راهنما | عامل | نادرست | مرحله زمانی بهروزرسانی: خاموش ( نادرست )، دستی ، یا فاکتور زمانی که timemethodexp روی ab3loc تنظیم شده است . |
updtstep | نادرست | کتابچه راهنمای | نادرست | مرحله زمانی بهروزرسانی: خاموش ( نادرست ) یا دستی ، زمانی که tstepping برای روشهای timemethodexp erk و ab3 به طور واضح تنظیم شده است . |
به طور پیش فرض می توانید فرآیند حل معادلات خطی یا غیرخطی را در هر مرحله زمانی به صورت دستی کنترل کنید. برای یک مشکل جفت شده، این کار از طریق ویژگی های Damp ، Dtech ، Hnlin ، Initstep ، Jtech ، Maxiter ، Minstep و Rstep که در زیر femnlin فهرست شده اند انجام می شود . برای یک مشکل تفکیک شده، ویژگی های فهرست شده در زیر femstatic که مربوط به حل کننده تفکیک شده هستند، در دسترس هستند. هنگامی که Timemethod روی bdf تنظیم می شود ، می توان از حل کننده غیرخطی داخلی انتگرالگر زمان استفاده کرد. این را می توان با تنظیم به دست آوردNlsolver به خودکار .
خواص atol ، atolmethod ، atolglobal ، atolglobalfactor ، atolglobalmethod ، atolglobalvaluemethod ، atoludot و atoludotactive نیاز به توضیح بیشتری دارد. مقدار پیشفرض تلورانس مطلق برای همه فیلدها با ویژگی atolglobalfactor یا atolglobal بسته به تنظیم برای atolglobalvaluemethod داده میشود . اصلاح کننده atolglobalmethod مشخص می کند که آیا مقدار داده شده atolglobal یا خیرباید برای متغیرهای مقیاس شده یا بدون مقیاس اعمال شود. برای متغیرهایی که مقیاس خودکار کار خوبی انجام می دهد، یا جایی که از مقیاس دستی استفاده شده است، تعیین تلورانس مطلق در متغیرهای مقیاس شده بسیار آسان تر است. اگر برای یک یا چند متغیر قدر مطلق یا روش مقیاسبندی متفاوتی نسبت به روش atolglobal و atolglobal میخواهید، میتوانید از ویژگیهای atol و atolmethod استفاده کنید . atol را بهعنوان رشتهای با فاصله با نامهای فیلد و تلورانسهای متناوب وارد کنید (مثلاً «u 1e-3 v 1e-6») . روش atolmethod را به عنوان رشته ای با فاصله با نام فیلدهای متناوب و یکی از رشته های سراسری وارد کنید .مقیاسپذیر یا بدون مقیاس (مثلاً « u unscaled v scaled »). به طور پیش فرض اتل متد برابر با global برای همه فیلدها است. لیست های atol و atolmethod لازم نیست شامل همه فیلدها باشد. مواردی که حضور ندارند تلورانس مطلقی را دریافت می کنند که توسط روش atolglobal و atolglobal مشخص شده است . هنگام حل معادلات نوع موج، مشتقات زمانی همه میدانها نیز به عنوان مجهول در نظر گرفته میشوند و بنابراین باید تلورانسهای مطلق نیز برای این مؤلفهها مشخص شود. به طور پیش فرض این تلورانس ها به طور خودکار انتخاب می شوند. در برخی شرایط ممکن است لازم باشد که آنها را به صورت دستی با ویژگی های atoludot مشخص کنیدو آتولودوتاکتیو . برای روشن کردن مشخصات دستی، مثلاً، دو فیلد u و v ، ویژگی atoludotactive را روی رشته «u on v on» قرار دهید . اگر atoludot مشخص نشده باشد، این دو مشتق زمانی تلورانس مطلق پیش فرض 1e-3 را دریافت می کنند . برای تعیین سایر تلورانس های مطلق، atoludot را به عنوان مثال روی رشته “u 1e-4 v 1e-7” تنظیم کنید . روش تحمل مطلق برای تمام مشتقات زمانی مانند روشی است که برای خود فیلد مشخص شده است.
حداکثر خطای نسبی مجاز در هر مرحله زمانی (خطای محلی) با استفاده از rtol مشخص می شود . با این حال، برای مولفههای کوچک بردار راهحل U ، الگوریتم فقط سعی میکند تا خطای محلی مطلق در U را کمتر از تحمل مطلق داده شده کاهش دهد.
هیچ تضمینی وجود ندارد که تحمل خطا به شدت رعایت شود. یعنی برای مشکلات سخت می توان از آنها فراتر رفت.
برای پارامتر تلورانس در معیار همگرایی برای سیستم های خطی، از حداکثر اعداد rtol و itol استفاده می شود.
اگر اعداد مختلط در فرآیند حل رخ می دهد از complex=on استفاده کنید .
ویژگی Consistent مقدار دهی اولیه یک سیستم معادله جبری دیفرانسیل (DAE) را کنترل می کند. مقدار Consistent=off به این معنی است که مقادیر اولیه سازگار هستند (این به ندرت اتفاق می افتد زیرا مقدار اولیه مشتق زمانی 0 است). در غیر این صورت، حل کننده سعی می کند مقادیر اولیه را به گونه ای تغییر دهد که سازگار شوند. مقدار consent=on را می توان (وقتی timemethod=bdf و nlsolver=automatic ) برای DAE های index-1 استفاده کرد. سپس حل کننده مقادیر DOF های دیفرانسیل را ثابت می کند و مقادیر اولیه DOF های جبری و مشتق زمانی DOF های دیفرانسیل را حل می کند. مقدار Consistent=bweulerمی تواند برای هر دو DAE index-1 و index-2 استفاده شود. سپس حل کننده با برداشتن یک گام اویلر به عقب، مقادیر اولیه همه DOFها را مختل می کند.
برای یک سیستم DAE، اگر Estrat=exclude ، DOFهای جبری از هنجار خطای خطای گسسته سازی زمانی حذف می شوند.
زمانی که timemethod روی bdf مشخصه initialstepdopri5 وقتی timemethod روی dopri5 تنظیم میشود و خاصیت initialstepgenalpha وقتی timemethod روی genalpha تنظیم میشود ، میتوانید اندازه مرحله زمانی اولیه را پیشنهاد کنید . همچنین برای فعال بودن مرحله اولیه مشخص شده، باید یکی از ویژگی های initialstepbdfactive ، initialstepdopri5active ، یا initialstepgenalphaactive را روی on قرار دهید.
به طور پیشفرض، حلکننده با بررسی ردیفها یا ستونهای صفر در ماتریس جرم، تعیین میکند که آیا سیستم دیفرانسیل-جبری است. اگر یک DAE دارید که در آن ماتریس جرم هیچ سطر یا ستون صفر ندارد، masssingular=yes را قرار دهید .
maxorder ویژگی حداکثر درجه چند جمله ای درون یابی را در روش BDF نشان می دهد (وقتی timemethod=bdf ).
اگر timemethod = bdf و maxstepconstraintbdf = const ، آنگاه ویژگی maxstepbdf یک حد بالایی برای اندازه گام زمانی قرار می دهد (این ویژگی زمانی که tstepsbdf = دستی مجاز نیست ). اگر در عوض maxstepconstraintbdf = expr ، ویژگی maxstepexpressionbdf حداکثر اندازه گام را از طریق عبارتی که هنگام حل ارزیابی میشود، کنترل میکند. همین امر برای ویژگی های maxstep مرتبط صادق است اگر timemethod = genalpha یا timemethod = rk و rkmethod= dopri5 .
از ویژگی odesolver برای انتخاب روش مرحله زمانی برای ODE ها استفاده می شود:
• | در حالت ضمنی (پیشفرض)، یک روش گامبهگام زمانی ضمنی مانند BDF استفاده میشود. |
• | در روش صریح ، از روش گامهای زمانی صریح، مانند روشهای صریح خانواده Runge-Kutta استفاده میشود. |
ویژگی timemethod برای انتخاب روش گام زمانی ضمنی استفاده می شود:
• | با timemethod=bdf ، حل کننده IDA (که از یک فرمول تمایز متغییر مرتبه برگشتی استفاده می کند) استفاده می شود. |
• | با timemethod=genalpha ، از روش α تعمیم یافته استفاده می شود. با α تعمیمیافته ، میرایی عددی را میتوان با دادن مقداری کنترل کرد، ![]() |
• | هنگامی که timemethod برای شروع تنظیم می شود ، حل کننده مقادیر اولیه ثابت (برای زمان شروع، همانطور که توسط ویژگی tlist تعریف شده است ) را برای سیستم محاسبه می کند و سپس متوقف می شود. مشتقات زمانی متغیرهای جبری و توابع نشانگر ممکن است هنوز پس از این عملیات از حالت اولیه خارج شوند. چنین کمیت های بدون مقدار اولیه با NaN (نه یک عدد) در شی محلول نشان داده می شوند. |
از ویژگی timemethodexp برای انتخاب روش مرحله زمانی صریح استفاده می شود:
• | با timemethodexp=rk ، از روش Runge–Kutta صریح استفاده میشود: RK34، Cash–Karp 5، یا Dormand–Prince 5 (مشخص شده با استفاده از ویژگی rkmethod ). |
• | با timemethodexp=erk ، یک روش کلاسیک Runge–Kutta از مرتبه 1-4 استفاده میشود. |
• | با timemethodexp=ab3 ، از روش صریح Adams–Bashforth 3 استفاده میشود. |
• | با timemethodexp=ab3loc ، یک نسخه محلی مرحلهای از Adams–Bashforth 3 استفاده میشود (برای رابط PDE شکل موج). |
ویژگی reacf محاسبه و ذخیره نیروی واکنش محدودیت را کنترل می کند. مقدار reacf=on (پیشفرض) به این معنی است که حلکننده بردار L باقیمانده FEM را در شی حل ذخیره میکند. از آنجا که L = N F Λ برای یک راه حل همگرا، باقیمانده همان نیروی محدودیت است. فقط اجزای L که مربوط به ردیف های غیر صفر N F هستند ذخیره می شوند. برای هر بار که راه حل درخواست می شود، یک مجموعه برداری باقیمانده اضافی انجام می شود. مقدار reacf=off هیچ محاسبه یا ذخیره سازی نیروی واکنش را ارائه نمی دهد و بنابراین می تواند مقداری در زمان محاسباتی صرفه جویی کند.
tlist ویژگی باید بردار کاملاً یکنواختی از اعداد واقعی باشد. معمولاً بردار از یک زمان شروع و یک زمان توقف تشکیل شده است. اگر بیش از دو عدد داده شود، زمانهای میانی را میتوان به عنوان زمانهای خروجی یا برای کنترل اندازه مراحل زمانی استفاده کرد (به زیر مراجعه کنید). اگر فقط یک عدد داده شود، نشان دهنده زمان توقف است و زمان شروع 0 است.
ویژگی tout زمان هایی را که در خروجی رخ می دهد تعیین می کند. اگر tout=tsteps ، خروجی شامل هر N امین گام زمانی است (که در آن N با استفاده از ویژگی tstepsstore مشخص می شود ؛ پیش فرض: 1) که توسط حل کننده گرفته شده است. اگر tout=tlist ، خروجی حاوی راه حل های درون یابی شده برای زمان ها در ویژگی tlist است . اگر tout=tstepscloses t. پیش فرض tout=tlist است .
ویژگیهای tstepsbdf (قابل اجرا زمانی که timemethod=bdf )، tstepsdopri5 (قابل اجرا زمانی که timemethod=dopri5 است )، و tstepsgenalpha (قابل اجرا زمانی که timemethod=genalpha است ) انتخاب مراحل زمانی را کنترل میکنند. اگر هر یک از این ویژگی ها آزاد باشد ، حل کننده مراحل زمانی را بر اساس منطق خود انتخاب می کند و زمان های میانی در بردار tlist را نادیده می گیرد . اگر هر یک از ویژگی ها روی سخت تنظیم شود ، مراحل زمانی که توسط حل کننده برداشته می شود حاوی زمان های موجود در tlist است . اگر هر یک از ویژگی ها روی متوسط تنظیم شود، پس حداقل یک مرحله زمانی در هر بازه از بردار tlist وجود دارد . اگر tstepsgenalpha روی دستی تنظیم شده باشد ، حل کننده مرحله زمانی مشخص شده در ویژگی timestepgenalpha را دنبال می کند . اگر timestepgenalpha یک مقدار اسکالر باشد، این مرحله زمانی در کل شبیه سازی گرفته می شود. هنگامی که timestepgenalpha یک بردار عددی (کاملاً یکنواخت) است، حلکننده جواب را در زمانهای بردار محاسبه میکند. زمان شروع و زمان توقف همچنان از tlist بدست می آید . بردار داده شده در timestepgenalphaبا استفاده از اولین و آخرین مرحله زمانی در بردار کوتاه یا بزرگ می شود تا زمان شروع و زمان توقف با مقادیر موجود در tlist مطابقت داشته باشد . در نهایت، عبارتی با استفاده از متغیرهایی با دامنه جهانی که منجر به اسکالر می شود، می تواند به عنوان timestepgenalpha استفاده شود .
برای مسائل از نوع موج، منطقی که توسط آن حل کننده مرحله زمانی را انتخاب می کند، گاهی اوقات می تواند منجر به یک مرحله زمانی شود که به شیوه ای ناکارآمد نوسان می کند. زمانی که timemethod=genalpha (حلکنندهای که معمولاً برای مسائل نوع موج استفاده میشود)، میتوانید با استفاده از ویژگیهای incrdelay و incrdelayactive از چنین نوسانهایی در مرحله زمانی اجتناب کنید . هنگامی که incrdelayactive=on ، شمارنده تعداد گام های زمانی متوالی را که برای آنها افزایش گام زمانی تضمین شده است، پیگیری می کند. هنگامی که این شمارنده از عدد داده شده در ویژگی incrdelay بیشتر شود ، گام زمانی افزایش یافته و شمارنده روی صفر تنظیم می شود.
ترتیب روش Runge–Kutta را می توان با خاصیت erkorder تنظیم کرد . اندازه گام زمانی از طریق ویژگی rktimestep کنترل می شود و می تواند به عنوان یک مقدار اسکالر منفرد، یک بردار عددی (کاملاً یکنواخت) یا یک عبارت با استفاده از متغیرهایی با دامنه جهانی ارائه شود که منجر به یک اسکالر می شود. برای Adams–Bashforth 3، فقط یک مقدار ثابت اسکالر از گام زمانی مجاز است. گام برداشتن از عبارات tstepping=elemexprs برای رابط PDE شکل موج مفید است. یک نسخه محلی مرحلهای از Adams–Bashforth 3 برای واسط Wave Form PDE توسط timemethodexp = ab3loc در دسترس است .
سازگاری
حل کننده TimeExplicit time-explicit از Model Builder در نسخه 6.1 حذف شده است اما همچنان می تواند در API استفاده شود. در عوض، از حلکننده زمان با یک پله زمانی صریح استفاده کنید. اگر مدلهای ایجاد شده در نسخههای قبل از 6.1 را که شامل حلکننده TimeExplicit هستند باز کنید ، حلکننده TimeExplicit همچنان در Model Builder موجود است .
![]() | در مدلهای مکانیک سازه، جابجاییها اغلب بسیار کوچک هستند، و بسیار مهم است که یک مقدار تحمل مطلق تعریفشده توسط کاربر کوچکتر از جابجاییهای واقعی انتخاب شود. |
![]() | برای اطلاعات بیشتر در مورد حلگر وابسته به زمان؛ حل وابسته به زمان را در کتابچه راهنمای مرجع Multiphysics COMSOL ببینید . |