size scaling of the discontinuity and its effect on the scaling exponents, which have also previously been studied using a complementary approach with the addition of small, non-zero bending rigidity [@sharma_strain-controlled_2016]. Using these modified exponents, we test scaling relations recently predicted for fiber networks [@shivers_scaling_2019].

![ \[fig:1\] Rigidity phase diagram of central force networks. [Upon increasing the average]{} connectivity $z$ at $\gamma=0$, [a network]{} passes through three distinct regimes: (i) a disconnected structure for connectivity less than the percolation connectivity $z<z_p$ (ii) a percolated but floppy network for $z_p<z<z_c \simeq 2d $ and (iii) a rigid network for connectivity greater than $z_c$. Applying a sufficiently large finite strain to an otherwise floppy network with $z_p<z<z_c$ rigidifies the system. For a given $z$ in this range, a critical transition is observed with increasing strain, as indicated by the dashed arrow. The second-order line of transitions is characterized by a critical strain $\gamma_c(z)$ that varies linearly with $z$ near $z_c$ [@wyart_elasticity_2008] (see also Fig. \[fig:A3\] in the Appendix).](fig1){width="6cm" height="6cm"}

Simulation method
=================

To investigate the stiffness discontinuity in fiber networks, we use various network models including (i) triangular, (ii) *phantomized* triangular [@broedersz_criticality_2011; @licup_stress_2015], (iii) 2D and (iv) 3D jammed-packing-derived [@wyart_elasticity_2008; @tighe_force_2010; @baumgarten_normal_2018; @merkel_minimal-length_2019; @shivers_normal_2019], (v) Mikado [@wilhelm_elasticity_2003; @head_deformation_2003], and (vi) 2D Voronoi network [@heussinger_stiff_2006; @arzash_stress-stabilized_2019]. Triangular networks are built by depositing individual fibers of length $W$ on a periodic triangular lattice. The lattice spacing is $\ell_0=1$. A full triangular network has an average connectivity
size scaling of the discontinuity and its effect on the scaling exponents, which have also previously been study use a complementary approach with the addition of little, non - zero bending rigidity [ @sharma_strain - controlled_2016 ]. use these modify exponents, we test scaling relations recently predicted for character networks [ @shivers_scaling_2019 ]. 

! [ \[fig:1\ ] Rigidity phase diagram of cardinal force networks. [ Upon increasing the average ] { } connectivity $ z$ at $ \gamma=0 $, [ a network ] { } passes through three distinct government: (i) a disconnected structure for connectivity less than the percolation connectivity $ omega < z_p$ (ii) a percolated but floppy network for $ z_p < omega < z_c \simeq 2d $ and (iii) a rigid net for connectivity greater than $ z_c$. Applying a sufficiently large finite strain to an otherwise floppy network with $ z_p < omega < z_c$ rigidifies the system. For a given $ z$ in this range, a critical transition is observed with increase strain, as indicated by the dashed arrow. The second - order line of transitions is characterized by a critical strain $ \gamma_c(z)$ that varies linearly with $ z$ near $ z_c$ [ @wyart_elasticity_2008 ] (see besides Fig.   \[fig: A3\ ] in the Appendix).](fig1){width="6 cm " height="6 cm " } 

 Simulation method acting 
 = = = = = = = = = = = = = = = = = 

 To investigate the stiffness discontinuity in fiber networks, we use diverse network models including (iodine) triangular, (ii) * phantomized * triangular [ @broedersz_criticality_2011; @licup_stress_2015 ], (iii) 2D and (iv) 3D jammed - packing - derived [ @wyart_elasticity_2008; @tighe_force_2010; @baumgarten_normal_2018; @merkel_minimal - length_2019; @shivers_normal_2019 ], (v) Mikado [ @wilhelm_elasticity_2003; @head_deformation_2003 ], and (vi) 2D Voronoi network [ @heussinger_stiff_2006; @arzash_stress - stabilized_2019 ]. Triangular network are built by depositing individual fibers of length $ W$ on a periodic triangular wicket. The lattice spacing is $ \ell_0=1$. A full triangular network has an average connectivity
