مقالة منهجية

نهج نمذجة العناصر المحدودة القائمة على الصور بالموجات فوق الصوتية داخل الأوعية الدموية لقياس الخصائص الميكانيكية في الجسم الحي للشريان التاجي البشري

1.4K مشاهدة

DOI:

10.3791/67573

ديسمبر 6, 2024

* These authors contributed equally

في هذه المقالة

ملخص

في الجسم الحي ، تظهر صور الموجات فوق الصوتية داخل الأوعية الدموية حركة المقطع العرضي التاجي المقابلة لظروف تحميل الضغط المختلفة. استنادا إلى نموذج العناصر المحدودة ، تم استخدام مخطط تكراري لتحديد الخصائص الميكانيكية الخاصة بالمريض للشرايين التاجية في الجسم الحي من خلال مطابقة حركة الشريان التاجي من النموذج الحسابي والصور الطبية.

الملخص

يمكن أن يوفر تحديد الخصائص الميكانيكية لجدران الشرايين التاجية معلومات مفيدة لتشخيص أمراض الشريان التاجي وإدارتها وعلاجها. نظرا لعدم توفر عينات الشريان التاجي الخاصة بالمريض للمرضى الذين يحتاجون إلى مراقبة مستمرة ، يصبح الاختبار التجريبي المباشر لخصائص مادة الأوعية الدموية مستحيلا. عادة ما تستخدم نماذج الشريان التاجي الحالية معلمات المواد من الأدبيات المتاحة ، مما يؤدي إلى أخطاء كبيرة في حساب الإجهاد / الإجهاد الميكانيكي. هنا ، سنقدم نهج التحديث القائم على نموذج العناصر المحدودة (FEMBUA) لتحديد خصائص المواد الخاصة بالمريض في الجسم الحي للشرايين التاجية بناء على الصور الطبية. في الجسم الحي ، تم الحصول على صور الموجات فوق الصوتية داخل الأوعية الدموية (IVUS) والأنسجة الافتراضية (VH) - IVUS للشرايين التاجية من مريض مصاب بمرض الشريان التاجي. تم تقسيم صور Cine IVUS التي تظهر حركة الأوعية الدموية خلال دورة قلبية واحدة ، وتم اختيار إطارين IVUS بمحيط التجويف الأقصى والحد الأدنى لتمثيل هندسة الشريان التاجي في ظل ظروف الضغط الانقباضي والانبساطي ، على التوالي. تم أيضا تقسيم صورة VH-IVUS للحصول على ملامح الأوعية الدموية ، وأضيف سمك طبقة يبلغ 0.05 سم إلى ملامح VH-IVUS لإعادة بناء هندسة الشريان التاجي. تم إنشاء نموذج عنصر محدود حسابي باستخدام نموذج مادة موني ريفلين متباين الخواص يستخدم لوصف الخصائص الميكانيكية للوعاء وظروف ضغط الدم النابض الموصوفة على السطح اللمعي التاجي لجعله ينقبض ويتوسع. بعد ذلك ، تم استخدام نهج التحديث التكراري لتحديد المعلمات المادية لنموذج Mooney-Rivlin متباين الخواص من خلال مطابقة محيط التجويف الأدنى والحد الأقصى من نموذج العناصر المحدودة الحسابية مع تلك الموجودة في صور IVUS السينمائية. يمكن توسيع نهج التحديث القائم على نموذج العناصر المحدودة المستند إلى الصور بنجاح لتحديد الخصائص المادية لجدران الشرايين في أسرة الأوعية الدموية المختلفة ويحمل إمكانية تقييم مخاطر أمراض القلب والأوعية الدموية.

المقدمة

يعد مرض الشريان التاجي (CAD) أحد الأسباب الرئيسية للوفيات والمراضة ، حيث يمثل أكثر من 9.14 مليون حالة وفاة في عام 2019 على مستوىالعالم 1،2. غالبا ما يكون تطور أمراض الشريان التاجي ، مثل تصلب الشرايين والتضيق ، مصحوبا بتغيرات في القوى الميكانيكية وتغيرات في خصائص مواد جدار الأوعيةالدموية 3. الخصائص المادية للشرايين التاجية ليست فقط حجر الزاوية لتحديد استجابتها الميكانيكية للحمل الفسيولوجي ولكن أيضا العناصر الرئيسية لمحاكاة السلوك الميكانيكي للأوعية الدموية ، والتنبؤ بتطور آفات تصلب الشرايين ، وتقييم التأثير العلاجي لمختلف الأجهزة الطبية4،5. وبالتالي ، فإن الفهم العميق والقياس الكمي الدقيق لخصائص المواد التاجية لها قيمة قصوى للتشخيص المبكر للمرض والطب الدقيق وتقييم التشخيص6.

التجارب الميكانيكية للأنسجة التاجية المعزولة ، مثل الاختبار ثنائي المحور المستوي ، واختبار المسافة البادئة ، وتمديد التضخم ، واختبار التمدد أحادي المحور ، هي طرق شائعة لتحديد الخواص الميكانيكية لجدران الأوعية التاجية خارج الجسم الحي7،8،9. من هذه الأساليب ، تم الحصول على عينات من الشريان التاجي من المرضى أو التجارب. تم إجراء اختبار ميكانيكي لتحديد استجابات الإجهاد لجدار الوعاء تحت ظروف الإجهاد المختلفة ، ثم تم تحديد معلمات المواد من خلال تركيب البيانات التجريبية10. أظهرت الدراسات السابقة أن خصائص الشريان التاجي غير خطية للغاية ومتباينةالخواص 11. على الرغم من أن التجارب خارج الجسم الحي يمكن أن توفر بيانات دقيقة عن خصائص المواد ، إلا أن هناك قيودا كبيرة أيضا ، وهي كما يلي: أولا ، سيكون السلوك الميكانيكي للعينة بعد إخراجها من الأشخاص الأحياء مختلفا عن السلوك الموجود في الظروف في الجسم الحي ، مما قد يؤثر على دقة نتائج الاختبار. ثانيا ، بسبب القيود الأخلاقية والعملية ، من الصعب الحصول على مجموعة كبيرة من الأنسجة الطبيعية أو المرضية للشرايين التاجية لإجراء الاختبار الميكانيكي.

للتغلب على هذه القيود ، استكشف الباحثون تقنيات جديدة للقياس الكمي في الجسم الحي ، وفي الوقت الفعلي ، والخاص بالمريض لخصائص المواد التاجية. من بينها ، يبشر نهج التحديث القائم على نموذج العناصر المحدودة (FEMBUA) القائم على الصورة الطبية بمعالجة هذه القضايا الصعبة. يستخدم هذا النهج تقنيات التصوير المتقدمة مثل الموجات فوق الصوتية داخل الأوعية الدموية (IVUS) وعلم الأنسجة الافتراضي (VH) -IVUS لالتقاط هندسة الشريان التاجي التفصيلية وتركيبات الأنسجة وحركتها12. من خلال إنشاء نماذج العناصر المحدودة ثلاثية الأبعاد (FE) ودمج ظروف ضغط الدم الفسيولوجية الخاصة بالمريض ، يمكن استعادة السلوك الوعائي الديناميكي أثناء دورات القلب عن طريق تحسين معلمات المواد لمطابقة بيانات الصورة من أجل القياس الكمي السريع والدقيق لخصائص المواد التاجية13. تشمل مزايا نهج تحديث FE في الجسم الحي على تجارب خارج الجسم الحي التقييم في الجسم الحي دون استئصال الأنسجة ، وتسهيل التقييمات واسعة النطاق ، ومحاكاة ديناميكيات الأوعية الدموية في ظل ظروف معقدة للمساعدة في فهم الفيزيولوجيا المرضية لأمراض الشريان التاجي.

