版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

yiran909

至尊木虫 (知名作家)

孤独的跑者

应助: 5 (幼儿园)
金币: 19418.5
散金: 4687
红花: 105
沙发: 65
帖子: 9935
在线: 532.5小时
虫号: 608646
注册: 2008-09-21
专业: 零件成形制造

[交流] 请问哪位能够提供有关hill 48 各向异性屈服准则的详细资料已有3人参与

大家好，我最近在看abaqus documentation 中的各向异性屈服准则的内容，
里边写的很简单就几个公式，也看不懂，
而且查了一些文献，但是那个屈服准则和函数的公式表达也不完全一样，
感觉看不懂，想推导推导，请求哪位能提供详细的关于Hill 48的屈服准则和函数的资料，
不胜感激。

回复此楼

» 猜你喜欢

三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

» 本主题相关价值贴推荐，对您同样有帮助:

ansys怎么模拟型材到达弯曲极限后破裂呢？各向异性硬化模型怎么定义破坏准则啊已经有7人回复

1楼2013-06-27 15:55:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

linrn

至尊木虫 (著名写手)

应助: 26 (小学生)
金币: 15056.6
散金: 5
红花: 2
帖子: 2171
在线: 283.3小时
虫号: 492936
注册: 2008-01-09
性别: GG
专业: 固体力学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

Google学术搜索里面一堆文献。
自己动手，丰衣足食。

[ 发自手机版 http://muchong.com/3g ]

赞一下(1人)

回复此楼

2楼2013-06-28 12:55:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yiran909

至尊木虫 (知名作家)

孤独的跑者

应助: 5 (幼儿园)
金币: 19418.5
散金: 4687
红花: 105
沙发: 65
帖子: 9935
在线: 532.5小时
虫号: 608646
注册: 2008-09-21
专业: 零件成形制造

引用回帖:

2楼: Originally posted by linrn at 2013-06-28 12:55:47
Google学术搜索里面一堆文献。
自己动手，丰衣足食。

哥，关键是我看不懂英文啊，
想找点中文的，像教材集中讲的
不是文献里边哪种，简单提起一下
活着是力学推导，都看不懂啊

赞一下

回复此楼

3楼2013-06-28 17:06:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

linrn

至尊木虫 (著名写手)

应助: 26 (小学生)
金币: 15056.6
散金: 5
红花: 2
帖子: 2171
在线: 283.3小时
虫号: 492936
注册: 2008-01-09
性别: GG
专业: 固体力学

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

看不懂英文，最好不学什么屈服方程。

[ 发自手机版 http://muchong.com/3g ]

赞一下(1人)

回复此楼

4楼2013-06-28 22:42:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

吕泽华

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 651.2
散金: 51
红花: 1
帖子: 225
在线: 54.1小时
虫号: 9181567
注册: 2018-06-24
专业: 金属基复合材料

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

***********************************************************************************
**  UMAT for Abaqus/Standard, large deformation formulation for plane stress    **
**  elements S4R. Implicit integration with continuum Jacobian.                **
***********************************************************************************
C==================================================================================
C                                                                               ==
C *  Copyright:  2007 by Zang Shun lai.  All Rights Reserved.                ==
C *             Date Aug, 24, 2007                                           ==
C *  Address: Xi'an Jiaotong University, Xianning Road, Xi'an Shaanxi,    ==
C *             P.R.China                                                    ==
C *                                                                         ==
C *  E-mail:    shunlai.zang@gmail.com                                     ==
C *                                                                         ==
C==================================================================================
C==================================================================================
C    NOTE:  UMAT FOR ABAQUS
C==================================================================================
C    (1) The statement: include 'aba_param.inc' must be added as below.
C
C    (2) This is current as of version 6.5 of abaqus.
C
C    (3) The current yield function is Hill1948 with three material parameters.
C    The anisotropy of initial yield stress and Lankford values can be described
C    by Hill1948 respectively.
C
C    (4) The current hardening model is isotropic hardening, which can %NOT% model
C    the Bauschinger effect, transient behavior in cyclic loading.
C
C    (5) The default type of variables,whose prefix is i-n��is integer. The name of
C    variables must be cared.
C
C    (6) The engineering shear strain is used in UMAT subroutine interface.
C
C    (7) For different unit system, the TOLER should be adjusted.
C
C==================================================================================
C
*USER SUBROUTINE
   SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,
   1  DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,
   2  CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
   3  CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)

      INCLUDE 'ABA_PARAM.INC'
      CHARACTER*8 CMNAME

      DIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),
   1  DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),TIME(2),
   2  PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),
   3  DFGRD0(3,3),DFGRD1(3,3)
