


{"id":449823,"date":"2026-01-06T14:39:07","date_gmt":"2026-01-06T14:39:07","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/nonlinear-adaptive-meshing-for-archard-contact-wear-simulation-is-not-working\/"},"modified":"2026-01-06T14:39:07","modified_gmt":"2026-01-06T14:39:07","slug":"nonlinear-adaptive-meshing-for-archard-contact-wear-simulation-is-not-working","status":"publish","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/nonlinear-adaptive-meshing-for-archard-contact-wear-simulation-is-not-working\/","title":{"rendered":"Nonlinear adaptive meshing for Archard contact wear simulation is not working"},"content":{"rendered":"<p>&lt;p&gt;I am following this example here:&nbsp;<a href=\"https:\/\/discuss.ansys.com\/home\/leaving?allowTrusted=1&amp;target=https%3A%2F%2Fansyshelp.ansys.com%2Fpublic%2Faccount%2Fsecured%3Freturnurl%3D%2F%2F%2FViews%2FSecured%2Fcorp%2Fv242%2Fen%2Fwb_wbtec%2FWBtecwear.html\" rel=\"nofollow\">(https:\/\/ansyshelp.ansys.com\/public\/account\/secured?returnurl=\/\/\/Views\/Secured\/corp\/v242\/en\/wb_wbtec\/WBtecwear.html)<\/a>. I am trying to simulate fretting wear for a polymer pin and a metal flat using a 2D Archard wear model. I have scaled the wear coefficient to simulate wear for 10 million cycles in 1s. Below are the commands for the asymmetric frictional contact (only pin wears, COF 0.35) and non-linear adaptive meshing. I was able to generate a stress contour when the wear is off, but the stress contour disappears for the step where wear is ON. The number of contacting elements stays the same, as no real movement is simulated, just the wear coefficient. The solution converges, but no change in mesh happens, so no change occurs in the volume of the bodies. Why can I not produce wear? Here is a graphic for the model:&nbsp;<a href=\"https:\/\/discuss.ansys.com\/home\/leaving?allowTrusted=1&amp;target=https%3A%2F%2Fimgur.com%2Fa%2FXBmTeSH\" rel=\"nofollow\">https:\/\/imgur.com\/a\/XBmTeSH<\/a>&lt;\/p&gt;&lt;p&gt;<img decoding=\"async\" src=\"https:\/\/innovationspace.ansys.com\/forum\/wp-content\/uploads\/sites\/2\/2026\/01\/06-01-2026-1767710169-Screenshot 2026-01-05 212754.jpg\" alt=\"\" \/>&lt;\/p&gt;&lt;p&gt;keyo,cid,5,1 !KEYOPT(5) = 1 is set to close the gap with an auto contact surface offset (CNOF)&lt;br&gt;keyo,cid,10,2 !KEYOPT(10) = 2 is set to perform a contact stiffness update each iteration&lt;\/p&gt;&lt;p&gt;! Wear model setup&lt;br&gt;TB,WEAR,PIN,,,ARCD&lt;br&gt;TBFIELD,TIME,0&lt;br&gt;TBDATA,PIN,0,1,1,0,0 ! Wear OFF initially&lt;br&gt;TBFIELD,TIME,1&lt;br&gt;TBDATA,PIN,0,1,1,0,0 ! Wear OFF&lt;br&gt;TBFIELD,TIME,1.01&lt;br&gt;TBDATA,PIN,K_ansys,H,1,1,0 ! Wear ON with scaled K&lt;br&gt;TBFIELD,TIME,2&lt;br&gt;TBDATA,PIN,K_ansys,H,1,1,0 ! Wear ON with scaled K&lt;\/p&gt;&lt;p&gt;ANTYPE,STATIC&lt;br&gt;NLGEOM,ON&lt;br&gt;KBC,0 !ramped load instead of stepped (KBC=1)&lt;\/p&gt;&lt;p&gt;!ESEL,S,TYPE,,PIN&lt;br&gt;!CM,wear,ELEM&lt;br&gt;!*GET,numelem,CM,cont_wear,ELEM ! Returns number of elements in component&lt;\/p&gt;&lt;p&gt;NLADAPTIVE,all,add,contact,wear,0.50 !wear at any contact point exceeds 50% of the average height of the solid element underlying the contact element&lt;br&gt;NLADAPTIVE,all,on,all,all,1,,2 !activate adaptivity&lt;br&gt;NLADAPTIVE,all,list,all,all !list adaptivity settings&lt;br&gt;&lt;br&gt;Listing the parameter just in case:&nbsp;&lt;\/p&gt;&lt;p&gt;! Parameters&lt;br&gt;AMP = 0.0001 ! 100 &micro;m amplitude&lt;br&gt;FREQ = 50 ! 50 Hz&lt;br&gt;CYCLES = 1E7 ! 10 million cycles&lt;\/p&gt;&lt;p&gt;! Calculations&lt;br&gt;S_per_cycle = 2<em>AMP ! 0.0002 m\/cycle&lt;br&gt;S_total = CYCLES*<\/em>S_per_cycle ! 2,000 m&lt;br&gt;T_exp = CYCLES \/ FREQ ! 200,000 seconds&lt;br&gt;T_sim_wear = 1 ! 1 second simulation wear time&lt;br&gt;SCALE = T_exp \/ T_sim_wear ! 200,000&lt;br&gt;V_avg = S_total \/ T_exp ! 0.01 m\/s&lt;\/p&gt;&lt;p&gt;! Material properties&lt;br&gt;K_base = 4.56E-7 ! Experimental dimensionless wear coefficient&lt;br&gt;H = 120E6 ! Hardness (Pa)&lt;\/p&gt;&lt;p&gt;! Scaled wear coefficient&lt;br&gt;K_ansys = K_base*<em>V_avg*<\/em>SCALE ! m^3\/N&middot;s&lt;\/p&gt;&lt;p&gt;*MSG,INFO&lt;br&gt;CALCULATIONS:&lt;br&gt;*MSG,INFO&lt;br&gt;AMP=%G m, S_total=%G m, V_avg=%G m\/s&lt;br&gt;*VWRITE,AMP,S_total,V_avg&lt;br&gt;(E10.3,E12.3,E12.5)&lt;br&gt;*MSG,INFO&lt;br&gt;Scale=%G, K_ansys=%G m&sup2;\/N&lt;br&gt;*VWRITE,SCALE,K_ansys&lt;br&gt;(E12.3,E12.3)&lt;\/p&gt;&lt;p&gt;!TBDATA,mat_id, K, H, m, n&lt;br&gt;! K= Wear coefficient, dimensionless&lt;br&gt;! H= hardness, N\/m^2&lt;br&gt;! m=1 !exponent on contact pressure&lt;br&gt;! n=1 !exponent on relative surface velocity&lt;\/p&gt;<\/p>\n","protected":false},"template":"","class_list":["post-449823","topic","type-topic","status-publish","hentry","topic-tag-adaptivemesh-2","topic-tag-archard-wear-1","topic-tag-frictional-contact"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_forum_id":["27791"],"_bbp_topic_id":["449823"],"_bbp_author_ip":["129.2.89.236"],"_bbp_last_reply_id":["450498"],"_bbp_last_active_id":["450498"],"_bbp_last_active_time":["2026-01-16 00:12:08"],"_bbp_reply_count":["12"],"_bbp_reply_count_hidden":["0"],"_bbp_voice_count":["3"],"_bbp_engagement":["636825","58821","240"],"_btv_view_count":["161"],"_bbp_topic_status":["unanswered"],"_edit_lock":["1767837454:636825"],"_bbp_subscription":["58821","636825","240"]},"test":"somaroyumd-edu"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/449823","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/449823\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=449823"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}