في هذه الورقة ، يتم تقديم الخطوات الرئيسية لنهج التحديث القائم على نموذج العناصر المحدودة ، والتي تشمل تجزئة ومعالجة مفصلة لصورة IVUS و VH-IVUS السينمائية ، وإعادة بناء نموذج هيكل الطبقة الرقيقة الحسابي فقط ، وتنفيذ المخطط التكراري للبحث عن معلمات المواد المثلى لأنسجة الشرايين التاجية. الهدف من هذا البروتوكول هو تحديد الخصائص المادية للشريان التاجي من عينة مريضة مصابة بمرض الشريان التاجي باستخدام طريقة FEMBUA كعرض توضيحي ، وخاصة توضيح الطرق خطوة بخطوة. اختتمنا بمناقشة أهمية هذه الطريقة والجوانب الأخرى لهذه الطريقة في الجسم الحي .

المشاركة المختارة هي أنثى تبلغ من العمر 64 عاما وليس لها تاريخ سريري سابق لمرض الشريان التاجي. تم تشخيص هذا المريض بمرض الشريان التاجي بعد ظهور أعراض ألم في الصدر. تم إجراء تصوير الأوعية التاجية وفحص IVUS لتأكيد التشخيص. تم العثور على آفة لويحية مع تضيق بنسبة 60٪ في منتصف الشريان الهابط الأمامي الأيسر. بعد التقييم ، تم اعتماد العلاج الطبي الأمثل لعلاج المريض.

البروتوكول

تم الحصول على البيانات السريرية غير المحددة ، بما في ذلك صور IVUS في الجسم الحي وبيانات ضغط الدم ، من مريض مصاب بمرض الشريان التاجي في مستشفى Zhongda ، جامعة ساوث إيست ، مع الحصول على موافقة مستنيرة. تم اختيار عينة المريض من مجموعة المرضى في دراسة سريرية حول آفات تصلب الشرايين التاجية المتوسطة لتوضيح طريقة تحديد الخصائص المادية للأوعية التاجية الخاصة بالمريض14. أجريت الدراسة وفقا للبروتوكول المعتمد من قبل لجنة أخلاقيات البحث السريري في مستشفى تشونغدا ، جامعة ساوث إيست (رقم الموافقة: 2017ZDSYLL023-p01).

1. الحصول على البيانات ومعالجتها

  1. الحصول على صور Cine IVUS و VH-IVUS
    1. ضع قسطرة IVUS البعيدة عن آفة تصلب الشرايين بتوجيه من تصوير الأوعية التاجية واسحبها مرة أخرى إلى الجانب القريب. قم بإنشاء صور IVUS ذات تدرج رمادي على طول الطريق لتصور المقطع العرضي للوعاء التاجي.
    2. استخدم نظام التصوير IVUS المجهز بقسطرة بلاتينية عالية الأداء 20 ميجاهرتز و 2.9F للحصول على صور IVUS للجزء الشرياني التاجي مع لويحات تصلب الشرايين من المريض المصاب بمرض الشريان التاجي (الشكل 1).
    3. أثناء الحصول على صورة IVUS ، أوقف القسطرة مؤقتا في موقع الآفة المحدد مسبقا لمدة ثانية تقريبا للحصول على سلسلة من صور IVUS تسمى cine IVUS. أظهرت صور IVUS السينمائية بوضوح التغيرات المقطعية الديناميكية في موقع اللويحات المحددة خلال الدورة القلبية.
    4. قم بإنشاء صور VH-IVUS بناء على إطارات IVUS المكتسبة في لحظة ذروة R على تخطيط كهربية القلب لتصور مكونات البلاك في شكل مرمز بالألوان باستخدام نظام التصوير.
      ملاحظة: توفر صور VH-IVUS خريطة ألوان بديهية لأربعة مكونات رئيسية للبلاك في لويحات تصلب الشرايين: اللب النخر الغني بالدهون (الدهون) باللون الأحمر ، والتكلس باللون الأبيض ، والأنسجة الليفية باللون الأخضر الداكن ، والأنسجة الدهنية الليفية باللون الأخضر الفاتح.
    5. احفظ صور VH-IVUS و cine IVUS بتنسيق DICOM للتحليل دون اتصال بالإنترنت.
  2. تجزئة الصور ومعالجتها
    1. افتح ملفات DICOM باستخدام العارض، وانقر نقرا مزدوجا فوق اسم التسلسل المقابل لفتح الصورة وانقر فوق تصدير > تصدير الصور لحفظ كل إطار سينمائي IVUS أو إطار VH-IVUS كصورة فردية بتنسيق BMP. تحتوي كل صورة BMP على 500 × 500 بكسل ، كما هو موضح في ملف DICOM الأصلي.
    2. افحص صورة cine IVUS إطارا تلو الآخر للعثور على الإطارات المتتالية التي تم الحصول عليها في موقع البلاك المحدد مسبقا خلال دورة قلبية واحدة. كان هناك 26 إطارا للأنابيب السينمائية IVUS تم إنشاؤها خلال دورة قلبية واحدة لموقع البلاك العينة هذا.
    3. افحص جميع صور VH-IVUS التي تم إنشاؤها للعثور على صورة VH-IVUS التي تم الحصول عليها في موقع البلاك المحدد. تم إنشاء صور VH-IVUS المستخدمة هنا باستخدام إطار IVUS واحد من إطارات IVUS السينمائية في دورة قلبية واحدة.
    4. صور VH-IVUS و cine IVUS باستخدام برنامج ImageJ للحصول على ملامح حدود السفينة وحدود مكونات البلاك (انظر الشكل 1C).
      1. حدد علامة التبويب Straight > Segmented Line وحدد يدويا ملامح التجويف والحدود الخارجية للوعاء التاجي ومكونات البلاك على صور سينمائية IVUS و VH-IVUS. بالنسبة لصورة IVUS السينمائية ، قم فقط بتقسيم ملامح التجويف وحدود الوعاء الخارجي ، بينما بالنسبة لصورة VH-IVUS ، قم بتقسيم ملامح التجويف وحدود الوعاء الخارجي وحدود مكونات البلاك.
      2. من أجل البساطة ، احتفظ فقط بمكونات الدهون الكبيرة لإنشاء نموذج عنصر محدود وتجاهل المكونات الدهنية الصغيرة المعزولة. كان هناك دهون واحدة فقط في عينة البلاك هذه. تراكب الخطوط المحددة على الصور الأصلية باستخدام علامة التبويب Image > Overlay > Add Selection (إضافة تحديد ).
      3. انتقل إلى قائمة إلى مدير العائد على الاستثمار لإدارة الخطوط وضبط الخصائص وتعيين الألوان وعرض الخط على القيم المناسبة لتحسين المرئيات. حدد علامة التبويب خصائص ، وقم بتعيين لون الحد إلى لون مختلف، واملأ عرض الخط عند العرض. هنا ، قم بتعيين ألوان الخط على أنها أخضر وأزرق وأحمر لخطوط التجويف والحدود الخارجية والدهون ، على التوالي ، وعرض الخط على 3.
      4. قم بتنعيم الخطوط باستخدام Edit > Selection > Fit Spline من شريط الأوامر بعد تحديد محيط معين للحصول على محيط سلس. تستخدم هذه العملية تقنيات تركيب منحنى الشريحة لتنعيم الخطوط تلقائيا.
      5. انقر فوق ملف > حفظ باسم > علامة التبويب إحداثيات XY لحفظ إحداثيات النقاط لكل محيط مثل التجويف والحد الخارجي ومكون اللوحة في ملف txt منفصل. يحتوي هذا الملف على قيم الإحداثيات x و y للنقاط التي تشكل الكفاف ، مع وحدات البكسل كوحدة.
    5. سجل الحجم المادي الحقيقي لكل بكسل في صور cine IVUS و VH-IVUS (يشار إليها بحجم البكسل) من ملف DICOM. المسافة الحقيقية لبكسل واحد في بيانات IVUS المستخدمة هنا هي 0.002 سم. سيتم استخدام هذه المعلومات لتحويل إحداثيات النقطة مع البكسل كوحدة إلى المسافة الحقيقية مع سم كوحدة.
  3. معالجة البيانات الكنتورية
    1. معالجة بيانات كفاف Cine IVUS
      1. اقرأ ملفات txt لملامح التجويف من جميع صور IVUS السينمائية في دورة قلبية واحدة باستخدام MATLAB.
      2. اضرب جميع خطوط التجويف بحجم البكسل للحصول على الحجم الفعلي لخطوط التجويف.
      3. احسب محيطات التجويف لجميع خطوط التجويف وحدد إطارات IVUS بمحيط التجويف الأقصى (Cmax) والحد الأدنى (Cmin) ، والتي تمثل الظروف الانبساطي والانقباضي ، على التوالي.
    2. معالجة بيانات كنتور VH-IVUS
      1. اقرأ ملفات txt لملامح مكونات التجويف والحدود الخارجية والبلاك من صورة VH-IVUS باستخدام MATLAB.
      2. اضرب جميع الخطوط بحجم البكسل للحصول على الحجم الفعلي لجميع الخطوط العريضة.
      3. أعد تقسيم كل كفاف إلى 100 نقطة متباعدة بشكل متساو وقم بتنعيم 2D للحصول على بيانات كفاف VH-IVUS جديدة لتحل محل البيانات القديمة.

