Konuya özgü kas-iskelet modeli kemik zorlanma sırasında dinamik hareket eğitimi

Özet

İniş sırasında alt vücut kemikleri büyük mekanik yükler deneyim ve deforme. Etkileri ile ilişkili kemik stres yaralanma mekanizmaları daha iyi anlamak için kemik deformasyonu ölçmek için önemlidir. Konuya özgü kas-iskelet modelleme ve Sonlu elemanlar analizi entegre yeni bir yaklaşımın tibial zorlanma sırasında dinamik hareketleri ölçmek için kullanılır.

Özet

Kemik stres kırıklarının spor ve askeri eğitimleri yaygındır. Tekrarlayan büyük kara etkisi Kuvvetleri eğitim sırasında neden olabilir. Kemik stres yaralanma mekanizmaları daha iyi anlamak için alt-vücut kemik deformasyon üzerinde yüksek zemin etkisi etkisini zorlar belirlemek için önemlidir. Geleneksel ağırlık ölçme esneklik detektörler ölçüm içinde vivo tibia deformasyon çalışma kullanılmaya başlanmıştır. Bu yöntem invasiveness yordam, kaç insan denekler ve küçük kemik yüzey alanları sınırlı zorlanma verileri katılımı da dahil olmak üzere sınırlamalar ile ilişkilidir. Çalışmada tibia kemik baskı koşulları yükleme yüksek etki altında çalışmaya yeni bir yaklaşımın tanıtmak niyetinde. Bir konuya özgü kas-iskelet modeli sağlıklı erkek (19 yıl, 80 kg, 1800 mm) temsil etmek için oluşturuldu. Esnek Sonlu elemanlar tibia modeli bir bilgisayarlı tomografi (CT) inceden inceye gözden geçirmek konunun doğru Tibia üzerinde tabanlı oluşturuldu. Laboratuvar hareket yakalama damla-iniş kinematik ve zemin tepki kuvvetlerinin farklı yükseklikte (26, 39, 52 cm) elde etmek için gerçekleştirildi. Esnek tibia kalıcı bir analizi ile kombine multibody dinamik bilgisayar simülasyonları tibia zorlanma damla-iniş sırasında ölçmek için yapıldı. Hesaplanan tibia zorlanma veri önceki in vivo çalışmalar ile iyi anlaşma vardı. Non-invaziv bu yaklaşım daha iyi tibia stres kırıkları yaralanma mekanizmasının anlaşılması için yol açacak büyük bir kohort için yüksek darbe faaliyetleri sırasında tibia kemik suşu çalışmaya uygulanabilir belirgindir.

Giriş

Kemik stres, stres kırıkları gibi uzun süre kurtarma gerektiren ve önemli tıbbi maliyet^1,2üstlenmeden şiddetli aşırı yaralanmaları yaralar. Stres kırıkları ortak atletik ve askeri nüfus her ikisi de. Arasında tüm spor yaralanmaları, stres kırıkları hesap için toplam³yüzde 10'u ile ilgili. Özellikle, sporcular yüzü daha yüksek bir yaralanma oranı % 20⁴izlemek. Askerler de stres kırıkları yüksek orandaki deneyim. Örneğin, %6 yaralanma oranı için ABD Ordusu¹ bildirildi ve % 31 yaralanma oranı İsrail Ordusu⁵' te bildirildi. Tüm bildirilen stres kırıkları arasında tibia stres kırığı en yaygın bir^6,7,8' dir.

Spor ve fiziksel eğitimler tibia stres kırığı daha yüksek bir risk ile normalde yüksek zemin etkileri (örneğin, atlama, açılış ve kesme) ile ilişkilidir. Ayak yere kurduğunda hareket sırasında vücut için bir yere darbe kuvveti uygulanır. Bu darbe kuvveti kas-iskelet sistemi ve ayakkabı tarafından harcanmış. İskelet sistemi kolları kas Kuvvetleri zemin etkisi⁹emmek için uygulamak izin veren bir dizi olarak hizmet vermektedir. Bacak kasları yeterli zemin etkisi azaltmak olamaz zaman, alt-vücut kemikleri kalan gücü absorbe gerekir. Kemik yapısı bu işlemi sırasında deformasyon yaşayacaksınız. Artık darbe kuvveti tekrarlayan emilimini birikir ve stres kırıkları haline kemik microdamages içinde neden olabilir. Bugüne kadar kemik için ilgili bilgiler dış kara etkisi Kuvvetleri tepki sınırlıdır. Tibia kemik dinamik hareketleri sırasında yüksek etki güçleri tarafından tanıtılan mekanik yük nasıl yanıt vereceğini çalışma önemlidir. Tibia kemik deformasyonu yüksek darbe etkinlikleri sırasında inceleyerek tibia stres kırığı mekanizmasının daha iyi anlamak için neden olabilir.

Kemik deformasyonu vivo ölçmek için kullanılan geleneksel teknikleri Araçlı gerilim ölçerler^{10,11,12,13,14,15}tarihinde güveniyor. Cerrahi işlemler kemik yüzeyindeki gerginlik ölçme implant için ihtiyaç vardır. İnvaziv doğası nedeniyle, in vivo çalışmalar gönüllü küçük bir örnek ile sınırlıdır. Buna ek olarak, ağırlık ölçme esneklik detektörler sadece kemik yüzeyine küçük bir bölgenin izleyebilirsiniz. Son zamanlarda, non-invaziv bir yöntem kemik zorlanma çözümlemek için bilgisayar simülasyonu kullanan tanıtılan^16,17yaşındaydım. Bu metodoloji insan hareketi sırasında kemik suşu çalışmaya kas-iskelet modelleme ve bilgisayar simülasyonları birleştirme olanağını sağlar.

