Open AccessArticle

Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat

Min Wu

Li Huang

^2,*,

Xiaoyu Zhang

²,

Jianzhong Chen

² and

Yong Lv

School of Science, Wuhan University of Technology, Wuhan 430070, China

Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanics, Wuhan University of Technology, Wuhan 430070, China

Author to whom correspondence should be addressed.

Polymers 2020, 12(12), 3072; https://doi.org/10.3390/polym12123072

Submission received: 19 November 2020 / Revised: 16 December 2020 / Accepted: 16 December 2020 / Published: 21 December 2020

(This article belongs to the Special Issue Reinforced Polymer Composites II)

Download

Browse Figures

Graphical abstract
"> Figure 1
Schematic diagram of electrode layout. "> Figure 2
Front and back of the specimen: (a) polymer-matrix carbon mat attachment position on the front; (b) strain gauge attachment position on the back. "> Figure 3
Device diagram of loading: (a) loading along the 1st-direction; (b) loading along the 2nd-direction. "> Figure 4
Test curve of sample No. 1: (a) loading along the 1st-direction; (b) loading along the 2nd-direction. "> Figure 5
Test curve of sample No. 2: (a) loading along the 1st-direction; (b) loading along the 2nd-direction. "> Figure 6
Schematic diagram of the conductive network of the carbon fiber mat: (a) before loading; (b) after loading. "> Figure 7
Test curve of sample No. 3: (a) loading along the 1st-direction; (b) loading along the 2nd-direction. ">

Versions Notes

Abstract

Based on the piezoresistive effect, the piezoresistive constitutive relation of a carbon fiber mat under orthogonal strain was deduced. Considering the Poisson effect, the piezoresistive responses and measurement of the sensitivity factor of a polymer-matrix carbon fiber mat under bidirectional strain were studied by a two-times uniaxial tension loading method in different directions, which was pasted in the center area of a cruciform aluminum substrate. The relations between the resistance change rate and the orthogonal strains were established, the reasonability of which was confirmed by comparison with the experimental results. The results show that the longitudinal piezoresistive sensitivity factor C₁₁ is 21.55, and the lateral piezoresistive sensitivity factor C₁₂ is 24.15. Using these factors, the resistance change rate of another polymer-matrix carbon mat was predicted, which was made by the same technique, and the error between the predicted and the experimental results was 1.3% in the longitudinal direction and 6.1% in the lateral direction.

Keywords:

polymer-matrix; carbon fiber mat; orthogonal strain; piezoresistive sensitivity factor; piezoresistive effect

Graphical Abstract

1. Introduction

Carbon fiber materials are widely used in military, aerospace and other fields due to their high strength, light weight, fatigue resistance, corrosion resistance and other excellent properties. Since the piezoresistive responses of carbon fiber-reinforced plastics (CFRP) were discovered in 1980s, the piezoresistive responses characteristics of carbon fiber material made by different methods (e.g., carbon fiber bundle [1], single carbon fiber epoxy composites [2], continuous carbon fiber prepreg [3]), and under different loading modes, for instance, unidirectional cyclic loading [4,5], uniaxial tension [6] and unidirectional compression [7], have been studied by the majority of scholars. Utilizing the piezoresistive responses characteristics of carbon fiber materials, when carbon fiber material was embedded in concrete it could be used to detect internal defects and damage [8,9]; if it was laid on the surface of a concrete beam, called a carbon fiber smart layer, it could be used for surface crack detection or beam deformation monitoring [10,11]. In order to establish the mapping relationship between strain and resistance of a carbon fiber smart layer, the existing research mainly focuses on the behaviors of a carbon fiber monofilament [12], carbon fiber bundle [13] and polymer-matrix carbon fiber mat [14].

Regarding the piezoresistive constitutive research of a polymer-matrix carbon fiber mat, the joint effect of multi-directional strain has not been considered in existing reports. The relationship between resistance change (rate) and unidirectional stress or strain was established, and did not involve complex strain states. However, in the actual project, forced components are often in a complex stress state, and the piezoresistive constitutive derived from the unidirectional strain field is no longer unconvincing. Therefore, it is of great significance to study the piezoresistive response of a polymer-matrix carbon fiber mat under complex stress or strain, which lays a foundation for the application of a carbon fiber mat in the actual structure detection.

