A nonlinear static analysis is performed in two load steps. Geometric nonlinearity is included in the analysis, and automatic time increments are used. The unsymmetric solver is used due to the large friction coefficient. For the first load step, when wear is not active, a larger time step is used.

/solu
antype,0                         ! static analysis
nlgeom,on                        ! Turn on Large Deformation Effects
_thickRatio=  0                  ! Ratio of thick parts in the model
eqsl,sparse,,,,,1
nropt,unsym          
! Use full Newton-Raphson with unsymmetric matrices as a large !friction coefficient (mu 0.8)
cntr,print,1                      ! print out contact info and also make no initial contact an error
dmpoption,emat,no                 ! Don't combine emat file for DANSYS
dmpoption,esav,no                 ! Don't combine esav file for DANSYS
nldiag,cont,subs                  ! print out contact info each substeps
rescontrol,define,last,last,,dele ! Program Controlled
/com,****************************************************
/com,******************* SOLVE FOR LS 1 OF 2 ****************
!Applying fluctuating stress field
esel,s,type,,7
nsle
sf,all,pres,%_loadvari55%
nsel,all
esel,all
!Applying compressive force
f,9213,fx,0.
f,9213,fy,fapp
!Applying occilatory motion
d,9213,ux,%_loadvari53ux%
d,9213,rotz,%_loadvari53rotz%
autots,on                         ! User turned on automatic time stepping
deltim,0.1,0.01,0.2,OFF
time,1.
outres,erase
outres,all,none
outres,nsol,all,
outres,rsol,all
outres,eangl,all
outres,etmp,all
outres,veng,all
outres,strs,all,
outres,epel,all,
outres,eppl,all,
outres,cont,all,
nldiag,cont,subs
cutcontrol,dsplimit,0.1
solve
/com *************** Write FE CONNECTORS *********
CEWRITE,file,ce,,INTE
/com ****************************************************
/com *************** FINISHED SOLVE FOR LS 1 *************
/com ****************************************************
/com ******************* SOLVE FOR LS 2 OF 2 ****************
allsel,all
autots,on                        ! User turned on automatic time stepping
deltim,0.002,0.0001,0.002,OFF
time,4.
frqout= -300                     !50 data points for each cycle
outres,erase
outres,all,none
outres,nsol,frqout,
outres,rsol,frqout,
outres,eangl,frqout,
outres,etmp,frqout
outres,veng,frqout
outres,strs,frqout,
outres,epel,frqout,
outres,eppl,frqout,
outres,cont,frqout,
solve

The repositioning of contact nodes during wear can result in changing contact status. If the wear increment is too large, all contact elements may go from a closed to an open status, resulting in rigid body motion. To prevent this, a very small time increment (or scale factor) is used so that the wear increment is also small and contact status changes are minimized. Contact data is written to the .cnd file at each substep for the tabulating scaled time. To avoid large output result files, the result is written at 300 equally spaced data points for the 6 simulation cycles, yielding 50 data points for each simulation cycle.