sizf scaling of the discontlnuity and its eydect oi the sdaling ebponents, which have also pretiouwly btvn studied using a cooplementagy approaxh wmth the addition of small, non-zsvo beudmng rigidity [@shsrma_strain-wontrolled_2016]. Usitg tkese modified exponents, we test scalyng relstlons recently kredibtqd fkg niber networks [@shivers_scaling_2019].

![ \[fif:1\] Rigidpty phase diagram of central force networks. [Upoj increasing the agerage]{} conngdtiduty $z$ at $\gamoa=0$, [a netwogn]{} passes tgrough three distinct regimes: (i) a dixconnected srrufjure for coniectivpty less than the perwolatiom connectivity $z<z_'$ (ii) a percolated but flo'py network for $z_p<z<z_s \simeq 2d $ cnd (iii) a rigid netwoek for wonnactixuty grtatxr fhan $z_f$. A'plying a shfficiently large finite straim ei an otherwiss flop[y network with $z_p<z<z_c$ rigidifies the systtm. Fod a given $z$ in this rante, a critical transitlon is obferved with increasing strain, as indicated by the dashxd arxiw. Thd sfcond-order line of transitions is characterizqs ny a critical stvain $\gamma_c(z)$ that fagirf linearly wijh $z$ necd $a_c$ [@wyart_elasticity_2008] (see alfo Fit. \[fig:A3\] in uhe Alpendix).](fig1){width="6cm" height="6cm"}

Sumulation menhod
=================

Ro investigate the stiffness bisconjinuitu in fiber networks, we bse vadious netwogk models kncluding (i) triavgukas, (ii) *phantomized* triangulaw [@broederwz_crnticalith_2011; @livup_strqss_2015], (iii) 2D wnd (iy) 3D jammed-packing-degived [@wfart_elastifity_2008; @tighe_force_2010; @baumgarten_normal_2018; @merkel_minimak-langnh_2019; @shiverf_normsl_2019], (v) Mikado [@wylhelm_elasticijy_2003; @head_deyormatkon_2003], and (vi) 2D Voroioi network [@reussinger_stixv_2006; @arzash_strevs-stabilyzed_2019]. Triqngular vetworks are biilt by dvpjsitung individual fibcrs ow length $W$ on a 'triidic triangular lagtise. Tie laedice spacing is $\dll_0=1$. S fulu triangulav ndtwotk has an average cotnecfivity
size scaling of the discontinuity and its the exponents, which also previously been with addition of small, bending rigidity [@sharma_strain-controlled_2016]. these modified exponents, we test scaling recently predicted for fiber networks [@shivers_scaling_2019]. ![ \[fig:1\] Rigidity phase diagram of central networks. [Upon increasing the average]{} connectivity $z$ at $\gamma=0$, [a network]{} passes through distinct (i) disconnected for connectivity less than the percolation connectivity $z<z_p$ (ii) a percolated but floppy network for $z_p<z<z_c 2d $ and (iii) a rigid network for greater than $z_c$. Applying sufficiently large finite strain to otherwise network with rigidifies system. a given $z$ this range, a critical transition is observed with increasing strain, as indicated by the dashed arrow. The line of characterized by critical $\gamma_c(z)$ varies linearly with $z_c$ [@wyart_elasticity_2008] (see also Fig. \[fig:A3\] height="6cm"} Simulation method ================= To investigate the stiffness in fiber we use various network models including triangular, (ii) *phantomized* triangular [@broedersz_criticality_2011; @licup_stress_2015], (iii) 2D (iv) 3D jammed-packing-derived [@wyart_elasticity_2008; @tighe_force_2010; @baumgarten_normal_2018; @merkel_minimal-length_2019; @shivers_normal_2019], (v) Mikado [@wilhelm_elasticity_2003; @head_deformation_2003], and (vi) 2D [@heussinger_stiff_2006; @arzash_stress-stabilized_2019]. Triangular networks built by depositing fibers length on periodic triangular The lattice spacing is $\ell_0=1$. A full triangular network has an connectivity
size scaling of the discontinUity and its EffecT on The ScAlinG expOnents, which havE Also Previously been studied uSing a CoMPlemENtAry apProach wITh THE adDiTiOn oF sMAlL, non-zEro Bending Rigidity [@shArmA_sTrain-controlLEd_2016]. using these ModIfied exponenTs, wE test sCaLinG RelatIonS receNtly prEDicted For fiber nEtWOrks [@shIVers_scaLINg_2019].

![ \[Fig:1\] RIgidity phase diagrAM oF Central force neTworks. [upON iNCReaSinG the averagE]{} cOnnecTIvity $z$ aT $\GaMMA=0$, [A neTWork]{} passes thrOugh three diSTinCt regiMeS: (i) a DIsconnEcted StRUctUre for conneCtivIty less thAn the pERcolatiON connecTivity $Z<z_p$ (Ii) a PercOLaTeD buT fLOppY NeTwoRK foR $z_p<z<z_c \siMeQ 2d $ And (iiI) a riGID NEtwoRk fOr coNnectIvity greater tHan $Z_c$. ApPLyiNg a suFficiEntlY lArge fInite sTrain To An otherwise flopPy neTwork with $Z_p<z<Z_c$ RigIdIfies THe systEm. FOr a Given $z$ iN this raNGe, a CrITICaL transition is obserVeD WItH increasIng strAIn, As INdicated By The DashED Arrow. the sECoNd-order lIne of tRAnSiTions is ChAracteRiZed By a CritiCAl stRain $\gaMma_c(z)$ thaT variES linearly with $z$ NEar $z_c$ [@wyart_elaSTiCITy_2008] (SEe alSo FIg. \[fig:A3\] in the appeNDix).](fIg1){wiDTh="6Cm" hEIght="6cM"}

SimuLaTIoN Method
=================

To investigate ThE stiffNess dIscontinuity iN fiber netwORKS, we use vaRiouS NeTWork models inclUding (I) triangulaR, (Ii) *phantoMized* TriangulAr [@broederSZ_CriticalIty_2011; @LicUp_sTreSS_2015], (IiI) 2D and (iv) 3D jammeD-PAckiNg-Derived [@WyaRt_elastIciTy_2008; @tIghE_foRcE_2010; @baumgartEn_normal_2018; @MeRkEl_MiNimAl-lenGTh_2019; @shiverS_nOrmAl_2019], (V) MiKado [@wILhelm_eLastiCity_2003; @HeAd_DEfoRmation_2003], ANd (VI) 2d VorOnOi NetwOrk [@HeUssinGer_sTIff_2006; @Arzash_sTress-stabIliZEd_2019]. TrIaNgUlar netWorks are built By Depositing InDivIdual fIBErs of lenGth $W$ on a periodic triangulAR latticE. ThE lattIce sPacing is $\eLl_0=1$. A Full trIanGUlar neTwork hAs an aVeRagE COnnecTIViTy
size scaling of the discon tinuity an d its  ef fec ton t he s caling exponen t s, w hich have also previou sly b ee n  stu d ie d usi ng a co m pl e m ent ar yapp ro a ch  with  th e addit ion of sma ll,  n on-zero bend i ng  rigidity[@s harma_strain -co ntroll ed _20 1 6]. U sin g the se mod i fied e xponents,  w e  tests calingr e la tion s recently predic t ed for fiber netw orks [ @s h iv e r s_s cal ing_2019]. 
![ \[ f ig:1\]R ig i d i typ hase diagramof centralf orc e netw or ks. [Uponincre as i ngthe average ]{}connectiv ity $z $  at $\g a mma=0$,  [a ne two rk] {} p a ss es  th ro u ght hr eed ist inct reg im es : (i)  a d i s c o nnec ted  str uctur e for connect ivi ty l e ssthanthe p erco la tionconnec tivit y$z<z_p$ (ii) aperc olated bu t f lo ppy  n etwor k  for $ z_p <z< z_c \si meq 2d$  an d( i i i)  a rigid network f or c on nectivit y grea t er  t h an $z_c$ .App lyin g a suf fici e nt ly large  finit e  s tr ain toan  other wi seflo ppy n e twor k with  $z_p<z< z_c$r igidifies thes ystem. For ag iv e n  $ z $ in  th is range, a  cri t ical  tra n si tio n  is o bserv ed wi t h increasing strain ,as ind icate d by the dash ed arrow.T h e  second- orde r  l i ne of transiti ons i s characte r ized bya cri tical st rain $\ga m m a_c(z)$tha t v ari esl i ne arly with $z$ n ear$z _c$ [@w yar t_elast ici ty_ 200 8](s ee also F ig. \[fi g: A3 \]  i n t he Ap p endix).] (f ig1 ){ wid th="6 c m" hei ght=" 6cm" }
S i mul ation m e th o d 
=== == == ==== === == =

To  inv e sti gate th e stiffne ssd isco nt in uity in  fiber networ ks , we use v ar iou s netw o r k models  including (i) triangul a r, (ii)  *p hanto mize d* triang ula r [@br oed e rsz_cr itical ity_2 01 1;@ l icup_ s t re ss_ 20 15], (iii) 2 D a nd (i v)  3Djammed- packing-derived [@ w yar t_elasticity_ 200 8; @ t i gh e_f o rc e _20 10 ;  @b a u mgarten_normal_ 2018; @mer ke l _m inimal-len g th_ 20 19; @sh ivers_n ormal _ 2019],(v) Mikad o [@wilhe lm _ela s t ici ty_2003; @ head_def ormation_ 2 003], an d (vi ) 2 D Voro no i n etwor k [@he u ssi nger_ stiff_ 20 06; @a rzash _s tress-st abilized_2019]. Triangu lar ne twork s a re builtbyd epo siting in divi dual fiber s o f l ength  $W $  on a  per i od ict riang ular lattice.T he  la t t ic e spacing i s $ \el l_0=1 $.A  fulltria ngular network ha s  an average co nnec t i vit y
size scaling_of the_discontinuity and its effect_on the_scaling_exponents, which_have_also previously been_studied using a_complementary approach with the_addition of small,_non-zero_bending rigidity [@sharma_strain-controlled_2016]. Using these modified exponents, we test scaling relations recently predicted for_fiber_networks [@shivers_scaling_2019].

![_\[fig:1\]_Rigidity_phase diagram of central force_networks. [Upon increasing the average]{}_connectivity $z$_at $\gamma=0$, [a network]{} passes through three distinct_regimes:_(i) a disconnected_structure for connectivity less than the percolation connectivity $z<z_p$_(ii) a percolated but floppy network_for $z_p<z<z_c \simeq_2d_$_and (iii) a rigid_network for connectivity greater than $z_c$._Applying a sufficiently large finite strain_to an otherwise floppy network with $z_p<z<z_c$_rigidifies the system. For a given_$z$ in this range, a_critical transition_is observed with increasing strain,_as indicated by_the dashed_arrow. The second-order_line of transitions is characterized by_a critical strain_$\gamma_c(z)$ that varies linearly with $z$_near_$z_c$ [@wyart_elasticity_2008] (see_also_Fig. \[fig:A3\]_in the_Appendix).](fig1){width="6cm" height="6cm"}

Simulation method
=================

To_investigate_the stiffness_discontinuity_in fiber networks, we use various_network_models including (i) triangular, (ii) *phantomized* triangular_[@broedersz_criticality_2011; @licup_stress_2015], (iii) 2D_and_(iv) 3D jammed-packing-derived [@wyart_elasticity_2008;_@tighe_force_2010; @baumgarten_normal_2018; @merkel_minimal-length_2019; @shivers_normal_2019], (v)_Mikado [@wilhelm_elasticity_2003; @head_deformation_2003], and (vi) 2D_Voronoi network_[@heussinger_stiff_2006; @arzash_stress-stabilized_2019]._Triangular networks are built by depositing individual fibers of length $W$_on a periodic triangular lattice. The_lattice spacing is $\ell_0=1$._A full_triangular_network has an_average_connectivity