Name: Intravascular Ultrasound Image-Based Finite Element Modeling Approach for Quantifying In Vivo Mechanical Properties of Human Coronary Artery (Video) | JoVE Journal
Uploaded: 2024-12-06T14:20:00.000+00:00
Duration: 6 min 18 s
Description: Discover finite element modeling using intravascular ultrasound images to quantify coronary artery mechanical properties in vivo. This computational approach enables patient-specific cardiovascular diagnostics.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

This research was supported in part by Shandong Province Medical Health Science and Technology Project (No. 202425020256, and 202403010254), National Natural Science Foundation of China grants 11972117 and 11802060, the Natural Science Foundation of Jiangsu Province under grant number BK20180352, and the Natural Science Foundation of Shandong Province under grant number ZR2024QA110.

De-identified clinical data, including in vivo IVUS images and blood pressure data, were acquired from a patient with CAD at Zhongda Hospital, Southeast University, with informed consent obtained. The sample patient was selected from the patient pool of a clinical study on intermediate coronary atherosclerotic lesions to demonstrate the method for quantifying the material properties of patient-specific coronary vessels14. The study was conducted following the protocol approved by the Clinical Research Ethics Committee of Zhongda Hospital, Southeast University (approval number: 2017ZDSYLL023-p01).
1. Data acquisition and processing
<ol>
	<li>Cine IVUS and VH-IVUS image acquisition
	<ol>
		<li>Place the IVUS catheter distal to the atherosclerotic lesion with the guidance of a coronary angiogram and pull it back to the proximal side. Generate grayscale IVUS images along the way to visualize the cross-section of the coronary vessel.</li>
		<li>Use the IVUS imaging system equipped with a high-performance 20 MHz, 2.9F platinum catheter to acquire IVUS images of the coronary arterial segment with atherosclerotic plaque from the patient with CAD (Figure 1).</li>
		<li>During the IVUS image acquisition, pause the catheter at the pre-selected lesion site for approximately 2 s to obtain a sequence of IVUS images called cine IVUS. The cine IVUS images clearly showed the dynamic cross-sectional changes at the given plaque site over the cardiac cycle.</li>
		<li>Generate VH-IVUS images based on IVUS frames acquired at the instant of R peak on electrocardiograms to visualize plaque components in color-coded form using the imaging system. 
		NOTE: VH-IVUS images provide an intuitive color map for four key plaque components in the atherosclerotic plaque: lipid-rich necrotic core (lipid) in red, calcification in white, fibrous tissue in dark green, and fibrous adipose tissue in light green.</li>
		<li>Save VH-IVUS and cine IVUS images in DICOM format for offline analysis.</li>
	</ol>
	</li>
	<li>Image segmentation and processing
	<ol>
		<li>Open DICOM files using the viewer, double-click the Corresponding Sequence Name to open the image and click Export &gt; Export Images to save each cine IVUS frame or VH-IVUS frame as an individual image in BMP format. Each BMP image contains 500 x 500 pixels, as indicated by the original DICOM file.</li>
		<li>Examine the cine IVUS image frame by frame to find the consecutive frames obtained at the pre-selected plaque site during one cardiac cycle. There were 26 cine IVUS frames generated during one cardiac cycle for this sample plaque site.</li>
		<li>Examine all the generated VH-IVUS images to find the VH-IVUS image obtained at the given plaque site. The VH-IVUS images used here were created using one IVUS frame from the cine IVUS frames in one cardiac cycle.</li>
		<li>Segment VH-IVUS and cine IVUS images using ImageJ software to obtain the contours of vessel boundaries and plaque component boundaries (See Figure 1C).
		<ol>
			<li>Select Straight &gt; Segmented Line tab and manually delineate the contours of the lumen, outer boundary of the coronary vessel, and plaque components on cine IVUS and VH-IVUS images. For the cine IVUS image, only segment the contours of the lumen and outer vessel boundary, while for the VH-IVUS image, segment the contours of the lumen, outer vessel boundary, and plaque component boundaries.</li>
			<li>For simplicity, keep only large lipid components for generating a finite element model and ignore small, isolated lipid components. Only one lipid was present in this plaque sample. Overlay the delineated contours onto the original images using Image &gt; Overlay &gt; Add Selection tab.</li>
			<li>Navigate to the To ROI Manager menu to manage contours, adjust properties, and set colors and line width to proper values for better visualization. Select the Properties tab, set the stroke color to a different color, and fill in the line width at width. Here, set the line colors as green, blue, and red for lumen, outer-boundary, and lipid contours, respectively, and the line width as 3.</li>
			<li>Smooth contours using Edit &gt; Selection &gt; Fit Spline from the command bar after selecting a specific contour to get a smooth contour. This operation utilizes spline curve fitting techniques to automatically smooth the contours.</li>
			<li>Click File &gt; Save As &gt; XY Coordinates tab to save point coordinates of each contour like lumen, outer-boundary, and plaque component in a separate txt file. This file contains the x and y coordinate values of the points constituting the contour, with pixels as the unit.</li>
		</ol>
		</li>
		<li>Record the real physical size of each pixel in cine IVUS and VH-IVUS images (denoted as pixel size) from the DICOM file. The real distance for one pixel in the IVUS data used here is 0.002 cm. This information would be used to convert the point coordinates with the pixel as the unit to real distance with cm as the unit.</li>
	</ol>
	</li>
	<li>Contour data processing
	<ol>
		<li>Cine IVUS contour data processing
		<ol>
			<li>Read txt files of lumen contours from all cine IVUS images in one cardiac cycle with MATLAB.</li>
			<li>Multiply all lumen contours by pixel size to get the actual size of the lumen contours.</li>
			<li>Calculate lumen circumferences for all lumen contours and identify the IVUS frames with maximum (Cmax) and minimum (Cmin) lumen circumferences, representing diastolic and systolic conditions, respectively.</li>
		</ol>
		</li>
		<li>VH-IVUS contour data processing
		<ol>
			<li>Read txt files of the contours of the lumen, outer-boundary, and plaque components from the VH-IVUS image with MATLAB.</li>
			<li>Multiply all contours by pixel size to get the actual size of all contours.</li>
			<li>Re-divide each contour into 100 equally spaced points and perform 2D smoothing to get new VH-IVUS contour data to replace the old ones.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
