创新的邯郸网站建设个人备案购物网站

创新的邯郸网站建设,个人备案购物网站,网站制作 苏州,wordpress主题ajax修正的考虑剪切和静水压力的各向异性和各向同性 GTN 模型 Abaqus_Vumat 子程序以及跑通的案例#xff0c;包括 三个子程序 for 文件#xff08;各向同性#xff0c;自己推导修改以及各向异性#xff09; 、cae 文件和 inp 文件等#xff0c;还有一些论文资料。 注#x…修正的考虑剪切和静水压力的各向异性和各向同性 GTN 模型 Abaqus_Vumat 子程序以及跑通的案例包括 三个子程序 for 文件各向同性自己推导修改以及各向异性 、cae 文件和 inp 文件等还有一些论文资料。 注1. 各项异性模型为 hill48三维应力状态 2. abaqus 版本为 2020更低版本可以导入 inp 文件查看 及求解 3. 子程序是在一个各向同性子程序基础上改的推导了相关公式修改了感觉不太对的地方 子程序都会提供可以自己推导。 4.基质屈服应力考虑了温度和应变率 另外程序中提供了两种常用的本构。搞金属成型仿真的老司机应该都踩过GTN模型的坑——这玩意儿在剪切和静水压力下的预测总是飘忽不定。最近为了项目需求硬刚了半个月的VUMAT子程序总算把各向同性和各向异性的修正版GTN模型都跑通了。这里分享些实战心得附上三个版本的子程序文件各向同性原版、各向同性修正版、Hill48各向异性版和验证案例。一、屈服函数的暴力改造原版各向同性模型最坑的地方在于q2参数的处理。在静水压力主导的工况下屈服函数中的指数项会直接崩掉。我们团队推导时加入了修正系数η! 原版屈服函数 Phi (q/sigmabar)**2 2.0*q1*fstar*cosh(1.5*q2*p/sigmabar) - (1.0 q3*fstar**2) ! 修正后加入压力修正项 eta 1.0/(1.0 alpha*(p/sigmabar)**beta) Phi (q/sigmabar)**2 2.0*q1*fstar*cosh(1.5*q2*eta*p/sigmabar) - (1.0 q3*fstar**2)这个eta函数在高压状态下能有效抑制cosh项爆炸参数α和β需要通过单轴压缩试验标定。调试时发现当β取0.5时304不锈钢的预测误差从23%降到了7%以内。二、Hill48各向异性的嫁接手术各向异性改造的核心是把等效应力替换成Hill48准则。这里有个巨坑Hill的应力分量顺序和Abaqus内置顺序不一致! Hill48等效应力计算注意应力分量顺序 F 0.5*(HillF*(S(2)-S(3))**2 HillG*(S(3)-S(1))**2 HillH*(S(1)-S(2))**2) F F HillL*S(4)**2 HillM*S(5)**2 HillN*S(6)**2 seq sqrt(F)这里S(1)~S(6)对应Abaqus的应力分量顺序是σ11,σ22,σ33,σ12,σ13,σ23。当初没注意这个直接套用论文公式导致预测结果完全反常识浪费了两天时间查错...三、温度应变率的双杀处理基质材料的屈服应力采用修正Johnson-Cook模型支持两种模式切换! 模式1: 标准J-C模型 IF(MODEL_TYPE 1) THEN sigmay (A B*epbar**n) * (1.0 C*log(epdot/epdot0)) * (1.0 - (T-Troom)/(Tmelt-Troom)**m) ! 模式2: 考虑动态回复的修正版 ELSE sigmay A*(1.0 epbar/ep0)**B * exp(-C*(T-Troom)) * (1.0 (epdot/epdot0)**D) ENDIF在高温高速工况下比如热冲压模式2比传统模型精度提升约15%。但要注意D参数不能为负否则会在低应变率时出现应力骤降的异常现象。四、验证案例翻车实录提供的案例包中包含三个典型工况单轴拉伸验证损伤演化平面应变压缩测试压力修正项扭转剪切各向异性响应跑案例时遇到最诡异的问题各向异性模型在扭转工况下损伤值比各向同性模型还高后来发现是Hill参数标定错误F参数多输了个小数点...修正后结果终于符合理论预期。资源清单/VUMATISOOriginal.f (原始各向同性版本)/VUMATISOModified.f (压力修正版)/VUMATAnisoHill48.f (Hill48各向异性版)/Case_Uniaxial.inp (单轴拉伸案例)/MaterialData.pdf (参数标定指南)建议用Abaqus2020直接提交Job低版本用户可能会遇到材料定义不兼容问题。遇到收敛困难时尝试把最大增量步从0.1调到0.05亲测有效