C
C LOCAL ARRAYS
C ===================================================================
C EELAS  - ELASTIC STRAINS
C EPLAS  - PLASTIC STRAINS
C ===================================================================
C
   PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,FOUR=4.D0,
1       FIVE=5.D0,SIX=6.D0,SEVEN=7.D0,EIGHT=8.D0,
   2       ENUMAX=.4999D0,NEWTON=60,TOLER=1.0D-9,
   3       E=2.718281828459045235D0,THRESHOLD=50.D0,KKSTEP=5)
   DIMENSION EELAS(NTENS+1),EPLAS(NTENS+1),
1       RPOTENTIAL(KKSTEP),
2       DEELAS(NTENS),DEPLAS(NTENS),
3       DSIGMA(NTENS),RSIGMA(NTENS),
   4       STRESSOLD(NTENS),EELASOLD(NTENS+1),EPLASOLD(NTENS+1),
5       STRESS_B(NTENS),STRESSNEW(NTENS),
6       PFPSIGMA(NTENS), PAPSIGMA(NTENS,NTENS),
   7       O(NTENS,NTENS),OREVERSE(NTENS,NTENS),ATOREVERSE(NTENS),
8       CEA(NTENS),T1(NTENS),
   9       R(NTENS,NTENS),RA(NTENS),ATR(NTENS),RAATR(NTENS,NTENS),
   1       XIDEN(NTENS,NTENS),STR(NTENS,NTENS),IS(NTENS),JS(NTENS)
C
EMOD = PROPS(1)
   ENU    = MIN(PROPS(2),ENUMAX)
   EG    = EMOD/(TWO*(ONE+ENU))
   ELAM = (EMOD*ENU)/((ONE+ENU)*(ONE-TWO*ENU))
EG2    = TWO*EG
ELAMBAR = (EG2*ELAM)/(ELAM+EG2)
C
C ELASTIC STIFFNESS
C
   DO K1=1,NTENS
DO K2=1,NTENS
   DDSDDE(K1,K2) = ZERO
END DO
END DO

DO K1=1,NDI
      DO K2=1,NDI
      DDSDDE(K2,K1)=ELAMBAR
      END DO
      DDSDDE(K1,K1)=EG2+ELAMBAR
   END DO
   DO K1=NDI+1,NTENS
      DDSDDE(K1,K1)=EG
   END DO
C

C       CALCULATE NEW EFFECTIVE STRESS EFSTRESS
C
      CALL KCALES(STRESSNEW,YR0,YR45,YR90,EFSTRESS)

DELTA_LAMBDA = DELTA_LAMBDA + DDELTA_LAMBDA
DEQPLAS = DELTA_LAMBDA
EQPLAS=EQPLASOLD+DEQPLAS
C
C
C check yield criterion with updated constitutive values
C
SYIELD= SYIEL0 + EQ*(ONE-E**(-EB*EQPLAS))
RHS = EFSTRESS-SYIELD
IF(DABS(RHS) .LT. (TOLER*SYIEL0)) GOTO 11

   END DO !ENDDO FOR DO KEWTON=1,NEWTON
C
C       WRITE WARNING MESSAGE TO THE .MSG FILE
C
C       WRITE(7,2) KEWTON, RHS, EFSTRESS, EQPLAS

11 CONTINUE

   END IF !END FOR IF(ZY .GT. THRESHOLD)
C
C
C    UPDATE STRESS, ELASTIC AND PLASTIC STRAINS AND
C    EQUIVALENT PLASTIC STRAIN
C
   CALL CALDR(STRESSNEW,YR0,YR45,YR90,PFPSIGMA)
   DO K1=1,NTENS
      STRESS(K1)=STRESSNEW(K1)
      DEELAS(K1)=DSTRAN(K1)-DELTA_LAMBDA*PFPSIGMA(K1)
      DEPLAS(K1)=DELTA_LAMBDA*PFPSIGMA(K1)
      END DO
C
C    FORMULATE THE JACOBIAN (MATERIAL TANGENT)
C
C    CALCULATE PFPSIGMA(NTENS)
C
      CALL CALDR(STRESS,YR0,YR45,YR90,PFPSIGMA)
      CALL KCALES(STRESS,YR0,YR45,YR90,EFSTRESS)