2. Finite element model
<ol>
	<li>Coronary vessel geometry reconstruction
	<ol>
		<li>Create one layer of contours in three-dimensional space by adding z coordinate value for all points of VH-IVUS contours, including lumen, outer-boundary, and lipid, and set z = 0 for all points (Figure 2).</li>
		<li>Create another layer of contours by adding z coordinate value for all points of VH-IVUS contours and reset z = 0.05 cm for all points. 
		NOTE: These two layers of contours reconstruct the 3D coronary vessel geometry for the thin-layer structure-only model by adding a fixed 0.05 cm layer thickness to VH-IVUS contours (Figure 2).</li>
	</ol>
	</li>
	<li>Finite element mesh generation
	<ol>
		<li>Create two auxiliary contours by linearly interpolating lumen and outer-boundary contours with weights 1/3 and 2/3 (Figure 3A) for each layer.</li>
		<li>Divide vessel area into 8 circumferential parts and 3 radial parts (see Figure 3B) by connecting the lumen/outer boundary to the closest point on the lipid contour (e.g., points A and B in Figure 3B) or two auxiliary contours with radial lines.</li>
		<li>Connect all points between layers with straight lines, forming a 3D structure with 3 x 8 volumes (Figure 3B). Divide each volume using hexahedral elements to generate the finite element mesh (Figure 3C) and different material groups (Figure 3D).</li>
		<li>Perform a mesh analysis by refining mesh density by 10% until changes in solutions &lt; 5%.</li>
	</ol>
	</li>
	<li>Material property definition
	<ol>
		<li>Use a modified anisotropic Mooney-Rivlin material model to describe the material properties of the coronary vessel wall. Coronary vessels and plaque components were assumed to be hyperelastic, anisotropic, nearly incompressible, and homogeneous materials, and the strain energy density function of the modified anisotropic Mooney-Rivlin material model is: 
		<img alt="Static equilibrium formula W=W&lt;sub&gt;iso&lt;/sub&gt;+W&lt;sub&gt;aniso&lt;/sub&gt;, equation, educational use." src="/files/ftp_upload/67573/67573eq01.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 1">     (1) 
		<img alt="Isothermal process equation W_iso for material science modeling." src="/files/ftp_upload/67573/67573eq02v2.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 2">   (2) 
		<img alt="Anisotropic energy equation W_aniso formula used in material science research analysis." src="/files/ftp_upload/67573/67573eq03.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 3">    (3) 
		where I1 and I2 are the first and second invariants of the right Cauchy-Green deformation tensor C defined as c = [cij] = XTX, X = [Xij] = [<img alt="Partial derivative ∂xi/∂aj, mathematical equation, symbolic representation, calculus concept." src="/files/ftp_upload/67573/67573eq06.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 6">], (Xi) was the current position (aj) was the original position, I4 = cij(nc)i(nc)j, nc was the unit vector in the circumferential direction of the vessel. c1, c2, D1, D2, K1 and K2 were patient-specific material parameters.</li>
		<li>Assign the initial values of material parameters for a patient-specific coronary vessel according to the ex vivo biaxial testing results, that is, c1 = −1,312.9 kPa, c2 = 114.7 kPa, D1 = 629.7 kPa, D2 = 2.0, K1 = 35.9 kPa and K2 = 23.5 (Figure 4A-B)13,15.</li>
		<li>Assign the material parameters for the plaque component if present. More specifically, for lipids, c1=0.5 kPa, c2=0, D1=0.5 kPa, and D2=1.5; for calcification, we used c1=920 kPa, c2=0, D1=360 kPa, and D2=2.0 (Figure 4B)16. 
		NOTE: Plaque components (lipid and calcification) were assumed to be hyperelastic, isotropic, and nearly incompressible, and their mechanical properties were described by isotropic Mooney-Rivlin material model with the strain energy density functions given in formula (2).</li>
	</ol>
	</li>
	<li>Governing equations and boundary conditions setting
	<ol>
		<li>Define governing equations for the thin-layer structure-only model, which includes motion equation, nonlinear Cauchy-Green strain-displacement relation, and coronary vessel material model11.</li>
		<li>Prescribe patient-specific blood pressure waveforms on the lumen surface to simulate real physiological conditions (Figure 4C). To obtain patient-specific blood pressure waveforms, scale a typical aortic pressure waveform with systolic and diastolic pressure values measured by arm cuff (Figure 4D).</li>
	</ol>
	</li>
