企业官网建设流程全解析

朝阳区网站建设,网推拉新怎么做,计算机学院网站建设系统可行性分析,苏醒主题wordpress晶体塑性有限元显示动力学cpfem_vumat子注意此程序为\界面调用程序\#xff0c;用于显示动力学。最近在搞晶体塑性有限元仿真的时候#xff0c;发现显式动力学框架下的材料子程序开发真是个磨人的小妖精。特别是那个cpfem_vumat#xff0c;说是界面调用程…晶体塑性有限元显示动力学cpfem_vumat子注意此程序为\界面调用程序\用于显示动力学。最近在搞晶体塑性有限元仿真的时候发现显式动力学框架下的材料子程序开发真是个磨人的小妖精。特别是那个cpfem_vumat说是界面调用程序听起来像是个中间商实际上它决定了材料模型能不能和求解器愉快地玩耍。先看个刺激的——这是VUMAT入口的骨架结构SUBROUTINE VUMAT( ! 这里省略ABAQUS官方定义的几十个参数 ) ! 材料参数读取 nstatev ... ! 从props数组里抠出状态变量总数 dt ... ! 当前时间步长 ! 材料计算核心 DO km 1, nblock ! 把变形梯度F从数组里提出来 F_new(1:3,1:3) reshape( (...), (/3,3/) ) ! 调用晶体塑性本构 call crystal_plasticity( F_new, stress, statev ) END DO END SUBROUTINE这段代码最要命的是那个nblock循环——显式算法的并行特性要求同时处理多个材料点。我上次忘了把状态变量数组的索引和km关联结果所有积分点的状态都串了仿真出来像毕加索的画作。晶体塑性核心算法里有个魔鬼细节滑移系 Schmid 因子计算。看看这个让人头皮发麻的片段do is1,nSlip schmid 0.0 do i1,3 do j1,3 schmid schmid dir_0(i,is)*nor_0(j,is)*Fc(i,j) enddo enddo tau(is) schmid * stress enddo这里dir0和nor0存储着初始滑移方向/法向。有次我把张量双点乘写成单点乘结果剪应力计算全错但程序居然没报错后来发现应力-应变曲线像心电图才意识到是这里埋了雷。说到时间积分显式算法对步长敏感得像初恋少女。这个稳定性判断我栽过跟头dt_stable 0.8 * element_size / sqrt( youngs_modulus / density ) if (dt dt_stable) then call xplb_abort() ! 直接崩掉比给出错误结果强 endif0.8这个安全系数别随便改。有次为了省时间调到0.85模型在20%应变时突然像爆米花一样炸开。后来用特征值分析才发现某个积分点的雅可比矩阵已经鬼畜了。最后给新手们一个忠告多存几套状态变量备份。有次我在迭代计算时没做好中间状态保存导致在显式框架下的隐式迭代出现了时间旅行悖论——当前步的结果竟然依赖下一步的状态调试的时候感觉自己像在对付量子计算机。代码框里的缩进可能显示异常实际开发建议用等宽字体查看