Bir kas-iskelet modeli bir iskelet ve iskelet kasları tarafından temsil edilir. İskelet rijit veya deforme olmayan organları olan kemik parçadan oluşur. İskelet kasları ilerici-integral-türev (PID) algoritması kullanılarak denetleyicileri olarak modellenmiştir. Üç dönemlik PID kontrol hataları tahmininde çıkış doğruluğu¹⁸geliştirmek için kullanır. Özünde, PID denetleyicileri kasları temsil eden kasların uzunluğu değişiklikler zamanla üretmek için gerekli kuvvetler geliştirerek vücut hareketleri oluşturmaya çalışın. PID denetleyicisi hata uzunluğu/saat tonundaki hareketi çoğaltılması için kuvvet değiştirmek için kullanır. Bu simülasyon işlemi tüm kas iskelet hareket ve vücut hareketleri üretmek amacıyla birlikte çalışmak için koordine etmek için uygun bir çözüm oluşturur.

İskelet kas-iskelet modeli bir veya daha fazla segmentlerinde deformasyon ölçümü izin vermek için esnek organları modellenebilir. Örneğin, tibia kemik öğeleri ve düğümleri binlerce oluşan bir sonlu sayıda öğe, bölünebilir. Esnek tibia mekanik yükleme etkisi Sonlu elemanlar (FE) analizi ile incelenebilir. FE analiz tek tek öğeleri yükleme tepki zaman içinde hesaplar. Kemik öğeleri ve düğümleri artış sayısı olarak, FE analizin hesaplama saat önemli ölçüde artacaktır.

Esnek organları deformasyon doğru değerlendirilmesi ile Hesaplamalı maliyetini azaltmak için kalıcı FE analiz geliştirilmiş ve içinde otomotiv ve Havacılık ve uzay sanayi^19,20kullanılır. Bireysel FE elements yanıt zaman etki alanındaki mekanik yük için analiz yerine, bu yordamı farklı titreşim frekansları frekans etki alanındaki temel nesnenin mekanik yanıt değerlendiriyor. Bu yöntem deformasyon²⁰doğru ölçüm sağlarken hesaplama zaman önemli bir azalma olur. Her ne kadar kalıcı FE analizi yaygın mekanik yorgunluk otomotiv ve Havacılık alanlarında eğitim için kullanılmıştır, bu yöntemin uygulama çok sınırlı olmuştur insan hareketi bilim. Al Nazer vd., insan yürüyüş sırasında tibial deformasyon incelemek için kalıcı bir FE analiz kullanılan ve sonuçları^16,17teşvik bildirdi. Ancak, bu yöntem büyük ölçüde bilgisayar simülasyonları sürücü için sadece bir deney sınırlı kinematik verileri kullanarak etkilendi; Gerçek olduğunu zemin simülasyonlar yardımcı olmak için kullanılan etki güçleri. Bu yaklaşım yürüyüş gibi düşük etkisi yavaş hareketleri eğitimi için makul olabilir ama bu yüksek zemin etkisi hareketleri çalışmaya uygun bir çözüm değil. Böylece, dinamik yüksek darbe faaliyetleri sırasında alt vücut kemik reaksiyonları incelemek için daha önce bildirilen yöntemiyle ilişkili sınırlamaları gidermek için yenilikçi bir yaklaşım geliştirmek esastır. Özellikle, bir yöntem kullanmak doğru deneysel kinematik veri ve gerçek kara etkisi Kuvvetleri geliştirilen gerekir. Bu nedenle, bu çalışmanın amacı multibody dinamik simülasyonlar tibial zorlanma yüksek darbe faaliyetleri sırasında incelemek için kalıcı FE analizi ile gerçekleştirmek için bir konuya özgü kas-iskelet modeli geliştirmekti. Damla-iniş tarafından farklı yükseklikte temsil bir dinamik yüksek darbe hareketi yöntemini sınamanızı sağlamak üzere seçildi.

Protokol

Helsinki Bildirgesi altında yapılan. Veri toplama önce konu incelendi ve çalışmaya katılan önce üniversite kurumsal inceleme Kurulu tarafından onaylanan izin formu imzaladı.

1. CT görüntüleme Protokolü

1. Katılımcı bir tesise nerede CT tarayıcı yer alır al. CT taraması önce CT makine aşağıdaki parametrelerle yapılandırmak: CT dilim kalınlığı 0.625 mm, CT görüş alanı 15 cm x 15 cm ve auto için en yüksek kilo-gerilim (kVp) ve miliamper-saniye (mAs) makine algoritması kullanarak parametrelerini ayarlama.
2. CT tarayıcı ringde içine slaytlar bir tablo üzerinde yalan katılımcı sormak. Çok hala CT tarama sırasında kalmak için Katılımcı sormak. Femur distal sonuna calcaneus ayrı ayrı her bacak inceden inceye gözden geçirmek.
3. CT taraması tamamlandığında, CT görüntüleri dijital görüntüleme ve iletişim Tıp (DICOM) formatında dışa aktarın. Bir görüntü boyutu 512 x 512 piksel (gri tonlama) seçin.
   Not: CT görüntüleme Protokolü normalde az 1 h sürer. Radyasyon dozu en az düzeydedir. Bundan normal X-ray tıbbi işlemler sırasında karşılaşılan büyük risk sunar.