</ol>
3. Finite element model-based updating approach for patient-specific coronary artery material properties
NOTE: The iterative process to determine patient-specific coronary material properties is illustrated in Figure 5.
<ol>
	<li>Determine the no-load geometry corresponding to the zero-pressure condition as the initial geometry for the computational model by shrinking the coronary geometry reconstructed from the VH-IVUS image axially with a fixed shrinkage rate of 95% and circumferentially with circumferential shrinkage (denoted as S) initially set as 98%. 
	NOTE: Since the coronary geometry reconstructed from the VH-IVUS image was under in vivo conditions with blood pressure prescribed on the lumen and axial stretch from tethered distal and proximal coronary arterial segments, the in vivo coronary geometry should shrink circumferentially and axially to obtain the zero-pressure geometry.</li>
	<li>Fix the axial shrinkage rate at 95% and update the circumferential shrinkage during the following steps.</li>
	<li>Define the material ratio (denoted as k) to assign the patient-specific material properties of the coronary vessel as: that is, c1 = k*(−1,312.9) kPa, c2 = k*114.7 kPa, D1 = k*629.7 kPa, K1 = k*35.9 kPa, and fix D2 = 2.0 and K2 = 23.5. 
	NOTE: Since only two data points (minimum and maximum lumen circumferences corresponding to diastolic and systolic pressures) were obtained to determine the unknown parameters (circumferential shrinkage rate S and material parameters of Mooney-Rivlin model), we reduced the number of unknown parameters by assuming that the in vivo patient-specific material properties of coronary vessel were proportional to the initial guess with material ratio denoted as k: that is, c1 = k*(−1,312.9) kPa, c2 = k*114.7 kPa, D1 = k*629.7 kPa, K1 = k*35.9 kPa, while D2 = 2.0 and K2 = 23.5 were fixed.</li>
	<li>Update the k value set to an initial value k of 1 along with circumferential shrinkage rate S during the following iterative procedure.</li>
	<li>Run software to solve the computational model to obtain the numerical results.
	<ol>
		<li>Write all commands for creating the thin-layer structure-only model into a batch file (Supplementary File 1) using MATLAB.</li>
		<li>Load this batch file using the advanced user interface (AUI) to generate the model (Figure 6A). Solve the thin-layer structure-only model by clicking Data File/Solution and save it as a .dat file (Figure 6C). Simulate three cardiac cycles and adopt the solution in the last cycle to present numerical results.</li>
		<li>Export results of node coordinates to a txt file by navigating List &gt; Value List &gt; Zone and selecting X-POSITION, Y-POSITION, and Z-POSITION in Variables to List under Coordinate. Click Apply and Export to export the coordinate results.</li>
		<li>Save lumen contour data corresponding to diastolic and systolic pressure conditions to .txt files for lumen circumference calculations.</li>
	</ol>
	</li>
	<li>Compare lumen circumferences calculated by the FE model (thin-slice layer structure-only model) at diastolic pressure condition with in vivo cine IVUS data (Cmin) and check if the relative error was &lt;1%. If the condition was met, then go to the next step, or otherwise update the material ratio k using the secant method and go to step 3.3 to re-run again17,18. 
	NOTE: In the first iterative, the Newton method was used to update the material ratio instead of the secant method.</li>
	<li>Compare lumen circumferences calculated by the FE model at systolic pressure condition with in vivo cine IVUS data (Cmax) and check if the relative error was &lt;1%. If yes, then stop the iterative procedure, or otherwise update the circumferential shrinkage rate S and redirect to step 3.4 to re-run again. 
	NOTE: In the first iterative, the Newton method was used to update the circumferential shrinkage rate instead of the secant method.</li>
	<li>Record optimal S and k values and calculate the corresponding material parameters of the Mooney-Rivlin material model.</li>
	<li>Plot the circumferential and axial stress-stretch ratio curves of the coronary vessel (Figure 7), which can be derived as follows: 
	<img alt="Static equilibrium equation; stress-strain relationship for materials; variables in deformation analysis." src="/files/ftp_upload/67573/67573eq04.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 4">    (4) 
	where σ represents Cauchy stress, λ represents stretch ratio, i = c, a represents the circumferential and axial directions.
	<ol>
		<li>To draw a material curve in a particular direction, fix the stretch ratio in the other direction to 1. Calculate effective Young's modulus in the circumferential and axial (labeled YMc and YMa, respectively) as the slope of the scale function of the material curve at the stretch ratio interval [1.0, 1.1] to reflect the general material stiffness of the coronary artery13: 
		<img alt="Elastic stress formula, sigma_i=YM(lambda_i-1); equation for mechanical analysis." src="/files/ftp_upload/67573/67573eq05.jpg" style="margin: auto;" data-altupdatedat="1788304762000" data-alt="Equation 5">    (5)</li>
	</ol>
	</li>
	<li>Extract plaque stress/strain distributions at any time and record node distribution and maximum stress values during systolic and diastolic phases (Figure 8).</li>
