if ((elt_m .ge. 1.0).and.(sigma_N.ge.0)) then !---------matrix tensile cracking---------
                        Gnc=Gnct
                       
                          ftype=1.0
                          T0_eff=T_eff/sqrt(elt_m)
                          D0_eff=D_eff/sqrt(elt_m)
                        Gc=Gnc+(Gsc-Gnc)*(Gs/(Gn+Gs))**eta
                          Df_eff=D0_eff
                        Df_eff=statev(6)
                        E_mn=0.5*D0_eff*T0_eff
                        do while(E_mn.lt.Gc)
                            Df_eff=Df_eff+0.001
                                statev(6)=Df_eff
                                E_mn=E_mn+0.001*exp(T0_eff*D0_eff/(Gc-0.5*T0_eff*D0_eff)*(1.0-Df_eff/D0_eff))*T0_eff
                        end do
                        Df_eff=statev(6)+0.001
                        if(Df_eff.le.D0_eff) then
                                dm=dgdf2
                                ftype=1.5
                        endif

这是包含循环部分的程序，请帮我看看do while部分有没有错误，谢谢了！ 返回小木虫查看更多