هدف مقاله حاضر توسعه يک فرمول بندی محاسباتی المان محدود نموی- انطباقی بر مبنای الگوريتمی توانمند به منظور محاسبه فشارهای شبه استاتيک تماسی در سازه های ويسکوالاستيک می باشد. از مدل تعميم يافته ماکسول برای مدل سازی توابع وارهيدگی معادلات متشکله ويسکوالاستيک بهره گرفته شده که در قالب دو تابع وارهيدگی کلی در اتساع (پاسخ حجمی) و برش (پاسخ انحرافی) نشان داده می شود. در اين فرمول بندی، توابع وارهيدگی به وسيله حاصل جمع سری های توابع نمايی کاهشی با زمان بيان می شوند. بر اساس اعمال اصل کار مجازی، يک فرمول بندی مؤثر المان های محدود با يک فرآيند نموسازی در معادلات وارهيدگی ويسکوالاستيک بدست آمده است. به منظور تحليل مسائل تماسی ويسکوالاستيک نيز، بر اساس شرايط سينماتيکی و هندسی اجزای در حال تماس، از رهيافت لاگرانژ الحاقی انطباقی استفاده گرديده است. همگرايی مناسب نتايج حاصل از مثالهای عددی با نتايج حل تحليلی نشانگر قابليتهای کاربردی بالای مدل المانهای محدود ارائه شده میباشد.