*heading MATERIAL COMBINATIONS: MAT1: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=NEWTONIAN MAT2: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=POWER LAW MAT3: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=CARREAU-YASUDA MAT4: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=CROSS MAT5: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=HERSCHEL-BULKLEY MAT6: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=ELLIS-METER MAT7: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=POWELL-EYRING MAT8: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=TABULAR MAT9: *EOS,TYPE=USUP+*VISCOSITY,DEFINITION=USER *preprint, model=yes *part, name=part-1 *node,nset=a1 1, 0., 0., 1. 2, 1., 0., 1. 3, 1., 0., 0. 4, 0., 0., 0. 5, 0., 1., 1. 6, 1., 1., 1. 7, 1., 1., 0. 8, 0., 1., 0. 11, 0., 2. 12, 1., 2. 13, 1., 3. 14, 0., 3. *node,nset=a2 101, 2., 0., 1. 102, 3., 0., 1. 103, 3., 0., 0. 104, 2., 0., 0. 105, 2., 1., 1. 106, 3., 1., 1. 107, 3., 1., 0. 108, 2., 1., 0. 111, 2., 2. 112, 3., 2. 113, 3., 3. 114, 2., 3. *node,nset=a3 201, 4., 0., 1. 202, 5., 0., 1. 203, 5., 0., 0. 204, 4., 0., 0. 205, 4., 1., 1. 206, 5., 1., 1. 207, 5., 1., 0. 208, 4., 1., 0. 211, 4., 2. 212, 5., 2. 213, 5., 3. 214, 4., 3. *node,nset=a4 301, 6., 0., 1. 302, 7., 0., 1. 303, 7., 0., 0. 304, 6., 0., 0. 305, 6., 1., 1. 306, 7., 1., 1. 307, 7., 1., 0. 308, 6., 1., 0. 311, 6., 2. 312, 7., 2. 313, 7., 3. 314, 6., 3. *node,nset=a5 401, 8., 0., 1. 402, 9., 0., 1. 403, 9., 0., 0. 404, 8., 0., 0. 405, 8., 1., 1. 406, 9., 1., 1. 407, 9., 1., 0. 408, 8., 1., 0. 411, 8., 2. 412, 9., 2. 413, 9., 3. 414, 8., 3. *node,nset=a6 501, 10., 0., 1. 502, 11., 0., 1. 503, 11., 0., 0. 504, 10., 0., 0. 505, 10., 1., 1. 506, 11., 1., 1. 507, 11., 1., 0. 508, 10., 1., 0. 511, 10., 2. 512, 11., 2. 513, 11., 3. 514, 10., 3. *node,nset=a7 601, 12., 0., 1. 602, 13., 0., 1. 603, 13., 0., 0. 604, 12., 0., 0. 605, 12., 1., 1. 606, 13., 1., 1. 607, 13., 1., 0. 608, 12., 1., 0. 611, 12., 2. 612, 13., 2. 613, 13., 3. 614, 12., 3. *node,nset=a8 701, 14., 0., 1. 702, 15., 0., 1. 703, 15., 0., 0. 704, 14., 0., 0. 705, 14., 1., 1. 706, 15., 1., 1. 707, 15., 1., 0. 708, 14., 1., 0. 711, 14., 2. 712, 15., 2. 713, 15., 3. 714, 14., 3. *node,nset=a9 801, 16., 0., 1. 802, 17., 0., 1. 803, 17., 0., 0. 804, 16., 0., 0. 805, 16., 1., 1. 806, 17., 1., 1. 807, 17., 1., 0. 808, 16., 1., 0. 811, 16., 2. 812, 17., 2. 813, 17., 3. 814, 16., 3. ** ** 3d solid element *element, type=c3d8r, elset=c3d1 1, 1, 2, 3, 4, 5, 6, 7, 8 *element, type=c3d8r, elset=c3d2 101, 101, 102, 103, 104, 105, 106, 107, 108 *element, type=c3d8r, elset=c3d3 201, 201, 202, 203, 204, 205, 206, 207, 208 *element, type=c3d8r, elset=c3d4 301, 301, 302, 303, 304, 305, 306, 307, 308 *element, type=c3d8r, elset=c3d5 401, 401, 402, 403, 404, 405, 406, 407, 408 *element, type=c3d8r, elset=c3d6 501, 501, 502, 503, 504, 505, 506, 507, 508 *element, type=c3d8r, elset=c3d7 601, 601, 602, 603, 604, 605, 606, 607, 608 *element, type=c3d8r, elset=c3d8 701, 701, 702, 703, 704, 705, 706, 707, 708 *element, type=c3d8r, elset=c3d9 801, 801, 802, 803, 804, 805, 806, 807, 808 ** ** plane strain element *element, type=cpe4r, elset=pe2d1 2, 11, 12, 13, 14 *element, type=cpe4r, elset=pe2d2 102, 111, 112, 113, 114 *element, type=cpe4r, elset=pe2d3 202, 211, 212, 213, 214 *element, type=cpe4r, elset=pe2d4 302, 311, 312, 313, 314 *element, type=cpe4r, elset=pe2d5 402, 411, 412, 413, 414 *element, type=cpe4r, elset=pe2d6 502, 511, 512, 513, 514 *element, type=cpe4r, elset=pe2d7 602, 611, 612, 613, 614 *element, type=cpe4r, elset=pe2d8 702, 711, 712, 713, 714 *element, type=cpe4r, elset=pe2d9 802, 811, 812, 813, 814 ** *solid section,elset=c3d1,material=mat1 *solid section,elset=pe2d1,material=mat1 ** *solid section,elset=c3d2,material=mat2 *solid section,elset=pe2d2,material=mat2 ** *solid section,elset=c3d3,material=mat3 *solid section,elset=pe2d3,material=mat3 ** *solid section,elset=c3d4,material=mat4 *solid section,elset=pe2d4,material=mat4 ** *solid section,elset=c3d5,material=mat5 *solid section,elset=pe2d5,material=mat5 ** *solid section,elset=c3d6,material=mat6 *solid section,elset=pe2d6,material=mat6 ** *solid section,elset=c3d7,material=mat7 *solid section,elset=pe2d7,material=mat7 ** *solid section,elset=c3d8,material=mat8 *solid section,elset=pe2d8,material=mat8 ** *solid section,elset=c3d9,material=mat9 *solid section,elset=pe2d9,material=mat9 ** *nset,nset=nall a1,a2,a3,a4,a5,a6,a7,a8,a9 *nset,nset=bot 1,2,3,4,11,12 101,102,103,104,111,112 201,202,203,204,211,212 301,302,303,304,311,312 401,402,403,404,411,412 501,502,503,504,511,512 601,602,603,604,611,612 701,702,703,704,711,712 801,802,803,804,811,812 *nset,nset=left 1,4,5,8,11,14 101,104,105,108,111,114 201,204,205,208,211,214 301,304,305,308,311,314 401,404,405,408,411,414 501,504,505,508,511,514 601,604,605,608,611,614 701,704,705,708,711,714 801,804,805,808,811,814 *nset,nset=top 5,6,7,8,13,14 105,106,107,108,113,114 205,206,207,208,213,214 305,306,307,308,313,314 405,406,407,408,413,414 505,506,507,508,513,514 605,606,607,608,613,614 705,706,707,708,713,714 805,806,807,808,813,814 *nset,nset=face,generate 1,801,100 2,802,100 5,805,100 6,806,100 *elset,elset=eall,generate 1,801,100 2,802,100 *elset,elset=qa_test eall, *end part ** *assembly, name=assembly-1 *instance, name=part-1-1, part=part-1 *end instance *end assembly ** *parameter c0 = 1480 s = 1.93 Gamma0 = 0.880 rho0 = 2070.0 *parameter eta_0 = 1.0 eta_inf = 0.1 a = 0.644 n = 0.392 lamda = 0.11 ** eta_max = eta_0 eta_min = 0.1 # k = (eta_inf+(eta_0-eta_inf)*2^((n-1)/a))* lamda^(n-1.0) k = (eta_inf+(eta_0-eta_inf)*exp(log(2)*(n-1)/a))*exp(log(lamda)*(n-1.0)) ** tau_0 = 3.59 ** # shrRate_half = ((0.5)^(a/(n-1.0))-1.0)^(1.0/a) / lamda shrRate_half = exp(log(exp(log(0.5)*(a/(n-1.0)))-1.0)*(1.0/a)) / lamda tau_half = 0.5*(eta_0+eta_inf)*shrRate_half ** t0 = 308.1 R = 8.31434 E0 = 109100 ** *physical constants,universal gas constant=, absolute zero=0.0 *material,name=mat1 *density *eos, type=usup , , *viscosity, definition=Newtonian *trs, definition=arrhenius , ** *material,name=mat2 *density *eos, type=usup , , *viscosity, definition=Power Law , , , *trs, definition=arrhenius , ** *material,name=mat3 *density *eos, type=usup , , *viscosity, definition=Carreau-Yasuda , , , , *trs, definition=arrhenius , ** *material,name=mat4 *density *eos, type=usup , , *viscosity, definition=Cross , , , *trs, definition=arrhenius , ** *material,name=mat5 *density *eos, type=usup , , *viscosity, definition=Herschel-Bulkley , , , *trs, definition=arrhenius , ** *material,name=mat6 *density *eos, type=usup , , *viscosity, definition=Ellis-Meter , , , *trs, definition=arrhenius , ** *material,name=mat7 *density *eos, type=usup , , *viscosity, definition=Powell-Eyring , , *trs, definition=arrhenius , ** *material,name=mat8 **,strain rate regularization=linear *density *eos, type=usup , , *viscosity, definition=Tabular 1.00000, 0.0 0.83383, 1.0 0.76532, 2.0 0.71776, 3.0 0.68112, 4.0 0.65134, 5.0 0.62631, 6.0 0.60477, 7.0 0.58593, 8.0 0.56921, 9.0 0.55422, 10.0 0.54066, 11.0 0.52830, 12.0 0.51697, 13.0 0.50652, 14.0 0.49684, 15.0 *trs, definition=arrhenius , ** *material,name=mat9 *density *eos, type=usup , , *viscosity, definition=User, Properties=3 , , *trs, definition=arrhenius , ** *initial conditions, type=temperature assembly-1.part-1-1.nall, 330.0 ** *boundary assembly-1.part-1-1.bot, 1, 3 assembly-1.part-1-1.top, 2, 3 ** *restart,timemarks=no,write,number=10 ** *amplitude, name=ramp 0.0, 0.0, 1.0, 15.0 ** *step *dynamic,explicit ,1.0 ** shear *boundary, type=velocity, amplitude=ramp assembly-1.part-1-1.top, 1, 1, 1.0 ** *output,field *element output,elset=assembly-1.part-1-1.qa_test s *node output u *output,history *element output,elset=assembly-1.part-1-1.qa_test s12, le12, temp *energy output, var=preselect ** *end step