C
C    CALCULATE A_iso
C
   H = EB*EQ*(E**(-EB*EQPLAS))
      A_iso = H
C
C    CALCULATE THE PAPSIGMA \FRAC{\PARTIAL A}{\PARTIAL \SIGMA}
C
      CALL CAL2DR(STRESS,YR0,YR45,YR90,PAPSIGMA)
C
C    CALCULATE THE TEMP VARIABLE O
C
   DO K1=1,NTENS
   DO K2=1,NTENS
      O(K1,K2)=ZERO
   END DO
   END DO

   DO K1=1,NTENS
   DO K2=1,NTENS
      DO K3=1,NTENS
      O(K1,K2)=O(K1,K2)+DDSDDE(K1,K3)*PAPSIGMA(K3,K2)
      END DO
   END DO
   END DO

       DO K1=1,NTENS
   DO K2=1,NTENS
      O(K1,K2)=XIDEN(K1,K2)+DELTA_LAMBDA*O(K1,K2)
   END DO
   END DO
C
C    CALCULATE OREVERSE \Theta^{-1}
C
   DO K1=1,NTENS
   DO K2=1,NTENS
      OREVERSE(K1,K2) = O(K1,K2)
   END DO
   END DO

      CALL BRINV(OREVERSE,NTENS,L,IS,JS)

   IF (L .EQ. 0) THEN
C       WRITE(7,15) KEWTON
15    FORMAT(//,30X,'***WARNING - THE MATRIX \Theta12356 DID NOT ',
   1                'CONVERGE AFTER ',I3,' ITERATIONS')
   END IF
C
C    CALCULATE R(NTENS,NTENS)=OREVERSE * DDSDDE
C
      DO K1=1,NTENS
   DO K2=1,NTENS
      R(K1,K2) = 0.D0
   END DO
   END DO

   DO K1=1,NTENS
   DO K2=1,NTENS
      DO K3=1,NTENS
      R(K1,K2) = R(K1,K2) + OREVERSE(K1,K3)*DDSDDE(K3,K2)
      END DO
   END DO
   END DO
C
C    CALCULATE RA(NTENS)
C
      DO K1=1,NTENS
      RA(K1)=0.D0
   END DO

      DO K1=1,NTENS
   DO K2=1,NTENS
      RA(K1)=RA(K1)+R(K1,K2)*PFPSIGMA(K2)
   END DO
   END DO
C
C    CALCULAT ATR(NTENS)
C
      DO K1=1,NTENS
      ATR(K1)=0.D0
   END DO

   DO K1=1,NTENS
   DO K2=1,NTENS
      ATR(K1)=ATR(K1)+PFPSIGMA(K2)*R(K2,K1)
   END DO
   END DO
C
C    CALCULAT RAATR(NTENS,NTENS)
C
      DO K1=1,NTENS
   DO K2=1,NTENS
      RAATR(K1,K2)=RA(K1)*ATR(K2)
   END DO
   END DO
C
C    CALCULAT ATRA
C
      ATRA=0.D0
   DO K1=1,NTENS
      ATRA= ATRA + PFPSIGMA(K1)*RA(K1)
   END DO
C
C    COPY TO DDSDDE(NTENS,NTENS)
C
      DO K1=1,NTENS
   DO K2=1,NTENS
      DDSDDE(K1,K2)=ZERO
   END DO
   END DO

   DO K1=1,NTENS
      DO K2=1,NTENS
      DDSDDE(K1,K2)=R(K1,K2)-RAATR(K1,K2)/(ATRA+A_iso)
      END DO
   END DO
END IF !END IF FOR IF(ZY .GT. 0.D0)
*********************************************************************
C STORE ELASTIC AND (EQUIVALENT) PLASTIC STRAINS
C IN STATE VARIABLE ARRAY
*********************************************************************
C
C    SAVE 11 AND 22 COMPONENT
C
DO K1=1,NDI
      STATEV(K1)=EELASOLD(K1)+DEELAS(K1)
      STATEV(K1+(NTENS+1))=EPLASOLD(K1)+DEPLAS(K1)
   END DO
C
C    SAVE 33 COMPONENT
C
   STATEV(3)=EELASOLD(3)-(ENU*(DEELAS(1)+DEELAS(2)))/(ONE-ENU)
STATEV(7)=EPLASOLD(3)-DEPLAS(1)-DEPLAS(2)
C
C    SAVE 12 COMPONENT
C
STATEV(4)=EELASOLD(4)+DEELAS(3)
STATEV(8)=EPLASOLD(4)+DEPLAS(3)
C
C    SAVE EQUIVALENT PLASTIC STRAIN
C
   STATEV(9)=EQPLAS
   RETURN
   END

**
*********************************************************************
**          UTILITY SUBROUTINES                               *
*********************************************************************
** CALCULATE YOUNG'S MODULUS, THE CHANGE OF YOUNG'S MODULUS WITH
** EQUIVALENT PLASTIC STRAIN IS CONSIDERED
*********************************************************************
*USER SUBROUTINE
   SUBROUTINE XYMOD(EQPLAS,EMOD1,ESTATMOD,EEPLAS,EMOD)
C
   INCLUDE 'ABA_PARAM.INC'
C
   IF (EQPLAS .LE. 0.04) THEN
      EMOD=EMOD1-((EMOD1-ESTATMOD)*EQPLAS)/EEPLAS
ELSE
   EMOD=ESTATMOD
ENDIF
RETURN
END
**
*********************************************************************
** INVERSE MATRIX                                              *
*********************************************************************
*USER SUBROUTINE
   SUBROUTINE BRINV(A,N,L,IS,JS)
DIMENSION A(N,N),IS(N),JS(N)
DOUBLE PRECISION A,T,D
L=1
DO 100 K=1,N
   D=0.0
   DO 10 I=K,N
   DO 10 J=K,N
      IF (ABS(A(I,J)).GT.D) THEN
      D=ABS(A(I,J))
      IS(K)=I
      JS(K)=J
      END IF
10    CONTINUE
   IF (D+1.0.EQ.1.0) THEN
      L=0
C       WRITE(*,20)
      RETURN
   END IF
20    FORMAT(1X,'ERR**NOT INV')
   DO 30 J=1,N
      T=A(K,J)
      A(K,J)=A(IS(K),J)
      A(IS(K),J)=T
30    CONTINUE
   DO 40 I=1,N
      T=A(I,K)
      A(I,K)=A(I,JS(K))
      A(I,JS(K))=T
40    CONTINUE
   A(K,K)=1/A(K,K)
   DO 50 J=1,N
      IF (J.NE.K) THEN
      A(K,J)=A(K,J)*A(K,K)
      END IF
50    CONTINUE
   DO 70 I=1,N
      IF (I.NE.K) THEN
      DO 60 J=1,N
         IF (J.NE.K) THEN
            A(I,J)=A(I,J)-A(I,K)*A(K,J)
         END IF
60       CONTINUE
      END IF
70    CONTINUE
   DO 80 I=1,N
      IF (I.NE.K) THEN
      A(I,K)=-A(I,K)*A(K,K)
      END IF
80    CONTINUE
100 CONTINUE
DO 130 K=N,1,-1
   DO 110 J=1,N
      T=A(K,J)
      A(K,J)=A(JS(K),J)
      A(JS(K),J)=T
110    CONTINUE
   DO 120 I=1,N
      T=A(I,K)
      A(I,K)=A(I,IS(K))
      A(I,IS(K))=T
120    CONTINUE
130 CONTINUE
RETURN
END
**
*********************************************************************
** CALCULATE EFFECTIVE STRESS FOR HILL48                         *
*********************************************************************
*USER SUBROUTINE
   SUBROUTINE KCALES(SIGMA,XR0,XR45,XR90,EES)

INCLUDE 'ABA_PARAM.INC'
PARAMETER (M=3)
DIMENSION SIGMA(M),P(M,M),PSIGMA(M)

BETA12=XR0/(1.D0+XR0)
BETA66=(XR0+XR90)*(1.D0+2.D0*XR45)/(XR90*(1.D0+XR0))
BETA23=XR0/(XR90*(1.D0+XR0))
BETA22=BETA12+BETA23

   DO K1=1,M
DO K2=1,M
   P(K1,K2) = 0.D0
   PSIGMA(K1) = 0.D0
END DO
END DO

P(1,1)=2.D0*1.D0
P(1,2)=-2.D0*BETA12
P(2,1)=-2.D0*BETA12
P(2,2)=2.D0*BETA22
P(3,3)=2.D0*BETA66

DO K1=1,M
DO K2=1,M
      PSIGMA(K1)=PSIGMA(K1)+P(K1,K2)*SIGMA(K2)
   END DO
