應用有限元素分析於纖維化肺組織之力學研究

Application of Finite Element Analysis for Mechanics Study on Pulmonary with Fibrosis Diseases

廖晏晨(Yen-Chen Liao)

指導教授 ： 王建凱

淡江大學/工學院/土木工程學系碩士班/碩士(2019年)

https://doi.org/10.6846/TKU.2019.00447

摘要

有限元素法為學術研究與先進工程仰賴之計算方法，能夠針對具有幾何外型特徵之連續固體或結構系統，精準求得系統受外力刺激所造成的力學及物理反應的場域分佈，近年來除應用於固體力學與結構工程學門外，更受到生物醫學工程領域的重視。
本論文利用有限元素分析軟體FEBio結合美國國家衛生研究院(NIH)所提供的肺臟STL模型，模擬肺纖維化前後再進行吸氣時所受到的力學狀態，在解讀程式經計算後所得出的結果。研究主要可分為三個階段，分別為MATLAB程式撰寫、模型與分析條件設置、歸納分析之成果。
首先先介紹MATLAB程式撰寫，本研究肺臟模型所希望之元素設定為板殼元素，但由於STL模型匯入FEBio中無法設定為板殼元素，所以利用MATLAB開發轉檔程式，藉由讀取設定分析參數的文字輸入檔和STL模型檔，將參數與模型資料寫成FEBio專有的FEB檔案格式儲存，便可解決板殼元素設定之問題；二為模型與分析條件設置，根據參考文獻設計纖維化區域為散佈於肺葉上與集中在支氣管上，另外也尋找相關資料設計模型的邊界條件參數；最後階段為歸納分析結果，在同樣的分析條件下，從兩個肺纖維化模型的分析結果去比較，並找出相關的規則加以歸納之。

關鍵字

肺纖維化 ； 生物計算力學 ； 有限元素法 ； 板殼元素 ； FEBio software suite ； 程式語言設計

並列摘要

The finite element method is a calculation method based on academic research and advanced engineering. It can accurately determine the field distribution of mechanical and physical reactions caused by external force stimulation for continuous solid or structural systems with geometric appearance characteristics. Applied to solid mechanics and structural engineering, it is also valued in the field of biomedical engineering.
In this thesis, the finite element analysis software FEBio is combined with the lung STL model provided by the National Institutes of Health (NIH) to simulate the mechanical state of inhalation before and after pulmonary fibrosis. The results obtained after the calculation of the interpretation program are calculated. The research can be divided into three stages, which are the results of MATLAB programming, model and analysis condition setting, and inductive analysis.
First, the MATLAB programming is introduced. The desired element of the lung model is set as the shell element. However, since the STL model cannot be set as the shell element in FEBio, the MATLAB development of the conversion program is used to analyze the reading configuration. The text input file and STL model file of the parameter, the parameters and model data are written into FEBio's proprietary FEB file format storage, which can solve the problem of setting the plate and shell elements. Secondly, the model and analysis conditions are set, and the fiberized area is designed according to the reference. In order to spread on the lungs and concentrate on the bronchus, the boundary condition parameters of the relevant data design model are also sought. The final stage is the inductive analysis result. Under the same analysis conditions, the analysis results of the two pulmonary fibrosis models are compared. And find out the relevant rules to summarize them.

並列關鍵字

Pulmonary fibrosis ； biocomputational mechanics ； finite element method ； plate and shell elements ； FEBio software suite ； programming language design

參考文獻

[1] 吳常瑋、張皓鈞、陳天華、黃信端、劉家榮、王鶴健(2018)，「特發性肺纖維化」，內科學誌 2018：29：283-291，國立台灣大學醫學院附設醫院 內科部

[2] 林承志(2013)，「病因不明性肺纖維化診斷與治療新知」，內科學誌 2013：24：446-452，國軍高雄總醫院左營分院 內科部

[3] 謝英恒、蘇浤傑(2017)，「台灣懸浮微粒與人類季節性流感相關性之探討」，中國醫藥大學

[4] Steve Maas, Dave Rawlins, Dr. Jeffrey Weiss, Dr. Gerard Ateshian (2018). “FEBio User’s Manual Version 2.8”, Musculoskeletal Research Laboratories, University of Utah.

[5] Steve Maas, Dave Rawlins, Dr. Jeffrey Weiss, Dr. Gerard Ateshian (2018). “FEBio Theory Manual Version 2.7”, Musculoskeletal Research Laboratories, University of Utah.