In this paper, based on piezoresistive theory, the constitutive relations of piezoresistive under bidirectional strain was derived. In order to measure the resistance response of a carbon fiber mat under orthogonal strain, and to avoid the errors caused by the differences between different samples, it was measured by two-times unidirectional tensile experiments in different directions, which were attached to a cruciform specimen of aluminum. According to the experimental results, the resistance sensitive factor was solved, and the constitutive relation of the polymer-matrix carbon mat was established and verified.

2. Theory of Piezoresistive Model

The change of the resistivity of solid material induced by the load is called the piezoresistive effect. For three-dimensional materials, the piezoresistive coefficient matrix is a square matrix of 6 × 6, and Equation (1) is the resistivity change rate caused by stresses [15].

[\begin{array}{l} {(\frac{Δ ρ}{ρ_{0}})}_{1} \\ {(\frac{Δ ρ}{ρ_{0}})}_{2} \\ {(\frac{Δ ρ}{ρ_{0}})}_{3} \\ {(\frac{Δ ρ}{ρ_{0}})}_{4} \\ {(\frac{Δ ρ}{ρ_{0}})}_{5} \\ {(\frac{Δ ρ}{ρ_{0}})}_{6} \end{array}] = [\begin{array}{l} π_{11} π_{12} π_{13} π_{14} π_{15} π_{16} \\ π_{21} π_{22} π_{23} π_{24} π_{25} π_{26} \\ π_{31} π_{32} π_{33} π_{34} π_{35} π_{36} \\ π_{41} π_{42} π_{43} π_{44} π_{45} π_{46} \\ π_{51} π_{52} π_{53} π_{34} π_{55} π_{56} \\ π_{61} π_{62} π_{63} π_{64} π_{65} π_{66} \end{array}] [\begin{array}{l} σ_{1} \\ σ_{2} \\ σ_{3} \\ σ_{4} \\ σ_{5} \\ σ_{6} \end{array}]

(1)

where,

σ_{1}

σ_{2}

σ_{3}

σ_{4}

σ_{5}

σ_{6}

are the stress components describing the stress state of solid material;

σ_{1}

σ_{2}

σ_{3}

are independent normal stresses and

σ_{4}

σ_{5}

σ_{6}

are independent shear stresses.

When the three-dimensional solid material is an isotropic material,

π_{11} = π_{22} = π_{33}

;

π_{12} = π_{13} = π_{23} = π_{21} = π_{31} = π_{32}

;

π_{44} = π_{55} = π_{66}

. meanwhile, the piezoresistive effect has the following characteristics:

(1): It is impossible for shear stress to produce positive piezoresistive effect,

$π_{14} = π_{15} = π_{16} = π_{24} = π_{25} = π_{26} = π_{34} = π_{35} = π_{36} = 0 .$
(2): Normal stress cannot produce shear piezoresistive effect,

$π_{41} = π_{51} = π_{61} = π_{42} = π_{52} = π_{62} = π_{43} = π_{53} = π_{63} = 0 .$
(3): Shear stress can only produce piezoresistive effect in the plane which acts on,

$π_{45} = π_{46} = π_{34} = π_{56} = π_{64} = π_{65} = 0 .$

Furthermore, for the plane stress state,

σ_{3} = 0

σ_{5} = 0

σ_{6} = 0

, thus the piezoresistive equation can be expressed by Equation (2).