END DO

EES=0.D0
DO K1=1,M
   EES = EES + SIGMA(K1)*PSIGMA(K1)
END DO
EES = DSQRT(EES/2.D0)
   RETURN
   END
**
*********************************************************************
** CALCULATE FIRST ORDER DERIVATIVE FOR HILL48                   *
*********************************************************************
*USER SUBROUTINE
   SUBROUTINE CALDR(SIGMA,XR0,XR45,XR90,PFPS)
C
   INCLUDE 'ABA_PARAM.INC'
   PARAMETER(M=3)
   DIMENSION SIGMA(M),PFPS(M),P(M,M),PSIGMA(M)

BETA12=XR0/(1.D0+XR0)
BETA66=(XR0+XR90)*(1.D0+2.D0*XR45)/(XR90*(1.D0+XR0))
BETA23=XR0/(XR90*(1.D0+XR0))
BETA22=BETA12+BETA23

   DO K1=1,M
DO K2=1,M
   P(K1,K2) = 0.D0
   PSIGMA(K1) = 0.D0
END DO
END DO

P(1,1)=2.D0*1.D0
P(1,2)=-2.D0*BETA12
P(2,1)=-2.D0*BETA12
P(2,2)=2.D0*BETA22
P(3,3)=2.D0*BETA66

DO K1=1,M
DO K2=1,M
      PSIGMA(K1)=PSIGMA(K1)+P(K1,K2)*SIGMA(K2)
   END DO
END DO

EES=0.D0
DO K1=1,M
   EES = EES + SIGMA(K1)*PSIGMA(K1)
END DO

EES = DSQRT(EES/2.D0)

DO K1=1,M
   PFPS(K1)=PSIGMA(K1)/(2.D0*EES)
END DO

   RETURN
   END
**
*********************************************************************
** CALCULATE SECOND DERIVATIVE FOR HILL48                      *
*********************************************************************
*USER SUBROUTINE
   SUBROUTINE CAL2DR(SIGMA,XR0,XR45,XR90,PAPS)
C
   INCLUDE 'ABA_PARAM.INC'
   PARAMETER(M=3)
   DIMENSION SIGMA(M),PAPS(M,M),P(M,M),PSIGMA(M),PFPS(M)

BETA12=XR0/(1.D0+XR0)
BETA66=(XR0+XR90)*(1.D0+2.D0*XR45)/(XR90*(1.D0+XR0))
BETA23=XR0/(XR90*(1.D0+XR0))
BETA22=BETA12+BETA23

   DO K1=1,M
DO K2=1,M
   P(K1,K2) = 0.D0
   PSIGMA(K1) = 0.D0
END DO
END DO

P(1,1)=2.D0*1.D0
P(1,2)=-2.D0*BETA12
P(2,1)=-2.D0*BETA12
P(2,2)=2.D0*BETA22
P(3,3)=2.D0*BETA66

DO K1=1,M
DO K2=1,M
      PSIGMA(K1)=PSIGMA(K1)+P(K1,K2)*SIGMA(K2)
   END DO
END DO

EES=0.D0
DO K1=1,M
   EES = EES + SIGMA(K1)*PSIGMA(K1)
END DO

EES = DSQRT(EES/2.D0)

DO K1=1,M
   PFPS(K1)=PSIGMA(K1)/(2.D0*EES)
END DO

DO K1=1,M
DO K2=1,M
   PAPS(K1,K2)=P(K1,K2)/(2.D0*EES)-PFPS(K1)*PFPS(K2)/EES
END DO
END DO

   RETURN
   END

赞一下(4人)

回复此楼

» 本帖已获得的红花（最新10朵）

自为胜利而来

5楼2019-04-14 17:42:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

自为胜利而来

至尊木虫 (著名写手)

应助: 6 (幼儿园)
金币: 13420.8
红花: 2
帖子: 1944
在线: 143.4小时
虫号: 3637225
注册: 2015-01-09
性别: GG
专业: 金属材料的力学行为

送红花一朵

引用回帖:

5楼: Originally posted by 吕泽华 at 2019-04-14 17:42:36
***********************************************************************************
** UMAT for Abaqus/Standard, large deformation formulation for plane stress **
** elements S4R. Implicit ...

非常感谢，学习。

回复此楼

道系青年

6楼2022-05-18 17:19:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主不惑可乐的主题更新

返回列表