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