Strain Rate and Anisotropic Microstructure Dependent Mechanical Behaviors of Silkworm Cocoon Shells
Strain Rate and Anisotropic Microstructure Dependent Mechanical Behaviors of Silkworm Cocoon Shells
Jun Xu 0 1
Wen Zhang 0 1
Xiang Gao 0 1
Wanlin Meng 0 1
Juan Guan 0 1
0 1 Department of Automotive Engineering, School of Transportation Science and Engineering, Beihang University , Beijing , China , 2 Advanced Vehicle Research Center, Beihang University , Beijing , China , 3 Beijing Key Laboratory for High-efficient Power Transmission and System Control of New Energy Resource Vehicle, Beihang University , Beijing , China , 4 Department of Transportation Engineering, School of Transportation Science and Engineering, Beihang University , Beijing , China , 5 Department of Civil Engineering, School of Transportation Science and Engineering,Beihang University , Beijing , China , 6 School of Material Science and Engineering, Beihang University , Beijing , China
1 Editor: Jie Zheng, University of Akron, UNITED STATES
Silkworm cocoons are multi-layered composite structures comprised of high strength silk fiber and sericin, and their mechanical properties have been naturally selected to protect pupas during metamorphosis from various types of external attacks. The present study attempts to gain a comprehensive understanding of the mechanical properties of cocoon shell materials from wild silkworm species Antheraea pernyi under dynamic loading rates. Five dynamic strain rates from 0.00625 s-1 to 12.5 s-1 are tested to show the strain rate sensitivity of the cocoon shell material. In the meantime, the anisotropy of the cocoon shell is considered and the cocoon shell specimens are cut along 0°, 45° and 90° orientation to the short axis of cocoons. Typical mechanical properties including Young's modulus, yield strength, ultimate strength and ultimate strain are extracted and analyzed from the stressstrain curves. Furthermore, the fracture morphologies of the cocoon shell specimens are observed under scanning electron microscopy to help understand the relationship between the mechanical properties and the microstructures of the cocoon material. A discussion on the dynamic strain rate effect on the mechanical properties of cocoon shell material is followed by fitting our experimental results to two previous models, and the effect could be well explained. We also compare natural and dried cocoon materials for the dynamic strain rate effect and interestingly the dried cocoon shells show better overall mechanical properties. This study provides a different perspective on the mechanical properties of cocoon material as a composite material, and provides some insight for bio-inspired engineering materials.
Data Availability Statement; All relevant data are within the paper
Funding: This work is financially supported by
Fundamental Research Funds for the Central
Universities, Beihang University (JX), Startup fund for
“Zhuoyue 100” titled professorship, Beihang
University (JX), "FengRu" Research Funds of School
of Transportation Science and Engineering, Beihang
University (WZ) and Research Funds for National
Student Innovation Training Program (WZ). JG would
like to thank Beihang University for the associate
professorship of “Zhuoyue 100” program and the
startup fund. The funders had no role in study design,
Silkworms and a range of arthropods and insects including spiders and bees have evolved to
spin protein fibers that nicely combine high strength and toughness for multiple functions.
data collection and analysis, decision to publish, or
preparation of the manuscript.
Natural silk fibers and cocoons are “ideal” defense and survival tools for the animals [
on the other hand the silk fibers from domesticated Bombyx mori silkworms have long been
great attractions to the textile industry [
Cocoons, consisting of twining silk fibers of a continuous length of approximately 1000 m
and bonded together by sericin glue, are a typical natural fiber composite with hierarchical
1, 9, 10
]. The widely studied silk fibers from the Bombyx mori silkworm cocoons are used
not only as textile materials, but also as multifunctional materials [
]. To realize the
biological function of cocoons as a shield for the pupas, the fiber component in cocoons possessing a
semicrystalline structure offers balanced stiffness and extensibility, and the binder component
sericin which is an amorphous protein helps to maintain the network structure of cocoons.
Up to date, much research focus have been put on the mechanical behaviors of the silk fibers
unravelled from cocoons, their structure-property relations and the biodegradability. Some
researchers have focused on studying different kinds of cocoons to comprehensively
understand the interplay among the multilayer structures of cocoons, the mechanical properties of
cocoon fibers, and the “optimized” mechanical properties of cocoons [
9, 13, 14
]. More recently,
cocoons with characteristic mechanical properties resulting from its constituent silk fiber and
the corresponding multi-layered structure have been reported to be ideal candidates for
bioinspired engineering design models [
]. Another study showed a strong dependency of the
mechanical properties of different cocoon layers on the microstructures regardless of the same
composition among these layers through SEM (scanning electron microscopy) images and
mechanical testing . The typical mechanical properties, such as Young’s modulus, ultimate
strength and thermo-mechanical parameters, were experimentally measured to be different
from innermost pelade to the outermost floss in the thickness direction of Bombyx mori
]. In particular, an interesting step was made by Blossmen and et al. [
reported that Bombyx mori cocoons were proven to be able to conserve water for the pupa and
acting as a humidity trap while offering mechanical protection.
Nevertheless, the most important biological role for cocoon is to protect the pupa from
predators’ attacks (e.g. birds’ pecking) which in essence is an engineering problem of impact
and dynamic loading events. Previous studies draw major attention on the quasi-static
mechanical behavior of the cocoon material while the mechanical behaviors under dynamic
conditions with both inertia and strain rate effect need to be evaluated. Mortimer and Drodge
et al [
] studied the bullet impact effect on silk fibers and Chen et al [
] reported a
quasistatic “impact” experiment and analysis on cocoons. However, a direct measurement of cocoon
shells under dynamic loading conditions between quasi-static and impact rates is still lacking.
Moreover, more research is required to study the structure-property relationships of the wild
cocoon instead of the domesticated ones from a biomimetic point of view. Thus, in the present
study, a focused experiment investigation on the mechanical properties of cocoon shells from
wild silkworm species Antheraea pernyi of three different orientations under various dynamic
strain rates is conducted. Furthermore, the microstructural changes of cocoon shells during
deformation are observed using scanning electron microscopy to evaluate the effect of
microstructure on mechanical behaviors of cocoons. Finally, simple constitutive models are
referenced to provide fundamental understanding in the mechanical properties of cocoon shells
under dynamic loading conditions.
The Chinese tussah silkworm Antheraea pernyi cocoon we chose for the present study is one of
the most common wild cocoons in northern China and no endangered species were involved
2 / 16
in this work. All the Antheraea pernyi cocoons used in this study were purchased from a
professional silk rearing farmer named Desheng Song in Xianrendong Town, Dalian City,
Liaoning Province, China. To keep the consistency of all the cocoon samples, they were stored
under the same environmental conditions in the laboratory (a temperature below 15°C and
low humidity to prevent the live pupas from hatching). In order to study the effect of moisture
and the processing of drying which is used in silk industries on the mechanical properties of
cocoon shells, approximately half of the overall cocoon shell specimens were dried at a
temperature of 120°C in a laboratory oven for 48 hours under full ventilation. Cocoon shell specimens
from cocoons as received with live pupas inside serve as a comparison group of the natural
untreated ones in material testing.
A PD-153 (standard) vernier caliper (Pro’sKit) and a micrometer were used to measure the
geometry of the cocoons. For preparing cocoon shell specimens for tensile testing, cocoon shells
were first carefully cut open along its long axis to remove the pupas inside to prevent
contamination from the body fluid of the pupa. Then as Fig 1 shows, dog-bone shaped specimens were cut
from three directions, i.e. 0°, 45°, 90° respectively to the short axis of the cocoon shell and the
middle section of the dog-bone specimens were 3 mm in width and 16 mm in length. The
thicknesses of the samples were measured through a micrometer, which varied from 0.3 mm to 0.6
mm. Note that the possible delamination of layers and damage were carefully controlled during
sample preparation such that its influence on the mechanical behaviors of cocoon material was
considered trivial. Considering biological materials are highly variable in both mechanical
properties and morphologies, the effect of individual cocoon structural difference is ruled out by
Fig 1. Dog-bone shaped specimens cut from the (a) 0°, (b) 45°, (c) 90° direction of the same contact cocoon and (d) The size of designed dog-bone shaped
3 / 16
comparing only the specimens of three directions from the same cocoons. Each dog-bone shaped
specimen was marked and placed in plastic bags before and after the experiment.
Quasi-static uniaxial tensile tests were performed at room temperature of 14.8°C and relative
humidity of 22% on a MTS materials testing machine with a load cell of 0–200 N with a precision
grade of 0.5% to accurately record the deformation extension and force of the cocoon specimens.
In the meantime, dynamic-rate uniaxial tensile tests were performed utilizing a computer
controlled Instron E10000 test machine and the load cell was 0–250 N with a precision of 0.0125 N.
The dynamic loadings ranging from 10−4 m/s to 0.2 m/s, converted to nominal strain rate of
0.00625 s-1 to 12.5 s-1, with the consideration of the cocoon anisotropic morphology were tested.
The nominal stress and strain was calculated from the load force and cross-head displacement
divided by the cross-sectional area and the effective length of the middle section of the dog-bone
specimen. The thickness and width were measured at least at two different positions of each
dogbone specimen for its averaged cross-sectional area. Eight repeated tests were conducted for each
dynamic rate and each anisotropic specimen. The properties were analyzed and an average and
standard deviation for each property was calculated. The fractured specimens were recycled for
microstructural observations using a SEM (JSM 6010).
The general profiles of mechanical properties of cocoon materials in three different directions are
highly reproducible from a number of repeated material tests (n = 8). To represent the
characteristic stress-strain curve of the cocoon shell specimen without missing the discussion on
variability, the representative stress-strain curves are selected to match the averaged properties of each
group of stress-strain curves. For example, all the stress-strain curves (representative curve
highlighted in red) of dried dog-bone shaped specimens in the 0° direction under quasi-static
condition are presented in Fig 2A. Fig 2B illustrates the different stages of the stress-strain
behavior of the cocoon shell specimen using a representative stress-strain curve. It appears that the
trend of the stress-strain curve of cocoon materials is quite similar to that of steels [
although it is a type of fiber reinforced polymer material, where four stages can be clearly
observed, including elastic region (oa), yielding (bc), plastic stage (cd) and fracture (de).
For mechanical property analysis average values and their standard deviation of typical
properties, i.e. Young’s modulus E, yield strength σy, ultimate strength σu and ultimate tensile strain εu, are
derived from the stress-strain curves and listed in Table 1 (natural dog-bone specimens) and Table 2
(dried dog-bone specimens). Following the trend of the stress-strain curve, firstly, we see the stress
linearly increases with strain and the slope decreases gradually at a yield point of approximately
50.9 MPa. Then cocoon shell specimen enters its “plastic” stage, where a post-yield modulus can be
obtained. Finally, when stress reaches its maximum value with an ultimate strain higher than 20%,
the curve starts to fall rapidly when the cocoon material fails. Evidently, some of the cocoon
specimens in the experiment showed progressive and multi-stage failure behavior as shown in Fig 2A.
The microstructure of fractured specimen in the insert of Fig 2B suggests that there is no necking
but individual fiber break-up along an irregular fracture path through the width of the specimen.
4.1 Anisotropy in cocoon shell properties
Material tensile tests under five different nominal strain rates ranging from 0.00625 s-1 to 12.5
s-1 were performed for both dried and natural dog-bone shaped specimens along the 0°, 45°,
4 / 16
Fig 2. Representative stress-strain curves (a) stress-strain curves of dried dog-bone shaped specimens in
the 0-degree direction under quasi-static condition (b) different stages in typical stress-strain curve of dried
dog-bone shaped specimen in the 0-degree direction under a nominal strain rate of 0.00625 s-1.
5 / 16
90° directions from cocoons. Some representative stress-strain curves for each experimental
condition and each group of anisotropic specimens are shown in Fig 3. The SEM
microstructural images of different strain rate and orientation are inserted for the natural cocoon shell
specimens in Fig 4. It is found that all the analyzed mechanical properties including E and σy,
σu and εu are the highest in the 0° direction (along the minor axis of the ellipsoid shaped
cocoon), followed by those in 45° direction and samples and in 90° direction. This interesting
phenomenon of anisotropic morphology manifestation on the mechanical behaviors of
cocoons cannot be accidental, but may have essential meanings for cocoons’ biological
functions to protect pupae. As shown in Fig 4, it appears that the fracture of the silk fibers in the
specimens in the 0° direction concentrates more along the tensile direction, while along the
other two directions fibers tend to distribute in a more random way. Actually, previous studies
Fig 3. Representative stress-strain curves of dog-bone shaped specimens along different directions with a nominal strain rate of (a) 0.00625 s-1 (b) 0.0625
s1 (c) 0.625 s-1 (d) 6.25 s-1 (e) 12.5 s-1.
have shown that the Antheraea pernyi silk fiber have powerful mechanical properties with a
tensile strength of 0.3–1.0 GPa and a Young’s modulus of approximately 10 GPa, which is
7 / 16
Fig 4. SEM graphs of natural dog-bone shaped specimens along different directions with a nominal strain
rate of (a) 0.00625 s-1 (b) 0.0625 s-1 (c) 0.625 s-1 (d) 6.25 s-1 (e) 12.5 s-1.
8 / 16
about ten times of that of cocoon shells [
], indicating that the strength of silk fiber
contributes greatly to the strength of the cocoon shell. Therefore, due to the preferential
orientation of the silk fibers along the 0° direction there is a tendency that cocoon specimen in the 0°
direction show enhanced mechanical properties compared to two other directions. In other
words, the preferential orientation of silk fibers in cocoon shells is attributed to the differences
in their directional mechanical properties. This designing principle could be used for spatial
structures with optimized mechanical integrity given the same limited amount of raw materials
based on qualitative speculations in Ref . Another probable explanation is that the 0°
direction of the cocoon shell possess a higher density of silk fibers. This may be supported by
observations on silkworm’s spinning process, during which the silkworm appears to spin more silk
along the short axis circle of the cocoon ellipsoid. This interesting behavior may be driven by
the “defense” function of the cocoon and suggests that this particular position of the cocoon
may face more attacks than the other places on the cocoon surface. However, further evidence
is needed to give conclusive explanation.
Porter and Vollrath [
] put forward a connectivity model to quantitatively describe the
stress-strain behaviors of cocoons and other particulate composite materials, where the
stressstrain behavior of Antheraea pernyi cocoon was also modelled. We found that the anisotropic
effect of the cocoon stress-strain behaviors can be well modelled based on the connectivity
model, which will be discussed in section 4.3.
4.2 The influence of strain rate
The table listed data of averaged typical properties (E, σy, σu, εu) of cocoon materials with the
increase of nominal strain rate ranging from 0.00625 s-1 to 12.5 s-1 are also presented and
compared in Fig 5. These results show cocoon shell specimens of different orientations have different
strain rate sensitivity. Increased strain rate results in significantly increased Young’s modulus,
yield strength and decreased ultimate strain for the cocoon specimens along the 0° direction; but
for the other two directions, the increased strain rate results in the same trend for Young’s
modulus and ultimate strain, but only negligible difference in the yield and ultimate strength.
For polymer materials increasing strain rate has an equivalent effect as decreasing
temperature to result in a “hardened” mechanical response [
]. In other words, higher strain rate does
not allow polymer chains or segments to relax, hence a higher modulus is often observed
accompanying the higher strain rate. The decrease in the ultimate strain can also be explained
by the lack of polymer molecular relaxation due to the increased strain rate and shorter
response times. Drodge and et al. have studied high strain-rate constitutive behavior of Bombyx
mori silk fibers[
]. Other semi-crystalline polymers have been reported to have strain-rate
sensitive (dynamic rate range) stress-strain behaviors [
], although this effect on silkworm
cocoons was not studied in detail. Coming also from metal materials background, we see
cocoons behave similar to other strain rate dependent materials including high strength steels
21, 22, 30
], whose yield stress and ultimate strength tend to increase to keep pace with the
increasing strain rate. Bearing in mind that the strain-rate dependence of cocoon shells could
be intrinsically determined by the strain-rate dependence of the constituent silk fibers and
binder, we propose to understand the dynamic strain rate dependence of cocoon shells using
well-established models for structural materials and fiber composites.
4.3 Modelling the anisotropy and strain-rate dependence of mechanical
properties of cocoon shells
Chen et al. constructed a connectivity model (strain-activated), which is widely used for
predicting non-woven composite properties, for understanding the cocoon shell stress-strain
9 / 16
Fig 5. Typical parameters of cocoon materials with the increase of loading rate: the trends of (a) Young’s modulus (b) yield strength (c) ultimate strength and
(d) ultimate strain of dried cocoon materials in different directions with increasing loading rate.
mechanical properties using a typical fiber-binder bonded morphology. Eq (1) shows a
modulus or stiffness decreases as the numbers of broken bonds are activated by increasing strain
which can be quantified by an Arrhenius activation function.
where D is a parameter showing strain rate sensitivity of cocoon materials. Here we use only
one strain activation mechanism (e.g. bonds lost through the yield of sericin binder), and the
strain rate dependency part on the right in Eq (2) is adapted from the Johnson-cook model
which will be discussed next. The model parameters of f, εa as shown in Table 3 are derived
from the experimental data at quasi-static rate and then applied to predict the mechanical
properties of cocoon materials at dynamic strain rates. The 90o direction for both dried and
natural cocoons gives a greater “f” and a slightly smaller “εa”, which suggests that the sericin
bonding connectivity in the 90o direction contributes more to the overall stress-strain behavior.
“D” shows a clear trend of decreasing with the increase of the orientation degrees for the dried
cocoons, which suggests that the 0o direction has a greater strain-rate sensitivity.
Fig 6 shows the modelled and experimental stress-strain curves of the dried cocoon shell
specimens, which suggests a good agreement between the experiments and the improved
model and can be used for predicting the stress-strain behaviors of cocoon shells under varied
dynamic loading rates. Nevertheless, at higher strain rate, e.g. 12.5 s-1, the model curves start to
deviate from the experimental curves which may require a further understanding of the
combinational effect of the strain rate on both the non-woven structure and the material properties
Johnson-Cook model is one of the most widely used constitutive models to describe strain
rate dependency in metals [
] and it is interesting to find out the tensile behavior of
cocoon materials has similar profiles as metals under quasi-static and dynamic strain rates.
Although materials properties and structural changing mechanisms are distinct between metals
and silk, the form of Johnson-Cook model could be used to mathematically describe the strain
rate dependence of mechanical properties of cocoon materials. The expression of the
Johnsoncook model is written as
s ¼ ðA þ BεpnÞð1 þ C ln ε_ 0Þ
Fig 6. Comparison between experimental and the improved model under different strain rates.
where σ is the stress, εp is the plastic strain and ε_ =ε_ 0 is dimensionless relative strain rate, with
ε_ 0 ¼ 0:005 s 1 as a reference strain rate, A, B, C and exponential term n are four material
parameters. Note that this model is used for predicting the post-yield plastic mechanical
behaviors and the strain rate dependence of metals. A is the yield stress of cocoon materials at a strain
rate of 0.005 s-1, B and n reflect the nonlinearity of stress-strain behaviors of the cocoon
materials after yield, which can be obtained by curve fitting at a strain rate of 0.005 s-1, while C is a
parameter showing strain rate sensitivity of cocoon materials, determined by plotting different
yield stresses against different relative strain rates and calculating the slope. Table 4
summarized parameters of this model which could describe the post-yield mechanical behaviors of
cocoon materials of three orientations. The values of the strain-rate parameter C are directly
used as D in the improved connectivity model.
The comparison between experimental and the Johnson-Cook model is presented in Fig 7.
Take 0° direction of the dried cocoon materials for example. The Johnson-Cook model could
nicely reflect the trend of modulus decreasing above the yield for cocoon shell structures.
To summarize the section, a modified connectivity model with the introduction of the
strain-rate dependence parameter could well describe the overall stress-strain behaviors and
the strain rate dependency of the cocoon shell materials, and the strain-rate sensitivity of
cocoon shell materials may use the mathematical description from Johnson-cook model, which
well-explains the post-yield tensile mechanical behaviors of cocoon shells.
4.4 The differences between dried and natural cocoon materials
The general variation trends of the main mechanical property indexes of both dried and
natural cocoon materials with the increase of loading rate are shown in Fig 8. There is a tendency
Fig 7. Comparison between experimental and the Johnson-Cook model under different strain rates.
12 / 16
Fig 8. Comparison of some typical parameters ((a) Young’s modulus (b) yield strength (c) ultimate strength (d) ultimate strain) between dried cocoon
materials and natural cocoon materials.
that E, σy and σu of dried cocoon materials are higher than those of natural ones. In addition to
that, dried cocoon materials appear to have lower ultimate strain when compared with natural
cocoon materials as well. In other words, natural cocoons with water lubrication tend to
possess lower Young’s modulus but higher ultimate strain as a trade-off which may result in an
overall larger deformation before the cocoon fails. Previous study [
] has shown that silk
fibers from natural cocoons becomes stiffer after being heated at the temperature of 120°C.
Furthermore, since silkworms spin protein polymers into multifunctional silk fiber as an
essential building block of cocoons, cocoon materials dried at a temperature of 120°C in a laboratory
oven for 48 hours under full ventilation are expected to have changed physicochemical
properties and protein structures of silks due to the loss of water molecules, resulting in enhanced
mechanical properties such as higher modulus and strength in an engineering sense.
Cocoons that naturally evolve to act as impact-resistant armor for silkworm species may serve
as an inspiration for impact protection material designs in engineering applications. In this
paper, we focused on the anisotropic shape effect and the dynamic-rate effect on tensile
mechanical properties of the cocoon shell from silkworm species Antheraea pernyi. Mechanical
13 / 16
properties including Young’s modulus, yield strength, ultimate strength and ultimate strain of
the cocoon shell material have been shown to depend on the microstructure of the cocoon
ellipsoid, which may be further attributed to the silk fiber orientation relative to the short axis
of the cocoon ellipsoid. It is interesting that the cocoon material specimens in the cocoon
short-axis direction showed enhanced mechanical behaviors compared to the specimens in the
other directions. Furthermore, we investigated the strain rate dependence of the tensile
mechanical behaviors of the cocoon shell material. Young’s modulus and ultimate strain
showed more obvious sensitivity to strain rate, and higher dynamic rates resulted in
significantly greater modulus and suppressed strain. Connectivity model and Johnson-cook model
were referenced and modified to describe the mechanical properties of the cocoon materials
upon the strain rate dependency. Additionally, the comparison between dried and natural
cocoon materials showed that dried cocoon materials possessed enhanced mechanical
properties, which suggests thermal treatments through the removal of water and the structural
changes in the silk fibers may be used to tailor the overall mechanical performances of the
This research serves as a preliminary step toward understanding the strain rate dependent
mechanical behaviors of cocoons and other natural materials, and hence may inspire us to
understand the mechanisms of biological defense and further to facilitate the design of
impactresistant artificial composites in the future.
This work is financially supported by Fundamental Research Funds for the Central
Universities, Beihang University, Startup fund for “Zhuoyue 100” titled professorship, Beihang
University, "FengRu" Research Funds of School of Transportation Science and Engineering, Beihang
University and Research Funds for National Student Innovation Training Program. GJ would
like to thank Beihang University for the associate professorship of “Zhuoyue 100” program and
the startup fund, and the National Natural Science Foundation of China (51503009).
Conceived and designed the experiments: JX GJ. Performed the experiments: WZ XG WM.
Analyzed the data: WZ XG WM JX. Contributed reagents/materials/analysis tools: JG. Wrote
the paper: JX JG WZ.
14 / 16
15 / 16
1. Shah DU , Vollrath F , Porter D. Silk cocoons as natural macro-balloon fillers in novel polyurethanebased syntactic foams . Polymer . 2015 ; 56 : 93 - 101 . doi: 10 .1016/j.polymer. 2014 . 09 .021
2. Vollrath F , Porter D . Silks as ancient models for modern polymers . Polymer . 2009 ; 50 ( 24 ): 5623 - 32 . doi: 10 .1016/j.polymer. 2009 . 09 .068
3. Cranford SW , Tarakanova A , Pugno NM , Buehler MJ . Nonlinear material behaviour of spider silk yields robust webs . Nature . 2012 ; 482 ( 7383 ): 72 - 6 . doi: 10 .1038/nature10739 PMID: 22297972 .
4. Trancik JE , Czernuszka JT , Bell FI , Viney C . Nanostructural features of a spider dragline silk as revealed by electron and X-ray diffraction studies . Polymer . 2006 ; 47 ( 15 ): 5633 - 42 . doi: 10 .1016/j. polymer. 2005 . 01 .110
5. Shah DU , Porter D , Vollrath F. Can silk become an effective reinforcing fibre? A property comparison with flax and glass reinforced composites . Composites Science and Technology . 2014 ; 101 : 173 - 83 . doi: 10 .1016/j.compscitech. 2014 . 07 .015
6. Chung da E , Kim HH , Kim MK , Lee KH , Park YH , Um IC . Effects of different Bombyx mori silkworm varieties on the structural characteristics and properties of silk . International journal of biological macromolecules . 2015 ; 79 : 943 - 51 . doi: 10 .1016/j.ijbiomac. 2015 . 06 .012 PMID: 26072984 .
7. Jo YN , Um IC . Effects of solvent on the solution properties, structural characteristics and properties of silk sericin . International journal of biological macromolecules . 2015 ; 78 : 287 - 95 . doi: 10 .1016/j. ijbiomac. 2015 . 04 .004 PMID: 25869308 .
8. Tulachan B , Meena SK , Rai RK , Mallick C , Kusurkar TS , Teotia AK , et al. Electricity from the silk cocoon membrane . Scientific reports . 2014 ; 4 : 5434 . doi: 10 .1038/srep05434 PMID: 24961354; PubMed Central PMCID : PMC4069722 .
9. Chen F , Porter D , Vollrath F. Structure and physical properties of silkworm cocoons . Journal of the Royal Society, Interface / the Royal Society . 2012 ; 9 ( 74 ): 2299 - 308 . doi: 10 .1098/rsif. 2011 .0887 PMID: 22552916; PubMed Central PMCID : PMC3405738 .
10. Chen F , Porter D , Vollrath F. Silkworm cocoons inspire models for random fiber and particulate composites . Physical Review E . 2010 ; 82 ( 4 ): 041911 .
11. Wang F , Xu H , Wang Y , Wang R , Yuan L , Ding H , et al. Advanced silk material spun by a transgenic silkworm promotes cell proliferation for biomedical application . Acta biomaterialia . 2014 ; 10 ( 12 ): 4947 - 55 . doi: 10 .1016/j.actbio. 2014 . 06 .031 PMID: 24980060
12. Chen F , Porter D , Vollrath F. Morphology and structure of silkworm cocoons . Materials Science and Engineering: C. 2012 ; 32 ( 4 ): 772 - 8 . doi: 10 .1016/j.msec. 2012 . 01 .023
13. Zhang J , Kaur J , Rajkhowa R , Li JL , Liu XY , Wang XG . Mechanical properties and structure of silkworm cocoons: A comparative study of Bombyx mori, Antheraea assamensis, Antheraea pernyi and Antheraea mylitta silkworm cocoons . Materials Science and Engineering: C. 2013 ; 33 ( 6 ): 3206 - 13 . doi: 10 . 1016/j.msec. 2013 . 03 .051
14. Reed EJ , Viney C . Mechanical properties of Bombyx mori silkworm silk subjected to microwave radiation . Journal of Materials Research . 2014 ; 29 ( 07 ): 833 - 42 . doi: 10 .1557/jmr. 2014 .54
15. Chen F , Hesselberg T , Porter D , Vollrath F. The impact behaviour of silk cocoons . The Journal of experimental biology . 2013 ; 216 (Pt 14): 2648 - 57 . doi: 10 .1242/jeb.082545 PMID: 23804671 .
16. Zhao H-P , Feng X-Q , Cui W-Z , Zou F-Z . Mechanical properties of silkworm cocoon pelades . Engineering Fracture Mechanics . 2007 ; 74 ( 12 ): 1953 - 62 . doi: 10 .1016/j.engfracmech. 2006 . 06 .010
17. Zhao H-P , Feng X-Q , Yu S-W , Cui W-Z , Zou F-Z. Mechanical properties of silkworm cocoons . Polymer . 2005 ; 46 ( 21 ): 9192 - 201 . doi: 10 .1016/j.polymer. 2005 . 07 .004
18. Blossman-Myer B , Burggren WW . The silk cocoon of the silkworm, Bombyx mori: macro structure and its influence on transmural diffusion of oxygen and water vapor . Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology . 2010 ; 155 ( 2 ): 259 - 63 .
19. Drodge DR , Mortimer B , Holland C , Siviour CR . Ballistic impact to access the high-rate behaviour of individual silk fibres . Journal of the Mechanics and Physics of Solids . 2012 ; 60 ( 10 ): 1710 - 21 . doi: 10 . 1016/j.jmps. 2012 . 06 .007 PMID: WOS: 000307694800002 .
20. Huh H , Kim S-B , Song J-H , Lim J-H. Dynamic tensile characteristics of TRIP-type and DP-type steel sheets for an auto-body . International Journal of Mechanical Sciences . 2008 ; 50 ( 5 ): 918 - 31 . doi: 10 . 1016/j.ijmecsci. 2007 . 09 .004
21. Yu H , Guo Y , Lai X. Rate-dependent behavior and constitutive model of DP600 steel at strain rate from 10−4 to 103s−1 . Materials & Design. 2009 ; 30 ( 7 ): 2501 - 5 . doi: 10 .1016/j.matdes. 2008 . 10 .001
22. Singh NK , Cadoni E , Singha MK , Gupta NK . Dynamic tensile behavior of multi phase high yield strength steel . Materials & Design . 2011 ; 32 ( 10 ): 5091 - 8 . doi: 10 .1016/j.matdes. 2011 . 06 .027
23. Fu C , Porter D , Chen X , Vollrath F , Shao Z. Understanding the Mechanical Properties of Antheraea Pernyi Silk-From Primary Structure to Condensed Structure of the Protein . Advanced Functional Materials . 2011 ; 21 ( 4 ): 729 - 37 . doi: 10 .1002/adfm.201001046
24. Rajkhowa R , Kaur J , Wang X , Batchelor W. Intrinsic tensile properties of cocoon silk fibres can be estimated by removing flaws through repeated tensile tests . Journal of The Royal Society Interface . 2015 ; 12 ( 107 ): 20150177 -. doi: 10 .1098/rsif. 2015 .0177
25. Du S , Li J , Zhang J , Wang X . Microstructure and mechanical properties of silk from different components of the Antheraea pernyi cocoon . Materials & Design . 2015 ; 65 : 766 - 71 . doi: 10 .1016/j.matdes. 2014 . 09 .066
26. Chen F , Porter D , Vollrath F. Silkworm cocoons inspire models for random fiber and particulate composites . Physical Review E . 2010 ; 82 ( 4 ). doi: 10 .1103/PhysRevE.82.041911 PMID: WOS: 000282946500008 .
27. Porter D , Gould PJ . Predictive nonlinear constitutive relations in polymers through loss history . International Journal of Solids and Structures . 2009 ; 46 ( 9 ): 1981 - 93 . doi: 10 .1016/j.ijsolstr. 2009 . 01 .012 PMID: WOS: 000264974200016 .
28. Dusunceli N , Colak OU . High density polyethylene (HDPE): Experiments and modeling . Mechanics of Time-Dependent Materials . 2006 ; 10 ( 4 ): 331 - 45 . doi: 10 .1007/s11043-007-9026-5 PMID : WOS: 000248006900004 .
29. Bahners T , Schlosser U , Gutmann J . Characterization of the Mechanical Properties of Technical Fibers at Extreme Strain Rates. Macromolecular Materials and Engineering . 2012 ; 297 ( 6 ): 550 - 8 . doi: 10 . 1002/mame.201100312 PMID: WOS: 000304711500008 .
30. Curtze S , Kuokkala VT , Hokka M , Peura P . Deformation behavior of TRIP and DP steels in tension at different temperatures over a wide range of strain rates . Materials Science and Engineering: A . 2009 ; 507 ( 1-2 ): 124 - 31 . doi: 10 .1016/j.msea. 2008 . 11 .050
31. Umbrello D. Finite element simulation of conventional and high speed machining of Ti6Al4V alloy . Journal of Materials Processing Technology . 2008 ; 196 ( 1-3 ): 79 - 87 . doi: 10 .1016/j.jmatprotec. 2007 . 05 .007
32. Clausen AH , Børvik T , Hopperstad OS , Benallal A . Flow and fracture characteristics of aluminium alloy AA5083-H116 as function of strain rate, temperature and triaxiality . Materials Science and Engineering: A. 2004 ; 364 ( 1-2 ): 260 - 72 . doi: 10 .1016/j.msea. 2003 . 08 .027
33. Smerd R , Winkler S , Salisbury C , Worswick M , Lloyd D , Finn M. High strain rate tensile testing of automotive aluminum alloy sheet . International Journal of Impact Engineering . 2005 ; 32 ( 1-4 ): 541 - 60 . doi: 10 .1016/j.ijimpeng. 2005 . 04 .013
34. Khan AS , Sung Suh Y , Kazmi R . Quasi-static and dynamic loading responses and constitutive modeling of titanium alloys . International Journal of Plasticity . 2004 ; 20 ( 12 ): 2233 - 48 . doi: 10 .1016/j.ijplas. 2003 . 06 .005
35. Dey S , Børvik T , Hopperstad OS , Leinum JR , Langseth M. The effect of target strength on the perforation of steel plates using three different projectile nose shapes . International Journal of Impact Engineering . 2004 ; 30 ( 8-9 ): 1005 - 38 . doi: 10 .1016/j.ijimpeng. 2004 . 06 .004
36. Guan J , Porter D , Vollrath F. Thermally induced changes in dynamic mechanical properties of native silks . Biomacromolecules . 2013 ; 14 ( 3 ): 930 - 7 . doi: 10 .1021/bm400012k PMID: 23405856 .