</ol>

 IVUS-Based Computational Method to Quantify Coronary Mechanical Properties

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

The authors declare no conflict of interest.

Critical steps in the protocol 
The most critical step in the finite element model-based updating approach lies in the iterative procedure. In the approach, the finite element model should accurately recover the coronary vessel motion on the vascular cross-section from in vivo cine IVUS images. To this purpose, minimizing the lumen circumference difference between the finite element model and in vivo images was adopted in this study to find the proper material properties. There were other key steps in the protocols, including image segmentation and processing, finite element model construction, and iterative schemes to quantify coronary vessel material properties. An improvement was to implement the automatic delineation of coronary plaque components to save time.
Significance of finite element modeling-based updating approach 
Classically, mechanical experiments such as uniaxial/biaxial tensile testing, indentation testing, and pressure inflation testing were conducted to quantify the mechanical behavior of the cardiovascular tissues ex vivo8,16. The clinical application of these approaches was limited due to the following reasons: 1) Ex vivo coronary arterial tissue samples are often not available. It is almost impossible to harvest the normal coronary samples in the clinical setting; 2) Material properties of arterial tissues might alter when taken out of living subjects; 3) it is not suitable for continuous monitoring in a patient-specific setting for personal management and precision medicine. Fortunately, FEMBUA provides another way to determine patient-specific tissue properties in vivo. The in vivo method could be easily modified to successfully apply to other biological tissues, such as aortic tissues19 and cardiac tissue20, and even non-biological materials like metal21. Prior studies have demonstrated that patient-specific in vivo tissue material properties had a significant influence on cardiovascular biomechanics, especially in strain calculation, compared to ex vivo material properties11. Therefore, patient-specific in vivo material properties are desirable for personalized treatment.
Even though in vivo and ex vivo methods are different, they could be integrated to inspire other approaches to quantify the mechanical properties of cardiovascular tissues. A hybrid approach that combines FEMBUA to match stress-strain data from biaxial/uniaxial experiments has been proposed for the aortic aneurysmal tissue22.
Potential clinical applications of Finite Element Modeling Based Updating Approach 
FEMBUA based on IVUS image is of great importance in cardiovascular materials science, medical image analysis, and personalized medical device design. This proposed method does not need to cause dramatic damage to the investigated coronary vessel wall as in a situation like taking the tissue out of the living human body, so it is suitable for continuous patient monitoring, such as investigating the effect of in vivo mechanical properties on patient prognosis, which cannot be done by a classical mechanical experimental approach. In addition, the proposed FEMBUA method plays a key role in the treatment optimization of coronary stents and the diagnosis of cardiovascular diseases. Finely analyzing the structure and lesion characteristics of blood vessel walls could guide the selection of stent size and location, thus maintaining the structural and mechanical stability within tissues to a higher degree, reducing complications, and improving patient prognosis23. In summary, the FEMBUA method has extensive and far-reaching application potential in the cardiovascular field.
Comparison of results from FEMBUA and ex vivo experimental approaches 
Comparison analysis between FEMBUA and classical ex vivo experimental approaches has been conducted to assess the accuracy and efficacy of the novel in vivo approach. For simplification, coronary tissue stiffness from the mechanical approach was compared to those from in vivo studies using FEMBUA, and tissue stiffness from both methods was generally in the same magnitude range11. The consistency was also confirmed in other vascular beds like aorta and carotid artery24. Of note that the limited number of ex vivo studies could influence the conclusions abovementioned, since the tissue stiffness variability across different individuals is also pronounced. Nevertheless, these conclusions suggested that FEMBUA is an accurate and effective approach to quantifying the material properties of arterial walls.
Validation and robustness of FEMBUA method were also investigated by performing both the in vivo FEMBUA method and ex vivo experiment on the same arterial tissue for validation purposes25,26,27. Liu et al. conducted in vivo and ex vivo experimental approaches on aortic tissues. Their findings revealed a close correlation between the material behavior curves generated by both methods, with an average mean absolute percentage error smaller than 5%25. Further, Cosentino et al. performed a similar comparative analysis based on a larger sample cohort (n=10), with similar conclusions obtained27. These findings collectively showed that FEMBUA could yield minimal discrepancies in biomechanical outcomes for identical samples. A further study also showed that the proposed method is reproducible and robust, as the study demonstrated that variations in mechanical properties had minimal impact on the simulated biomechanical outcomes in coronary artery models through finite element analysis27.
Modeling assumptions and limitations 
There are some assumptions involved in FEM for in vivo identification of the material properties of coronary arterial walls, which would impact the results from FEMBUA as described here. The axial shrinkage rate was assumed to be 95% since the actual axial shrinkage could not be obtained under in vivo conditions. The impact of axial stretch on material properties was investigated in prior studies, and the results indicated that smaller axial stretch led to greater slice shrinkage and softer material stiffness estimation11,28,29,30. Therefore, patient-specific axial stretch data should be used when available. The quality evaluation of finite element mesh involves several key indicators, among which the shape of the mesh element is the main evaluation criterion in this paper, which is directly related to the stability of the numerical solution and iterative convergence. In practical applications, due to the complexity of the problem, it is necessary to weigh the indicators to meet the analysis needs, and if necessary, the mesh quality can be optimized by designing new segmentation methods and selecting suitable mesh types31. Special attention should be paid to small lipids at key sites during actual segmentation32,33,34. In this implementation, small lipids in key sites (fibrous cap and shoulder) are retained, while lipids in other sites which have little influence on the stress/strain conditions. Blood pressure measured by arm cuff was used as a surrogate for on-site intracoronary pressure since invasive intracoronary pressure was not available for the patient. Active stress in the coronary artery was not considered since early evidence suggested that its contribution to the elastic properties of the living blood vessel was very small35. Residual stress information was not available and thus was not included in this model36,37. Structure-only models rather than more complex fluid-structure interaction models were used in the FEMBUA method, given it is more computationally efficient as iterative procedures typically need to solve the computational models several times to find the material properties constant.
In summary, image-based FEMBUA, different from other classical ex vivo experimental approaches, could be used to effectively determine the patient-specific material properties of coronary vessels in vivo. Since arterial tissue stiffness has already been employed in clinical settings as a risk factor for cardiovascular diseases, FEMBUA makes it possible to continuously monitor the coronary arterial mechanical properties under in vivo conditions and thus holds the potential in clinical applications for personalized treatment and precision medicine.

