شبیه‌سازی رشد عضله‌ دوسررانی‌ بر اثر کشش‌ با بکارگیری یک‌ مدل چندمقیاسه‌

شناخت‌ روند رشد ماهیچه‌ و تعیین‌ نواحی‌ بحرانی‌ تحت‌ آسیب‌ و یا پارگی‌ ابزاری برای تشخیص‌ روش صحیح‌ درمان برای متخصصان طب‌ فیزیکی‌ و توان بخشی‌ و ارتوپدی می‌باشند. هدف این‌ مقاله‌، بررسی‌ نحوه رشد سلول عضلانی‌اسکلتی‌ و هم‌چنین‌ تعیین‌ نواحی‌ آسیب‌ عضله‌ دوسررانی‌ تحت‌ کشش‌های غیر فعال وارده بر آن است‌. با تجزیه‌ تانسور گرادیان تغییر شکل‌ به‌ دو بخش‌ الاستیک‌ و رشد، روابط‌ رشد محدود برای ماهیچه‌ با رفتار ماده هایپرالاستیک‌ تعیین‌ شدهاند. روابط‌ مکانیک‌ محیط‌ پیوسته‌ با معادله‌ تکامل‌ رشد تلفیق‌ و معادلات دیفرانسیل‌ بیولوژیکی‌ و مکانیکی‌ حاصل‌ شدند. برای حل‌ آنها از روش اجزای محدود در نرمافزار آباکوس با نوشتن‌ زیربرنامه‌ای برای رفتار ماده در زبان فرترن استفاده شد. عضله‌ی دوسررانی‌ قسمت‌ سر بلند به‌ شکل‌ یک‌ استوانه‌ فرض شده و شبیه‌سازی آن برای یک‌ دوره ٦ هفته‌ای و به‌میزان ١٠% افزایش‌ طول اولیه‌ انجام شد. نتایج‌ نشان می‌دهند که‌ ماهیچه‌ به‌طور ناهمگن‌ رشد می‌کند و بیشینه‌ کشیدگی‌ها در رویه‌ها اتفاق می‌افتند و نه‌ در داخل‌ عضله‌ که‌ در ناحیه‌ بالا در رویه‌ بیرونی‌ معادل1/045 ودر ناحیه‌ پایین‌ در رویه‌ درونی‌ معادل 1/06 می‌باشند. به‌علاوه، نتایج‌ می‌تواند به‌نحوه درمان صحیح‌ و بهینه‌ و توان بخشی‌ بیماران و جراحی‌های ارتوپدی کمک‌ کند.

Simulation of Biceps Femoris Muscle Growth Based on Stretch Using a Multiscale Model

Understanding the process of muscle tissue growth is important to professionals who are involved in curing musculoskeletal disorders, physical medicine and rehabilitation specialists and orthopedic surgeons. This article investigates the development of a musculoskeletal cell and also determining the vulnerable areas of biceps femoris muscles due to passive strains applied on it. By decomposing the deformation gradient tensor to two parts, the elastic and growth, the finite growth relations have been applied for an isotropic hyperelastic muscle material behavior. Consequently, the continuum relations were combined with the growth evolution equation whrer a series of mechanobiological relations were obtained. To solve them, a FORTRAN userdefined material subroutine (UMAT) for the finite element Abaqus software was written and executed. The biceps femoris – long head muscle was simulated based on a 6week period assuming as a cylinder with 10% increase in initial length. Results of the simulation indicate that maximum strains occur in the surfaces, not inside the muscle. They reach 1.045 near the proximal muscletendon junction in the posterior layer and 1.06 in distal musclejunction in interior surface. Also, these results can help a correct and optimal treatment, patient’s rehabilitation and orthopedic surgeries.