2. نموذج العناصر المحدودة

  1. إعادة بناء هندسة الأوعية التاجية
    1. قم بإنشاء طبقة واحدة من الخطوط في الفضاء ثلاثي الأبعاد عن طريق إضافة قيمة إحداثيات z لجميع نقاط خطوط VH-IVUS ، بما في ذلك التجويف والحدود الخارجية والدهون ، وتعيين z = 0 لجميع النقاط (الشكل 2).
    2. قم بإنشاء طبقة أخرى من الخطوط عن طريق إضافة قيمة إحداثيات z لجميع نقاط خطوط VH-IVUS وإعادة تعيين z = 0.05 سم لجميع النقاط.
      ملاحظة: تعيد هاتان الطبقتان من الخطوط بناء هندسة الأوعية التاجية ثلاثية الأبعاد لنموذج هيكل الطبقة الرقيقة فقط عن طريق إضافة سماكة طبقة ثابتة 0.05 سم إلى ملامح VH-IVUS (الشكل 2).
  2. توليد شبكة العناصر المحدودة
    1. قم بإنشاء معالم مساعدة عن طريق الاستيفاء الخطي للتجويف وملامح الحدود الخارجية بأوزان 1/3 و 2/3 (الشكل 3 أ) لكل طبقة.
    2. قسم مساحة الوعاء إلى 8 أجزاء محيطية و 3 أجزاء شعاعية (انظر الشكل 3 ب) عن طريق توصيل التجويف / الحدود الخارجية بأقرب نقطة على محيط الدهون (على سبيل المثال ، النقطتان A و B في الشكل 3 ب) أو اثنين من الخطوط المساعدة بخطوط شعاعية.
    3. قم بتوصيل جميع النقاط بين الطبقات بخطوط مستقيمة ، لتشكيل هيكل ثلاثي الأبعاد بأحجام 3 × 8 (الشكل 3 ب). قسم كل حجم باستخدام عناصر سداسية السطوح لإنشاء شبكة العناصر المحدودة (الشكل 3C) ومجموعات المواد المختلفة (الشكل 3D).
    4. قم بإجراء تحليل شبكي عن طريق تحسين كثافة الشبكة بنسبة 10٪ حتى < التغييرات في المحاليل 5٪.
  3. تعريف الخاصية المادية
    1. استخدم نموذج مادة موني ريفلين متباين الخواص المعدل لوصف الخصائص المادية لجدار الأوعية الدموية التاجية. كان من المفترض أن تكون الأوعية التاجية ومكونات اللويحات شديدة المرونة ، ومتباينة الخواص ، وغير قابلة للانضغاط تقريبا ، ومتجانسة ، ووظيفة كثافة طاقة الإجهاد لنموذج مادة موني ريفلين المعدل متباين الخواص هي:
      figure-protocol-1(1)
      figure-protocol-2(2)
      figure-protocol-3(3)
      حيث I1 و I2 هما الثوابت الأولى والثانية لموتر تشوه Cauchy-Green الأيمن C المعرف بأنه c = [cij] = XTX ، X = [Xij] = [figure-protocol-4] ، (Xi) كان الموضع الحالي (aj) هو الموضع الأصلي ، I4 = cij (nc) i (nc) j ، nc كان متجه الوحدة في الاتجاه المحيطي للوعاء. ج1 ، (ج)2 و D1 و D2 و K1 و K2 هي معلمات مادية خاصة بالمريض.
    2. قم بتعيين القيم الأولية لمعلمات المواد لوعاء تاجي خاص بالمريض وفقا لنتائج الاختبار ثنائي المحور خارج الجسم الحي ، أي c1 = −1،312.9 kPa ، c2 = 114.7 kPa ، D1 = 629.7 kPa ، D2 = 2.0 ، K1 = 35.9 kPa و K2 = 23.5 (الشكل 4A-B) 13،15.
    3. قم بتعيين معلمات المواد لمكون البلاك إن وجد. وبشكل أكثر تحديدا ، بالنسبة للدهون ، ج1 = 0.5 كيلو باسكال ، ج2 = 0 ، د1 = 0.5 كيلو باسكال ، و D2 = 1.5 ؛ للتكلس ، استخدمنا c1 = 920 كيلو باسكال ، ج2 = 0 ، D1 = 360 كيلو باسكال ، و د2 = 2.0 (الشكل 4 ب) 16.
      ملاحظة: تم افتراض أن مكونات البلاك (الدهون والتكلس) شديدة المرونة ، متناحية، وغير قابلة للانضغاط تقريبا ، وتم وصف خصائصها الميكانيكية بواسطة نموذج مادة موني ريفلين متناحي الخواص مع وظائف كثافة طاقة الإجهاد الواردة في الصيغة (2).
  4. تحديد المعادلات الحاكمة والشروط الحدودية
    1. تحديد المعادلات الحاكمة لنموذج هيكل الطبقة الرقيقة فقط ، والذي يتضمن معادلة الحركة ، وعلاقة الإجهاد والإزاحة غير الخطية Cauchy-Green ، ونموذج مادة الأوعية التاجية11.
    2. وصف أشكال موجات ضغط الدم الخاصة بالمريض على سطح التجويف لمحاكاة الظروف الفسيولوجية الحقيقية (الشكل 4 ج). للحصول على أشكال موجة ضغط الدم الخاصة بالمريض ، قم بقياس شكل موجة ضغط الأبهر النموذجي بقيم ضغط الانقباضي والانبساطي المقاسة بواسطة صفعة الذراع (الشكل 4 د).