2. antropometrik ölçüm Protokolü

1. Hareket yakalama önce laboratuvar ziyaret sırasında katılımcının vücut kitle (kg), vücut yüksekliği (mm), anterior superior İliak spines (ASISs) (mm), bacak uzunluğu (mm), diz eklem genişliği (mm) ve ayak bileği eklem genişliği (mm) arasındaki uzaklığı ölçmek.
2. ASIS mesafe ölçümleri Inter: bir kumpas sol ASIS ve doğru ASIS arasındaki doğrusal mesafe ölçmek için kullanın.
3. Bacak uzunluk ölçümü: bir kaset ölçü birimi ölçü doğrusal mesafe ASIS ve bileşenleri kullanmak için her iki bacak.
4. Diz eklem genişliği ölçümü: bir kumpas her iki diz için femur lateral ve medial epicondyles arasındaki doğrusal mesafe ölçmek için kullanın.
5. Ayak bileği eklem genişliği ölçümü: bir kumpas lateral ve medial malleoli için her iki bacak arasında doğrusal mesafe ölçmek için kullanın.
   Not: Inter ASIS mesafe, bacak uzunluk, diz ve ayak bileği genişliği biyomekanik yazılımındaki bir konu model oluşturmak için kullanılır (bkz. Tablo reçetesi) kinematik ve Kinetik hesaplamaları gerçekleştirmek için.

3. hareket yakalama iletişim kuralı

Not: Tüm yazılım ve kullanılan araçlar için Malzemeler tablo bakın.

1. Yansıtıcı işaretleri yerleşimini
   1. 14-mm yansıtıcı işaretleri aşağıdaki anatomik kemikli yerler katılımcının vücuduna yerleştirin: acromion süreçleri, sternoclavicular eklem, göğüs kemiğinin tabanı, 10^th torasik vertebra, ASISs, posterior-Superior İliak posterior işlemi dikenleri (PSISs), yanal diz eklem hatları, medial diz eklem hatları üzerinde 1,5 cm yukarıda 1,5 cm malleoli, medial malleoli, posterior topuklar, ikinci ayaktarakları üsleri ve beşinci ayaktarakları üsleri lateral.
   2. 4-marker kümeleri ile yarı sert plastik tabaklar uyluk ve shanks, sırasıyla yerleştirin.
      Not: bir daha iyi hareket yakalama sonuç için Katılımcı yalınayak ve deri giyim giymek tavsiye edilir. Buna ek olarak, işaretleyici yerleştirme yordam değiştirilmiş "Plug-in-yürüyüş" protokolü²¹izler. Toplam 39 yansıtıcı işaretleri hareket yakalama için kullanılır ve 34 tanesi alt beden bağlıdırlar.
2. 5 min için otomatik olarak seçilen bir hızda bir motorlu koşu bandı, yürüme ısınmak için Katılımcı talimat.
3. Kalibrasyon için hareket yakalama yordamı Oda alanı
   1. Hareket yakalama sistemi (12 yüksek hızlı kızılötesi kameralar) ve iki tabak zorlamak. Bir hareket yakalama yazılım programı açın. Ana program penceresinde, 'Kaynaklar' bölmesini açın. "Sistem" sekmesini Yapılandır 200 Hz kamera frekansta ve 2.000 Hz frekansta plaka zorlamak.
   2. Ana program penceresinde, 'Araçlar' bölmesini açın. "Sistem Hazırlama" düğmesini tıklatın. "Kameralar kalibre"'ı tıklatın. "Başlat". Dalga Dinamik Kalibrasyon damla-Açılış hareketleri gerçekleştirilecek nerede Oda alanı içinde gerçekleştirmek için bir standart 5-marker kalibrasyon değnek için bir araştırma personeli sor. 5 sonra "Dur"'ı tıklatın s değnek veri satın aldı.
   3. Kalibrasyon değnek bir köşesinde (origin) kalibre edilmiş alanı için bir başvuru konumu belirtme amacıyla bir kuvvet plaka ile hizalamak için düz zemine yerleştirin. "Ayarla birim Origin" içinde 'Sistem Hazırlama' araçları Bölmesi'ni tıklatın.
4. Hareket yakalama yazılım programı katılımcı hazırlık
   1. Ana program penceresinde, 'Kaynaklar' bölmesini açın. "Yeni konu bir etiketleme iskelet oluşturmak" düğmesini tıklayın "Konu" sekmesini tıklatın. Etiketleme şablon şablon dosyalarını sağlanan listeden seçin.
   2. 'Özellikler' penceresinde ilgilinin adı ve vücut kitle (kg), vücut yüksekliği (mm), Inter-ASIS mesafe (mm), sol ve sağ ayak Uzunluk (mm), sol ve sağ diz genişliği (mm) ve sol ve sağ ayak bileği genişliği (mm) değerleri girin. 'Konu kaynaklar' bölmesinde özne adını sağ tıklatın ve "Konu Kaydet'i" tıklatın.