Coronary artery disease (CAD) is one of the leading causes of mortality and morbidity, accounting for more than 9.14 million deaths in 2019 globally1,2. The development of coronary artery diseases, such as atherosclerosis and stenosis, is often accompanied by alterations in mechanical forces and changes in vascular wall material properties3. The material properties of coronary arteries are not only the cornerstone to determine their mechanical response to the physiological loading but also the key elements to simulate the mechanical behavior of blood vessels, predict the development of atherosclerotic lesions, and evaluate the therapeutic effect of various medical devices4,5. Consequently, a profound understanding and accurate quantification of coronary material properties hold paramount value for early disease diagnosis, precision medicine, and prognosis assessment6.
Mechanical experiments of isolated coronary tissues, such as planar biaxial testing, indentation testing, inflation-extension, and uniaxial extension testing, are common approaches to quantify the mechanical properties of coronary vessel walls ex vivo7,8,9. From these approaches, coronary artery samples were obtained from patients or experimental animals. Mechanical testing was carried out to determine the strain responses of the vessel wall under different stress conditions, and then the material parameters were determined by fitting the experimental data10. Prior studies have shown that coronary properties are highly nonlinear and anisotropic11. Although ex vivo experiments can provide accurate material properties data, significant limitations also exist, which are as follows: First, the mechanical behavior of the sample after taking out from the living subjects would be different from that under in vivo conditions, which may affect the accuracy of testing results. Second, due to ethical and practical constraints, it is difficult to obtain a large collection of normal or pathological tissues of coronary arteries to perform the mechanical testing.
To overcome these limitations, researchers have explored novel techniques for in vivo, real-time, and patient-specific quantification of coronary material properties. Among them, the finite element model based updating approach (FEMBUA) based on medical image holds the promise to address these challenging issues. This approach makes use of advanced imaging techniques like intravascular ultrasound (IVUS) and virtual histology (VH)-IVUS to capture detailed coronary geometry, tissue compositions, and its movement12. By constructing 3D finite element (FE) models and incorporating patient-specific physiological blood pressure conditions, dynamic vascular behavior during cardiac cycles could be recovered by optimizing material parameters to match image data for rapid and accurate quantification of coronary material properties13. The advantages of the in vivo FE updating approach over ex vivo experiments include in vivo assessment without tissue excision, facilitating large-scale evaluations, and simulating vascular dynamics under complex conditions to aid the pathophysiology understanding of coronary diseases.
In this paper, key steps of the finite element model-based updating approach are introduced, which include a detailed segmentation and processing of cine IVUS and VH-IVUS image, reconstruction of computational thin-layer structure-only model, execution of the iterative scheme to search optimal material parameters for coronary arterial tissues. The aim of this protocol is to quantify the material properties of the coronary artery from a sample patient with CAD using the FEMBUA method as a demonstration, especially the illustration of step-by-step methods. We concluded by discussing the significance and other aspects of this in vivo method.
The selected participant is a 64-year-old female with no previous clinical history of coronary artery disease. This patient was diagnosed with coronary artery disease after having symptoms of chest pain. The coronary angiogram and IVUS scan were performed to confirm the diagnosis. A plaque lesion with 60% stenosis was found in the middle of the left anterior descending artery. After assessment, optimal medical therapy was adopted to treat the patient.