3. نهج التحديث القائم على نموذج العناصر المحدودة لخصائص مادة الشريان التاجي الخاصة بالمريض

ملاحظة: يتم توضيح العملية التكرارية لتحديد خصائص المواد التاجية الخاصة بالمريض في الشكل 5.

  1. حدد هندسة عدم الحمل المقابلة لحالة الضغط الصفري كهندسة أولية للنموذج الحسابي عن طريق تقليص الهندسة التاجية المعاد بناؤها من صورة VH-IVUS محوريا بمعدل انكماش ثابت بنسبة 95٪ ومحيطيا مع انكماش محيطي (يشار إليه باسم S) تم تعيينه في البداية على أنه 98٪.
    ملاحظة: نظرا لأن الهندسة التاجية التي أعيد بناؤها من صورة VH-IVUS كانت في ظروف في الجسم الحي مع ضغط الدم الموصوف على التجويف والتمدد المحوري من الأجزاء الشريانية التاجية البعيدة والقريبة المربوطة ، يجب أن تتقلص الهندسة التاجية في الجسم الحي محيطيا ومحوريا للحصول على هندسة الضغط الصفر.
  2. قم بإصلاح معدل الانكماش المحوري عند 95٪ وقم بتحديث الانكماش المحيطي خلال الخطوات التالية.
  3. حدد نسبة المواد (المشار إليها ب k) لتعيين خصائص المواد الخاصة بالمريض للوعاء التاجي على النحو التالي: c1 = k * (−1،312.9) kPa ، c2 = k * 114.7 kPa ، D1 = k * 629.7 kPa ، K1 = k * 35.9 kPa ، وإصلاح D2 = 2.0 و K2 = 23.5.
    ملاحظة: نظرا لأنه تم الحصول على نقطتي بيانات فقط (الحد الأدنى والحد الأقصى لمحيط التجويف المقابل للضغوط الانبساطي والانقباضي) لتحديد المعلمات غير المعروفة (معدل الانكماش المحيطي S ومعلمات المواد لنموذج Mooney-Rivlin) ، فقد قمنا بتقليل عدد المعلمات غير المعروفة بافتراض أن الخصائص المادية الخاصة بالمريض في الجسم الحي للوعاء التاجي كانت متناسبة مع التخمين الأولي مع نسبة المواد المشار إليها باسم k: أي ، c1 = k * (−1،312.9) كيلو باسكال ، ج2 = ك * 114.7 كيلو باسكال ، د1 = ك * 629.7 كيلو باسكال ، ك1 = ك * 35.9 كيلو باسكال ، بينما D2 = 2.0 و K2 = 23.5 تم إصلاحها.
  4. قم بتحديث قيمة k المعينة إلى قيمة أولية k من 1 جنبا إلى جنب مع معدل الانكماش المحيطي S أثناء الإجراء التكراري التالي.
  5. قم بتشغيل برنامج لحل النموذج الحسابي للحصول على النتائج العددية.
    1. اكتب جميع الأوامر لإنشاء نموذج بنية الطبقة الرقيقة فقط في ملف دفعي (ملف تكميلي 1) باستخدام MATLAB.
    2. قم بتحميل هذا الملف الدفعي باستخدام واجهة المستخدم المتقدمة (AUI) لإنشاء النموذج (الشكل 6 أ). قم بحل نموذج هيكل الطبقة الرقيقة فقط بالنقر فوق ملف البيانات / الحل واحفظه كملف .dat (الشكل 6C). محاكاة ثلاث دورات قلبية واعتماد الحل في الدورة الأخيرة لتقديم النتائج العددية.
    3. قم بتصدير نتائج إحداثيات العقدة إلى ملف txt عن طريق التنقل في قائمة > قائمة القيم > المنطقة وتحديد X-POSITION وY-POSITION وZ-POSITION في المتغيرات المراد إدراجها ضمن الإحداثيات. انقر فوق تطبيق وتصدير لتصدير نتائج الإحداثيات.
    4. احفظ بيانات محيط التجويف المقابلة لظروف الضغط الانبساطي والانقباضي .txt الملفات لحسابات محيط التجويف.
  6. قارن محيط التجويف المحسوبة بواسطة نموذج FE (نموذج هيكل الطبقة الرقيقة فقط) في حالة الضغط الانبساطي مع بيانات IVUS السينمائية في الجسم الحي (Cmin) وتحقق مما إذا كان الخطأ النسبي <1٪). إذا تم استيفاء الشرط ، فانتقل إلى الخطوة التالية ، أو قم بتحديث نسبة المواد k باستخدام طريقة secant وانتقل إلى الخطوة 3.3 لإعادة التشغيل مرة أخرى17،18.
    ملاحظة: في التكرار الأول ، تم استخدام طريقة نيوتن لتحديث نسبة المواد بدلا من طريقة secant.
  7. قارن محيط التجويف المحسوبة بواسطة نموذج FE في حالة الضغط الانقباضي مع بيانات IVUS السينمائية في الجسم الحي (Cmax) وتحقق مما إذا كان الخطأ النسبي <1٪. إذا كانت الإجابة بنعم، فقم بإيقاف الإجراء التكراري، أو قم بتحديث معدل الانكماش المحيطي S وأعد التوجيه إلى الخطوة 3.4 لإعادة التشغيل مرة أخرى.
    ملاحظة: في التكرار الأول ، تم استخدام طريقة نيوتن لتحديث معدل الانكماش المحيطي بدلا من طريقة secant.
  8. سجل قيم S و k المثلى واحسب معلمات المواد المقابلة لنموذج مادة موني ريفلين.
  9. ارسم منحنيات نسبة الإجهاد والتمدد المحيطية والمحورية للوعاء التاجي (الشكل 7) ، والتي يمكن اشتقاقها على النحو التالي:
    figure-protocol-5(4)
    حيث يمثل σ إجهاد كوشي ، يمثل λ نسبة التمدد ، أنا = ج ، يمثل a الاتجاهات المحيطية والمحورية.
    1. لرسم منحنى مادي في اتجاه معين ، قم بإصلاح نسبة التمدد في الاتجاه الآخر إلى 1. احسب معامل يونغ الفعال في المحيط والمحوري (المسمى YMc و YMa ، على التوالي) كميل دالة المقياس لمنحنى المادة عند فاصل نسبة التمدد [1.0 ، 1.1] لتعكس صلابة المواد العامة للشريان التاجي13:
      figure-protocol-6(5)
  10. استخراج توزيعات الإجهاد / الإجهاد البلاك في أي وقت وتسجيل توزيع العقدة وقيم الإجهاد القصوى خلال المرحلتين الانقباضية والانبساطي (الشكل 8).

النتائج

نصف بالتفصيل طريقة FEMBUA ، والتي تتيح تحليل سريع لمواد البلاك والإجهاد للويحات التاجية بعد تصوير IVUS في الوقت الفعلي ويمكنها تحديد خصائص المواد في الجسم الحي والنتائج الميكانيكية الحيوية للويحات. يتم توفير معلمات المواد في الجسم الحي لنموذج مادة Mooney-Rivlin لهذا الوعاء التاجي في الجدول 1. تم رسم نتائج المحاكاة لنموذج العناصر المحدودة ، بما في ذلك توزيعات الإجهاد / الإجهاد في الوعاء التاجي ، في الشكل 8. وفيما يلي تحليل مفصل للنتائج.

يستغرق الأمر حوالي 2 ساعة لإنهاء العملية برمتها والحصول على معلمات المواد في الجسم الحي لهذا الوعاء التاجي ، بما في ذلك تجزئة الصور السينمائية IVUS و VH-IVUS ، وبناء نموذج FE الحسابي ، وتنفيذ الإجراء التكراري لتحديد خصائص المواد. يمكن أن تلبي هذه الطريقة متطلبات الفحص السريري السريع مع حسن التوقيت

بالنسبة لهذا الوعاء التاجي ، تم تحديد الكميتين باتباع طريقة FEMBUA على أنها k = 1.125 و S = 98٪. يتم توفير معلمات المواد الأخرى في الجسم الحي لنموذج مادة Mooney-Rivlin لهذا الوعاء التاجي في الجدول 1. يتم رسم منحنيات المواد في كل من الاتجاهين المحيطي والمحوري في الشكل 7. أظهر كلا الاتجاهين منحنى مادة على شكل حرف J ، مما يشير إلى علاقة الإجهاد والتمدد غير الخطية مع زيادة تصلب الأنسجة عندما ترتفع نسبة التمدد. كانت الخطوط المستقيمة هي وظائف المقياس لتناسب منحنيات المواد غير الخطية هذه ، وكانت منحدراتها هي وحدات يونغ الفعالة. YMc و YMa لهذا الوعاء التاجي هو 1055.41 و 1835.77 كيلو باسكال ، على التوالي.

تم رسم نتائج المحاكاة لنموذج العناصر المحدودة ، بما في ذلك توزيعات الإجهاد / الإجهاد في الوعاء التاجي ، في الشكل 8. توزيع الإجهاد المقابل لظروف الضغط الانقباضي أعلى من توزيع الضغط الانبساطي. أيضا ، يمكن ملاحظة تركيز الإجهاد الموضعي في المناطق ذات سمك الأوعية الرقيقة أو مناطق الغطاء الليفي. بالنظر إلى أن تمزق اللويحات غالبا ما يحدث في الغطاء الليفي ومناطق الكتف ، يتم الإبلاغ عن قيم ذروة الإجهاد والإجهاد فوق هذه المناطق الحرجة في الجدول 2.

أجرينا تحليل الشبكة عن طريق زيادة كثافة الشبكة تدريجيا ، وسيتوقف تحسين الشبكة إذا كانت الشبكة الدقيقة يمكن أن تؤدي إلى تغييرات ضئيلة في المحاليل. يتم عرض توزيع العنصر / الشبكة ونتائج الإجهاد لهذا النموذج في الجدول التكميلي 1. يعطي الجدول نتائج أكثر مقارنة وتفصيلا ، مما يبين اختلافات ضئيلة بين النماذج ذات أرقام العناصر المختلفة. تظهر النتائج أنه بعد تنقية كثافة الشبكة أربع مرات ، تكون الأخطاء النسبية للضغط الأقصى والحد الأدنى في حدود 5٪. نظرا لأن طريقة FEMBUA تتضمن عملية تكرارية ، فقد تم اعتماد كثافة شبكة منخفضة نسبيا لتوفير التكاليف الحسابية. وفقا للنتائج ، يمكن ملاحظة أن كمية الشبكة الحالية (رقم العنصر = 900) تلبي المتطلبات ولها تأثير ضئيل على نتائج الإجهاد.

figure-results-1
الشكل 1: الاستحواذ والاستخراج الكنتوري للسينما IVUS و VH-IVUS. (أ) اتجاه سحب القسطرة IVUS وموقع القسطرة المتوقفة مؤقتا (مميز بنجمة) للحصول على صور Cine IVUS. (ب) Cine IVUS مع محيط تجويف Cmin و Cmax وصورة VH-IVUS المقابلة. (ج) استخراج ملامح IVUS. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-2
الشكل 2: هندسة الأوعية التاجية ثلاثية الأبعاد المعاد بناؤها بسمك 0.05 سم. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-3
الشكل 3: عملية توليد شبكة العناصر المحدودة لنموذج هيكل الشريحة الرقيقة فقط. (أ) نقاط كفاف على مستوى ثنائي الأبعاد. (ب) الخطوط الهندسية لتوليد الشبكات. (ج) شبكة العناصر المحدودة. (د) أنشأت مجموعات مادية مختلفة. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-4
الشكل 4: نموذج هيكل الشريحة الرقيقة الحسابي فقط للوعاء التاجي. يتضمن النموذج تعريف خاصية المواد وتحديد حالة الحدود. (أ) وضع مجموعات مادية مختلفة داخل الشريان التاجي. (ب) منحنى مواد الأوعية التاجية خارج الجسم الحي . (ج) تحميل حالة الحدود. (د) منحنى الضغط. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-5
الشكل 5: مخطط تكراري. المخطط التكراري لتحديد خصائص مادة البلاك في الجسم الحي للأوعية التاجية بناء على صور IVUS. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-6
الشكل 6: استخدم الواجهة. واجهة المستخدم المتقدمة لبرنامج ADINA والعمليات ذات الصلة. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-7
الشكل 7: منحنيات المواد. منحنيات المواد في الاتجاهات المحيطية والمحورية للوعاء التاجي وتقدير وحدات يونغ الفعالة. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

figure-results-8
الشكل 8: نموذج هيكل شريحة رقيقة فقط لهندسة الأوعية التاجية ونتائج محاكاة. يتضمن الشكل خريطة سحابة توزيع الإجهاد في ظروف الضغط الانبساطي والانقباضي. (أ) العقد المعروضة. (ب) الإجهاد الانبساطي. (ج) الإجهاد الانقباضي. يشير شريط الألوان إلى قيمة الإجهاد. الرجاء النقر هنا لعرض نسخة أكبر من هذا الرقم.

معلمات الموادS(%)ج1 (كيلو باسكال)ج2 (كيلو باسكال)د1 (كيلو باسكال)K1 (كيلو باسكال)YMa (كيلو باسكال)YMc (كيلو باسكال)
98%-147.7612.91708.6740.41055.411835.77

الجدول 1: القيم المرجعية لمعلمات المواد لنموذج العينة.

الوقت في الضغط الانبساطيالوقت في الضغط الانقباضي
PWS (كيلو باسكال)PWSnPWS (كيلو باسكال)PWSn
56.530.0540264.510.05402

الجدول 2: قيم ذروة الإجهاد / الإجهاد على المناطق الحرجة من النموذج الحسابي. الاختصارات: ذروة إجهاد الجدار = PWS; ذروة إجهاد الجدار = PWSn. وحدة الإجهاد: كيلو باسكال.

الجدول التكميلي 1: ينتج عن إجهاد المقارنة كميات شبكية مختلفة في نفس النموذج.> الرجاء النقر هنا لتنزيل هذا الملف.

الملف التكميلي 1: رموز MATLAB. الرجاء النقر هنا لتنزيل هذا الملف.

المناقشة

الخطوات الحاسمة في البروتوكول
تكمن الخطوة الأكثر أهمية في نهج التحديث المستند إلى نموذج العناصر المحدودة في الإجراء التكراري. في هذا النهج ، يجب أن يستعيد نموذج العناصر المحدودة بدقة حركة الأوعية التاجية على المقطع العرضي للأوعية الدموية من صور IVUS السينمائية في الجسم الحي . لهذا الغرض ، تم اعتماد تقليل فرق محيط التجويف بين نموذج العناصر المحدودة والصور في الجسم الحي في هذه الدراسة للعثور على خصائص المواد المناسبة. كانت هناك خطوات رئيسية أخرى في البروتوكولات ، بما في ذلك تجزئة الصور ومعالجتها ، وبناء نموذج العناصر المحدودة ، والمخططات التكرارية لتحديد خصائص مادة الأوعية التاجية. كان التحسين هو تنفيذ التحديد التلقائي لمكونات البلاك التاجي لتوفير الوقت.

أهمية نهج التحديث القائم على نمذجة العناصر المحدودة
بشكل كلاسيكي ، تم إجراء تجارب ميكانيكية مثل اختبار الشد أحادي المحور / ثنائي المحور ، واختبار المسافة البادئة ، واختبار تضخم الضغط لتحديد السلوك الميكانيكي لأنسجة القلب والأوعية الدموية خارج الجسم الحي8،16. كان التطبيق السريري لهذه الأساليب محدودا للأسباب التالية: 1) غالبا ما تكون عينات الأنسجة الشريانية التاجية خارج الجسم الحي غير متوفرة. يكاد يكون من المستحيل حصاد عينات الشريان التاجي الطبيعية في البيئة السريرية. 2) قد تتغير الخصائص المادية للأنسجة الشريانية عند إخراجها من الأشخاص الأحياء. 3) إنه غير مناسب للمراقبة المستمرة في بيئة خاصة بالمريض للإدارة الشخصية والطب الدقيق. لحسن الحظ ، يوفر FEMBUA طريقة أخرى لتحديد خصائص الأنسجة الخاصة بالمريض في الجسم الحي. يمكن تعديل طريقة in vivo بسهولة لتطبيقها بنجاح على الأنسجة البيولوجية الأخرى ، مثل أنسجة الأبهر19 وأنسجة القلب20 ، وحتى المواد غير البيولوجية مثل المعدن21. أظهرت الدراسات السابقة أن خصائص مواد الأنسجة الخاصة بالمريض في الجسم الحي كان لها تأثير كبير على الميكانيكا الحيوية للقلب والأوعية الدموية ، خاصة في حساب الإجهاد ، مقارنة بخصائص المواد خارج الجسم الحي 11. لذلك ، فإن خصائص المواد الخاصة بالمريض في الجسم الحي مرغوبة للعلاج الشخصي.

على الرغم من اختلاف طرق in vivo و ex vivo ، إلا أنه يمكن دمجها لإلهام مناهج أخرى لتحديد الخصائص الميكانيكية لأنسجة القلب والأوعية الدموية. تم اقتراح نهج هجين يجمع بين FEMBUA لمطابقة بيانات الإجهاد والإجهاد من التجارب ثنائية المحور / أحادية المحور لأنسجة تمدد الأوعية الدموية الأبهري22.

التطبيقات السريرية المحتملة لنهج التحديث القائم على نمذجة العناصر المحدودة
تعتبر FEMBUA القائمة على صورة IVUS ذات أهمية كبيرة في علوم مواد القلب والأوعية الدموية ، وتحليل الصور الطبية ، وتصميم الأجهزة الطبية الشخصية. لا تحتاج هذه الطريقة المقترحة إلى إحداث أضرار جسيمة لجدار الأوعية الدموية التاجية التي تم فحصها كما هو الحال في حالة مثل إخراج الأنسجة من جسم الإنسان الحي ، لذا فهي مناسبة للمراقبة المستمرة للمريض ، مثل التحقيق في تأثير الخصائص الميكانيكية في الجسم الحي على تشخيص المريض ، والذي لا يمكن القيام به من خلال نهج تجريبي ميكانيكي كلاسيكي. بالإضافة إلى ذلك ، تلعب طريقة FEMBUA المقترحة دورا رئيسيا في تحسين علاج دعامات الشريان التاجي وتشخيص أمراض القلب والأوعية الدموية. يمكن أن يؤدي التحليل الدقيق لخصائص الهيكل والآفة لجدران الأوعية الدموية إلى توجيه اختيار حجم الدعامة وموقعها ، وبالتالي الحفاظ على الاستقرار الهيكلي والميكانيكي داخل الأنسجة بدرجة أعلى ، وتقليل المضاعفات ، وتحسين تشخيص المريض23. باختصار ، تتمتع طريقة FEMBUA بإمكانيات تطبيق واسعة النطاق وبعيدة المدى في مجال القلب والأوعية الدموية.

مقارنة النتائج من FEMBUA والنهج التجريبية خارج الجسم الحي
تم إجراء تحليل مقارنة بين FEMBUA والأساليب التجريبية الكلاسيكية خارج الجسم الحي لتقييم دقة وفعالية النهج الجديد في الجسم الحي . من أجل التبسيط ، تمت مقارنة تصلب الأنسجة التاجية من النهج الميكانيكي بتلك التي أجريت في الدراسات في الجسم الحي باستخدام FEMBUA ، وكانت تصلب الأنسجة من كلتا الطريقتين بشكل عام في نفس نطاقالحجم 11. تم تأكيد الاتساق أيضا في أسرة الأوعية الدموية الأخرى مثل الشريان الأورطي والشريان السباتي24. تجدر الإشارة إلى أن العدد المحدود من الدراسات خارج الجسم الحي يمكن أن يؤثر على الاستنتاجات المذكورة أعلاه ، حيث أن تباين تصلب الأنسجة بين الأفراد المختلفين واضح أيضا. ومع ذلك ، أشارت هذه الاستنتاجات إلى أن FEMBUA هو نهج دقيق وفعال لقياس الخصائص المادية لجدران الشرايين.

كما تم التحقيق في التحقق من صحة ومتانة طريقة FEMBUA من خلال إجراء كل من طريقة FEMBUA في الجسم الحي وتجربة خارج الجسم الحي على نفس الأنسجة الشريانية لأغراض التحققمن الصحة 25،26،27. أجرى Liu et al. في الجسم الحي والخارج في الجسم الحي مناهج تجريبية على أنسجة الأبهر. كشفت النتائج التي توصلوا إليها عن وجود علاقة وثيقة بين منحنيات السلوك المادي الناتجة عن كلتا الطريقتين ، بمتوسط النسبة المئوية المطلقة للخطأ أقل من 5٪ 25. علاوة على ذلك ، Cosentino et al. أجرى تحليلا مقارنا مشابها بناء على مجموعة عينة أكبر (ن = 10) ، مع استنتاجات مماثلة تم الحصولعليها 27. أظهرت هذه النتائج مجتمعة أن FEMBUA يمكن أن ينتج عنه الحد الأدنى من التناقضات في النتائج الميكانيكية الحيوية للعينات المتطابقة. أظهرت دراسة أخرى أيضا أن الطريقة المقترحة قابلة للتكرار وقوية ، حيث أظهرت الدراسة أن الاختلافات في الخواص الميكانيكية كان لها تأثير ضئيل على النتائج الميكانيكية الحيوية المحاكاة في نماذج الشريان التاجي من خلال تحليل العناصرالمحدودة 27.

نمذجة الافتراضات والقيود
هناك بعض الافتراضات المتضمنة في FEM لتحديد الخصائص المادية لجدران الشرايين التاجية في الجسم الحي ، والتي من شأنها أن تؤثر على نتائج FEMBUA كما هو موضح هنا. كان من المفترض أن يكون معدل الانكماش المحوري 95٪ حيث لا يمكن الحصول على الانكماش المحوري الفعلي في ظل الظروف في الجسم الحي . تم التحقيق في تأثير التمدد المحوري على خصائص المواد في دراسات سابقة ، وأشارت النتائج إلى أن التمدد المحوري الأصغر أدى إلى انكماش أكبر للشريحة وتقدير صلابة المواد الأكثر ليونة11،28،29،30. لذلك ، يجب استخدام بيانات التمدد المحوري الخاصة بالمريض عند توفرها. يتضمن تقييم جودة شبكة العناصر المحدودة عدة مؤشرات رئيسية ، من بينها شكل عنصر الشبكة هو معيار التقييم الرئيسي في هذه الورقة ، والذي يرتبط ارتباطا مباشرا باستقرار الحل العددي والتقارب التكراري. في التطبيقات العملية ، نظرا لتعقيد المشكلة ، من الضروري وزن المؤشرات لتلبية احتياجات التحليل ، وإذا لزم الأمر ، يمكن تحسين جودة الشبكة من خلال تصميم طرق تجزئة جديدة واختيار أنواع الشبكات المناسبة31. يجب إيلاء اهتمام خاص للدهون الصغيرة في المواقع الرئيسية أثناء التجزئة الفعلية32،33،34. في هذا التنفيذ ، يتم الاحتفاظ بالدهون الصغيرة في المواقع الرئيسية (الغطاء الليفي والكتف) ، بينما يتم الاحتفاظ بالدهون في المواقع الأخرى التي لها تأثير ضئيل على ظروف الإجهاد / الإجهاد. تم استخدام ضغط الدم المقاس بواسطة صفعة الذراع كبديل للضغط داخل الشريان التاجي في الموقع حيث لم يكن الضغط الغازي داخل الشريان التاجي متاحا للمريض. لم يتم النظر في الإجهاد النشط في الشريان التاجي لأن الأدلة المبكرة تشير إلى أن مساهمته في الخصائص المرنة للأوعية الدموية الحية كانت صغيرةجدا 35. لم تكن المعلومات عن الإجهاد المتبقي متاحة وبالتالي لم يتم تضمينها في هذا النموذج36،37. تم استخدام نماذج الهيكل فقط بدلا من نماذج تفاعل هيكل السوائل الأكثر تعقيدا في طريقة FEMBUA ، نظرا لأنها أكثر كفاءة من الناحية الحسابية حيث تحتاج الإجراءات التكرارية عادة إلى حل النماذج الحسابية عدة مرات للعثور على خصائص المواد ثابتة.

باختصار ، يمكن استخدام FEMBUA المستندة إلى الصور ، والتي تختلف عن الأساليب التجريبية الكلاسيكية الأخرى خارج الجسم الحي ، لتحديد الخصائص المادية الخاصة بالمريض للأوعية التاجية في الجسم الحي بشكل فعال. نظرا لأن تصلب الأنسجة الشريانية قد تم استخدامه بالفعل في البيئات السريرية كعامل خطر لأمراض القلب والأوعية الدموية ، فإن FEMBUA يجعل من الممكن المراقبة المستمرة للخصائص الميكانيكية للشريان التاجي في ظل ظروف في الجسم الحي ، وبالتالي يحمل الإمكانات في التطبيقات السريرية للعلاج الشخصي والطب الدقيق.

الإفصاحات

ويعلن أصحاب البلاغ عدم وجود تضارب في المصالح.

شكر وتقدير

تم دعم هذا البحث جزئيا من قبل مشروع علوم وتكنولوجيا الصحة الطبية في مقاطعة شاندونغ (رقم 202425020256 و 202403010254) ، ومنح المؤسسة الوطنية للعلوم الطبيعية في الصين 11972117 و 11802060 ، ومؤسسة العلوم الطبيعية لمقاطعة جيانغسو بموجب المنحة رقم BK20180352 ، ومؤسسة العلوم الطبيعية لمقاطعة شاندونغ بموجب المنحة رقم ZR2024QA110.

المواد

قائمة المواد المستخدمة في هذه المقالة
الاسمالشركةرقم فهرسيالتعليقات
  Bee DICOM Viewer SinoUnion للرعاية الصحيةالإصدار 3.5.1برنامج قارئ صور DICOM
ADINA Adina R & Dالإصدار 9.0حل العناصر المحدودة
ImageJ & nbsp ؛المعاهد الوطنية للصحةIVUS المجزأة MATLAB
MathWorks   الإصدار R2018aمنصة البرمجة التجارية
نظام التصوير Volcano s5Volcanoنظام التصوير بالموجات فوق الصوتية داخل الأوعية الدموية  
شركة شركة

المراجع

  1. Roth, G. A., et al. Global burden of cardiovascular diseases and risk factors, 1990-2019: update from the GBD 2019 study. J Am Coll Cardiol. 76 (25), 2982-3021 (2020).
  2. Kanwar, S. S., et al. Acute coronary syndromes without coronary plaque rupture. Nat Rev Cardiol. 13 (5), 257-265 (2016).
  3. Milzi, A., et al. Coronary plaque composition influences biomechanical stress and predicts plaque rupture in a morpho-mechanic OCT analysis. Elife. 10, e64020(2021).
  4. Stefanati, M., et al. Effect of variability of mechanical properties on the predictive capabilities of vulnerable coronary plaques. Comput Methods Programs Biomed. 254, 108271(2024).
  5. Laurent, S., et al. Expert consensus document on arterial stiffness: methodological issues and clinical applications. Eur Heart J. 27 (21), 2588-2605 (2006).
  6. Daisuke, K., et al. Coronary plaque phenotype associated with positive remodeling. J Cardiovasc Comput Tomogr. 18 (4), 401-407 (2024).
  7. Macrae, R. A., Miller, K., Doyle, B. J. Methods in mechanical testing of arterial tissue: A review. Strain. 52, 380-399 (2016).
  8. Hayashi, K. Experimental approaches on measuring the mechanical properties and constitutive laws of arterial walls. J Biomech Eng. 115 (4B), 481-488 (1993).
  9. Sacks, M. S., Sun, W. Multiaxial mechanical behavior of biological materials. Annu Rev Biomed Eng. 5, 251-284 (2003).
  10. Charis, C., et al. Impact of combined plaque structural stress and wall shear stress on coronary plaque progression, regression, and changes in composition. Eur Heart J. 40 (18), 1411-1422 (2019).
  11. Wang, L., et al. Quantifying patient-specific in vivo coronary plaque material properties for accurate stress/strain calculations: An IVUS-based multi-patient study. Front Physiol. 12, 721195(2021).
  12. Guo, X. Y., et al. Quantify patient-specific coronary material property and its impact on stress/strain calculations using in vivo IVUS data and 3D FSI models: a pilot study. Biomech Model Mechanobiol. 16 (1), 333-344 (2017).
  13. Wang, L., et al. Quantification of patient-specific coronary material properties and their correlations with plaque morphological characteristics: An in vivo IVUS study. Int J Cardiol. 371, 21-27 (2023).
  14. Lu, W., et al. Accurate identification of potential critical coronary lesions for the reduction of risk of cardiovascular events: study protocol for a randomized, open-label, active-controlled multi-center trial. Clin Trial Degenerat Dis. 3 (3), 106-110 (2018).
  15. Lv, R., et al. Using optical coherence tomography and intravascular ultrasound imaging to quantify coronary plaque cap stress/strain and progression: A follow-up study using 3D thin-layer models. Front Bioeng Biotechnol. 9, 713525(2021).
  16. Camasao, D. B., Mantovani, D. The mechanical characterization of blood vessels and their substitutes in the continuous quest for physiological-relevant performances. A critical review. Mater Today Bio. 10, 100106(2021).
  17. Holistic Numerical Methods Institute at University of South Florida. Secant Method. , Florida, FL, USA. (2003).
  18. Barnes, J. An algorithm for solving nonlinear equations based on the secant method. Comp J. 8 (8), 66-72 (1965).
  19. Liu, M., Liang, L., Sun, W. A new inverse method for estimation of in vivo mechanical properties of the aortic wall. J Mech Behav Biomed Mater. 72, 148-158 (2017).
  20. Yu, H., et al. Patient-specific in vivo right ventricle material parameter estimation for patients with tetralogy of Fallot using MRI-based models with different zero-load diastole and systole morphologies. Int J Cardiol. 276, 93-99 (2019).
  21. Meuwissen, M. H. H., Oomens, C. W. J., Baaijens, F. P. T., Petterson, R., Janssen, J. D. Determination of the elasto-plastic properties of aluminium using a mixed numerical-experimental method. J Mater Process Technol. 75, 204-211 (1998).
  22. Davis, F. M., Luo, Y., Avril, S., Duprey, A., Lu, J. Local mechanical properties of human ascending thoracic aneurysms. J Mech Behav Biomed Mater. 61, 235-249 (2016).
  23. Fan, L., Wang, H., Kassab, G. S., Lee, L. C. Review of cardiac-coronary interaction and insights from mathematical modeling. WIREs Mech Dis. 16 (3), e1642(2024).
  24. Andreas, W., et al. A finite element updating approach for identification of the anisotropic hyperelastic properties of normal and diseased aortic walls from 4D ultrasound strain imaging. J Mech Behav Biomed Mater. 58, 122-138 (2016).
  25. Liu, M., et al. Identification of in vivo nonlinear anisotropic mechanical properties of ascending thoracic aortic aneurysm from patient-specific CT scans. Sci Rep. 9, 12983(2019).
  26. Trabelsi, O., Duprey, A., Favre, J. P., Avril, S. Predictive models with patient specific material properties for the biomechanical behavior of ascending thoracic aneurysms. Ann Biomed Eng. 44, 84-98 (2016).
  27. Cosentino, F., et al. On the role of material properties in ascending thoracic aortic aneurysms. Comput Biol Med. 109, 70-78 (2019).
  28. Krauz, K., et al. The role of epicardial adipose tissue in acute coronary syndromes, post-infarct remodeling and cardiac regeneration. Int J Mol Sci. 25 (7), 3583(2024).
  29. McCracken, I. R., Smart, N. Control of coronary vascular cell fate in development and regeneration. Semin Cell Dev Biol. 155 (Pt C), 50-61 (2024).
  30. Holzapfel, G. A., et al. Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling. Am J Physiol Heart Circ Physiol. 289 (5), H2048-H2058 (2005).
  31. Han, Y., et al. Ventricle stress/strain comparisons between Tertalogy of Fallot patients and healthy using models with different zero-load diastole and systole morphologies. PLoS One. 14 (8), e0220328(2019).
  32. Ryo, K., et al. Role of the low-density lipoprotein-cholesterol/high-density lipoprotein-cholesterol ratio in predicting serial changes in the lipid component of coronary plaque. Circ J. 81 (10), 1439-1446 (2017).
  33. Thomas, T. W., et al. In vivo characterization and quantification of atherosclerotic carotid plaque components with multidetector computed tomography and histopathological correlation. Arterioscler Thromb Vasc Biol. 26 (10), 2366-2372 (2006).
  34. Masanori, K., et al. In vivo quantitative tissue characterization of human coronary arterial plaques by use of integrated backscatter intravascular ultrasound and comparison with angioscopic findings. Circulation. 105 (21), 2487-2492 (2002).
  35. Roach, M. R., Burton, A. C. The reason for the shape of the distensibility curves of arteries. Can J Biochem Physiol. 35, 681-690 (1957).
  36. Fung, Y. C., Liu, S. Q. Strain distribution in small blood vessel with zero-stress state taken into consideration. Am J Physiol. 262, 544-552 (1992).
  37. Ohayon, J., et al. Influence of residual stress/strain on the biomechanical stability of vulnerable coronary plaques: Potential impact for evaluating the risk of plaque rupture. Am J Physiol Heart Circ Physiol. 293, 1987-1996 (2007).

إعادة الطباعة والأذونات

الوسوم

VH IVUS Cine IVUS