5. Bir statik vücut kalibrasyon poz kaydetmek
   1. Ayaklarınızı omuz genişliği üst ekstremite yanal böylece vücuttaki tüm yansıtıcı işaretleri de kameralar için maruz genişletme sırasında apart ile kalibre edilmiş odanın ortasında hareketsiz durmak için Katılımcı sormak.
   2. Ana program penceresinde, araç bölmesini açın. "Konu hazırlık" sekmesini tıklatın. Konu yakalama bölümünde, statik kalibrasyon deneme olmak kayıt 3 sn hareket deneme "Başlat"'ı tıklatın.
6. Fonksiyonel ortak merkezleri belirleme yordamı
   1. Fonksiyonel Kalça eklemi Merkezi
      1. Bir ayak ile katılımcı sormak ve tamamen diğer bacağını hafifçe öne uzatın. Katılımcı aşağıdaki sırayla genişletilmiş bacak ve Kalça eklemi etrafında taşımak için talimat: anteriorly taşımak ve için tarafsız dönmek, anterior yanal hareket ve için tarafsız dönmek, yanal hareket ve için tarafsız dönmek, arka yanal hareket ve dönmek bağımsız, özafagusu taşımak ve tarafsız ve sirkumdiksiyon hareketi geri dönmek.
      2. Ana program penceresi içinde belgili tanımlık 'Alet' cam açın, "Yakalama" sekmesini tıklatın. Yakalama bölümünde, "her işlevsel kalça hareket için hareket deneme kaydetmek" Başlat.
   2. Fonksiyonel diz eklem Merkezi
      1. Katılımcı ile bir ayak ve bir 30 ° kalça hiper-eklenti diğer bacak korumak için sor. Katılımcı için 5 kez 45 ° diz fleksiyon sigara ağırlık taşıyan bacak ile gerçekleştirmek için talimat.
      2. Belgili tanımlık 'Alet' cam 'Yakalama' bölümünde, "her işlevsel diz hareket için hareket deneme kaydetmek" Başlat.
         Not: fonksiyonel ortak yordamı ilgili ayrıntılı bilgi için lütfen bkz: Schwarz, vd. ²²
7. Damla-açılış hareketlerinin hareket yakalama
   1. Üç farklı damla-açılış yükseklikte (26 cm, 39 cm ve 52 cm)¹⁴kullanarak sipariş rastgele.
   2. Yer yüksekliği 50 x 50 cm² yerdeki paspas tarafından kaplı tahta kutu üst yüzey alanı ile ayarlanabilir. Tahta kutu kuvvet plakaların kenarlarından 11 cm'dir. Kutusunu yüzey üzerinde durmak için Katılımcı sormak.
   3. Katılımcının doğrudan kutusunun önünde baskın onların ayak uzatmak ve onların ağırlığı ileri ve kutusundan geri çekil talimat. Her iki bacak ile ayrı kuvvet tabakta her ayak ile aynı anda yere iniş için Katılımcı sormak.
   4. Deneme hareket yakalama tamamlanana kadar ayakta kalmayı katılımcı sormak. Üç hareket denemeler her yükseklik için toplamak için hareket yakalama üç kez tekrarlayın.
8. Hareket yakalama veri işleme
   1. Bir hareket yakalama yazılım programı açın. Programın ana penceresi içinde 'İletişim' bölmesine git. Kaydedilen hareket denemeler birini seçin "Veri yönetimi" sekmesini tıklatın ve programda açın.
   2. 'Araçlar' Bölmesi'nde, "Boru hattı" sekmesini tıklatın. 'Geçerli boru hattı', "Yeniden" boru hattı seçin. Üç-boyut (3D) yörüngeleri yansıtıcı işaretleyicilerin görünümünü elde etmek için yeniden yapılandırma işlemini başlatmak için "Çalıştır" düğmesini tıklatın.
   3. 'Araçlar' Bölmesi'nde, "Etiket/düzen" sekmesini tıklatın. 'Manuel etiketleme' bölümünde, bireysel marker isimleri seçin ve ilgili 3D yörüngeleri etiket. "Düğme ne zaman etiketleme araç çubuğunun tamamlandı Kaydet" i tıklayın.
   4. 'Araçlar' Bölmesi'nde, "Boru hattı" sekmesini tıklatın. 'Kullanılabilir işlemleri' bölümünde, "Dosyayı Aktar" seçin. "İhracat C3D boru hattı" çift tıklatın. İşlenmiş hareket deneme bir koordinat üç-boyut (C3D) biçiminde bir dosyaya vermek için "Çalıştır" düğmesini tıklatın.