<table><tbody><tr><td>&amp;nbsp;Bee DICOM Viewer&amp;nbsp;</td><td>SinoUnion Healthcare Inc.</td><td>Version 3.5.1</td><td>A DICOM image reader software</td></tr><tr><td>ADINA&amp;nbsp;</td><td>Adina R &amp;amp; D</td><td>Version 9.0</td><td>Finite element solver</td></tr><tr><td>ImageJ&amp;nbsp;</td><td>National Institutes of Health</td><td></td><td>Segmented IVUS contours</td></tr><tr><td>MATLAB</td><td>MathWorks&amp;nbsp;</td><td>&amp;nbsp;Version R2018a</td><td>Commercial programming platform</td></tr><tr><td>Volcano s5 imaging system</td><td>Volcano Company</td><td></td><td>Intravascular ultrasound imaging system&amp;nbsp;</td></tr></tbody></table>

intravascular ultrasound image-based finite element modeling approach for quantifying in vivo mechanical properties of human coronary artery

Intravascular Ultrasound Image-Based Finite Element Modeling Approach for Quantifying In Vivo Mechanical Properties of Human Coronary Artery

Quantifying the mechanical properties of coronary arterial walls could provide meaningful information for the diagnosis, management, and treatment of coronary artery diseases. Since patient-specific coronary samples are not available for patients requiring continuous monitoring, direct experimental testing of vessel material properties becomes impossible. Current coronary models typically use material parameters from available literature, leading to significant mechanical stress/strain calculation errors. Here, we would introduce a finite element model-based updating approach (FEMBUA) to quantify patient-specific in vivo material properties of coronary arteries based on medical images. In vivo cine intravascular ultrasound (IVUS) and virtual histology (VH)-IVUS images of coronary arteries were acquired from a patient with coronary artery disease. Cine IVUS images showing the vascular movement over one cardiac cycle were segmented, and two IVUS frames with maximum and minimum lumen circumferences were selected to represent the coronary geometry under systolic and diastolic pressure conditions, respectively. VH-IVUS image was also segmented to obtain the vessel contours, and a layer thickness of 0.05 cm was added to the VH-IVUS contours to reconstruct the coronary geometry. A computational finite element model was created with an anisotropic Mooney-Rivlin material model used to describe the vessel's mechanical properties and pulsatile blood pressure conditions prescribed to the coronary luminal surface to make it contract and expand. Then, an iterative updating approach was employed to determine the material parameters of the anisotropic Mooney-Rivlin model by matching minimum and maximum lumen circumferences from the computational finite element model with those from cine IVUS images. This image-based finite element model-based updating approach could be successfully extended to determine the material properties of arterial walls in various vascular beds and holds the potential for risk assessment of cardiovascular diseases.

