بررسی عددی رفتار الاستیک- پلاستیک و آسیب استخوان کورتیکال با اعمال یک مدل آسیب جدید

به‌دلیل نیاز‌هایی که در جراحی‌های ارتوپدی وجود دارد، رفتار مکانیکی استخوان کورتیکال در بار‌گذاری سیکلی و نرخ کرنش فیزیولوژیک، مورد بررسی قرار گرفته است. تاکید بر‌روی توسعه قانون ساختاری است که بتواند رفتار حین بارگذاری، باربرداری و بارگذاری مجدد که در آزمایش‌ها مشاهده شده را ایجاد کند. این مدل‌ها با ترکیب المان‌های رئولوژیک و اصول انرژی فرمول‌بندی خواهند شد. ابتدا دو مدل یک‌بعدی مستقل از نرخ کرنش فرمول‌بندی شده، یکی دارای یک متغیر آسیب و دیگری دارای دو متغیر آسیب مجزا در کشش و فشار است، مورد بررسی قرار گرفته و با استفاده از داده‌های آزمایشگاهی، ضرایب هر مدل به‌دست آمده است. با مقایسه نتایج شبیه‌سازی و داده‌های آزمایشگاهی، اصلاح لازم بر‌روی مدل‌ها صورت گرفته است. درنهایت با ترکیب معیار تسلیم ناهمسانگرد برسلر- پیستر و مدل یک‌بعدی مستقل از نرخ همراه دو متغیر آسیب، مدل سه‌بعدی متناظر حاصل شد. این مدل سه‌بعدی در قالب روش اجزاء محدود صریح پیاده‌سازی شده و نتیجه آن با نتایج شبیه‌سازی مدل یک‌بعدی و داده‌های تجربی سازگاری قابل قبولی نشان داد. این مدل سه‌بعدی مناسب شبیه‌سازی هندسه‌های پیچیده خواهد بود. ضریب‌تعیین برای مدل‌های یک‌بعدیri ، ±ri، ri و  ±ri اصلاح شده و مدل سه بعدی به ترتیب مقادیر 0/882174 ، 0/965665، 0/995508 ، 0/996279 و  0/984866  حاصل شده است.

numerical investigation of elastoplastic and damage behavior of cortical bone by applying a new damage model

due to the need for orthopedic surgery, the mechanical behavior of the cortical bone in cyclic loading and physiological strain rate has been investigated. the emphasis is on developing a structural law that can establish the behavior during loading, unloading and reloading observed in experiments. these models will be formulated by combining rheological elements and energy principles. first, two one-dimensional models independent of the strain rate are formulated, one with one damage variable and the other with two different damage variables in tension and compression, are examined, and using laboratory data, the coefficients of each model are obtained. by comparing the simulation results and laboratory data, the necessary modifications have been made to the models. finally, by combining the bresler-pister anisotropic yield criterion and the one-dimensional model independent of the rate associated with the two damage variables, the corresponding three-dimensional model was obtained. this three-dimensional model was implemented in the form of an explicit finite element method and the result showed acceptable compatibility with the simulation results of the one-dimensional model and experimental data. this three-dimensional model will be suitable for simulating complex geometries. the coefficient of determination for one-dimensional modelsri ,ri± ,ri , ri±and  has been modified and the three-dimensional model has obtained values of 0.882174, 0.965665, 0.995508, 0.996279, and 0.984866, respectively.