隨著科技不斷地進步，人類在醫學生物體上的研究也日益更新，生物體的軟組織研究，如皮膚、眼睛、關節、心臟、肺部等等，這些都與人體的健康密不可分，而目前關於人體肺部之慢性疾病影響，例如特發性肺纖維化(IPF)、哮喘(Asthma)、慢性阻塞性肺部疾病(COPD)和2019年爆發的新型冠狀病毒肺炎(COVID-19)這些疾病為主要全球高齡化人口之主要死因。肺組織就像許多生物軟組織一樣，都是具有高度黏彈性的機械性質，且肺部的機械性質與上述的這些疾病有著相當大的關係。因此若能更進一步的了解這些組織特性，研究其性質，並將其數據應用於建立人工模型，則從實驗獲取數據的方式變得非常重要。
在本研究當中，將會對兩隻豬之新鮮豬肺以圓柱形平端壓頭與圓球形壓頭進行較淺與較深之壓痕試驗與應力鬆弛實驗，將數據以均方根誤差法擬合Generalized Maxwell黏彈性模型獲得參數，並進行建模擬合與實驗比較，最後在加入暫態擠壓彈液動潤滑系統探討是否有助於改善模型，使之更接近真實情況。
實驗顯示兩種壓頭在不同深度實驗上大致上都擁有相同之趨勢，而在第二次實驗改善初始接觸之問題能使結果更為真實且使再現性提升，並藉由擬合參數與建模方式來重現實驗之趨勢。而以本研究擬合中間深度之黏彈參數，去模擬不同深度並與實驗值比較發現，在壓入淺深度的區域中能與實驗值接近，而在較深的區域加入彈液動潤滑模型能使模擬結果更靠近實驗值。

With the advance of science and technology,research on medical organisms becoming more and more updated. The soft tissue research of organisms, such as skin、eyes、joints、heart、lungs、etc.. are inseparable from human health. The current impact of chronic diseases in the human lungs, such as Idiopathic pulmonary fibrosis (IPF)、Asthma、Chronic dual lung disease (COPD) and Coronavirus disease 2019 (COVID-19)。These diseases,like many biological soft tissues,have highly viscoelastic mechanical properties,for lung this property has a considerable relationship with the above-mentioned diseases.Therefore,if we can further understand these organizational characteristics,study their properties, and merge their data to build artificial models,the way to obtain data from experiments becomes very important.

In this study, the fresh pig lungs of two pigs will be subjected to shallower and deeper indentation tests and stress relaxation tests with a cylindrical flat-end indenter and a spherical indenter, and the data will be simulated by the root mean square error method. Integrate the Generalized Maxwell viscoelastic model to obtain the parameters, and carry out the modeling fitting and experiment comparison.Experiments show that the two types of indenters generally have the same trend in different depth experiments. In the second experiment, improving the initial contact problem can make the results more realistic and improve the reproducibility.

After fitting the simulated value of the experimental result of the intermediate depth, and adding the elastic fluid motion, the result shows that it can be closer to the real situation.

[1] MS. Wilson, TA,Wynn, “Pulmonary fibrosis: pathogenesis, etiology and regulation,” Mucosal Immunology, vol. 2, no. 2, pp. 103-121, 2009.
[2] JC. Hogg, W,Timens, “The Pathology of Chronic Obstructive Pulmonary Disease,” Annual Review of Pathology-Mechanisms of Disease, vol. 4, pp. 435-459, 2009.
[3] HS. Shi, XY. Han, “Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: a descriptive study,” Lancet Infectious Diseases, vol. 20, no. 4, pp. 425-434, 2020.
[4] P. Burney, D. Jarvis, “The global burden of chronic respiratory disease in adults,” International Journal of Tuberculosis And Lung Disease, vol. 19, no. 1, pp. 10-20, 2015.
[5] J. Jaffar,SH. Yang, “Greater cellular stiffness in fibroblasts from patients with idiopathic pulmonary fibrosis,” American Journal of Physiology-Lung Cellular And Molecular Physiology, vol. 315, no. 1, pp. 59-65, 2018.
[6] T. Ebihara, N, Venkatesan, “Changes in Extracellular Matrix and Tissue Viscoelasticity in Bleomycin-Induced Lung Fibrosis - Temporal Aspects,” American Journal of Respiratory and Critical Care Medicine,vol. 162, no. 4, pp. 1569-1576.
[7] M. Eskandari, AL. Arvayo, ME. Levenston, “Mechanical Properties of the Airway Tree: Heterogeneous and Anisotropic Pseudoelastic and Viscoelastic Tissue Responses,” Journal of Applid Physiology, vol. 123, no. 3, pp. 878-888.
[8] M. Griffin, Y. Premakumar, A. Seifalian, P. E. Butler, and M. Szarko,“Biomechanical Characterization of Human Soft Tissues Using Indentation andTensile Testing,” Journal of Visualized Experiments, no. 118, Dec.2016.
[9] R. M. Delaine-Smith, S. Burney, F. R. Balkwill, and M. Knight,“Experimental Validation of a Flat Punch Indentation Methodology Calibratedagainst Unconfined Compression Tests for Determination of Soft TissueBiomechanics,” Journal of the Mechanical Behavior of Biomedical Materials,vol. 60, pp. 401–415, 2016.
[10] M. H. Lu, W. Yu, Q. H. Huang, Y. P. Huang, and Y. P. Zheng, “A Hand-HeldIndentation System for the Assessment of Mechanical Properties of Soft Tissuesin Vivo,” IEEE Transactions on Instrumentation and Measurement, vol. 58, no.9, pp. 3079–3085, 2009.
[11] Y. Zheng, A. F. Mak, and B. Lue, “Objective Assessment of Limb TissueElasticity: Development of a Manual Indentation Procedure.,” Journal of Rehabilitation Research and Development, vol. 36, no. 2, pp. 71–85, 1999.
[12] S. Budday, R. Nay, R. deRooij, P. Steinmann, T. Wyrobek, T. C. Ovaert, and E.Kuhl, “Mechanical Properties of Gray and White Matter Brain Tissue byIndentation,” Journal of the Mechanical Behavior of Biomedical Materials, vol.46, pp. 318–330, 2015.
[13] L. Boltzmann , “Zur Theorie der elastischen Nachwirkung,1874.
[14] F. J. Lockett, “Nonlinear viscoelastic solids,”, 1972.
[15] F. Bueche, “Journal of Applied Physics, ” ,vol. 26(6), pp. 738-749, 1955.
[16] J. Hildebrandt, “Pressure-Volume Data of Cat Lung Interpreted by a Plastoelastic, Linear Viscoelastic Model, Journal of Applied Physiology, vol. 28, no. 3, pp. 365-372, 1970.
[17] D. T. V. Pawluk, “Dynamic Contact of the Human Fingerpad Against a Flat Surface, ” Journal of Biomechanical Engineering, vol. 121, pp. 605-611, 1999.
[18] M. Rodríguez-Nieto, “Viscoelastic properties of doxorubicin‑treated HT‑29 cancer cells by atomic force microscopy: the fractional Zener model as an optimal viscoelastic model for cells, ” Biomechanics And Modeling in Mechanobiology, vol. 7, no. 2-3, pp. 147-152, 2019.
[19] JHT. Bates, BS. Ma, “A Progressive Rupture Model of Soft Tissue Stress Relaxation, ” Annals of Biomedical Engineering, vol. 41, no. 6, pp. 1129-1138, 2013.
[20] B. Babaei, A. Davarian, “Efficient and optimized identification of generalized Maxwell viscoelastic relaxation spectra, ” Journal of the Mechanical Behavior of Biomedical Materials, vol. 55, no. 6, pp. 32-41, 2016.
[21] E. Bayliss, G. W. Robertson, “The Visco-Elastic Properties of the Lungs,” Quarterly Journal of Experimental Physiology and Cognate Medical Sciences,vol. 29, pp. 27-47.
[22] L. E. Mount, “The Ventilation Flow-Resistance and Compliance of Rat Lungs,” Journal of Physiology-London,vol. 127, no. 1, pp. 157-167, 1955.
[23] R. Marshall, J. G. Widdicombe, “Stress Relaxation of Human Lung,” Clinical Science,vol. 20, no. 1, pp. 19-31.
[24] F. G. Salerno, M. S. Ludwig, “Elastic Moduli of Excised Constricted Rat Lungs, ” Journal of Applied Physiology, vol. 86, no. 1, pp. 66-70, 1999.
[25] B. Weed, S. Patnaik, “Experimental Evidence of Mechanical Isotropy in Porcine Lung Parenchyma, ” Materials, vol. 8, no. 5, pp. 2454-2466, 2015.
[26] R. C. Tai, G. C. Lee, “Isotropy and Homogeneity of Lung Tissue Deformation, ” Journal of Biomechanics, vol. 14, no. 4, pp. 243-252, 1981.
[27] AM. Birzle, WA. Wall, “A Viscoelastic Nonlinear Compressible Material Model of Lung Parenchyma - Experiments and Numerical Identification, ” Journal of the Mechanical Behavior of Biomedical Materials, vol. 84, pp. 164-175, 2019.
[28] R. H. J. de Hilster, PK. Sharma, “Human Lung Extracellular Matrix Hydrogels Resemble the Stiffness and Viscoelasticity of Native Lung Tissue,” American Journal of Physiology-Lung Cellular And Molecular Physiology, vol. 318, no. 4, pp. 698-704, 2020.
[29] Z. Dai, Y. Peng, HA .Mansy, “A Model of Lung Parenchyma Stress Relaxation Using Fractional Viscoelasticity,” Medical Engineering And Physics, vol. 37, no. 8, pp. 752-758, 2015.
[30] C. Wex, S. Arndt, A. Stoll, C. Bruns, and Y. J. B. E. B. T. Kupriyanova, “Isotropic incompressible hyperelastic models for modelling the mechanical behaviour of biological tissues: a review," vol. 60, no. 6, pp. 577-592, 2015.
[31] C. Roelands, “Correlational aspects of the viscosity-temperature-pressure relationship of lubricating oils,” University of Groningen, Netherlands, 1966.
[32] D. Dowson and G. R. Higginson, “Elasto-hydrodynamic Lubrication: The Fundamentals of Roller And Gear Lubrication.” Pergamon Press, 1965.
[33] H. Christensen, “Elastohydrodynamic Theory of Spherical Bodies in Normal Approach, ”Journal of Lubrication Technology, vol. 92,no. 1,pp 145-153, 1970.
[34] H. Noori-Dokht, “Finite element analysis of elastohydrodynamic lubrication in an artificial hip joint under squeeze film motion using fluid–structure interaction method,” Proceedings of the Institution of Mechanical Engineers Part J-Journal of Engineering Tribology, vol. 231,no. 9,pp 1171-1181, 2017.
[35] HM. Chu, WL. Li, MD,Chen , “Elastohydrodynamic Lubrication of Circular Contacts at Pure Squeeze Motion with non-Newtonian Lubricants.” Tribology International, vol. 39, no. 9, pp.897-905, 2006.
[36] H. Hertz, “On the contact of elastic solids, ” Z. Reine Angew. Mathematik, vol. 92, pp. 156-171, 1881. S.
[37] Timoshenko, J. N. Goodier, “Theory of Elasticity,” 3rd , 1951.
[38] E. Riande, R. Diaz-Calleja, “Polymer Viscoelasticity: Stress and Strain in Practice,” CRC Press Inc, 1999.
[39] A. Drozdov, Mechanics of Viscoelastic Solids, Wiley,1998.
[40] Teodorovich, E. V. (1978). Sliding of a cylinder on a viscoelastic foundation. Journal of Applied Mathematics and Mechanics, 42(2), 384–389.
[41] Wang, X., Schoen, J. A., & Rentschler, M. E. (2013). A quantitative comparison of soft tissue compressive viscoelastic model accuracy. Journal of the Mechanical Behavior of Biomedical Materials, 20, 126–136.