*HEADING
Direct cyclic analysis of a cylinder head under cyclic loadings
*RESTART,WRITE,FREQ=25
*INCLUDE,INPUT=dircyccylinderhead_dcm_mesh.inp
*INCLUDE,INPUT=dircyccylinderhead_dcm_sets.inp
*INCLUDE,INPUT=dircyccylinderhead_dcm_eqc.inp
*SOLID SECTION,ELSET=ZONE1SOL,MATERIAL=AS7V
*SOLID SECTION,ELSET=ZONE2SOL,MATERIAL=AS7V
*SOLID SECTION,ELSET=ZONE3SOL,MATERIAL=AS7NV
*SOLID SECTION,ELSET=PT000002,MATERIAL=MT000002
*SOLID SECTION,ELSET=PT000003,MATERIAL=MT000002
*SOLID SECTION,ELSET=PT000004,MATERIAL=MT000004
*SOLID SECTION,ELSET=PT000006,MATERIAL=MT000006
*SOLID SECTION,ELSET=PT000007,MATERIAL=MT000006
*MATERIAL,NAME=AS7V
*ELASTIC
70000, .33,  20.
70000,.33,  100.
65000,.33,  150.
65000, .33,  200.
62000, .33,  250.
57000, .33,  300.
*PLASTIC, HARDENING=KINEMATIC
62,0,20
595.5,.25,20
45,0,100
478.5,.25,100
32,0,150
315.5,.25,150
28,0,200
253,.25,200
26,0,250
176,.25,250
17,0,300
142,.25,300
*VISCOUS,LAW=STRAIN
3.1085e-13 , 7.4, 0.,  0.2857  ,   20.
7.8288e-16 , 7.4, 0.,  0.2857  ,  100.
6.4675e-17 , 7.8, 0.,  0.3077  ,  150.
4.9974e-17 , 7.8, 0.,  0.3846  ,  200.
1.0136e-16 , 7.4, 0.,  0.5161  ,  250.
3.1855e-16 , 7.8, 0.,  0.5614  ,  300.
*EXPANSION,TYPE=ISO,ZERO=20.
 22.6E-06 ,   20.                                                               
 22.6E-06 ,  175.                                                               
 25.1E-06 ,  225.                                                               
 25.1E-06 ,  400.     
*MATERIAL,NAME=AS7NV
*ELASTIC
70000, .33,  20.
70000,.33,  100.
65000,.33,  150.
65000, .33,  200.
62000, .33,  250.
57000, .33,  300.
*PLASTIC, HARDENING=KINEMATIC
62,0,20
595.5,.25,20
45,0,100
478.5,.25,100
32,0,150
315.5,.25,150
28,0,200
253,.25,200
26,0,250
176,.25,250
17,0,300
142,.25,300
*VISCOUS,LAW=STRAIN
3.1085e-13 , 7.4, 0.,  0.2857  ,   20.
7.8288e-16 , 7.4, 0.,  0.2857  ,  100.
6.4675e-17 , 7.8, 0.,  0.3077  ,  150.
4.9974e-17 , 7.8, 0.,  0.3846  ,  200.
1.0136e-16 , 7.4, 0.,  0.5161  ,  250.
3.1855e-16 , 7.8, 0.,  0.5614  ,  300.
*EXPANSION,TYPE=ISO,ZERO=20.
 22.6E-06 ,   20.                                                               
 22.6E-06 ,  175.                                                               
 25.1E-06 ,  225.                                                               
 25.1E-06 ,  400.     
*MATERIAL,NAME=MT000002                                                                   
*EXPANSION                                                                      
 13E-06, 300.                                                                   
 14E-06, 500.                                                                   
*ELASTIC                                                                        
 190.E+03, 0.3,  20.                                                            
 170.E+03, 0.3, 300.                                                            
 155.E+03, 0.3, 500.                                                            
*MATERIAL,NAME=MT000004                                                                          
*EXPANSION                                                                      
 17.2E-06,                                                                      
*ELASTIC                                                                        
 200.E+03, 0.3
*MATERIAL,NAME=MT000006                                                                       
*EXPANSION                                                                      
 20.5E-06                                                                      
*ELASTIC                                                                        
 106.E+03, 0.35
*BOUNDARY
BLOCVIS, 1, 1,0
BLOCVIS, 2, 2,0
BLOCVIS, 3, 3,0
*INITIAL CONDITIONS,TYPE=TEMPERATURE,FILE=dircyccylinderhead_heat.odb,STEP=6,INC=20
*STEP,INC=100
step 1 fit valve seats into valve ports
*STATIC
 1.,    1.,0.1
*BOUNDARY
100000, 1, 1 ,.0725
100001, 1, 1 ,.0810
100002, 1, 1 ,.025
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*NODE FILE,FREQ=0
*EL FILE,FREQ=0
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT
e,s
*END STEP
*STEP,INC=500
Step 2 Monocylinder under cyclic loadings
*DIRECT CYCLIC
.735,30,,,20,20,,100
*controls,type=direct cyclic
1,5e-3,1e-4,5e-3,1e-4
*TEMPERATURE,FILE=dircyccylinderhead_heat.odb,BSTEP= 5,BINC= 1,ESTEP= 6,EINC= 20
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*NODE FILE,FREQ=0
*EL FILE,FREQ=0
*OUTPUT,HISTORY,FREQ=1
*ELEMENT OUTPUT,ELSET=POSTRAIT
S11,
E11,
PE11
VE11
sener,
pener,
vener,
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
s,e
*END STEP