We describe in detail the FEMBUA method, which enables rapid plaque material and stress analysis of coronary plaques after real-time IVUS imaging and can determine the in vivo material properties and biomechanical results of plaques. The in vivo material parameters of the Mooney-Rivlin material model for this coronary vessel are provided in Table 1. The simulation results of the finite element model, including the stress/strain distributions in the coronary vessel, are plotted in Figure 8. A detailed analysis of the results is as follows.
It takes about 2 h to finish the entire process and to obtain in vivo material parameters for this coronary vessel, including cine IVUS and VH-IVUS image segmentation, construction of the computational FE model, and execution of the iterative procedure to quantify the material properties. This method can meet the requirements of rapid clinical examination with high timeliness.
For this coronary vessel, the two quantities were determined following the FEMBUA method as k = 1.125 and S = 98%. Other in vivo material parameters of the Mooney-Rivlin material model for this coronary vessel are provided in Table 1. Its material curves in both circumferential and axial directions are plotted in Figure 7. Both directions exhibited a J-shape material curve, indicating the nonlinear stress-stretch relationship with increased tissue stiffness when the stretch ratio was elevated. The straight lines were the scale functions to fit these nonlinear material curves, and their slopes were the effective Young's moduli. The YMc and YMa for this coronary vessel are 1055.41 and 1835.77 kPa, respectively.
The simulation results of the finite element model, including the stress/strain distributions in the coronary vessel, are plotted in Figure 8. The stress distribution corresponding to systolic pressure conditions is higher than that of diastolic pressure conditions. Also, local stress concentration could be observed in regions with thin vessel thickness or fibrous cap regions. Given that plaque rupture often occurs at the fibrous cap and the shoulder regions, the peak stress and strain values over these critical regions are reported in Table 2.
We performed the mesh analysis by gradually increasing the mesh density, and the mesh refinement would stop if a finer mesh could yield insignificant changes in solutions. The element/grid distribution and stress results for this model are shown in Supplementary Table 1. The table gives more comparative and detailed results, showing insignificant differences between models with different element numbers. The results show that, after refining the mesh density four times, the relative errors of maximum and minimum stress are within 5%. Since the FEMBUA method involves an iterative process, a relatively low mesh density was adopted to save computational costs. According to the results, it can be seen that the current mesh quantity (element number = 900) meets the requirements and has little influence on the stress results.
<img alt="Intravascular ultrasound, vessel cross-sections, lipid analysis, coronary artery imaging." class="xfigimg" src="/files/ftp_upload/67573/67573fig01.jpg" data-altupdatedat="1788304762000" data-alt="Figure 1"> 
Figure 1: Acquisition and contour extraction of Cine IVUS and VH-IVUS. (A) IVUS catheter withdrawal direction and catheter paused site (marked with a star) to obtain Cine IVUS images. (B) Cine IVUS with lumen circumferences of Cmin and Cmax and corresponding VH-IVUS image. (C) IVUS contours extraction. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="3D geometry diagram showing a structure with thickness in XYZ axes for spatial analysis." class="xfigimg" src="/files/ftp_upload/67573/67573fig02.jpg" data-altupdatedat="1788304762000" data-alt="Figure 2"> 
Figure 2: Reconstructed 3D coronary vessel geometry with a 0.05 cm thickness. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Lipid structure diagram with auxiliary points and labeled lipid regions for analysis." class="xfigimg" src="/files/ftp_upload/67573/67573fig03.jpg" data-altupdatedat="1788304762000" data-alt="Figure 3"> 
Figure 3: Finite element mesh generation process for thin-slice structure-only model. (A) Contour points on 2D plane. (B) Geometric lines for mesh generation. (C) Finite element mesh. (D) Generated different material groups. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Stress-strain analysis, arterial model, pressure graph; coronary stress, pressure-time data analysis." class="xfigimg" src="/files/ftp_upload/67573/67573fig04.jpg" data-altupdatedat="1788304762000" data-alt="Figure 4"> 
Figure 4: Computational thin-slice structure-only model of the coronary vessel. The model includes material property definition and boundary condition setting. (A) Set different material groups within the coronary artery. (B) Ex vivo coronary vessel materials curve. (C) Boundary condition loading. (D) Pressure curve. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Coronary vessel material properties flowchart; steps include FE simulation, iteration, secant method." class="xfigimg" src="/files/ftp_upload/67573/67573fig05.jpg" data-altupdatedat="1788304762000" data-alt="Figure 5"> 
Figure 5: Iterative scheme. The iterative scheme to quantify in vivo plaque material properties of coronary vessels based on IVUS images. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig05large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Finite element analysis, ADINA software, structural dynamics, mesh model, simulation interface." class="xfigimg" src="/files/ftp_upload/67573/67573fig06.jpg" data-altupdatedat="1788304762000" data-alt="Figure 6"> 
Figure 6: Use interface. The advanced user interface of ADINA software and relevant operations. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig06large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Stress vs. stretch ratio graph; Young's Modulus comparison; data analysis; material mechanics." class="xfigimg" src="/files/ftp_upload/67573/67573fig07.jpg" data-altupdatedat="1788304762000" data-alt="Figure 7"> 
Figure 7: Material curves. The material curves in circumferential and axial directions of the coronary vessel and their effective Young's moduli estimation. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig07large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Finite element analysis; stress distribution diagram; structural stress results; pressure levels comparison." class="xfigimg" src="/files/ftp_upload/67573/67573fig08.jpg" data-altupdatedat="1788304762000" data-alt="Figure 8"> 
Figure 8: Thin-slice structure-only model of the coronary vessel geometry and simulated results. The figure includes stress distribution cloud map at diastolic and systolic pressure conditions. (A) Nodes Displayed. (B) Diastolic Stress. (C) Systolic Stress. The color bar indicates stress value. <a href="https://www.jove.com/files/ftp_upload/67573/67573fig08large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td rowspan="2">Material parameters</td>
			<td>S(%)</td>
			<td>C1(kPa)</td>
			<td>C2(kPa)</td>
			<td>D1(kPa)</td>
			<td>K1(kPa)</td>
			<td>YMa(kPa)</td>
			<td>YMc(kPa)</td>
		</tr>
		<tr>
			<td>98%</td>
			<td>-147.76</td>
			<td>12.91</td>
			<td>708.67</td>
			<td>40.4</td>
			<td>1055.41</td>
			<td>1835.77</td>
		</tr>
	</tbody>
</table>
Table 1: Reference values for material parameters of the sample model.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td colspan="2">Time at Diastolic Pressure</td>
			<td colspan="2">Time at Systolic Pressure</td>
		</tr>
		<tr>
			<td>PWS (kPa)</td>
			<td>PWSn</td>
			<td>PWS (kPa)</td>
			<td>PWSn</td>
		</tr>
		<tr>
			<td>56.53</td>
			<td>0.05402</td>
			<td>64.51</td>
			<td>0.05402</td>
		</tr>
	</tbody>
</table>
Table 2: Peak stress/strain values over the critical regions from the computational model. Abbreviations: Peak wall stress = PWS; Peak wall strain = PWSn. Stress unit: kPa.
Supplementary Table 1: The comparison stress results of different grid quantities in the same model.&gt; <a href="https://www.jove.com/files/ftp_upload/67573/Supplementary%20Table%201_RE.xlsx" target="_blank">Please click here to download this File.</a>
Supplementary File 1: MATLAB codes. <a href="https://www.jove.com/files/ftp_upload/67573/Supplementary%20File%20S1.zip" target="_blank">Please click here to download this File.</a>