9. Hareket yakalama veri biyomekanik analiz
   1. Daha fazla işlem hareket yakalama veri için biyomekanik yazılım programını açın. Üst menüden "Dosya" tıklayın ve "Aç/Ekle" düğmesini tıklatın. Biyomekanik yazılım programı almak için ham C3D dosyaları seçin.
   2. Üst menüden "Model"'i tıklatın. Tıklatın "oluşturmak (statik kalibrasyon dosyası eklemek)". Alt menüde "Hibrid modeli üzerinden C3DFile". Seçin ve statik kalibrasyon C3D dosyasını açın.
   3. Üst menüden "Model"'i tıklatın. Aşağı açılan listeden "Modeli şablonunu Uygula"'ı tıklatın. Seçin ve bir modeli şablon dosyası açın. "Modelleri" sekmesini araç çubuğunda'i tıklatın. "Konu veri / ölçümleri" sekmesini tıklatın. 'Konu veri' penceresi içinde 'Kitle' ve 'Yükseklik' modele konuya özgü değerleri değiştirin.
   4. "Modelleri" sekmesini araç çubuğunda'i tıklatın. Üst menü çubuğunun "Modeli Oluşturucu gelişmiş posta işleme" düğmesini tıklatın. "Modeli Oluşturucu gelişmiş posta işleme" açılan pencerede "Fonksiyonel eklem" sekmesini seçin "hareket dosya üzerinden çalışma alanı ekle"'ı tıklatın.
   5. Fonksiyonel ortak Merkezi C3D dosyaları seçin. İçe aktarılmış bir işlevsel ortak dosyanın vurgulayın. Dosya eşleştirme fonksiyonel esrar vurgulayın. "Ayarla başlatmak için geçerli kare" ve "Son kare için geçerli çerçeveyi ayarla" deneme hareket uygun bölümlerini seçmek için kullanın. "Kontrol yer hesaplamak" düğmesini tıklatın. İskelet modeli geliştirmek için diğer fonksiyonel ortak merkezleri hesaplamak için bu işlemi yineleyin.
   6. Üst menü çubuğundaki "Model" düğmesini tıklatın. "Modeli hareket dosyaları için ata"'yı seçin. "Atama modelleri için hareket verileri" açılan pencerede konuya özgü iskelet modeli tüm hareket denemeler için geçerlidir.
   7. Araç çubuğundaki "Boru hattı" düğmesini tıklatın. "Boru hattı Atölyesi" açılan pencerede "Açık boru hattı" düğmesini tıklatın. "Filtreleme hedefleri boru hattı" seçin. Kesme 10 Hz frekansta hareket yakalama denemelerin 3D yörüngeleri ile dördüncü dereceden alçak geçiren Butterworth filtre gerçekleştirmek için "Boru hattı yürütmek" düğmesini tıklatın.
   8. Araç çubuğundaki "Boru hattı" düğmesini tıklatın. "Boru hattı Atölyesi" açılan pencerede "Açık boru hattı" düğmesini tıklatın. "Filtreleme Kuvvetleri boru hattı" seçin. Kesme 60 Hz frekansta zemin mukabele Kuvvetleri hareket yakalama denemeler ile dördüncü dereceden alçak geçiren Butterworth filtre gerçekleştirmek için "Boru hattı yürütmek" düğmesini tıklatın.
   9. Üst menü çubuğunun "Ayarlar" düğmesini tıklatın. "Kullanım işleme analogları için zemin tepki kuvvet hesaplamaları" ve "İşlenen hedefleri için Model/Segment/LinkModelBased kullanım" yanındaki onay işareti koyun.
   10. Araç çubuğundaki "Boru hattı" düğmesini tıklatın. "Boru hattı Atölyesi" açılan pencerede "Açık boru hattı" düğmesini tıklatın. "Model alan hesaplama" boru hattı seçin. Alt vücut ortak kinematik ve Kinetik hesaplamaları gerçekleştirmek için "Boru hattı yürütmek" düğmesini tıklatın.
   11. Araç çubuğundaki "Boru hattı" düğmesini tıklatın. "Boru hattı Atölyesi" açılan pencerede "Açık boru hattı" düğmesini tıklatın. "İhracat C3D koordinatları" boru hattı seçin. Alt vücut görsel işaretleri C3D dosyasına işlenmiş 3D koordinatlarını vermek için "Boru hattı yürütmek" düğmesini tıklatın.
   12. Araç çubuğundaki "Boru hattı" düğmesini tıklatın. "Boru hattı Atölyesi" açılan pencerede "Açık boru hattı" düğmesini tıklatın. "İhracat zemin mukabele Kuvvetleri" boru hattı seçin. İşlenmiş 3D yere tepki vermek için "Boru hattı yürütmek" düğmesini zorlar bir ikili dosyada tıklatın (dosya uzantısı: MAT).
       Not:, ham toprak reaksiyonu kuvvet veri²³süzmek için bir kesme frekans 60 Hz yüksek darbe doruklarına iniş sırasında korumak için kullanılır.
10. Bilgisayar simülasyonları için hareket yakalama veri hazırlama
    1. Bir bilgisayar yazılım programlama açın. Filtre uygulanmış C3D veri dosyası ve MAT veri dosyasını içe aktarın.
    2. Alt vücut ortak Merkezi koordinatları içeren bir metin dosyasına verme. C3D veri dosyası ve MAT veri dosyasının metin dosyalarını dönüştürme (dosya uzantısı: slf) multibody dinamik simülasyon programı tarafından kullanılacak.

4. Konu özel yordamın modelleme

1. Alt vücut iskelet modeli oluşturma
   1. Dinamik simülasyon yazılım programı eklenti modelleme insan vücudu ile yüklü multibody açın. Bu işlem sırasında eklenti modülü modelleme insan vücudu otomatik olarak açılır. Giriş ekranı içinde Denetim Masası bina modeli açmak için "Yeni Model" simgesini çift tıklatın.
   2. Ana modelleme paneli içinde "Antropometrik veritabanı kitaplığı" bölümünde, genel vücut (GeBOD) aþaðý açýlan listeden seçin. Ana modelleme paneli içinde belirtin vücut kitle (kg), vücut yüksekliği (mm), cinsiyet ve yaş (ay).
   3. "Vücut yapılandırma" bölümünde, ana modelleme paneli içindeki "Alt beden" radyo düğmesini tıklatın. "Milimetre Kilogram Newton" "Birim" aþaðý açýlan listeden seçin. Ana modelleme paneli içinde Vücut ölçüleri kabul etmek için "Beden ölçüm tablo oluşturma" bölümünde "Uygula" düğmesini tıklatın. Bir alt-vücut iskelet temel model oluşturmak için "İnsan segmentleri oluşturmak" bölümündeki "Uygula" düğmesini tıklatmaya devam edin.
      Not: Bu model göre ölçeklendirilir bireyin yüksekliğine, kitle, yaş ve cinsiyet. Model yedi parçadan oluşur: bir leğen kemiği, iki uyluk, iki shanks ve ayakları (şekil 1). Tüm parçaların katı cisimlerin modellenmiştir.
2. Alt vücut eklem modelleme
   1. Ana modelleme paneli içinde "eklemler" ortak yapılandırma panelini açmak için ana menü aşağı açılan listeden seçin.
   2. "Ortak döndürme öğeleri" bölümünde, ortak yapılandırma paneli içinde belgili tanımlık "hazırlamak modeli ile kayıt ek" yanındaki düğmeyi tıklatın. "Bahar tamponlar ve ortak sınırları özellikleri" bölümünde, aşağıdaki parametreleri girin: 1 Nmm / ° Nominal eklem sertliği, Nominal ortak Damping 0.1 Nmm∙s / °, 3.38E7 durdurmak eklem sertliği Nmm / °. Radyo düğmeleri adlarının yanındaki kontrol ederek "Sol bacak" ve "Sağ bacak" seçmeye devam. Ortak yapılandırmalar kabul etmek için "Apply" tuşuna tıklayın.
   3. Ana modelleme paneli içinde "İş akışı" ana menü aşağı açılan listeden seçin. "Yürüyüş" ve "Ayarla" alt menü aşağı açılan listeden seçin. "Ortak merkezi veri" bölümünde katılımcının alt vücut ortak Merkezi dosya girin.
   4. Ortak merkezleri konumlarını değiştirmek için verileri almak için "Yükle" düğmesini tıklatın. "Yük statik deneme" bölümünde, statik kalibrasyon hareket yakalama deneme girin (slf dosya biçiminde üretimi açıklanan adımları 3.8-3.10). Alt vücut iskelet modeli parametreleştirmek için dosyasını içe aktarmak için "Yükle" düğmesini tıklatın.
      Not: varsayılan olarak, kalça eklemleri özgürlük üç derece ile küresel eklem olarak yapılandırılmış, diz eklemleri revolute eklem özgürlük bir derece ile yapılandırılır ve ayak bileği eklem evrensel eklem ile özgürlük iki derece olarak yapılandırılır.
3. İskelet kasları modelleme
   1. Ana modelleme paneli içinde "Yumuşak dokular" ana menü aşağı açılan listeden seçin. "Create Bankası doku küme" alt menü aşağı açılan listeden seçin. "Kas CONTRACTILE öğeleri" bölümünde, "Kayıt kas öğeler ile modeli hazırlamak"'ı tıklatın.
   2. "Küresel kayıt öğesi kas özellikleri" bölümünde "Updated 45 kas grubu" radyo düğmesini tıklatın.
   3. "Küresel kayıt öğesi kas özellikleri" bölümünde, kas özellikleri için aşağıdaki varsayılan ayarları kabul etmek: 0.4448 N/mm pasif sertliği, pasif Damping 1.75 E-2 Ns/mm, kas istirahat sürü 0.4448 i. onay "Sol bacak" radyo düğmeleri ve "Bacak kas atamaları için sağ". Yapılandırmaları kabul etmek için "Apply" tuşuna tıklayın.
      Not: 45 bacak kas kümesi aşağıdaki kas içerir: Adductor Brevis, Adductor Longus, Adductor Magnus (üç grup), pazı kemiği uzun kafa, pazı kemiği kısa kafa, ekstansiyon Digitorum, ekstansiyon Hallucis, fleksör Digitorum, fleksör Hallucis, Gastroknemius, Gemellus, Gluteus Maximus (üç grup), Gluteus Medias (üç grup), Gluteus Minimis (üç grup), Gracilis, Hamstring, Iliacus, Lateral gastroknemius, Medial gastroknemius, Pectineus, Peroneus Brevis, Peroneus Longus, Peroneus Tertius, Piriformis, Psoas, kuadriseps kemiği, rektus kemiği, Sartorius, Semimembranosus, Semitendinosus, Soleus, Tensor silmeli arşitrav Latae, Tibialis Anterior, Tibialis Posterior, Vastus Intermedius, Vastus kasılmasına, Vastus Medialis.

5. çok vücut dinamiği simülasyonları

1. Ters kinematik simülasyon gerçekleştirme
   1. Ana modelleme paneli içinde "İş akışı" ana menü aşağı açılan listeden seçin. "Yürüyüş" ve "Deneme" alt menü aşağı açılan listeden seçin. "Dinamik deneme veri" bölümünde, dinamik hareket yakalama deneme (slf dosya biçiminde) dosyanın adını girin ve verileri almak için "Yükle" düğmesini tıklatın. Veri (slf dosya biçiminde) dosya ve verileri almak için "Yükle" düğmesini tıklatın ilgili yere müdahale gücü girmek devam edin.
   2. Ana modelleme paneli içinde "_Analyze" ana menü aşağı açılan listeden seçin. Başlangıçta dinamik davanın duruş eşleşecek şekilde modeli duruş ayarlamak için Reparameterize analizi yap.
   3. Simülasyon Masası'nı açın. Yerçekimi ve zemin mukabele Kuvvetleri etkilerini devre dışı bırakın. Tüm hareket deneme uzunluğu simülasyon olarak seçin.
   4. 100 adım/s. hareket yakalama veri tarafından yönlendirilen bir ters kinematik benzetim çalıştırın bir simülasyon zaman adım belirtin. Ters kinematik simülasyon çözümlemeyi kaydedin.
2. Bir hareket İzleyicisi Aracısı'nı oluşturma
   1. Hareket izci Aracısı oluşturma panelini açın. Varsayılan izci adı kabul edin: MA_Track.
   2. Translasyonel sertlik ve dönme sertlik 10 N/mm ve 1000 Nmm / °, anılan sıraya göre ayarlayın. Translasyonel Damping ve dönme Damping 10 Ns/mm ve 1.000 Nmms / °, anılan sıraya göre ayarlayın. Tüm translasyonel hem de rotasyonel serbestlik Driven ayarlayın.
   3. Dikkat edin. Yalnızca alt gövdeli modeli ileri dinamik simülasyon için kullanıldığı gibi bir hareket izleyici'yi hesabınıza istikrarsızlık üst vücut hareket eksikliği nedeniyle gereklidir.
3. Eğitim bacak kasları
   1. Yumuşak doku yapılandırma Masası'nı açın. Kapalı-döngü basit kas modeli seçin. Kas modeli için aşağıdaki parametreleri ayarlayın: orantılı kazanç, 1.0E6, ayrılmaz kazanç, 1.0E6 ve türev kazanç, 1.0E4.
   2. Ters kinematik simülasyon analiz kas eğitim hedefi olması için seçin. Kas eğitimi uygulanır.
4. Esnek bir tibia içe aktarma
   1. Esnek vücut alma panelini açın. Hizalama eşleme üç bilinen üreticileri ile ve onların karşılık gelen düğümleri esnek tibia yüzeyinde gerçekleştirin.
   2. Katı tibia esnek tibia tarafından değiştirilmesi için seçin. Esnek tibia temsil eden MNF dosyayı seçin. Bacak kaslarını esnek Tibia reattaching için kas eki eşleme dosyasını seçin. Esnek tibia için kas-iskelet modeli alın.
5. İleri dinamik simülasyon esnek tibia ile yerinde gerçekleştirme
   1. Simülasyon Masası'nı açın. Yerçekimi ve zemin mukabele Kuvvetleri etkileri etkinleştirin. Hareket ajanların etkileri devre dışı bırakın.
   2. Tüm hareket uzunluğu için simülasyon deneme çalıştırmak seçin. Simülasyon zaman adım 100 adım/s. Çalıştır / ileri bir dinamik simülasyon eğitimli kasları tarafından tahrik ayarlayın. İleri dinamik çözümlemeyi kaydedin.

6. bir esnek Tibia modeli oluşturma

1. Bir 3B yüzey kafes model oluşturma
   1. Bir görüntü işleme programı açın. DICOM formatında alma CT dilimleri. Bölge büyüyen yöntemi çevre yumuşak dokulardan kemik dokusu ayırarak bir maskesi oluşturun.
   2. Tibia ve fibula nereye bağlı CT dilimler için arama. Tibia ve fibula maskeyi boyunca iki kemik birlikte silerek ayırın.
   3. Sadece tibia kemik eklemek için yöntem büyüyen bölge kullanarak ikinci bir maskesi oluşturun. Tibia maskede mevcut boşluklar ortaya çıkarmak için CT dilimler geçmesi. Maskenin içindeki boşlukları doldurun. Tibia maskesine göre bir 3D tibia nesnesi oluşturun. 3D tibia nesneyi çizim Değişim Biçimi (DXF) dosyası olarak verin.
2. Bir Sonlu elemanlar tibia modeli oluşturma
   1. FE analiz yazılım programını açın. DXF uzantılı 3D tibia modeli dosyasını alın.
   2. Yinelenen öğeleri ve düğüm kaldırmak için süpürme komutunu gerçekleştirir. Tüm öğeler için aşağıdaki malzeme özellikleri 3 mm x 3 mm x 3 mm. ata altıgen unsurları ile FE tibia model oluşturmak için ses Mesh komutunu gerçekleştirir: Young katsayısı 17 not ortalaması, 1.9E yoğunluğu ve Poisson oranı 0,3-6 Kg/cm³.
      Not: Malzeme özellikleri her öğesine kemik dokusu izotropik baskı sırasında dinamik hareketleri^24,25,26kemik tarafından deneyimli aralıklarda varsayımı ile atanır.
3. Esnek tibia modeli oluşturma
   1. Ana kontrol paneli içinde "Geometri & kafes" sekmesini seçin "geometri & kafes"'ı tıklatın. "& Geometri Mesh" açılır penceresinde, "Mesh" bölümünde, "düğümleri Ekle" yi tıklayın diz ve ayak bileği eklemleri merkezlerinden temsil etmek için iki yeni düğüm oluşturmak için.
   2. Ana Denetim Masası'ndaki "Bağlantılar" sekmesini seçin RBE2'ın ' ı tıklatın. RBE2'ın açılan pencerede, diz ve ayak bileği yüzeylerde öğesinin türü 2 katı vücut (RBE2) ortak düğümleri ve yüzey düğümler arasındaki bağlantılar oluşturun.
   3. Ana Denetim Masası'ndaki "Sınır koşullar" sekmesini tıklatın. "Sınır koşulları" bölümünde bulunan "Yeni" düğmesini tıklatın. "DOF_Set düğümleri" seçin. "Sınır koşulu özellikleri" açılan pencerede özgürlük altı derece her iki RBE2 ortak düğüm atayarak bir sınır koşulu oluşturun.
   4. Ana Denetim Masası'nda "Loadcases" sekmesini tıklatın. "Loadcases" bölümünde, "Yeni"'ı tıklatın, "Adams Craig-Bampton"¹⁹seçin. "Loadcase özellikleri" açılan pencerede "DOF-küme düğümleri"'ı tıklatın. Yukarıdaki adımda oluşturduğunuz dofset_nodes seçin.
   5. Ana Denetim Masası'ndaki "İşleri" sekmesini tıklatın. "İşleri" bölümünde, "Yeni" seçeneğini tıklatın. "Yapısal" seçin. "İş özellikleri" açılır penceresinde, önceki adımda oluşturduğunuz loadcase seçin. "İş sonuçları" düğmesini tıklatın. "Sonuçlar" açılan pencerede, "Stres" ve "Soy" seçin. Ayrıca uzunluk ve zaman için "ikinci" için "Kilogram" kitle için güç için "Newton", "Mm" seçin. "Çalıştır" düğmesini tıklatın.
   6. "İşin yürütülmesini" açılan pencerede FE simülasyon için işini göndermek için ve tibia¹⁶kalıcı tarafsız dosyası (MNF) oluşturmak için "Gönder" düğmesini tıklatın.

7. zorlanma veri analizleri

1. Kemik zorlanma veri verme
   1. Yazı işlemci multibody simülasyon programı açın. Dayanıklılık eklenti programını yükleyin.
   2. Simülasyon ile esnek tibia simülasyon adını tıklatarak açın. Asıl suşları maksimum ve minimum ve maksimum kesme zorlanma tibial diaphysis antero-medial yönünü temsil eden düğüm dışa aktarın.
2. Ham baskı veri işleme
   1. Bir bilgisayar veri işleme için yazılım programlama açın. Ham baskı veri. 15 Hz kesme frekans ile ham veri için dördüncü sırada alçak-geçiren Butterworth filtre uygulayın.

Sonuçlar

Çalışma için gönüllü oldum bir sağlıklı kafkas erkek (19 yaş, yükseklik 1.800 mm, kitle 80 kg). Veri toplama önce konu incelendi ve çalışmaya katılan önce üniversite kurumsal inceleme Kurulu tarafından onaylanan izin formu imzaladı. Helsinki Bildirgesi altında yapılan. Deneme temel alınarak aşağıdaki protokol gerçekleştirildi.

İleri dinamik simülasyon doğruluğunu denetlemek için alt-vücut simülasyon ortak açılardan biyomekanik a...

Tartışma

Bu çalışmada tibia deformasyon yüksek darbe etkinlikleri sırasında belirlemek için non-invaziv bir yöntem geliştirmekti. Tibia zorlanma nedeniyle etkisi yükleme miktarının daha iyi tibia stres kırığı anlamak için yol açacaktır. Bu çalışmada bir konuya özgü kas-iskelet modeli geliştirilmiştir ve bilgisayar simülasyonları bir laboratuvar ortamında gerçekleştirilen damla-Açılış hareketleri çoğaltmak için işletilmiştir. Tibial zorlanma damla-açılış yükseklik etkisi incelenmiştir. Bu çalışmada, biz bu bırakma-açılış olarak gö...

Malzemeler

Bu makalede kullanılan malzemelerin listesi

Ad	Şirket	Katalog numarası	Yorumlar
CT Tarayıcı	GE Tıbbi Sistem	N/A	Işık Hızı VCT. Tibia BT taraması yapmak için.
Hareket Yakalama Sistemi	Vicon Inc	N / A	Vicon FX40 yüksek hızlı kameralar. 3D hareket yakalama gerçekleştirmek için.
Kuvvet plakaları	AMTI Inc	N/A	3D zemin reaksiyon kuvvetlerinin toplanması
Vicon Nexus	Vicon Inc	N/A	Hareket yakalama yazılım programı. Görsel işaretleyici yörünge verilerini işlemek için.
Visual 3D	C-Motion Inc	N/A Biyomekanik	analiz yazılımı. İnsan hareketlerinin 3B kinematiğini ve kinetiğini hesaplamak için.
MATLAB	Mathworks Inc	N/A	Bilgisayar programlama yazılımı. Ham veri filtreleme, veri dönüştürme ve veri işleme gerçekleştirmek için.
ADAMS 2012	MSC Software Inc	N/A	Çok gövdeli dinamik bilgisayar simülasyon programı.
LifeMOD	, Lifemodeler, Inc	N/A	ADAMS'ta bir yazılım eklentisi. İnsan vücudu kas-iskelet modelleri oluşturmak için.
MIMICS 13	Materialise Inc	N/A	Görüntü işleme programı. Görüntüleme verilerini işlemek için bir 3B modelleme aracı. CT taramalarından 3D tibia modeli oluşturmak için.
MARC 2012	MSC Software Inc	N/A	Sonlu elemanlar analiz yazılımı. Hacimsel ağ oluşturma gerçekleştirmek, tibia FE modeli oluşturmak ve modal FE analizini çalıştırmak için.
SPSS 19	IBM Inc	N/A	İstatistiksel analiz yazılımı.