[ABAQUS应用] 一个inp文件中调用多个mat的方法

leejack

　　在ABAUQS帮助文档UMAT或VUMAT中提到可以同时调用不同材料模型：

1. 　　IF (CMNAME(1:4) .EQ. 'MAT1') THEN
   
2. 
   
3. 　　CALL UMAT_MAT1(argument_list) .or. CALL VUMAT_MAT1(argument_list)
   
4. 
   
5. 　　ELSE IF(CMNAME(1:4) .EQ. 'MAT2') THEN
   
6. 
   
7. 　　CALL UMAT_MAT2(argument_list) .or. CALL VUMAT_MAT2(argument_list)
   
8. 
   
9. 　　END IF
   

复制代码

　　这里以VUMAT为例，
　　subroutine vumat(
　　C Read only -

1. 　　1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
   
2. 
   
3. 　　2 stepTime, totalTime, dt, cmname, coordMp, charLength,
   
4. 
   
5. 　　3 props, density, strainInc, relSpinInc,
   
6. 
   
7. 　　4 tempOld, stretchOld, defgradOld, fieldOld,
   
8. 
   
9. 　　3 stressOld, stateOld, enerInternOld, enerInelasOld,
   
10. 
    
11. 　　6 tempNew, stretchNew, defgradNew, fieldNew,
    
12. 
    
13. 　　C Write only -
    
14. 
    
15. 　　5 stressNew, stateNew, enerInternNew, enerInelasNew )
    
16. 
    
17. 　　C
    
18. 
    
19. 　　include 'vaba_param.inc'
    
20. 
    
21. 　　dimension props(nprops), density(nblock),
    
22. 
    
23. 　　1 coordMp(nblock,*),
    
24. 
    
25. 　　2 charLength(*), strainInc(nblock,ndir+nshr),
    
26. 
    
27. 　　3 relSpinInc(*), tempOld(*),
    
28. 
    
29. 　　4 stretchOld(*), defgradOld(*),
    
30. 
    
31. 　　5 fieldOld(*), stressOld(nblock,ndir+nshr),
    
32. 
    
33. 　　6 stateOld(nblock,nstatev), enerInternOld(nblock),
    
34. 
    
35. 　　7 enerInelasOld(nblock), tempNew(*),
    
36. 
    
37. 　　8 stretchNew(*), defgradNew(*), fieldNew(*),
    
38. 
    
39. 　　9 stressNew(nblock,ndir+nshr), stateNew(nblock,nstatev),
    
40. 
    
41. 　　1 enerInternNew(nblock), enerInelasNew(nblock),
    
42. 
    
43. 　　2 ustran(2),ustress(2),udstran(2),uddsdde(2,2)
    
44. 
    
45. 　　C
    
46. 
    
47. 　　character*80 cmname
    
48. 
    
49. 　　if (cmname(1:4) .eq. 'MAT1') then
    
50. 
    
51. 　　call VUMAT_mat1(
    
52. 
    
53. 　　C Read only -
    
54. 
    
55. 　　1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
    
56. 
    
57. 　　2 stepTime, totalTime, dt, cmname, coordMp, charLength,
    
58. 
    
59. 　　3 props, density, strainInc, relSpinInc,
    
60. 
    
61. 　　4 tempOld, stretchOld, defgradOld, fieldOld,
    
62. 
    
63. 　　3 stressOld, stateOld, enerInternOld, enerInelasOld,
    
64. 
    
65. 　　6 tempNew, stretchNew, defgradNew, fieldNew,
    
66. 
    
67. 　　C Write only -
    
68. 
    
69. 　　5 stressNew, stateNew, enerInternNew, enerInelasNew )
    
70. 
    
71. 　　else if (cmname(1:4) .eq. 'MAT2') then
    
72. 
    
73. 　　call VUMAT_mat2(
    
74. 
    
75. 　　C Read only -
    
76. 
    
77. 　　1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal,
    
78. 
    
79. 　　2 stepTime, totalTime, dt, cmname, coordMp, charLength,
    
80. 
    
81. 　　3 props, density, strainInc, relSpinInc,
    
82. 
    
83. 　　4 tempOld, stretchOld, defgradOld, fieldOld,
    
84. 
    
85. 　　3 stressOld, stateOld, enerInternOld, enerInelasOld,
    
86. 
    
87. 　　6 tempNew, stretchNew, defgradNew, fieldNew,
    
88. 
    
89. 　　C Write only -
    
90. 
    
91. 　　5 stressNew, stateNew, enerInternNew, enerInelasNew )
    
92. 
    
93. 　　end if
    
94. 
    
95. 　　!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    
96. 
    
97. 　　RETURN
    
98. 
    
99. 　　END
    

复制代码

　　其他就是将MAT1和MAT2复制在这后面就可以了。
　　注意： “if (cmname(1:4) .eq. 'MAT1') then”中MAT1必须是大写的，如果是小写，将在子程序运行中出现：

　　ERROR：BAD MATERIAL DEFINITION。


转自：http://blog.sina.com.cn/s/blog_642b4b710100ibdd.html