根据tsaiwu屈服准则编写的，屈服函数大于0后，应力基本不变

星尘17

根据tsaiwu屈服准则编写的，用的是一个单元做的测试，屈服函数大于0后，应力基本不改变
每次判断屈服后进入循环，算得的塑性应变增加基本等于总应变增量导致的应力没变化
可是公式应该都没错，不知道问题出在哪，求各位指点一下，方便的话可以QQ或者邮箱指导嘛，万分感谢

Alex1478

楼主用的哪个版本的lsdyna

星尘17

Alex1478 发表于 2021-1-7 05:38
楼主用的哪个版本的lsdyna

用的5.1.1   其他的license搞不定

xiaoduan

楼主的程序我有两点不明白：
1）在没有给出历史变量hsv（1）-hsv（6）的等式下直接用是否不妥，初始下所有历史变量是0；
2）应力计算是基于塑性应变增加计算的，也就是应力增量等于应力偏量增量和净水压力增量之和得到的；
我在楼主的程序里没有看到这块的代码；

星尘17

xiaoduan 发表于 2021-1-20 11:43
楼主的程序我有两点不明白：
1）在没有给出历史变量hsv（1）-hsv（6）的等式下直接用是否不妥，初始下所有 ...

1）历史变量是需要先给定义嘛，初始时候历史变量为0，直接用是会产生什么问题的吗
2）我原先是用的下面的代码
      alt=cm(1)*cm(2)/(1.0+cm(2))/(1.0-2.0*cm(2))
      au=cm(1)/2.0/(1.0+cm(2))
      g2 =cm(1)/(1.0+cm(2))
      g  =0.5*g2
      davg=(-eps(1)-eps(2)-eps(3))/3.0
      p=-davg*cm(1)/((1.0-2.0*cm(2)))
      sig(1)=sig(1)+p+g2*(eps(1)+davg)
      sig(2)=sig(2)+p+g2*(eps(2)+davg)
      sig(3)=sig(3)+p+g2*(eps(3)+davg)
      sig(4)=sig(4)+g*eps(4)
      sig(5)=sig(5)+g*eps(5)
      sig(6)=sig(6)+g*eps(6)
      hsv(13)=-(sig(1)+sig(2)+sig(3))/3.0
      hsv(14)=1.0/6.0*((sig(1)-sig(2))**2+(sig(2)-sig(3))**2
     &        +(sig(3)-sig(1))**2+6.0*(sig(4)**2+sig(5)**2
     &        +sig(6)**2))
c

      f=hsv(14)-(cm(4)+cm(5)*hsv(13)+cm(6)*hsv(13)*hsv(13))
      k=0
111   if(f.gt.0.0) then
        if(k.le.100000)then
          fn11=(2.0/3.0-2.0/9.0*cm(6))*sig(1)-(1.0/3.0+2.0/9.0*cm(6))
     &          *(sig(2)+sig(3))+1.0/3.0*cm(5)
          fn22=(2.0/3.0-2.0/9.0*cm(6))*sig(2)-(1.0/3.0+2.0/9.0*cm(6))
     &          *(sig(1)+sig(3))+1.0/3.0*cm(5)
          fn33=(2.0/3.0-2.0/9.0*cm(6))*sig(3)-(1.0/3.0+2.0/9.0*cm(6))
     &          *(sig(1)+sig(2))+1.0/3.0*cm(5)
          fn12=2.0*sig(4)
          fn23=2.0*sig(5)
          fn13=2.0*sig(6)
          Fncfn=(alt+2.0*au)*(fn11**2+fn22**2+fn33**2)
     &          +au*(fn12**2+fn13**2+fn23**2)
     &          +2.0*alt*(fn11*fn22+fn22*fn33+fn11*fn33)
          adlt=f/Fncfn
          adsig11=adlt*((alt+2.0*au)*fn11+alt*fn22+alt*fn33)
          adsig22=adlt*((alt+2.0*au)*fn22+alt*fn11+alt*fn33)
          adsig33=adlt*((alt+2.0*au)*fn33+alt*fn22+alt*fn11)
          adsig12=adlt*au*fn12
          adsig23=adlt*au*fn23
          adsig13=adlt*au*fn13
          sig(1)=sig(1)-adsig11
          sig(2)=sig(2)-adsig22
          sig(3)=sig(3)-adsig33
          sig(4)=sig(4)-adsig12
          sig(5)=sig(5)-adsig23
          sig(6)=sig(6)-adsig13
          hsv(1)=hsv(1)+adlt*fn11
          hsv(2)=hsv(2)+adlt*fn22
          hsv(3)=hsv(3)+adlt*fn33
          hsv(4)=hsv(4)+adlt*fn12
          hsv(5)=hsv(5)+adlt*fn23
          hsv(6)=hsv(6)+adlt*fn13
          hsv(13)=-(sig(1)+sig(2)+sig(3))/3.0
          hsv(14)=1.0/6.0*((sig(1)-sig(2))**2+(sig(2)-sig(3))**2
     &            +(sig(3)-sig(1))**2+6.0*(sig(4)**2+sig(5)**2
     &            +sig(6)**2))
          f=hsv(14)-(cm(4)+cm(5)*hsv(13)+cm(6)*hsv(13)*hsv(13))
          k=k+1
          goto 111
        endif
      endif
c
      epsp=sqrt(2.0/9.0*((hsv(1)-hsv(2))**2+(hsv(2)-hsv(3))**2
     &     +(hsv(2)-hsv(3))**2+3.0/2.0*(hsv(4)**2+hsv(5)**2
     &     +hsv(6)**2)))
      hsv(15)=cm(7)+(hsv(13)/10.0**8-0.5)**2
c
      if(hsv(13).le.cm(3)) goto 222
      if(epsp.ge.hsv(15)) failel=.true.
      goto 333
c
222   sig(1)=-cm(3)
      sig(2)=-cm(3)
      sig(3)=-cm(3)
      sig(4)=0.0
      sig(5)=0.0
      sig(6)=0.0
      hsv(13)=cm(3)
c      
333   return
      end
      算起来还是有问题

xiaoduan

星尘17 发表于 2021-1-20 15:47
1）历史变量是需要先给定义嘛，初始时候历史变量为0，直接用是会产生什么问题的吗
2）我原先是用的下面的 ...

1）注意不发散就行，在某些情况下；
2）tsai-wu是复合材料的失效准则，前面求解sig分量怎么没体现正交各向异性的刚度矩阵呢，而采用各项同性不解？另外还采用J2屈服准则，而不是二次屈服面（可以退化为J2），没啥大问题。
3）111到222应该是tsai-wu准则部分，失效则单元应力置为净水压力，程序没有考虑属性加载吧。
楼主可否把111到222这段的公式贴出来，我可以看得明白些

星尘17

xiaoduan 发表于 2021-1-21 09:47
1）注意不发散就行，在某些情况下；
2）tsai-wu是复合材料的失效准则，前面求解sig分量怎么没体现正交各 ...

您方便用邮箱交流嘛  我的邮箱是19s030019@stu.hit.edu.cn
我参考材料的是两三篇文献 这里传文件挺慢的

阿寸的马甲

学习了。

CCT1996

请问楼主能发一下出自的参考资料吗？谢谢。邮箱2587399029@qq.com