In vivo cine intravascular ultrasound images show the coronary cross-sectional movement corresponding to different pressure loading conditions. Based on a finite element model, an iterative scheme was employed to determine the patient-specific mechanical properties of coronary arteries in vivo by matching coronary motion from the computational model and medical images.

Quantifying the mechanical properties of coronary arterial walls could provide meaningful information for the diagnosis, management, and treatment of ...

intravascular-ultrasound-image-based-finite-element-modeling-approach

Research

JoVE Journal

Engineering

This method outlines the use of Quantum Micro-Computed Tomography (MicroCT) to assess cardiac morphology, function, perfusion, metabolism and viability with iodinated contrast agent in mice with experimentally-induced myocardial ischemia. The technique can be applied for non-destructive high-throughput longitudinal in vivo imaging of various animal models of human heart disease.

In Vivo Quantitative Assessment of Myocardial Structure, Function, Perfusion and Viability Using Cardiac Micro-computed Tomography

Here, we describe a correlative workflow for the excision, pressurization, fixation, and imaging of the murine pulmonary valve to determine the gross conformation and local extracellular matrix structures.

Surgery and Sample Processing for Correlative Imaging of the Murine Pulmonary Valve

There is a need to determine which atherosclerotic lesions will progress in the coronary vasculature to guide intervention before myocardial infarction occurs. This article outlines the biomechanical modeling of arteries from Optical Coherence Tomography using fluid-structure interaction techniques in a commercial finite element solver to help predict this progression.

Optical Coherence Tomography Based Biomechanical Fluid-Structure Interaction Analysis of Coronary Atherosclerosis Progression

This protocol demonstrates the workflow of the SILICOFCM platform for automatically generating a parametric model of the left ventricle from patient-specific ultrasound images by applying a multi-scale electromechanical model of the heart. This platform enables in silico clinical trials intended to reduce real clinical trials and maximize positive therapeutic outcomes.

In Silico Clinical Trials for Cardiovascular Disease

Here, we describe a protocol to use in vivo four-dimensional ultrasound imaging and ex vivo mass spectrometry imaging to assess biomechanical and biomolecular alterations in the murine cardiovascular system. This technique is applied to analyze cardiac remodeling in surgically induced myocardial infarction and vascular changes in aging animals.

Multimodal Study of Murine Cardiovascular Remodeling: Four-Dimensional Ultrasound and Mass Spectrometry Imaging

Phase-shift nanoemulsions (PSNE) can be vaporized using high intensity focused ultrasound to enhance localized heating and improve thermal ablation in tumors. In this report, the preparation of stable PSNE with a narrow size distribution is described. Furthermore, the impact of vaporized PSNE on ultrasound-mediated ablation is demonstrated in tissue-mimicking phantoms.

Synthesis of Phase-shift Nanoemulsions with Narrow Size Distributions for Acoustic Droplet Vaporization and Bubble-enhanced Ultrasound-mediated Ablation

Organic photovoltaic (OPV) materials are inherently inhomogeneous at the nanometer scale. Nanoscale inhomogeneity of OPV materials affects performance of photovoltaic devices. In this paper, we describe a protocol for quantitative measurements of electrical and mechanical properties of OPV materials with sub-100 nm resolution.

Concurrent Quantitative Conductivity and Mechanical Properties Measurements of Organic Photovoltaic Materials using AFM

Ultrasound velocimetry is used to study mixing by fluid flow in liquid metal electrodes. The focus of this manuscript is to illustrate the methods used for making precise, spatially-resolved ultrasound measurements while limiting oxidation and controlling and monitoring temperature, applied current, and the heater power being supplied.

Ultrasound Velocity Measurement in a Liquid Metal Electrode

The fabrication of high contrast gratings as the parallel spectrum splitting dispersive element in a concentrated photovoltaic system is demonstrated. Fabrication processes including nanoimprint lithography, TiO2 sputtering and reactive ion etching are described. Reflectance measurement results are used to characterize the optical performance.

Fabrication of High Contrast Gratings for the Spectrum Splitting Dispersive Element in a Concentrated Photovoltaic System

This study describes a protocol to evaluate the targeting accuracy in the focal plane of an ultrasound-guided high-intensity focused ultrasound phased-array system.

Evaluating Targeting Accuracy in the Focal Plane for an Ultrasound-guided High-intensity Focused Ultrasound Phased-array System

Discover finite element modeling using intravascular ultrasound images to quantify coronary artery mechanical properties in vivo. This computational approach enables patient-specific cardiovascular diagnostics.

Intravascular Ultrasound Image-Based Finite Element Modeling Approach for Quantifying In Vivo Mechanical Properties of Human Coronary Artery (Video) | JoVE Journal

Intravascular Ultrasound Image-Based Finite Element Modeling Approach for Quantifying In Vivo Mechanical Properties of Human Coronary Artery

Explore More Videos