如何利用creep来编写分段函数,可以自动对比选择应力区间来,进而时时更新所需方程。 我想利用子程序来控制蠕变(应力松弛)过程,方程如下:
应力松弛曲线如下:
在编写过程中没有过多的考虑等效应力的具体计算问题,均用Misesstress。在文件中两个子程序采用同一个beam.inp文件,其中devap没有改变,这里不会有影响。我尝试了两种方式,第一是直接用creep来编写,利用if语句进行选择控制。见子程序“onlycreep.for”文件。第二是利用USDFLD和creep组合来编写程序,见“usdfldandcreep.for”子程序。这里我是想是不是可以利用USDFLD来时实提取应力然后用于creep中的判断。(这里产考了的帖子) 然而在编写过两种情况进行计算时发现,程序确实是可以选择在那个IF中计算,但是这只是根据最初的应力值判断,然后全部过程均用该等式,而不会向我预想的那样会随着应力的不同而不断的选择用哪个方程。下面附上单独用一种方程计算的对比结果compare.opj,如果你有修改调试了一下,应该与这三个曲线均不同,应该是其变相组合。 最后再次强调一下我的问题,如何利用creep来编写,可以自动对比选择应力区间来,进而时时更新所需方程。 读到这里,我诚恳的问一下阁下,有没有这方面的经验给指点一二。在此先谢过了!
|