abaqus修改inp直接建立工程
前面已经知道,通过修改以下inp的节点和单元编号,就可以新建模型,可是对于大的工程来说,逐个选取单元进行添加材料以及确定哪步进行填土仍是比较麻烦的(如果工程网格划分好并告知哪些单元好属于哪些材料,哪些单元号在哪步加载)。如果完全在abaqus中建模和计算就没必要看下面的介绍了。
首先建模:
1 *Heading 2 ** Job name: buildmodel Model name: Model-1 3 ** Generated by: Abaqus/CAE 6.14-1 4 *Preprint, echo=NO, model=NO, history=NO, contact=NO 5 ** 6 ** PARTS 7 ** 8 *Part, name=Part-1 9 *Node 10 1, 25., 15., 5. 11 2, 25., 10., 5. 12 3, 25., 5., 5. 13 4, 25., 0., 5. 14 5, 25., 15., 0. 15 6, 25., 10., 0. 16 7, 25., 5., 0. 17 8, 25., 0., 0. 18 9, 20., 15., 5. 19 10, 20., 10., 5. 20 11, 20., 5., 5. 21 12, 20., 0., 5. 22 13, 20., 15., 0. 23 14, 20., 10., 0. 24 15, 20., 5., 0. 25 16, 20., 0., 0. 26 17, 15., 15., 5. 27 18, 15., 10., 5. 28 19, 15., 5., 5. 29 20, 15., 0., 5. 30 21, 15., 15., 0. 31 22, 15., 10., 0. 32 23, 15., 5., 0. 33 24, 15., 0., 0. 34 25, 10., 15., 5. 35 26, 10., 10., 5. 36 27, 10., 5., 5. 37 28, 10., 0., 5. 38 29, 10., 15., 0. 39 30, 10., 10., 0. 40 31, 10., 5., 0. 41 32, 10., 0., 0. 42 33, 5., 15., 5. 43 34, 5., 10., 5. 44 35, 5., 5., 5. 45 36, 5., 0., 5. 46 37, 5., 15., 0. 47 38, 5., 10., 0. 48 39, 5., 5., 0. 49 40, 5., 0., 0. 50 41, 0., 15., 5. 51 42, 0., 10., 5. 52 43, 0., 5., 5. 53 44, 0., 0., 5. 54 45, 0., 15., 0. 55 46, 0., 10., 0. 56 47, 0., 5., 0. 57 48, 0., 0., 0. 58 *Element, type=C3D8R 59 1, 9, 10, 14, 13, 1, 2, 6, 5 60 2, 10, 11, 15, 14, 2, 3, 7, 6 61 3, 11, 12, 16, 15, 3, 4, 8, 7 62 4, 17, 18, 22, 21, 9, 10, 14, 13 63 5, 18, 19, 23, 22, 10, 11, 15, 14 64 6, 19, 20, 24, 23, 11, 12, 16, 15 65 7, 25, 26, 30, 29, 17, 18, 22, 21 66 8, 26, 27, 31, 30, 18, 19, 23, 22 67 9, 27, 28, 32, 31, 19, 20, 24, 23 68 10, 33, 34, 38, 37, 25, 26, 30, 29 69 11, 34, 35, 39, 38, 26, 27, 31, 30 70 12, 35, 36, 40, 39, 27, 28, 32, 31 71 13, 41, 42, 46, 45, 33, 34, 38, 37 72 14, 42, 43, 47, 46, 34, 35, 39, 38 73 15, 43, 44, 48, 47, 35, 36, 40, 39 74 *End Part 75 ** 76 ** 77 ** ASSEMBLY 78 ** 79 *Assembly, name=Assembly 80 ** 81 *Instance, name=Part-1-1, part=Part-1 82 *End Instance 83 ** 84 *End Assembly
看1好单元体:
*Element, type=C3D8R
1, 50, 14, 48, 82, 11, 1, 12, 46
在有限元中,节点编号呈逆时针排列,应该是看向z轴正方向时节点绕单元形心呈逆时针排列。
通过修改以上inp的节点和单元编号,就可以新建模型。
更进一步,给相应的单元以相应的材料:
在abaqus/cae中操作就是:
新建材料,并赋值参数:
然后建立截面,并对截面选择对应的材料属性或本构,截面是对结构的人设计的,比如梁的设计需要建梁并划分网格太麻烦,这个截面可以方便快速的把界面设计好。对于我们岩土来说这个就是走个流程,反倒是麻烦了一些。
下面将界面赋给对应的模型部分,对于我们岩土来说,相当于把对应的本构模型参数赋值给对应的模型,将来这部分模型划分网格后,这些网格就是这些材料了。
划分网格后,写出input文档,就可以看上面直接写inp和cae操作的对应了。
下面给个未完成的inp总结:
1 *Heading 2 ** Job name: exam Model name: Model-1 3 ** Generated by: Abaqus/CAE 6.14-5 4 *Preprint, echo=NO, model=NO, history=NO, contact=NO 5 ** 6 ** PARTS 7 ** 8 *Part, name=Part-1 9 *Node 10 1, 25., 0., 0. 11 2, 25., 15., 0. 12 3, 35., 0., 0. 13 4, 25., 0., 10. 14 5, 25., 15., 10. 15 6, 35., 0., 10. 16 7, 0., 0., 10. 17 8, 0., 0., 0. 18 9, 0., 15., 10. 19 10, 0., 15., 0. 20 11, 30., 0., 0. 21 12, 25., 5., 0. 22 13, 25., 10., 0. 23 14, 25., 0., 5. 24 15, 28.333334, 10., 0. 25 16, 31.666666, 5., 0. 26 17, 25., 15., 5. 27 18, 35., 0., 5. 28 19, 25., 5., 10. 29 20, 25., 10., 10. 30 21, 30., 0., 10. 31 22, 28.333334, 10., 10. 32 23, 31.666666, 5., 10. 33 24, 20., 0., 10. 34 25, 15., 0., 10. 35 26, 10., 0., 10. 36 27, 5., 0., 10. 37 28, 0., 0., 5. 38 29, 5., 0., 0. 39 30, 10., 0., 0. 40 31, 15., 0., 0. 41 32, 20., 0., 0. 42 33, 5., 15., 10. 43 34, 10., 15., 10. 44 35, 15., 15., 10. 45 36, 20., 15., 10. 46 37, 0., 15., 5. 47 38, 20., 15., 0. 48 39, 15., 15., 0. 49 40, 10., 15., 0. 50 41, 5., 15., 0. 51 42, 0., 5., 0. 52 43, 0., 10., 0. 53 44, 0., 10., 10. 54 45, 0., 5., 10. 55 46, 28.4895821, 4.53125, 0. 56 47, 27.2743053, 8.17708302, 0. 57 48, 25., 5., 5. 58 49, 25., 10., 5. 59 50, 30., 0., 5. 60 51, 28.333334, 10., 5. 61 52, 31.666666, 5., 5. 62 53, 28.489584, 4.53125, 10. 63 54, 27.2743073, 8.17708302, 10. 64 55, 20., 0., 5. 65 56, 15., 0., 5. 66 57, 10., 0., 5. 67 58, 5., 0., 5. 68 59, 20., 15., 5. 69 60, 15., 15., 5. 70 61, 10., 15., 5. 71 62, 5., 15., 5. 72 63, 20., 5., 0. 73 64, 15., 5., 0. 74 65, 10., 5., 0. 75 66, 5., 5., 0. 76 67, 20., 10., 0. 77 68, 15., 10., 0. 78 69, 10., 10., 0. 79 70, 5., 10., 0. 80 71, 20., 5., 10. 81 72, 15., 5., 10. 82 73, 10., 5., 10. 83 74, 5., 5., 10. 84 75, 20., 10., 10. 85 76, 15., 10., 10. 86 77, 10., 10., 10. 87 78, 5., 10., 10. 88 79, 0., 5., 5. 89 80, 0., 10., 5. 90 81, 27.2754631, 8.17361069, 5. 91 82, 28.4868832, 4.53935194, 5. 92 83, 20., 5., 5. 93 84, 15., 5., 5. 94 85, 10., 5., 5. 95 86, 5., 5., 5. 96 87, 20., 10., 5. 97 88, 15., 10., 5. 98 89, 10., 10., 5. 99 90, 5., 10., 5. 100 *Element, type=C3D8R 101 1, 50, 14, 48, 82, 11, 1, 12, 46 102 2, 82, 48, 49, 81, 46, 12, 13, 47 103 3, 49, 17, 51, 81, 13, 2, 15, 47 104 4, 51, 52, 82, 81, 15, 16, 46, 47 105 5, 52, 18, 50, 82, 16, 3, 11, 46 106 6, 50, 21, 53, 82, 14, 4, 19, 48 107 7, 82, 53, 54, 81, 48, 19, 20, 49 108 8, 54, 22, 51, 81, 20, 5, 17, 49 109 9, 54, 53, 82, 81, 22, 23, 52, 51 110 10, 53, 21, 50, 82, 23, 6, 18, 52 111 11, 12, 48, 83, 63, 1, 14, 55, 32 112 12, 48, 19, 71, 83, 14, 4, 24, 55 113 13, 63, 83, 84, 64, 32, 55, 56, 31 114 14, 83, 71, 72, 84, 55, 24, 25, 56 115 15, 64, 84, 85, 65, 31, 56, 57, 30 116 16, 84, 72, 73, 85, 56, 25, 26, 57 117 17, 65, 85, 86, 66, 30, 57, 58, 29 118 18, 85, 73, 74, 86, 57, 26, 27, 58 119 19, 66, 86, 79, 42, 29, 58, 28, 8 120 20, 86, 74, 45, 79, 58, 27, 7, 28 121 21, 13, 49, 87, 67, 12, 48, 83, 63 122 22, 49, 20, 75, 87, 48, 19, 71, 83 123 23, 67, 87, 88, 68, 63, 83, 84, 64 124 24, 87, 75, 76, 88, 83, 71, 72, 84 125 25, 68, 88, 89, 69, 64, 84, 85, 65 126 26, 88, 76, 77, 89, 84, 72, 73, 85 127 27, 69, 89, 90, 70, 65, 85, 86, 66 128 28, 89, 77, 78, 90, 85, 73, 74, 86 129 29, 70, 90, 80, 43, 66, 86, 79, 42 130 30, 90, 78, 44, 80, 86, 74, 45, 79 131 31, 2, 17, 59, 38, 13, 49, 87, 67 132 32, 17, 5, 36, 59, 49, 20, 75, 87 133 33, 38, 59, 60, 39, 67, 87, 88, 68 134 34, 59, 36, 35, 60, 87, 75, 76, 88 135 35, 39, 60, 61, 40, 68, 88, 89, 69 136 36, 60, 35, 34, 61, 88, 76, 77, 89 137 37, 40, 61, 62, 41, 69, 89, 90, 70 138 38, 61, 34, 33, 62, 89, 77, 78, 90 139 39, 41, 62, 37, 10, 70, 90, 80, 43 140 40, 62, 33, 9, 37, 90, 78, 44, 80 141 *Nset, nset=MatGiveSet-1 142 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 143 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 144 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 145 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 146 83, 84, 85, 86, 87, 88, 89, 90 147 *Elset, elset=MatGiveSet-1, generate 148 11, 40, 1 149 *Nset, nset=MatGiveSet-2 150 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 151 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 152 *Elset, elset=MatGiveSet-2, generate 153 1, 10, 1 154 ** Section: Section-2 155 *Solid Section, elset=MatGiveSet-2, material=E2 156 , 157 ** Section: Section-1 158 *Solid Section, elset=MatGiveSet-1, material=E1 159 , 160 *End Part 161 ** 162 ** 163 ** ASSEMBLY 164 ** 165 *Assembly, name=Assembly 166 ** 167 *Instance, name=Part-1-1, part=Part-1 168 *End Instance 169 ** 170 *End Assembly 171 ** 172 ** MATERIALS 173 ** 174 *Material, name=E1 175 *Elastic 176 210000., 0.3 177 *Material, name=E2 178 *Elastic 179 1e+06, 0.25 180 *Material, name=csuh1 181 *Depvar 182 20, 183 *User Material, constants=10 184 1.5, 0.3, 0.03, 0.125, 0.7, 1.2, 0., 0.6 185 0.4, 2.
下面接着进行,看堆载问题如何解决,工程中给出哪个计算步step堆载哪些单元了,那就把这些单元编号分组并写在对应的step中即可:
一般在abaqus中先建立step:
给激活和死掉的单元分组,并由step:
1 *Heading 2 ** Job name: exam Model name: Model-1 3 ** Generated by: Abaqus/CAE 6.14-5 4 *Preprint, echo=NO, model=NO, history=NO, contact=NO 5 ** 6 ** PARTS 7 ** 8 *Part, name=Part-1 9 *Node 10 1, 25., 0., 0. 11 2, 25., 15., 0. 12 3, 35., 0., 0. 13 4, 25., 0., 10. 14 5, 25., 15., 10. 15 6, 35., 0., 10. 16 7, 0., 0., 10. 17 8, 0., 0., 0. 18 9, 0., 15., 10. 19 10, 0., 15., 0. 20 11, 30., 0., 0. 21 12, 25., 5., 0. 22 13, 25., 10., 0. 23 14, 25., 0., 5. 24 15, 28.333334, 10., 0. 25 16, 31.666666, 5., 0. 26 17, 25., 15., 5. 27 18, 35., 0., 5. 28 19, 25., 5., 10. 29 20, 25., 10., 10. 30 21, 30., 0., 10. 31 22, 28.333334, 10., 10. 32 23, 31.666666, 5., 10. 33 24, 20., 0., 10. 34 25, 15., 0., 10. 35 26, 10., 0., 10. 36 27, 5., 0., 10. 37 28, 0., 0., 5. 38 29, 5., 0., 0. 39 30, 10., 0., 0. 40 31, 15., 0., 0. 41 32, 20., 0., 0. 42 33, 5., 15., 10. 43 34, 10., 15., 10. 44 35, 15., 15., 10. 45 36, 20., 15., 10. 46 37, 0., 15., 5. 47 38, 20., 15., 0. 48 39, 15., 15., 0. 49 40, 10., 15., 0. 50 41, 5., 15., 0. 51 42, 0., 5., 0. 52 43, 0., 10., 0. 53 44, 0., 10., 10. 54 45, 0., 5., 10. 55 46, 28.4895821, 4.53125, 0. 56 47, 27.2743053, 8.17708302, 0. 57 48, 25., 5., 5. 58 49, 25., 10., 5. 59 50, 30., 0., 5. 60 51, 28.333334, 10., 5. 61 52, 31.666666, 5., 5. 62 53, 28.489584, 4.53125, 10. 63 54, 27.2743073, 8.17708302, 10. 64 55, 20., 0., 5. 65 56, 15., 0., 5. 66 57, 10., 0., 5. 67 58, 5., 0., 5. 68 59, 20., 15., 5. 69 60, 15., 15., 5. 70 61, 10., 15., 5. 71 62, 5., 15., 5. 72 63, 20., 5., 0. 73 64, 15., 5., 0. 74 65, 10., 5., 0. 75 66, 5., 5., 0. 76 67, 20., 10., 0. 77 68, 15., 10., 0. 78 69, 10., 10., 0. 79 70, 5., 10., 0. 80 71, 20., 5., 10. 81 72, 15., 5., 10. 82 73, 10., 5., 10. 83 74, 5., 5., 10. 84 75, 20., 10., 10. 85 76, 15., 10., 10. 86 77, 10., 10., 10. 87 78, 5., 10., 10. 88 79, 0., 5., 5. 89 80, 0., 10., 5. 90 81, 27.2754631, 8.17361069, 5. 91 82, 28.4868832, 4.53935194, 5. 92 83, 20., 5., 5. 93 84, 15., 5., 5. 94 85, 10., 5., 5. 95 86, 5., 5., 5. 96 87, 20., 10., 5. 97 88, 15., 10., 5. 98 89, 10., 10., 5. 99 90, 5., 10., 5. 100 *Element, type=C3D8R 101 1, 50, 14, 48, 82, 11, 1, 12, 46 102 2, 82, 48, 49, 81, 46, 12, 13, 47 103 3, 49, 17, 51, 81, 13, 2, 15, 47 104 4, 51, 52, 82, 81, 15, 16, 46, 47 105 5, 52, 18, 50, 82, 16, 3, 11, 46 106 6, 50, 21, 53, 82, 14, 4, 19, 48 107 7, 82, 53, 54, 81, 48, 19, 20, 49 108 8, 54, 22, 51, 81, 20, 5, 17, 49 109 9, 54, 53, 82, 81, 22, 23, 52, 51 110 10, 53, 21, 50, 82, 23, 6, 18, 52 111 11, 12, 48, 83, 63, 1, 14, 55, 32 112 12, 48, 19, 71, 83, 14, 4, 24, 55 113 13, 63, 83, 84, 64, 32, 55, 56, 31 114 14, 83, 71, 72, 84, 55, 24, 25, 56 115 15, 64, 84, 85, 65, 31, 56, 57, 30 116 16, 84, 72, 73, 85, 56, 25, 26, 57 117 17, 65, 85, 86, 66, 30, 57, 58, 29 118 18, 85, 73, 74, 86, 57, 26, 27, 58 119 19, 66, 86, 79, 42, 29, 58, 28, 8 120 20, 86, 74, 45, 79, 58, 27, 7, 28 121 21, 13, 49, 87, 67, 12, 48, 83, 63 122 22, 49, 20, 75, 87, 48, 19, 71, 83 123 23, 67, 87, 88, 68, 63, 83, 84, 64 124 24, 87, 75, 76, 88, 83, 71, 72, 84 125 25, 68, 88, 89, 69, 64, 84, 85, 65 126 26, 88, 76, 77, 89, 84, 72, 73, 85 127 27, 69, 89, 90, 70, 65, 85, 86, 66 128 28, 89, 77, 78, 90, 85, 73, 74, 86 129 29, 70, 90, 80, 43, 66, 86, 79, 42 130 30, 90, 78, 44, 80, 86, 74, 45, 79 131 31, 2, 17, 59, 38, 13, 49, 87, 67 132 32, 17, 5, 36, 59, 49, 20, 75, 87 133 33, 38, 59, 60, 39, 67, 87, 88, 68 134 34, 59, 36, 35, 60, 87, 75, 76, 88 135 35, 39, 60, 61, 40, 68, 88, 89, 69 136 36, 60, 35, 34, 61, 88, 76, 77, 89 137 37, 40, 61, 62, 41, 69, 89, 90, 70 138 38, 61, 34, 33, 62, 89, 77, 78, 90 139 39, 41, 62, 37, 10, 70, 90, 80, 43 140 40, 62, 33, 9, 37, 90, 78, 44, 80 141 *Nset, nset=MatGiveSet-1 142 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 143 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 144 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 145 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 146 83, 84, 85, 86, 87, 88, 89, 90 147 *Elset, elset=MatGiveSet-1, generate 148 11, 40, 1 149 *Nset, nset=MatGiveSet-2 150 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 151 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 152 *Elset, elset=MatGiveSet-2, generate 153 1, 10, 1 154 ** Section: Section-2 155 *Solid Section, elset=MatGiveSet-2, material=E2 156 , 157 ** Section: Section-1 158 *Solid Section, elset=MatGiveSet-1, material=E1 159 , 160 *End Part 161 ** 162 ** 163 ** ASSEMBLY 164 ** 165 *Assembly, name=Assembly 166 ** 167 *Instance, name=Part-1-1, part=Part-1 168 *End Instance 169 ** 170 *Nset, nset=Set-2, instance=Part-1-1 171 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 172 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 173 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 174 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 175 83, 84, 85, 86, 87, 88, 89, 90 176 *Elset, elset=Set-2, instance=Part-1-1, generate 177 11, 40, 1 178 *Nset, nset=Set-3, instance=Part-1-1 179 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 180 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 181 *Elset, elset=Set-3, instance=Part-1-1, generate 182 1, 10, 1 183 *Nset, nset=allDeactSet-1, instance=Part-1-1, generate 184 1, 90, 1 185 *Elset, elset=allDeactSet-1, instance=Part-1-1, generate 186 1, 40, 1 187 *End Assembly 188 ** 189 ** MATERIALS 190 ** 191 *Material, name=E1 192 *Elastic 193 210000., 0.3 194 *Material, name=E2 195 *Elastic 196 1e+06, 0.25 197 *Material, name=csuh1 198 *Depvar 199 20, 200 *User Material, constants=10 201 1.5, 0.3, 0.03, 0.125, 0.7, 1.2, 0., 0.6 202 0.4, 2. 203 ** ---------------------------------------------------------------- 204 ** 205 ** STEP: Step-1 206 ** 207 *Step, name=Step-1, nlgeom=NO 208 *Static 209 1., 1., 1., 1. 210 ** 211 ** INTERACTIONS 212 ** 213 ** Interaction: deactInt-1 214 *Model Change, remove 215 allDeactSet-1, 216 ** 217 ** OUTPUT REQUESTS 218 ** 219 *Restart, write, frequency=0 220 ** 221 ** FIELD OUTPUT: F-Output-1 222 ** 223 *Output, field, variable=PRESELECT 224 ** 225 ** HISTORY OUTPUT: H-Output-1 226 ** 227 *Output, history, variable=PRESELECT 228 *End Step 229 ** ---------------------------------------------------------------- 230 ** 231 ** STEP: react1Step-2 232 ** 233 *Step, name=react1Step-2, nlgeom=NO 234 *Static 235 1., 1., 1e-05, 1. 236 ** 237 ** INTERACTIONS 238 ** 239 ** Interaction: react1Int-2 240 *Model Change, add 241 Set-2, 242 ** 243 ** OUTPUT REQUESTS 244 ** 245 *Restart, write, frequency=0 246 ** 247 ** FIELD OUTPUT: F-Output-1 248 ** 249 *Output, field, variable=PRESELECT 250 ** 251 ** HISTORY OUTPUT: H-Output-1 252 ** 253 *Output, history, variable=PRESELECT 254 *End Step 255 ** ---------------------------------------------------------------- 256 ** 257 ** STEP: react2Step-3 258 ** 259 *Step, name=react2Step-3, nlgeom=NO 260 *Static 261 1., 1., 1e-05, 1. 262 ** 263 ** INTERACTIONS 264 ** 265 ** Interaction: react2Int-3 266 *Model Change, add 267 Set-3, 268 ** 269 ** OUTPUT REQUESTS 270 ** 271 *Restart, write, frequency=0 272 ** 273 ** FIELD OUTPUT: F-Output-1 274 ** 275 *Output, field, variable=PRESELECT 276 ** 277 ** HISTORY OUTPUT: H-Output-1 278 ** 279 *Output, history, variable=PRESELECT 280 *End Step
添加体力和固定x=0的面:
施加边界条件后的阶段性inp:
1 *Heading 2 ** Job name: exam Model name: Model-1 3 ** Generated by: Abaqus/CAE 6.14-5 4 *Preprint, echo=NO, model=NO, history=NO, contact=NO 5 ** 6 ** PARTS 7 ** 8 *Part, name=Part-1 9 *Node 10 1, 25., 0., 0. 11 2, 25., 15., 0. 12 3, 35., 0., 0. 13 4, 25., 0., 10. 14 5, 25., 15., 10. 15 6, 35., 0., 10. 16 7, 0., 0., 10. 17 8, 0., 0., 0. 18 9, 0., 15., 10. 19 10, 0., 15., 0. 20 11, 30., 0., 0. 21 12, 25., 5., 0. 22 13, 25., 10., 0. 23 14, 25., 0., 5. 24 15, 28.333334, 10., 0. 25 16, 31.666666, 5., 0. 26 17, 25., 15., 5. 27 18, 35., 0., 5. 28 19, 25., 5., 10. 29 20, 25., 10., 10. 30 21, 30., 0., 10. 31 22, 28.333334, 10., 10. 32 23, 31.666666, 5., 10. 33 24, 20., 0., 10. 34 25, 15., 0., 10. 35 26, 10., 0., 10. 36 27, 5., 0., 10. 37 28, 0., 0., 5. 38 29, 5., 0., 0. 39 30, 10., 0., 0. 40 31, 15., 0., 0. 41 32, 20., 0., 0. 42 33, 5., 15., 10. 43 34, 10., 15., 10. 44 35, 15., 15., 10. 45 36, 20., 15., 10. 46 37, 0., 15., 5. 47 38, 20., 15., 0. 48 39, 15., 15., 0. 49 40, 10., 15., 0. 50 41, 5., 15., 0. 51 42, 0., 5., 0. 52 43, 0., 10., 0. 53 44, 0., 10., 10. 54 45, 0., 5., 10. 55 46, 28.4895821, 4.53125, 0. 56 47, 27.2743053, 8.17708302, 0. 57 48, 25., 5., 5. 58 49, 25., 10., 5. 59 50, 30., 0., 5. 60 51, 28.333334, 10., 5. 61 52, 31.666666, 5., 5. 62 53, 28.489584, 4.53125, 10. 63 54, 27.2743073, 8.17708302, 10. 64 55, 20., 0., 5. 65 56, 15., 0., 5. 66 57, 10., 0., 5. 67 58, 5., 0., 5. 68 59, 20., 15., 5. 69 60, 15., 15., 5. 70 61, 10., 15., 5. 71 62, 5., 15., 5. 72 63, 20., 5., 0. 73 64, 15., 5., 0. 74 65, 10., 5., 0. 75 66, 5., 5., 0. 76 67, 20., 10., 0. 77 68, 15., 10., 0. 78 69, 10., 10., 0. 79 70, 5., 10., 0. 80 71, 20., 5., 10. 81 72, 15., 5., 10. 82 73, 10., 5., 10. 83 74, 5., 5., 10. 84 75, 20., 10., 10. 85 76, 15., 10., 10. 86 77, 10., 10., 10. 87 78, 5., 10., 10. 88 79, 0., 5., 5. 89 80, 0., 10., 5. 90 81, 27.2754631, 8.17361069, 5. 91 82, 28.4868832, 4.53935194, 5. 92 83, 20., 5., 5. 93 84, 15., 5., 5. 94 85, 10., 5., 5. 95 86, 5., 5., 5. 96 87, 20., 10., 5. 97 88, 15., 10., 5. 98 89, 10., 10., 5. 99 90, 5., 10., 5. 100 *Element, type=C3D8R 101 1, 50, 14, 48, 82, 11, 1, 12, 46 102 2, 82, 48, 49, 81, 46, 12, 13, 47 103 3, 49, 17, 51, 81, 13, 2, 15, 47 104 4, 51, 52, 82, 81, 15, 16, 46, 47 105 5, 52, 18, 50, 82, 16, 3, 11, 46 106 6, 50, 21, 53, 82, 14, 4, 19, 48 107 7, 82, 53, 54, 81, 48, 19, 20, 49 108 8, 54, 22, 51, 81, 20, 5, 17, 49 109 9, 54, 53, 82, 81, 22, 23, 52, 51 110 10, 53, 21, 50, 82, 23, 6, 18, 52 111 11, 12, 48, 83, 63, 1, 14, 55, 32 112 12, 48, 19, 71, 83, 14, 4, 24, 55 113 13, 63, 83, 84, 64, 32, 55, 56, 31 114 14, 83, 71, 72, 84, 55, 24, 25, 56 115 15, 64, 84, 85, 65, 31, 56, 57, 30 116 16, 84, 72, 73, 85, 56, 25, 26, 57 117 17, 65, 85, 86, 66, 30, 57, 58, 29 118 18, 85, 73, 74, 86, 57, 26, 27, 58 119 19, 66, 86, 79, 42, 29, 58, 28, 8 120 20, 86, 74, 45, 79, 58, 27, 7, 28 121 21, 13, 49, 87, 67, 12, 48, 83, 63 122 22, 49, 20, 75, 87, 48, 19, 71, 83 123 23, 67, 87, 88, 68, 63, 83, 84, 64 124 24, 87, 75, 76, 88, 83, 71, 72, 84 125 25, 68, 88, 89, 69, 64, 84, 85, 65 126 26, 88, 76, 77, 89, 84, 72, 73, 85 127 27, 69, 89, 90, 70, 65, 85, 86, 66 128 28, 89, 77, 78, 90, 85, 73, 74, 86 129 29, 70, 90, 80, 43, 66, 86, 79, 42 130 30, 90, 78, 44, 80, 86, 74, 45, 79 131 31, 2, 17, 59, 38, 13, 49, 87, 67 132 32, 17, 5, 36, 59, 49, 20, 75, 87 133 33, 38, 59, 60, 39, 67, 87, 88, 68 134 34, 59, 36, 35, 60, 87, 75, 76, 88 135 35, 39, 60, 61, 40, 68, 88, 89, 69 136 36, 60, 35, 34, 61, 88, 76, 77, 89 137 37, 40, 61, 62, 41, 69, 89, 90, 70 138 38, 61, 34, 33, 62, 89, 77, 78, 90 139 39, 41, 62, 37, 10, 70, 90, 80, 43 140 40, 62, 33, 9, 37, 90, 78, 44, 80 141 *Nset, nset=MatGiveSet-1 142 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 143 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 144 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 145 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 146 83, 84, 85, 86, 87, 88, 89, 90 147 *Elset, elset=MatGiveSet-1, generate 148 11, 40, 1 149 *Nset, nset=MatGiveSet-2 150 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 151 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 152 *Elset, elset=MatGiveSet-2, generate 153 1, 10, 1 154 ** Section: Section-2 155 *Solid Section, elset=MatGiveSet-2, material=E2 156 , 157 ** Section: Section-1 158 *Solid Section, elset=MatGiveSet-1, material=E1 159 , 160 *End Part 161 ** 162 ** 163 ** ASSEMBLY 164 ** 165 *Assembly, name=Assembly 166 ** 167 *Instance, name=Part-1-1, part=Part-1 168 *End Instance 169 ** 170 *Nset, nset=Set-2, instance=Part-1-1 171 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 172 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 173 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 174 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 175 83, 84, 85, 86, 87, 88, 89, 90 176 *Elset, elset=Set-2, instance=Part-1-1, generate 177 11, 40, 1 178 *Nset, nset=Set-3, instance=Part-1-1 179 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 180 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 181 *Elset, elset=Set-3, instance=Part-1-1, generate 182 1, 10, 1 183 *Nset, nset=allDeactSet-1, instance=Part-1-1, generate 184 1, 90, 1 185 *Elset, elset=allDeactSet-1, instance=Part-1-1, generate 186 1, 40, 1 187 *Nset, nset=bodyforceSet-5, instance=Part-1-1, generate 188 1, 90, 1 189 *Elset, elset=bodyforceSet-5, instance=Part-1-1, generate 190 1, 40, 1 191 *Nset, nset=facex0Set-4, instance=Part-1-1 192 7, 8, 9, 10, 28, 37, 42, 43, 44, 45, 79, 80 193 *Elset, elset=facex0Set-4, instance=Part-1-1 194 19, 20, 29, 30, 39, 40 195 *End Assembly 196 ** 197 ** MATERIALS 198 ** 199 *Material, name=E1 200 *Elastic 201 210000., 0.3 202 *Material, name=E2 203 *Elastic 204 1e+06, 0.25 205 *Material, name=csuh1 206 *Depvar 207 20, 208 *User Material, constants=10 209 1.5, 0.3, 0.03, 0.125, 0.7, 1.2, 0., 0.6 210 0.4, 2. 211 ** 212 ** BOUNDARY CONDITIONS 213 ** 214 ** Name: xBC-1 Type: Displacement/Rotation 215 *Boundary 216 facex0Set-4, 1, 1 217 facex0Set-4, 2, 2 218 facex0Set-4, 3, 3 219 ** ---------------------------------------------------------------- 220 ** 221 ** STEP: Step-1 222 ** 223 *Step, name=Step-1, nlgeom=NO 224 *Static 225 1., 1., 1., 1. 226 ** 227 ** LOADS 228 ** 229 ** Name: Load-1 Type: Body force 230 *Dload 231 bodyforceSet-5, BX, -20. 232 ** 233 ** INTERACTIONS 234 ** 235 ** Interaction: deactInt-1 236 *Model Change, remove 237 allDeactSet-1, 238 ** 239 ** OUTPUT REQUESTS 240 ** 241 *Restart, write, frequency=0 242 ** 243 ** FIELD OUTPUT: F-Output-1 244 ** 245 *Output, field, variable=PRESELECT 246 ** 247 ** HISTORY OUTPUT: H-Output-1 248 ** 249 *Output, history, variable=PRESELECT 250 *End Step 251 ** ---------------------------------------------------------------- 252 ** 253 ** STEP: react1Step-2 254 ** 255 *Step, name=react1Step-2, nlgeom=NO 256 *Static 257 1., 1., 1e-05, 1. 258 ** 259 ** INTERACTIONS 260 ** 261 ** Interaction: react1Int-2 262 *Model Change, add 263 Set-2, 264 ** 265 ** OUTPUT REQUESTS 266 ** 267 *Restart, write, frequency=0 268 ** 269 ** FIELD OUTPUT: F-Output-1 270 ** 271 *Output, field, variable=PRESELECT 272 ** 273 ** HISTORY OUTPUT: H-Output-1 274 ** 275 *Output, history, variable=PRESELECT 276 *End Step 277 ** ---------------------------------------------------------------- 278 ** 279 ** STEP: react2Step-3 280 ** 281 *Step, name=react2Step-3, nlgeom=NO 282 *Static 283 1., 1., 1e-05, 1. 284 ** 285 ** INTERACTIONS 286 ** 287 ** Interaction: react2Int-3 288 *Model Change, add 289 Set-3, 290 ** 291 ** OUTPUT REQUESTS 292 ** 293 *Restart, write, frequency=0 294 ** 295 ** FIELD OUTPUT: F-Output-1 296 ** 297 *Output, field, variable=PRESELECT 298 ** 299 ** HISTORY OUTPUT: H-Output-1 300 ** 301 *Output, history, variable=PRESELECT 302 *End Step
上面的step步长调整可以计算后:
1 *Heading 2 ** Job name: exam Model name: Model-1 3 ** Generated by: Abaqus/CAE 6.14-5 4 *Preprint, echo=NO, model=NO, history=NO, contact=NO 5 ** 6 ** PARTS 7 ** 8 *Part, name=Part-1 9 *Node 10 1, 25., 0., 0. 11 2, 25., 15., 0. 12 3, 35., 0., 0. 13 4, 25., 0., 10. 14 5, 25., 15., 10. 15 6, 35., 0., 10. 16 7, 0., 0., 10. 17 8, 0., 0., 0. 18 9, 0., 15., 10. 19 10, 0., 15., 0. 20 11, 30., 0., 0. 21 12, 25., 5., 0. 22 13, 25., 10., 0. 23 14, 25., 0., 5. 24 15, 28.333334, 10., 0. 25 16, 31.666666, 5., 0. 26 17, 25., 15., 5. 27 18, 35., 0., 5. 28 19, 25., 5., 10. 29 20, 25., 10., 10. 30 21, 30., 0., 10. 31 22, 28.333334, 10., 10. 32 23, 31.666666, 5., 10. 33 24, 20., 0., 10. 34 25, 15., 0., 10. 35 26, 10., 0., 10. 36 27, 5., 0., 10. 37 28, 0., 0., 5. 38 29, 5., 0., 0. 39 30, 10., 0., 0. 40 31, 15., 0., 0. 41 32, 20., 0., 0. 42 33, 5., 15., 10. 43 34, 10., 15., 10. 44 35, 15., 15., 10. 45 36, 20., 15., 10. 46 37, 0., 15., 5. 47 38, 20., 15., 0. 48 39, 15., 15., 0. 49 40, 10., 15., 0. 50 41, 5., 15., 0. 51 42, 0., 5., 0. 52 43, 0., 10., 0. 53 44, 0., 10., 10. 54 45, 0., 5., 10. 55 46, 28.4895821, 4.53125, 0. 56 47, 27.2743053, 8.17708302, 0. 57 48, 25., 5., 5. 58 49, 25., 10., 5. 59 50, 30., 0., 5. 60 51, 28.333334, 10., 5. 61 52, 31.666666, 5., 5. 62 53, 28.489584, 4.53125, 10. 63 54, 27.2743073, 8.17708302, 10. 64 55, 20., 0., 5. 65 56, 15., 0., 5. 66 57, 10., 0., 5. 67 58, 5., 0., 5. 68 59, 20., 15., 5. 69 60, 15., 15., 5. 70 61, 10., 15., 5. 71 62, 5., 15., 5. 72 63, 20., 5., 0. 73 64, 15., 5., 0. 74 65, 10., 5., 0. 75 66, 5., 5., 0. 76 67, 20., 10., 0. 77 68, 15., 10., 0. 78 69, 10., 10., 0. 79 70, 5., 10., 0. 80 71, 20., 5., 10. 81 72, 15., 5., 10. 82 73, 10., 5., 10. 83 74, 5., 5., 10. 84 75, 20., 10., 10. 85 76, 15., 10., 10. 86 77, 10., 10., 10. 87 78, 5., 10., 10. 88 79, 0., 5., 5. 89 80, 0., 10., 5. 90 81, 27.2754631, 8.17361069, 5. 91 82, 28.4868832, 4.53935194, 5. 92 83, 20., 5., 5. 93 84, 15., 5., 5. 94 85, 10., 5., 5. 95 86, 5., 5., 5. 96 87, 20., 10., 5. 97 88, 15., 10., 5. 98 89, 10., 10., 5. 99 90, 5., 10., 5. 100 *Element, type=C3D8R 101 1, 50, 14, 48, 82, 11, 1, 12, 46 102 2, 82, 48, 49, 81, 46, 12, 13, 47 103 3, 49, 17, 51, 81, 13, 2, 15, 47 104 4, 51, 52, 82, 81, 15, 16, 46, 47 105 5, 52, 18, 50, 82, 16, 3, 11, 46 106 6, 50, 21, 53, 82, 14, 4, 19, 48 107 7, 82, 53, 54, 81, 48, 19, 20, 49 108 8, 54, 22, 51, 81, 20, 5, 17, 49 109 9, 54, 53, 82, 81, 22, 23, 52, 51 110 10, 53, 21, 50, 82, 23, 6, 18, 52 111 11, 12, 48, 83, 63, 1, 14, 55, 32 112 12, 48, 19, 71, 83, 14, 4, 24, 55 113 13, 63, 83, 84, 64, 32, 55, 56, 31 114 14, 83, 71, 72, 84, 55, 24, 25, 56 115 15, 64, 84, 85, 65, 31, 56, 57, 30 116 16, 84, 72, 73, 85, 56, 25, 26, 57 117 17, 65, 85, 86, 66, 30, 57, 58, 29 118 18, 85, 73, 74, 86, 57, 26, 27, 58 119 19, 66, 86, 79, 42, 29, 58, 28, 8 120 20, 86, 74, 45, 79, 58, 27, 7, 28 121 21, 13, 49, 87, 67, 12, 48, 83, 63 122 22, 49, 20, 75, 87, 48, 19, 71, 83 123 23, 67, 87, 88, 68, 63, 83, 84, 64 124 24, 87, 75, 76, 88, 83, 71, 72, 84 125 25, 68, 88, 89, 69, 64, 84, 85, 65 126 26, 88, 76, 77, 89, 84, 72, 73, 85 127 27, 69, 89, 90, 70, 65, 85, 86, 66 128 28, 89, 77, 78, 90, 85, 73, 74, 86 129 29, 70, 90, 80, 43, 66, 86, 79, 42 130 30, 90, 78, 44, 80, 86, 74, 45, 79 131 31, 2, 17, 59, 38, 13, 49, 87, 67 132 32, 17, 5, 36, 59, 49, 20, 75, 87 133 33, 38, 59, 60, 39, 67, 87, 88, 68 134 34, 59, 36, 35, 60, 87, 75, 76, 88 135 35, 39, 60, 61, 40, 68, 88, 89, 69 136 36, 60, 35, 34, 61, 88, 76, 77, 89 137 37, 40, 61, 62, 41, 69, 89, 90, 70 138 38, 61, 34, 33, 62, 89, 77, 78, 90 139 39, 41, 62, 37, 10, 70, 90, 80, 43 140 40, 62, 33, 9, 37, 90, 78, 44, 80 141 *Nset, nset=MatGiveSet-1 142 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 143 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 144 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 145 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 146 83, 84, 85, 86, 87, 88, 89, 90 147 *Elset, elset=MatGiveSet-1, generate 148 11, 40, 1 149 *Nset, nset=MatGiveSet-2 150 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 151 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 152 *Elset, elset=MatGiveSet-2, generate 153 1, 10, 1 154 ** Section: Section-2 155 *Solid Section, elset=MatGiveSet-2, material=E2 156 , 157 ** Section: Section-1 158 *Solid Section, elset=MatGiveSet-1, material=E1 159 , 160 *End Part 161 ** 162 ** 163 ** ASSEMBLY 164 ** 165 *Assembly, name=Assembly 166 ** 167 *Instance, name=Part-1-1, part=Part-1 168 *End Instance 169 ** 170 *Nset, nset=Set-2, instance=Part-1-1 171 1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 17, 19, 20, 24, 25 172 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41 173 42, 43, 44, 45, 48, 49, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64 174 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80 175 83, 84, 85, 86, 87, 88, 89, 90 176 *Elset, elset=Set-2, instance=Part-1-1, generate 177 11, 40, 1 178 *Nset, nset=Set-3, instance=Part-1-1 179 1, 2, 3, 4, 5, 6, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 180 21, 22, 23, 46, 47, 48, 49, 50, 51, 52, 53, 54, 81, 82 181 *Elset, elset=Set-3, instance=Part-1-1, generate 182 1, 10, 1 183 *Nset, nset=allDeactSet-1, instance=Part-1-1, generate 184 1, 90, 1 185 *Elset, elset=allDeactSet-1, instance=Part-1-1, generate 186 1, 40, 1 187 *Nset, nset=bodyforceSet-5, instance=Part-1-1, generate 188 1, 90, 1 189 *Elset, elset=bodyforceSet-5, instance=Part-1-1, generate 190 1, 40, 1 191 *Nset, nset=facex0Set-4, instance=Part-1-1 192 7, 8, 9, 10, 28, 37, 42, 43, 44, 45, 79, 80 193 *Elset, elset=facex0Set-4, instance=Part-1-1 194 19, 20, 29, 30, 39, 40 195 *End Assembly 196 ** 197 ** MATERIALS 198 ** 199 *Material, name=E1 200 *Elastic 201 210000., 0.3 202 *Material, name=E2 203 *Elastic 204 1e+06, 0.25 205 *Material, name=csuh1 206 *Depvar 207 20, 208 *User Material, constants=10 209 1.5, 0.3, 0.03, 0.125, 0.7, 1.2, 0., 0.6 210 0.4, 2. 211 ** 212 ** BOUNDARY CONDITIONS 213 ** 214 ** Name: xBC-1 Type: Displacement/Rotation 215 *Boundary 216 facex0Set-4, 1, 1 217 facex0Set-4, 2, 2 218 facex0Set-4, 3, 3 219 ** ---------------------------------------------------------------- 220 ** 221 ** STEP: Step-1 222 ** 223 *Step, name=Step-1, nlgeom=NO 224 *Static 225 0.01, 1., 1e-05, 0.1 226 ** 227 ** LOADS 228 ** 229 ** Name: Load-1 Type: Body force 230 *Dload 231 bodyforceSet-5, BX, -20. 232 ** 233 ** INTERACTIONS 234 ** 235 ** Interaction: deactInt-1 236 *Model Change, remove 237 allDeactSet-1, 238 ** 239 ** OUTPUT REQUESTS 240 ** 241 *Restart, write, frequency=0 242 ** 243 ** FIELD OUTPUT: F-Output-1 244 ** 245 *Output, field, variable=PRESELECT 246 ** 247 ** HISTORY OUTPUT: H-Output-1 248 ** 249 *Output, history, variable=PRESELECT 250 *End Step 251 ** ---------------------------------------------------------------- 252 ** 253 ** STEP: react1Step-2 254 ** 255 *Step, name=react1Step-2, nlgeom=NO 256 *Static 257 0.01, 1., 1e-05, 0.1 258 ** 259 ** INTERACTIONS 260 ** 261 ** Interaction: react1Int-2 262 *Model Change, add 263 Set-2, 264 ** 265 ** OUTPUT REQUESTS 266 ** 267 *Restart, write, frequency=0 268 ** 269 ** FIELD OUTPUT: F-Output-1 270 ** 271 *Output, field, variable=PRESELECT 272 ** 273 ** HISTORY OUTPUT: H-Output-1 274 ** 275 *Output, history, variable=PRESELECT 276 *End Step 277 ** ---------------------------------------------------------------- 278 ** 279 ** STEP: react2Step-3 280 ** 281 *Step, name=react2Step-3, nlgeom=NO 282 *Static 283 0.01, 1., 1e-05, 0.1 284 ** 285 ** INTERACTIONS 286 ** 287 ** Interaction: react2Int-3 288 *Model Change, add 289 Set-3, 290 ** 291 ** OUTPUT REQUESTS 292 ** 293 *Restart, write, frequency=0 294 ** 295 ** FIELD OUTPUT: F-Output-1 296 ** 297 *Output, field, variable=PRESELECT 298 ** 299 ** HISTORY OUTPUT: H-Output-1 300 ** 301 *Output, history, variable=PRESELECT 302 *End Step
由此可见,分布堆载时模型中间部分位移最大: