1
|
Yu L, He W, Qin W, Wang K, Guo W, Wang S. Noninvasive computed tomography derived fractional flow reserve simulation based on microvascular tree model reconstruction. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3643. [PMID: 36054275 DOI: 10.1002/cnm.3643] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 07/22/2022] [Accepted: 08/17/2022] [Indexed: 06/15/2023]
Abstract
To establish a novel method for noninvasive computed tomography derived fractional flow reserve (CT-FFR) simulation based on microvascular tree model reconstruction and to evaluate the feasibility and diagnostic performance of the novel method in coronary artery disease compared with invasive fractional flow reserve (FFR). Twenty patients (20 vessels) who underwent coronary computed tomography angiography (CCTA) and invasive FFR were retrospectively studied. The anatomic epicardial coronary artery model was reconstructed based on CCTA image, and the microvascular tree model was simulated based on patient-specific anatomical structures and physiological principles. Numerical simulation was subsequently performed using the CFD method with full consideration of the variation of viscosity in microvascular. Two patients with the FFR value of .80 were selected for adjusting the parameters of the model, while the remaining 18 patients were selected as a validation cohort. After simulation, CT-FFR was compared with invasive FFR with a threshold of .80. Eleven (55%) patients had an abnormal FFR that was less than or equal to .80. Compared with invasive FFR, sensitivity, specificity, positive predictive value, negative predictive value, and accuracy of CT-FFR with an optimal threshold of .80 were 100%, 77.8%, 81.8%, 100%, 88.89%, respectively. There were a good correlation and consistency between CT-FFR and invasive FFR. Time per patient of CT-FFR analysis was less than 15 min. CT-FFR based on microvascular tree model reconstruction is feasible with good diagnostic performance. It requires a short processing time with excellent accuracy. Large multicenter prospective studies are required for further demonstrating the diagnostic performance of this novel model in myocardium ischemia evaluation.
Collapse
Affiliation(s)
- Long Yu
- Department of aeronautics and astronautics, Fudan University, Shanghai, China
| | - Wei He
- Department of Vascular Surgery, Zhongshan Hospital, Fudan University, Shanghai, China
| | - Wang Qin
- Department of aeronautics and astronautics, Fudan University, Shanghai, China
| | - Keqiang Wang
- Institute of Panvascular Medicine, Zhongshan Hospital, Fudan University, Shanghai, China
| | - Weifeng Guo
- Department of Radiology, Zhongshan Hospital, Fudan University, Shanghai, China
| | - Shengzhang Wang
- Department of aeronautics and astronautics, Fudan University, Shanghai, China
- Institute of Biomedical Engineering Technology, Academy for Engineering and Technology, Fudan University, Shanghai, China
| |
Collapse
|
2
|
Wang J, Fan T, Zhang H, Ge Y, Lu W, Liu F, Rong D, Guo W. Aortic hemodynamic and morphological analysis before and after repair of thoracoabdominal aortic aneurysm using a G-Branch endograft. Front Physiol 2023; 14:1234989. [PMID: 37601633 PMCID: PMC10438984 DOI: 10.3389/fphys.2023.1234989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Accepted: 07/19/2023] [Indexed: 08/22/2023] Open
Abstract
Background and objective: The G-Branch endograft is a novel multibranched "off-the-shelf" device used to repair thoracoabdominal aortic aneurysms (TAAAs). This report describes the hemodynamic and morphological performance of the G-Branch endograft in a human patient with TAAA. Materials and methods: We retrospectively reviewed the computed tomography angiography scans and clinical data of a woman in whom TAAA was treated using a G-Branch endograft. Patient-specific three-dimensional models were reconstructed, and computational fluid dynamics and morphological and hemodynamic indicators were analyzed before and after implantation of the device. Results: From a morphological perspective, there was an increase in cross-sectional area in the G-Branch endograft and all bridging stent grafts over time. Blood flow was redistributed among the renovisceral arteries, with a decrease in flow rate in the celiac artery and an increase in the left renal artery. Laminar blood flow was smoother and more rapid after implantation of the G-Branch device and remained stable during follow-up. In the bridging stent grafts, flow recirculation zones were found in the bridging zones of the celiac artery and superior mesenteric artery as well as the distal sealing zones of both renal arteries. Furthermore, higher time-averaged wall shear stress and a lower oscillatory index and relative resident time were found in the G-Branch endograft and bridging stent grafts. Quantitative analysis showed obvious reduction in the surface area ratio of the elevated time-averaged wall shear stress area and surface area ratio of the relative resident time after G-branch implantation. Conclusion: The revascularization of branch vessels occurred following G-branch implantation, with improvements arising not only from morphological changes but also from hemodynamic alterations. The long-term performance of the G-Branch endograft needs further investigation and clinical validation.
Collapse
Affiliation(s)
- Jiabin Wang
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
- Medical School of Chinese PLA, Beijing, China
| | - Tingting Fan
- School of Biomedical Engineering, Capital Medical University, Beijing, China
| | - Hongpeng Zhang
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| | - Yangyang Ge
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| | - Weihang Lu
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| | - Feng Liu
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| | - Dan Rong
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| | - Wei Guo
- The First Medical Centre, Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China
| |
Collapse
|
3
|
Liu D, Wang X, Zhao D, Sun Z, Biekan J, Wen Z, Xu L, Liu J. Influence of MRI-based boundary conditions on type B aortic dissection simulations in false lumen with or without abdominal aorta involvement. Front Physiol 2022; 13:977275. [PMID: 36160847 PMCID: PMC9490059 DOI: 10.3389/fphys.2022.977275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Accepted: 08/15/2022] [Indexed: 11/13/2022] Open
Abstract
Most computational hemodynamic studies of aortic dissections rely on idealized or general boundary conditions. However, numerical simulations that ignore the characteristics of the abdominal branch arteries may not be conducive to accurately observing the hemodynamic changes below the branch arteries. In the present study, two men (M-I and M-II) with type B aortic dissection (TBAD) underwent arterial-phase computed tomography angiography and four-dimensional flow magnetic resonance imaging (MRI) before and after thoracic endovascular aortic repair (TEVAR). The finite element method was used to simulate the computational fluid dynamic parameters of TBAD [false lumen (FL) with or without visceral artery involvement] under MRI-specific and three idealized boundary conditions in one cardiac cycle. Compared to the results of zero pressure and outflow boundary conditions, the simulations with MRI boundary conditions were closer to the initial MRI data. The pressure difference between true lumen and FL after TEVAR under the other three boundary conditions was lower than that of the MRI-specific results. The results of the outflow boundary conditions could not characterize the effect of the increased wall pressure near the left renal artery caused by the impact of Tear-1, which raised concerns about the distal organ and limb perfused by FL. After TEVAR, the flow velocity and wall pressure in the FL and the distribution areas of high time average wall shear stress and oscillating shear index were reduced. The difference between the calculation results for different boundary conditions was lower in M-II, wherein FL did not involve the abdominal aorta branches than in M-I. The boundary conditions of the abdominal branch arteries from MRI data might be valuable in elucidating the hemodynamic changes of the descending aorta in TBAD patients before and after treatment, especially those with FL involving the branch arteries.
Collapse
Affiliation(s)
- Dongting Liu
- Department of Radiology, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| | - Xuan Wang
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Dongliang Zhao
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Zhonghua Sun
- Discipline of Medical Radiation Science, Curtin Medical School, Curtin University, Perth, WA, Australia
- *Correspondence: Jiayi Liu, ; Zhonghua Sun,
| | | | - Zhaoying Wen
- Department of Radiology, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| | - Lei Xu
- Department of Radiology, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| | - Jiayi Liu
- Department of Radiology, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
- *Correspondence: Jiayi Liu, ; Zhonghua Sun,
| |
Collapse
|
4
|
A dry immersion model of microgravity modulates platelet phenotype, miRNA signature, and circulating plasma protein biomarker profile. Sci Rep 2021; 11:21906. [PMID: 34753989 PMCID: PMC8578674 DOI: 10.1038/s41598-021-01335-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 10/26/2021] [Indexed: 11/08/2022] Open
Abstract
Ground based research modalities of microgravity have been proposed as innovative methods to investigate the aetiology of chronic age-related conditions such as cardiovascular disease. Dry Immersion (DI), has been effectively used to interrogate the sequelae of physical inactivity (PI) and microgravity on multiple physiological systems. Herein we look at the causa et effectus of 3-day DI on platelet phenotype, and correlate with both miRomic and circulating biomarker expression. The miRomic profile of platelets is reflective of phenotype, which itself is sensitive and malleable to the exposome, undergoing responsive transitions in order to fulfil platelets role in thrombosis and haemostasis. Heterogeneous platelet subpopulations circulate at any given time, with varying degrees of sensitivity to activation. Employing a DI model, we investigate the effect of acute PI on platelet function in 12 healthy males. 3-day DI resulted in a significant increase in platelet count, plateletcrit, platelet adhesion, aggregation, and a modest elevation of platelet reactivity index (PRI). We identified 15 protein biomarkers and 22 miRNA whose expression levels were altered after DI. A 3-day DI model of microgravity/physical inactivity induced a prothrombotic platelet phenotype with an unique platelet miRNA signature, increased platelet count and plateletcrit. This correlated with a unique circulating protein biomarker signature. Taken together, these findings highlight platelets as sensitive adaptive sentinels and functional biomarkers of epigenetic drift within the cardiovascular compartment.
Collapse
|
5
|
Feng Y, Liu J, Fan T, Zhang W, Yin X, E Y, Tan W, Huo Y. Vertebral Artery Stenoses Contribute to the Development of Diffuse Plaques in the Basilar Artery. Front Bioeng Biotechnol 2020; 8:168. [PMID: 32211395 PMCID: PMC7068210 DOI: 10.3389/fbioe.2020.00168] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2019] [Accepted: 02/19/2020] [Indexed: 12/19/2022] Open
Abstract
Vertebral artery (VA) stenosis is relevant to a high early risk of recurrent stroke and basilar artery (BA) is the most common intracranial site of atherosclerotic lesions. It is important to show predictive risk factors for transient ischemic attack (TIA) or posterior infarctions. The aim of the study is to investigate morphometry and hemodynamics in intracranial vertebral and basilar arteries of health and diseased patients to enhance the risk assessment. Based on the geometrical model reconstructed from CTA images in 343 patients, a transient three-dimensional computational model was used to determine the hemodynamics. Patients were classified in symmetric, asymmetric, hypoplastic, and stenotic groups while patients in the stenotic group were divided into unilateral, bilateral, bifurcation, and tandem stenotic sub-groups. Patients in bilateral, bifurcation, and tandem stenotic sub-groups had significantly lower basilar artery diameters than other groups. Patients in the stenotic group had significantly higher surface area ratio (SAR) of high time-averaged wall shear stress gradient (TAWSSG) and higher incidence of TIAs or posterior infarctions than other groups while patients in the tandem stenotic sub-group had the highest values (SAR-TAWSSG of 57 ± 22% and TIAs or posterior infarction incidence of 54%). The high SAR-TAWSSG is predisposed to induce TIAs or posterior infarction.
Collapse
Affiliation(s)
- Yundi Feng
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China.,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China
| | - Jian Liu
- Department of Cardiology, Peking University People's Hospital, Beijing, China
| | - Tingting Fan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Wenxi Zhang
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Xiaoping Yin
- Department of Radiology, Affiliated Hospital of Hebei University, Hebei University, Baoding, China
| | - Yajun E
- Department of Neurology, Affiliated Hospital of Hebei University, Hebei University, Baoding, China
| | - Wenchang Tan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China.,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China.,Shenzhen Graduate School, Peking University, Shenzhen, China
| | - Yunlong Huo
- PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China.,Institute of Mechanobiology & Medical Engineering, School of Life Sciences & Biotechnology, Shanghai Jiao Tong University, Shanghai, China
| |
Collapse
|
6
|
Fan T, Zhou Z, Fang W, Wang W, Xu L, Huo Y. Morphometry and hemodynamics of coronary artery aneurysms caused by atherosclerosis. Atherosclerosis 2019; 284:187-193. [DOI: 10.1016/j.atherosclerosis.2019.03.001] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/03/2018] [Revised: 02/18/2019] [Accepted: 02/22/2019] [Indexed: 12/20/2022]
|
7
|
Morphometry and hemodynamics of posterior communicating artery aneurysms: Ruptured versus unruptured. J Biomech 2018; 76:35-44. [DOI: 10.1016/j.jbiomech.2018.05.019] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2017] [Revised: 01/12/2018] [Accepted: 05/10/2018] [Indexed: 11/20/2022]
|
8
|
Fan T, Feng Y, Feng F, Yin Z, Luo D, Lu Y, Xu Y, Tan W, Huo Y. A comparison of postoperative morphometric and hemodynamic changes between saphenous vein and left internal mammary artery grafts. Physiol Rep 2018; 5:5/21/e13487. [PMID: 29122958 PMCID: PMC5688779 DOI: 10.14814/phy2.13487] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Revised: 09/24/2017] [Accepted: 10/05/2017] [Indexed: 11/24/2022] Open
Abstract
There is higher long‐term failure of the saphenous vein graft (SVG) compared with the left internal mammary artery (LIMA) graft, which is affected by the hemodynamic environment. A comprehensive analysis of postoperative structure‐function changes is important to study the atherogenesis in the SVG. A comparison of morphometric and hemodynamic parameters was carried out between LIMA grafts and SVGs and between different time points postoperatively. Various parameters were obtained from the image reconstruction and flow simulation in patients, who underwent CT exams for ~1 year, 5 and 10 years after revascularization. Morphometric data showed a decrease in lumen size in the entire SVG and anastomosis of different patients in a sequence of ~1 year, 5 and 10 years postoperatively despite negligible changes of LIMA size. Computational results indicated the fourfold increased surface area ratio (SAR) of low time‐averaged wall shear stress (TAWSS) in the SVG and anastomosis at postoperative 10 years than that at postoperative ~1 year. The SAR of high TAWSS gradient (TAWSSG) at the distal anastomosis between SVG and coronary arteries was significantly higher (14 ± 9% vs. 6 ± 8%) than that in the LIMA group at postoperative ~1 year. There were strong correlations between morphometric and hemodynamic parameters in the SVG and distal anastomosis at various time points postoperatively, which showed deterioration relevant to persistent diffuse diseases at postoperative ~10 years.
Collapse
Affiliation(s)
- Tingting Fan
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Yundi Feng
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Feng Feng
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Zhongjie Yin
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Dayou Luo
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Yuan Lu
- Department of Cardiology, Affiliated Hospital of Xuzhou Medical University, Xuzhou, China
| | - Yingjin Xu
- Department of Radiology, Affiliated Hospital of Hebei University Hebei University, Baoding, China
| | - Wenchang Tan
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China .,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China.,Shenzhen Graduate School, Peking University, Shenzhen, China
| | - Yunlong Huo
- Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China .,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China
| |
Collapse
|
9
|
Feng Y, Wang X, Fan T, Li L, Sun X, Zhang W, Cao M, Liu J, Li J, Huo Y. Bifurcation Asymmetry of Small Coronary Arteries in Juvenile and Adult Mice. Front Physiol 2018; 9:519. [PMID: 29867562 PMCID: PMC5962776 DOI: 10.3389/fphys.2018.00519] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2018] [Accepted: 04/23/2018] [Indexed: 11/13/2022] Open
Abstract
Background: Microvascular bifurcation asymmetry is of significance for regulation of coronary flow heterogeneity during juvenile and adult growth. The aim of the study is to investigate the morphometric and hemodynamic variation of coronary arterial bifurcations in mice of different ages. Methods: Pulsatile blood flows were computed from a Womersley-type model in the reconstructed left coronary arterial (LCA) trees from Micro-CT images in normal mice at ages of 3 weeks, 6 weeks, 12 weeks, 5-6 months, and >8 months. Diameter and flow ratios and bifurcation angles were determined in each bifurcation of the LCA trees. Results: The blood volume and inlet flow rate of LCA trees increase and decrease during juvenile and adult growth, respectively. As vessel diameters decrease, the increased ratios of small to large daughter vessel diameters (Ds/Dl) result in more uniform flows and lower velocities. There are significant structure-functional changes of LCA trees in mice of >8 months compared with mice of < 8 months. As Ds/Dl increases, the variation trend of bifurcation angle during juvenile growth is different from that during adult growth. Conclusions: Although inlet flows are different in adult vs. juvenile mice, the adult still have uniform flow and low velocity. This is accomplished through a decrease in diameter. The design ensures ordered dispersion of red cells through asymmetric branching patterns into the capillaries.
Collapse
Affiliation(s)
- Yundi Feng
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Xuan Wang
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Tingting Fan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Li Li
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Xiaotong Sun
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Wenxi Zhang
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Minglu Cao
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Jian Liu
- Department of Cardiology, Peking University People's Hospital, Beijing, China
| | - Jianping Li
- Department of Cardiology, Peking University First Hospital, Beijing, China
| | - Yunlong Huo
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| |
Collapse
|
10
|
Yin X, Huang X, Li Q, Li L, Niu P, Cao M, Guo F, Li X, Tan W, Huo Y. Hepatic Hemangiomas Alter Morphometry and Impair Hemodynamics of the Abdominal Aorta and Primary Branches From Computer Simulations. Front Physiol 2018; 9:334. [PMID: 29674973 PMCID: PMC5895747 DOI: 10.3389/fphys.2018.00334] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2017] [Accepted: 03/16/2018] [Indexed: 01/25/2023] Open
Abstract
Background: The formation of hepatic hemangiomas (HH) is associated with VEGF and IL-7 that alter conduit arteries and small arterioles. To our knowledge, there are no studies to investigate the effects of HH on the hemodynamics in conduit arteries. The aim of the study is to perform morphometric and hemodynamic analysis in abdominal conduit arteries and bifurcations of HH patients and controls. Methods: Based on morphometry reconstructed from CT images, geometrical models were meshed with prismatic elements for the near wall region and tetrahedral and hexahedral elements for the core region. Simulations were performed for computation of the non-Newtonian blood flow using the Carreau-Yasuda model, based on which multiple hemodynamic parameters were determined. Results: There was an increase of the lumen size, diameter ratio, and curvature in the abdominal arterial tree of HH patients as compared with controls. This significantly increased the surface area ratio of low time-averaged wall shear stress (i.e., SAR-TAWSS =Surface areaTAWSS≤4 dynes·cm−2Total surface area× 100%) (24.1 ± 7.9 vs. 5 ± 6%, 11.6 ± 12.8 vs. < 0.1%, and 44.5 ± 9.2 vs. 21 ± 24% at hepatic bifurcations, common hepatic arteries, and abdominal aortas, respectively, between HH and control patients). Conclusions: Morphometric changes caused by HH significantly deteriorated the hemodynamic environment in abdominal conduit arteries and bifurcations, which could be an important risk factor for the incidence and progression of vascular diseases.
Collapse
Affiliation(s)
- Xiaoping Yin
- Department of Radiology, Affiliated Hospital of Hebei University, Hebei University, Baoding, China
| | - Xu Huang
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Qiao Li
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Li Li
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Pei Niu
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Minglu Cao
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
| | - Fei Guo
- College of Medicine, Hebei University, Baoding, China
| | - Xuechao Li
- College of Medicine, Hebei University, Baoding, China
| | - Wenchang Tan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China.,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China.,Shenzhen Graduate School, Peking University, Shenzhen, China
| | - Yunlong Huo
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China.,College of Medicine, Hebei University, Baoding, China.,PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China
| |
Collapse
|