\{\begin{cases} {(\frac{Δ ρ}{ρ_{0}})}_{1} = π_{11} σ_{1} + π_{12} σ_{2} \\ {(\frac{Δ ρ}{ρ_{0}})}_{2} = π_{12} σ_{1} + π_{11} σ_{2} \\ {(\frac{Δ ρ}{ρ_{0}})}_{4} = π_{44} σ_{4} \end{cases}

(2)

Combined with Hooke’s law in the plane stress state,

[\begin{matrix} σ_{1} \\ σ_{2} \\ σ_{4} \end{matrix}] = \frac{E}{(1 - μ^{2})} [\begin{matrix} 1 & μ & 0 \\ μ & 1 & 0 \\ 0 & 0 & \frac{1 - μ}{2} \end{matrix}] [\begin{matrix} ε_{1} \\ ε_{2} \\ ε_{4} \end{matrix}]

(3)

Equation (2) can be expressed as the rate of resistivity change caused by strains, as shown in Equation (4).

\{\begin{cases} {(\frac{Δ ρ}{ρ_{0}})}_{1} = C_{11} ε_{1} + C_{12} ε_{2} \\ {(\frac{Δ ρ}{ρ_{0}})}_{2} = C_{12} ε_{1} + C_{11} ε_{2} \\ {(\frac{Δ ρ}{ρ_{0}})}_{4} = C_{44} ε_{4} \end{cases}

(4)

where,

C_{11} = \frac{E (π_{11} + μ π_{12})}{1 - μ^{2}}

C_{12} = \frac{E (μ π_{11} + π_{12})}{1 - μ^{2}}

C_{44} = \frac{E π_{44}}{21 + μ}

respectively represent the longitudinal piezoresistive sensitivity factor, the lateral piezoresistive sensitivity factor, and the shear piezoresistive sensitivity factor; µ is Poisson’s ratio and E is the elastic modulus.

3. Experimental Research

3.1. Raw Material

The PAN-based carbon fiber mat of 0.05 mm thickness was manufactured by Pengyuan fiberglass Products Co. Ltd. The raw material parameters of the mat are shown in Table 1. Conductive silver paint for electrodes adhering to the carbon fiber mat was manufactured by SPI Corporation. Other materials: 0.1 mm diameter copper core wire.

3.2. Specimens Preparation

The carbon fiber mat was cut to the size of 50 mm × 20 mm. In order to measure the resistance of the carbon fiber mat, the four-electrode method was used to eliminate the contact resistance between the electrodes and the carbon fiber mat. Four copper wires were bonded to the carbon fiber mat at predetermined positions with conductive silver glue to ensure that the distance between the actuation electrodes for current excitation was 40 mm, and the distance between the potential electrodes for voltage acquisition was 20 mm, as shown in Figure 1. It was assumed that the length direction of the carbon fiber mat was the 1st-direction, and the width direction was the 2nd-direction.

The process of polymer-matrix carbon fiber mat was as follows: the E-44 epoxy resin (product of Yueyang Petrochemical Complex), and 2-ethyl-4-methylimidazole were mixed at a mass ratio of 100:10, and then 4% glycidyl phenyl ether was added for dilution. The 3 mm-thickness insulated cruciform aluminum substrate with width of 20 mm was brushed the diluted epoxy resin. The carbon fiber mat with electrodes was pasted on the central position of the aluminum substrate, ensuring that the distance between the potential electrodes coincides with the width of the aluminum substrate. Then the diluted epoxy resin was brushed again on the upper surface of the mat. Finally, the specimen was covered with peel ply, and placed into a vacuum bag for pressure curing. After 7 days’ solidification at room temperature, the specimen was demoulded. Two perpendicular strain gauges were pasted on the back of the aluminum substrate, as illustrated in Figure 2.

3.3. The Experimental Process

Two-times uniaxial loading in different directions was carried out on an Instron 5882 material test machine. The loading speed was 2000 N/min and the maximum load was 5000 N under force control, as shown in Figure 3. During loading, Keithley 6221 was used to provide a constant 1 mA direct current (DC) to the actuation electrodes, and the voltage between the potential electrodes was collected by a Keithley 2700 with the acquisition frequency of 1 Hz. The resistance of the carbon fiber mat was calculated by Ohm’s law, that is, the voltage divided by current.

4. Results and Discussion

4.1. Experimental Results

Two groups of tests were conducted on the same batch of samples. The relationship between the resistance response of the polymer-matrix carbon fiber mat and the two orthogonal strains is as shown in Figure 4 and Figure 5, where the resistance change rate is the ratio of the resistance change to the original resistance.

The experimental results show that when the loading direction is consistent with the 1st-direction of the carbon fiber mat, the resistance change rate is greatly affected by the 1st-direction strain; when the loading direction is consistent with the 2nd-direction, due to the Poisson effect, strain in the 1st-direction is negative, and the resistance change rate still increases, which indicates the positive piezoresistive effect. This phenomenon is caused by the change of the gap between the fibers in the carbon fiber mat. During loading, the force caused the fibers to rotate or translate and made some overlapped fibers fall off, resulting in the change of the number of conductive channels. As shown in Figure 6, when the loading direction is consistent with the 2nd-direction, some fibers in the carbon fiber mat change before and after loading. The two conductive paths before loading form a parallel connection, and the conductive path is reduced to one after loading. According to the calculation theory of series and parallel resistance, the resistance will increase. This phenomenon demonstrates that the sensing characteristics of the carbon fiber mat are determined by strains in two directions, and the Poisson effect has a non-negligible influence on the sensing characteristics of the carbon fiber mat.

4.2. Solution of Piezoresistivity Sensitivity Factor

Because the carbon fiber mat is very thin and tightly pasted on the surface of the aluminum substrate, it can be considered that the strain of the two materials is the same. For the carbon fiber mat, the resistance expression is shown in Equation (5):

R = ρ \frac{L}{W}

(5)

where,

ρ

is the surface resistivity of the carbon fiber mat, L is the length and W is the width of the carbon fiber mat between the potential electrodes.

From Equation (5), the resistance change rate in the piezoresistivity can be derived, as shown in Equation (6):

\frac{Δ R}{R_{0}} = \frac{Δ ρ}{ρ_{0}} + \frac{Δ L}{L} - \frac{Δ W}{W} = \frac{Δ ρ}{ρ_{0}} + ε_{1} - ε_{2}

(6)

Combining Equation (6) with the first two formulas of Equation (4), the relationship between the change rate of resistance and strain of the carbon fiber mat can be obtained, as shown in Equation (7):

\{\begin{cases} {(\frac{Δ R}{R_{0}})}_{1} = (C_{11} + 1) ε_{1} + (C_{12} - 1) ε_{2} \\ {(\frac{Δ R}{R_{0}})}_{2} = (C_{12} - 1) ε_{1} + (C_{11} + 1) ε_{2} \end{cases}

(7)

According to the experimental results in Figure 4 and Figure 5, the change rate of resistance of the carbon fiber mat corresponding to the maximum strain is substituted into Equation (7), respectively, and the piezoresistive sensitivity factor of the carbon fiber mat under the action of bidirectional strain could be calculated, as shown in Table 2.

The piezoresistive constitutive of the polymer-matrix carbon fiber mat under the action of bidirectional strain can be obtained, as shown in Equation (8):

\{\begin{cases} {(\frac{Δ R}{R_{0}})}_{1} = 22 . 55 ε_{1} + 23 . 15 ε_{2} \\ {(\frac{Δ R}{R_{0}})}_{2} = 23 . 15 ε_{1} + 22 . 55 ε_{2} \end{cases}

(8)

4.3. Piezoresistive Constitutive Verification

In order to verify the accuracy of the piezoresistive constitutive relation, a similar force resistance experiment was performed on the No. 3 sample. The experimental results are shown in Figure 7:

The resistance change rate of the No. 3 polymer-matrix carbon fiber mat sample is predicted by Equation (8), and the specific results are shown in Table 3. Obviously, the predicted resistance change rate

{(\frac{Δ R}{R})}_{pre}

is very close to the experimental resistance change rate

{(\frac{Δ R}{R})}_{\exp}

. The results prove that this constitutive relation can better reflect the resistance response law of the carbon fiber mat under orthogonal strains.

5. Conclusions

(1): Based on piezoresistance theory, the constitutive model of the polymer-matrix carbon fiber mat under plane stress state was deduced in this paper, and the response law of its resistance under bidirectional strains was investigated. Combined with the experimental results, the piezoresistance sensitivity factor of the polymer-matrix carbon fiber mat was obtained, and the piezoresistance constitutive relationship in the orthogonal strain state was established. The rationality of the piezoresistance relationship was verified by comparing it with the experimental results.
(2): In this study, the shear piezoresistance sensitivity factor of the carbon fiber mat could not be measured. Therefore, further research is needed in the future.

Author Contributions

Conceptualization, M.W. and L.H.; validation, L.H., J.C. and X.Z.; investigation, M.W.; writing—original draft preparation, M.W.; writing—review and editing, M.W., and L.H.; supervision, X.Z.; project administration, J.C.; funding acquisition, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Composite Product Design and Development and Technical Services, grant number 20181h03461.

Acknowledgments

The first author would like to extend deep gratitude to all those who offered practical, cordial and selfless support in writing this thesis.

Conflicts of Interest

The authors declare no conflict of interest.

References

Kalashnyk, N.; Faulques, E.; Schjødt-Thomsen, J.; Jensen, L.R.; Rauhe, J.C.M.; Pyrz, R. Strain Sensing in Single Carbon Fiber Epoxy Composites by Simultaneous in-situ Raman and Piezoresistance Measurements. Carbon 2016, 109, 124–130. [Google Scholar] [CrossRef]
Kalashnyk, N.; Faulques, E.; Schjødt-Thomsen, J.; Jensen, L.R.; Rauhe, J.C.M.; Pyrz, R. Monitoring Self-sensing Damage of Multiple Carbon Fiber Composites Using Piezoresistivity. Synth. Met. 2017, 224, 56–62. [Google Scholar] [CrossRef]
Tang, J.; Zhu, S.R.; Li, Z.Q. Experimental Research on Piezoresistivity of Continuous Carbon Fiber Prepreg Sensitive Layer. J. Wuhan Univ. Technol. 2010, 32, 16–19. [Google Scholar]
Vavouliotis, A.; Paipetis, A.; Kostopoulos, V. On the fatigue life prediction of CFRP laminates using the Electrical Resistance Change method. Compos. Sci. Technol. 2011, 71, 630–642. [Google Scholar] [CrossRef] [Green Version]
Todoroki, A.; Haruyama, D.; Mizutani, Y.; Suzuki, Y.; Yasuoka, T. Electrical Resistance Change of Carbon/Epoxy Composite Laminates under Cyclic Loading under Damage Initiation Limit. J. Compos. Mater. 2014, 4, 22–31. [Google Scholar] [CrossRef] [Green Version]
Chen, J.F.; Morozov, E.V.; Shankar, K. Progressive failure analysis of perforated aluminium/CFRP fibre metal laminates using a combined elastoplastic damage model and including delamination effects. Compos. Struct. 2014, 114, 64–79. [Google Scholar] [CrossRef]
Takahashi, K.; Todoroki, A.; Matsuzaki, R. Simultaneous Measurement of Multiple Electrical Resistance Changes with Strain of CFRP. J. Solid Mech. Mater. Eng. 2010, 4, 557–567. [Google Scholar] [CrossRef]
Bontea, D.M.; Chung, D.D.L.; Lee, G.C. Damage in carbon fiber-reinforced concrete, monitored by electrical resistance measurement. Cem. Concr. Res. 2000, 30, 651–659. [Google Scholar] [CrossRef]
Zheng, L.X. Research on Electric Conductivity and Sensing Characteristic of Smart Concrete. Ph.D. Thesis, Wuhan University of Technology, Wuhan, China, 2011. [Google Scholar]
Wang, L.; Deng, Y.S.; Duan, B.Z. Crack Detection for Composite Plate under Cyclic Loading by Carbon Fiber Smart Layers. Technol. Eng. Sci. 2018, 18, 295–299. [Google Scholar]
Huang, J.J.; Liu, R.G.; Xu, Z.H.; Xie, G.H. Study on the Engineering Application of the Sensor Based on the Piezo-resistance Effect of CFRP. Fiber Reinf. Plast. Compos. 2017, 9, 46–51. [Google Scholar]
Liang, G.; Liu, R.G.; Xu, Z.H.; Xie, G.H.; Huang, J.J. Experimental Study on the Sensing Characteristics of CFRP Tendons Based on the Piezoresistivity. Fiber Reinf. Plast. Compos. 2017, 2, 60–63. [Google Scholar]
Zheng, H.S.; An, W.; Wu, J.W.; Zhao, Z. Piezoresistivity of polymer-matrix carbon fiber filament in plane stress state. Mater. Res. Express. 2019, 6, 085602. [Google Scholar] [CrossRef]
Li, J. Piezoresistive Response of Polymer Matrix Carbon Fiber Smart Layer in Thin Plate Bending. Guangzhou Chem. Ind. Technol. 2014, 42, 81–82. [Google Scholar]
Zhang, X.Y. Characteristic Research of Carbon Fiber Smart Layer and Its Field Monitoring. Ph.D. Thesis, Wuhan University of Technology, Wuhan, China, 2007. [Google Scholar]

Figure 1. Schematic diagram of electrode layout.

Figure 2. Front and back of the specimen: (a) polymer-matrix carbon mat attachment position on the front; (b) strain gauge attachment position on the back.

Figure 3. Device diagram of loading: (a) loading along the 1st-direction; (b) loading along the 2nd-direction.

Figure 4. Test curve of sample No. 1: (a) loading along the 1st-direction; (b) loading along the 2nd-direction.

Figure 5. Test curve of sample No. 2: (a) loading along the 1st-direction; (b) loading along the 2nd-direction.

Figure 6. Schematic diagram of the conductive network of the carbon fiber mat: (a) before loading; (b) after loading.

Figure 7. Test curve of sample No. 3: (a) loading along the 1st-direction; (b) loading along the 2nd-direction.

Table 1. Raw material parameters of PAN-based carbon fiber mat.

Name	Specification (g/m²)	Organic Matter Content (%)	Fiber Diameter (μm)	Water Content (%)
Parameters	10	6–12	7.02 ± 2	≤0.5

Table 2. Resistance sensitivity factor of carbon fiber mat under bidirectional strain.

Groups	Longitudinal Piezoresistive Sensitivity Factor/C₁₁	Lateral Piezoresistive Sensitivity Factor/C₁₂
No. 1	23.6	23.3
No. 2	19.5	25.0
Average value	21.55	24.15

Table 3. Predicted and experimental values of resistance change rate of polymer-matrix carbon mat.

Loading Direction	ε₁ × (10⁻⁶)	ε₂ × (10⁻⁶)	${(\frac{Δ R}{R})}_{\exp}$	${(\frac{Δ R}{R})}_{pre}$	Error
1-direction	1186.1	−618.1	0.0126036	0.0124375	1.3%
2-direction	−593.5	1071.7	0.0107288	0.0114264	6.1%

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Wu, M.; Huang, L.; Zhang, X.; Chen, J.; Lv, Y. Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat. Polymers 2020, 12, 3072. https://doi.org/10.3390/polym12123072

AMA Style

Wu M, Huang L, Zhang X, Chen J, Lv Y. Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat. Polymers. 2020; 12(12):3072. https://doi.org/10.3390/polym12123072

Chicago/Turabian Style

Wu, Min, Li Huang, Xiaoyu Zhang, Jianzhong Chen, and Yong Lv. 2020. "Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat" Polymers 12, no. 12: 3072. https://doi.org/10.3390/polym12123072

APA Style

Wu, M., Huang, L., Zhang, X., Chen, J., & Lv, Y. (2020). Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat. Polymers, 12(12), 3072. https://doi.org/10.3390/polym12123072

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Two-Dimensional Piezoresistive Response and Measurement of Sensitivity Factor of Polymer-Matrix Carbon Fiber Mat

Abstract

1. Introduction

2. Theory of Piezoresistive Model

3. Experimental Research

3.1. Raw Material

3.2. Specimens Preparation

3.3. The Experimental Process

4. Results and Discussion

4.1. Experimental Results

4.2. Solution of Piezoresistivity Sensitivity Factor

4.3. Piezoresistive Constitutive Verification

5. Conclusions

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI