Utafiti huu unaweka mkabala wa Modified Cam Clay, unaotumika kuiga tabia ya udongo wa mfinyanzi, theluji na nyenzo nyingine za punjepunje zinazobanika chini ya mabadiliko makubwa ya umbo, ndani ya mfumo wa plastisiti ya mabadiliko makubwa unaolingana wazi na kanuni za nishati na utawanyaji wa nishati za thermodynamics. Watafiti wametenganisha gradienti jumla ya deformation kwa kuzidisha katika sehemu za elastiki na plastiki; na wamefafanua nishati huru ya Helmholtz pamoja na kazi ya yield kama vipengele viwili vikuu vya kikatiba vya modeli. Milinganyo ya mtiririko wa plastiki na hardening imetolewa kwa kutumia kanuni ya maximum dissipation, hivyo kuonyesha wazi mahusiano ya nishati kati ya vipimo vya stress, deformation ya plastiki na kigezo cha hardening.

Katika formulasheni iliyopendekezwa, hardening haitokei kwa kusasisha moja kwa moja shinikizo la preconsolidation katika modeli ya kawaida ya Modified Cam Clay kwa sheria ya kimaabara ya kielelezo, bali kupitia kigezo cha ndani kilichofafanuliwa katika nishati huru na nguvu yake jumuishi ya thermodynamic. Chini ya mzigo wa shinikizo unaoongezeka, shinikizo la hydrostatic kuwa hasi zaidi kulingana na kanuni ya alama huongeza kasi ya hardening na kuleta maendeleo yanayofanana na tabia ya kawaida ya kubanika kwa kielelezo. Watafiti wanaonyesha kwamba mkabala huu unaendana na mfumo wa maximum dissipation na kwamba usemi wa dissipation ya jumla unaweza kuhesabiwa kwa uwazi.

Modeli imelinganishwa na majaribio ya oedometer kwa theluji na udongo wa Kaolin, pamoja na majaribio ya triaxial ya udongo katika shinikizo tatu tofauti za kuzunguka sampuli. Katika simulation ya theluji ilitumika elementi moja ya mwisho ya hexahedral yenye urefu wa upande wa 76 mm, na katika simulation ya oedometer ya udongo ilitumika elementi moja yenye upande wa 50 mm. Silinda ya jaribio la triaxial iliigwa kwa elementi 2.600 za mwisho za hexahedral za mstari. Mikunjo ya simulation ilionyesha kwa ujumla ulinganifu na matokeo ya majaribio ya ujazo, void ratio na equivalent shear stress; na katika majaribio ya triaxial ilibainishwa kwamba mpito wa elastic-plastic katika data ya majaribio ulikuwa laini zaidi kuliko katika modeli ya kihesabu.

Utafiti pia umechunguza athari za masharti ya mipaka. Dhana ya stress sawia, ambapo mwendo wa pembeni wa nyuso za juu na chini umeachwa huru, ililinganishwa na hali ya msuguano uliokithiri ambapo mienendo hiyo imezuiwa kabisa. Ilionekana kwamba tofauti kati ya mikunjo ya equivalent stress–axial strain ilikuwa ndogo, lakini usambazaji wa axial plastic strain ndani ya sampuli uliathiriwa na masharti ya mipaka. Hivyo, dhana ya takribani stress sawia inayotumiwa katika tathmini ya majaribio ya triaxial imeungwa mkono kwa njia ya kihesabu kwa masharti yaliyochunguzwa.

Kwa mtazamo wa Türkiye: Modeli iliyopendekezwa inaweza kutathminiwa kwa maendeleo ya modeli za constitutive katika uchanganuzi wa finite element wa matatizo ya geotechnical nchini Türkiye, kama vile udongo laini wa mfinyanzi, fills, slopes, foundations, tunnels, miili ya mabwawa na matatizo yanayoweza kuwa na displacements kubwa. Kabla ya matumizi, vigezo vinapaswa kukalibishwa kwa udongo wa mfinyanzi wa ndani nchini Türkiye kwa majaribio ya oedometer na triaxial; saturation, pore-water pressure, drainage, anisotropy, athari ya kasi na tabia ya cyclic loading zinapaswa kuthibitishwa kando. Kutokana na utafiti huu haiwezekani kutoa hitimisho la moja kwa moja maalumu kwa Türkiye kuhusu tabia ya udongo wakati wa tetemeko la ardhi, liquefaction, bearing capacity ya kiwango cha eneo, settlement au usalama wa landslide. Uthibitishaji uliowasilishwa kwa theluji pia unapaswa kujaribiwa upya kwa wiani wa theluji wa ndani, halijoto, muundo wa punje, damage na utegemezi wa kasi kabla ya kuhamishwa kwenda maeneo ya milimani ya Türkiye.

Modeli ya Modified Cam Clay inaeleza nini?

Modified Cam Clay ni modeli ya udongo inayotegemea critical state, inayotumiwa hasa kueleza tabia ya elastic-plastic ya udongo wa mfinyanzi ulioconsolidate kawaida au ulioconsolidate kupita kiasi kwa kiwango kidogo chini ya mizigo ya shinikizo na shear. Yield surface ya modeli katika stress space kwa kawaida huunda jiometri inayofanana na ellipse katika ndege ya mean stress na equivalent shear stress.

Katika hali za stress zilizo ndani ya yield surface, nyenzo kwa hali bora hujitokeza kwa tabia ya elastiki. Hali ya stress inapofikia uso huo, deformation ya plastiki huanza. Kadiri udongo unavyokuwa mnene na wenye nguvu zaidi kutokana na compression, yield surface hupanuka; mabadiliko haya huitwa hardening.

Katika matumizi ya kawaida ya Modified Cam Clay, hardening mara nyingi hutegemea uhusiano wa kielelezo uliobainishwa kati ya plastic volumetric strain na preconsolidation pressure. Ingawa mkabala huu ni muhimu kwa majaribio, ikiwa sheria ya elastiki iliyotumiwa haiwezi kutolewa kutoka kwenye potential ya nishati huru, uhifadhi wa nishati na ulinganifu wa thermodynamic havihakikishwi moja kwa moja.

Tatizo kuu ambalo utafiti unajaribu kutatua ni lipi?

Tatizo kuu linaloshughulikiwa na watafiti ni kutoa kwa ulinganifu milinganyo ya stress, nishati, mtiririko wa plastiki na hardening kutoka katika muundo mmoja wa thermodynamic wakati wa kutumia yield surface ya Modified Cam Clay ndani ya plastisiti ya multiplicative inayofaa kwa deformations kubwa.

Katika mkabala uliopendekezwa kuna kazi mbili za constitutive zinazobainisha tabia ya modeli:

• Nishati huru ya Helmholtz: Hufafanua nishati iliyohifadhiwa katika deformation ya elastiki na mchango wa nishati wa kigezo cha hardening.
• Kazi ya yield: Hubainisha mpaka wa eneo la elastiki, mwelekeo wa mtiririko wa plastiki na masharti ya mzigo ambayo hardening hukua.

Milinganyo ya ukuaji wa plastic deformation gradient na kigezo cha ndani cha hardening hutolewa kutoka kwa kazi hizi mbili kwa kutumia kanuni ya maximum dissipation. Kwa njia hii, sheria ya hardening inakuwa matokeo ya muundo wa nishati–yield wa modeli badala ya kuwa mlinganyo huru wa kimaabara ulioongezwa kutoka nje.

Kinematiki ya deformation kubwa imeundwaje?

Mwili unaodhaniwa mwanzoni kutokuwa na stress hufafanuliwa kwa material coordinates katika reference configuration. Kadiri muda unavyopita, nafasi ya sasa ya mwili hupatikana kupitia deformation map:

\[ \boldsymbol{x}=\boldsymbol{\varphi}(\boldsymbol{X},t) \]

Sehemu ya displacement ni:

\[ \boldsymbol{u}=\boldsymbol{x}-\boldsymbol{X} \]

Material gradient ya deformation map hutoa total deformation gradient:

\[ \boldsymbol{F}=\operatorname{Grad}(\boldsymbol{\varphi})=\frac{\partial\boldsymbol{\varphi}}{\partial\boldsymbol{X}} \]

Tensor F huhamisha vipande vidogo sana vya mstari kutoka reference configuration kwenda current configuration. Mabadiliko ya ujazo hupimwa kwa:

\[ J=\det(\boldsymbol{F}) \]

na hupimwa kwa usemi huu.

Dhana kuu ya finite deformation plasticity ni kugawanya total deformation gradient kwa kuzidisha katika sehemu za elastiki na plastiki:

\[ \boldsymbol{F}=\boldsymbol{F}^{e}\boldsymbol{F}^{p} \]

• F^e: Huwakilisha deformation ya elastiki inayozalisha stress na inayoweza kurudi mzigo unapoondolewa.
• F^p: Ni sehemu ya deformation ya plastiki ya kudumu ambayo haizalishi stress.

Elastic right Cauchy-Green tensor imefafanuliwa kama:

\[ \boldsymbol{C}^{e}=(\boldsymbol{F}^{e})^{T}\boldsymbol{F}^{e} \]

Mabadiliko ya uwezo wa baadaye wa kubeba mzigo kutokana na historia ya plastiki yanaelezwa kwa scalar hardening variable α.

Kwa nini multiplicative decomposition inatumika?

Katika nadharia ya small deformation, strain ya jumla inaweza kugawanywa kwa kujumlisha katika sehemu za elastiki na plastiki. Katika deformations kubwa, mgawanyo huu rahisi hautoshi kutokana na rotations na geometric nonlinearities.

Multiplicative decomposition huchukulia kwamba nyenzo kwanza huhamishwa kwa plastiki kwenda intermediate configuration na kisha kwa elastiki kwenda current configuration. Mkabala huu hutoa msingi wa kinematiki unaofaa zaidi kwa uchanganuzi wa udongo, theluji na nyenzo za punjepunje zenye compression, shear na rotations kubwa.

Vipimo vya stress vimetolewaje kutoka thermodynamics?

Katika utafiti, stress power imeandikwa kwa kutumia quantities zinazooana katika configurations tofauti:

\[ \mathcal{P}=\boldsymbol{P}:\dot{\boldsymbol{F}} =\boldsymbol{\Sigma}^{I}:\boldsymbol{L}^{I} \]

• P: Ni first Piola-Kirchhoff stress tensor.
• Σ^I: Ni Mandel stress katika intermediate configuration iliyofafanuliwa kwa elastic unloading.
• L^I: Ni velocity gradient katika intermediate configuration.

Velocity gradient katika intermediate configuration imegawanywa katika sehemu za elastiki na plastiki:

\[ \boldsymbol{L}^{I} =(\boldsymbol{F}^{e})^{-1}\dot{\boldsymbol{F}}^{e} +\dot{\boldsymbol{F}}^{p}(\boldsymbol{F}^{p})^{-1} =\boldsymbol{L}^{e}+\boldsymbol{L}^{p} \]

Total dissipation imefafanuliwa kwa kutoa mabadiliko ya nishati huru kwa muda kutoka external stress power:

\[ \mathcal{D}=\mathcal{P}-\dot{\psi} \]

Inapokubaliwa kwamba ni deformations za plastiki pekee zinazopaswa kuzalisha dissipation, expressions za stress hupatikana kutoka derivatives za nishati huru:

\[ \boldsymbol{P}=\frac{\partial\psi}{\partial\boldsymbol{F}} \]

\[ \boldsymbol{\Sigma}^{I} =2\boldsymbol{C}^{e}\frac{\partial\psi}{\partial\boldsymbol{C}^{e}} \]

Nguvu jumuishi ya thermodynamic ya hardening variable α ni:

\[ \beta=-\frac{\partial\psi}{\partial\alpha} \]

Kwa ufafanuzi huu, plastic dissipation hupata muundo ulio wazi:

\[ \mathcal{D} =\boldsymbol{\Sigma}^{I}:\boldsymbol{L}^{p} +\beta\dot{\alpha} \]

na hivyo hupata muundo ulio wazi.

Kanuni ya maximum dissipation ina maana gani?

Kanuni ya maximum dissipation huchukulia kwamba kati ya maendeleo yanayokubalika ya stress na variables za ndani, mtiririko unaoongeza matumizi ya nishati ya plastiki ya nyenzo hadi kiwango cha juu ndio hutokea. Katika matokeo ya optimization ambayo yield function inatumika kama constraint, plastic multiplier λ hutokea.

Kasi za plastic deformation na hardening:

\[ \dot{\boldsymbol{F}}^{p} =\dot{\lambda} \frac{\partial f_y}{\partial\boldsymbol{\Sigma}^{I}} \boldsymbol{F}^{p} \]

\[ \dot{\alpha} =\dot{\lambda}\frac{\partial f_y}{\partial\beta} \]

zimetolewa.

Karush-Kuhn-Tucker conditions hutumika kutenganisha hali za elastic loading, plastic loading na unloading:

\[ \dot{\lambda}\geq0,\qquad f_y\leq0,\qquad \dot{\lambda}f_y=0 \]

• Ikiwa f_y < 0, hali ya stress iko ndani ya yield surface na plastic multiplier ni sifuri.
• Ikiwa f_y = 0 na loading inaendelea kuelekea nje, plastic deformation inaweza kuendelea.
• Plastic multiplier hairuhusiwi kuwa hasi.

Nishati huru ya Helmholtz imechaguliwaje?

Total free energy imefafanuliwa kama jumla ya sehemu ya elastiki na sehemu inayotokana na hardening:

\[ \psi=\psi^{e}+\psi^{\alpha} \]

Watafiti wamechagua modeli ya aina ya modified Neo-Hooke kwa elastic response na nishati ya quadratic kwa hardening:

\[ \psi= \frac{\kappa}{2}\left(\ln J^{e}\right)^2 +\frac{\mu}{2}\left(\operatorname{tr}(\bar{\boldsymbol{C}}^{e})-3\right) +\frac{1}{2}H\alpha^2 \]

Hapa:

• κ: Ni bulk compression modulus.
• μ: Ni shear modulus.
• H: Ni hardening modulus.
• J^e = det(F^e): Ni elastic volume change.
• C̄^e = (J^e)^−2/3C^e: Ni elastic strain tensor iliyotenganishwa na sehemu ya volumetric.

Kutoka kwa quadratic hardening energy, thermodynamic hardening force hupatikana kama:

\[ \beta=-H\alpha \]

hupatikana.

Yield function iliyopendekezwa ya Modified Cam Clay ikoje?

Muundo wa kwanza wa yield function ni:

\[ f_y= \frac{1}{\Sigma_y} \left[ \frac{q^2}{M^2} +(p-p_t)\left(p-(p_c+\beta)\right) \right] \]

.

• p: Ni sehemu ya hydrostatic ya Mandel stress.
• q: Ni equivalent shear stress.
• M: Ni parameter inayobainisha slope ya critical state line.
• p_t: Ni threshold upande wa tension.
• p_c: Ni initial pressure threshold.
• β: Ni thermodynamic force ya hardening.
• Σ_y: Ni parameter ya dimensional consistency inayodumisha unit ya stress ya yield function.

Katika utafiti Σ_y imewekwa kiidadi kuwa 1. Sababu ya kuwepo kwa parameter hii ni kwamba katika mlinganyo wa plastic flow, yield function inapaswa kubeba unit ya stress.

Hydrostatic na equivalent stresses zimefafanuliwa kama:

\[ p=\frac{1}{3}\operatorname{tr}(\boldsymbol{\Sigma}^{I}) \]

\[ q=\sqrt{\frac{3}{2}}\, \left\|\operatorname{dev}(\boldsymbol{\Sigma}^{I})\right\| \]

. Katika sign convention ya utafiti, tensile stress ni chanya; kwa hiyo kadiri compression inavyoongezeka p huenda kwenye values hasi zaidi.

Uhusiano β = −Hα unapowekwa kwenye yield function hupatikana:

\[ f_y= \frac{1}{\Sigma_y} \left[ \frac{q^2}{M^2} +(p-p_t)\left(p-(p_c-H\alpha)\right) \right] \]

. Kadiri α inavyoongezeka, yield boundary upande wa compression hupanuka kuelekea compression stresses kubwa zaidi.

Sheria ya hardening inatokeaje?

Derivative ya yield function kwa β inapotumika, hardening rate hupatikana kama:

\[ \dot{\alpha} =\dot{\lambda}\frac{p_t-p}{\Sigma_y} \]

.

Kwa sababu yield surface haiwezi kuvuka threshold p_t kwenye hydrostatic axis, p ≤ p_t. Plastic multiplier pia haiwezi kuwa hasi, kwa hiyo mabadiliko ya α kwa muda si hasi.

Compression load inapoongezeka, p huwa hasi zaidi kwa sababu ya sign convention. Hivyo tofauti p_t − p na hardening rate huongezeka. α inayoongezeka hupanua upande wa compression wa yield surface na kuruhusu compression loads za juu zaidi. Mwingiliano huu huleta maendeleo ya positive feedback yanayofanana na exponential hardening behavior ya modeli ya kawaida ya Cam Clay.

Jambo muhimu katika mkabala huu ni kwamba tabia ya kielelezo haitokani na mlinganyo wa kielelezo uliofafanuliwa moja kwa moja kutoka nje, bali hutokea kutokana na uchaguzi wa pamoja wa energy potential na yield function.

Mlinganyo wa plastic flow una vipengele gani?

Derivative ya yield function kwa Mandel stress hutoa mwelekeo wa plastic deformation. Katika utafiti, maendeleo ya plastic deformation gradient yameandikwa kwa namna inayojumuisha kwa pamoja shear na volumetric components:

\[ \dot{\boldsymbol{F}}^{p} =\dot{\lambda} \left[ \frac{2q}{M^2} \sqrt{\frac{3}{2}} \frac{\operatorname{dev}(\boldsymbol{\Sigma}^{I})} {\left\|\operatorname{dev}(\boldsymbol{\Sigma}^{I})\right\|} +\frac{1}{3} \left(2p-p_c+H\alpha-p_t\right)\boldsymbol{I} \right] \frac{1}{\Sigma_y}\boldsymbol{F}^{p} \]

Term ya kwanza inawakilisha plastic flow katika mwelekeo wa deviatoric au shear, na term ya pili inawakilisha tendency ya volumetric plastic compression au expansion.

Kwa nini dissipation inaweza kuhesabiwa wazi?

Kwa sababu stress na hardening force zimetolewa kutoka free energy, huku plastic flow na hardening rate zikitolewa kutoka yield function, total plastic dissipation inaweza kuandikwa kama:

\[ \mathcal{D} =\boldsymbol{\Sigma}^{I}:\boldsymbol{L}^{p} -H\alpha\dot{\alpha} \]

. Katika makala, mlinganyo huu umeendelezwa hadi muundo wa wazi wa indices unaojumuisha components za elastic na plastic deformation gradients.

Watafiti wanaeleza kwamba nishati iliyohifadhiwa kwa hardening hupunguza kiasi cha total dissipation. Hili ni jambo linalotarajiwa; sehemu ya external work hutumika katika irreversible plastic processes, huku sehemu nyingine ikihifadhiwa katika internal energy inayohusishwa na hali ya hardening ya nyenzo. Modeli imeundwa ili ibaki ndani ya maximum dissipation framework chini ya masharti yanayofaa ya yield na loading.

Milinganyo ya finite element imeundwaje?

Ili modeli itumike katika uchanganuzi wa muundo au udongo, linear momentum balance imechukuliwa kuwa msingi:

\[ \operatorname{Div}(\boldsymbol{P}) +\rho_0\boldsymbol{b} =\rho_0\ddot{\boldsymbol{u}} \]

• ρ₀: Ni density katika reference configuration.
• b: Ni body acceleration au body force inayotenda kwa unit mass.
• ü: Ni derivative ya pili ya displacement field kwa muda.

Mlinganyo umebadilishwa kuwa nodal residual equations kwa kutumia Galerkin method, divergence theorem na finite-element shape functions. Mizani kati ya internal force, inertia na surface loads imeelezwa kwa residual vector.

Kutokana na geometric na material nonlinearities, nodal displacements hutatuliwa kwa iteration kwa kutumia Newton-Raphson method. Global stiffness matrix hutegemea consistent tangent operator inayojumuisha total derivative ya first Piola-Kirchhoff stress kwa deformation gradient.

Plastic internal variables zinasasishwaje kwa muda?

Kwa plastic deformation gradient, time integration ya aina ya exponential map imetumika:

\[ \boldsymbol{F}^{p}_{n+1} =\exp\left[ \lambda\frac{\partial f_y} {\partial\boldsymbol{\Sigma}^{I}} \right]\boldsymbol{F}^{p}_{n} \]

Hardening variable inasasishwa kwa implicit Euler method:

\[ \alpha_{n+1} =\alpha_n+ \lambda\frac{p_t-p}{\Sigma_y} \]

Unknowns katika time step mpya ni:

• Plastic deformation gradient F^p_n+1
• Hardening variable α_n+1
• Plastic consistency parameter λ

. Hutatuliwa kupitia local nonlinear equation system inayojumuisha plastic deformation residual, hardening residual na yield condition.

Kwa hiyo, computation ina viwango viwili vya iteration:

1. Local Newton-Raphson solution ya plastic internal variables katika kila integration point
2. Global Newton-Raphson equilibrium solution ya nodal displacements za muundo

Muundo huu ni muhimu kwa kupata consistent tangent stiffness na kwa finite-element solution kukaribia tabia ya quadratic convergence.

Jaribio la oedometer la theluji limeigwaje?

Kwa tabia ya kubanika ya theluji, data za majaribio ya oedometer zilizochapishwa hapo awali zimetumika. Katika modeli ya kihesabu, elementi moja ya hexahedral ya mstari yenye pande zote za 76 mm iliundwa. Lateral deformations zilizuiwa na displacement inayowakilisha loading-unloading history ya Kielelezo 1 ilitumika kwenye uso mmoja.

Loading history ina hatua kadhaa za compression, partial unloading na reloading. Matokeo ya simulation yalifuata kwa ujumla pointi za majaribio za kupungua kwa ujazo wa sampuli za theluji zenye initial volumes tofauti kadiri hydrostatic pressure ilivyoongezeka.

Katika makala, ulinganifu huu umetathminiwa kuwa “wa kutosha”, lakini hakuna error norm, mean absolute error, confidence interval au parameter uncertainty iliyoripotiwa. Kwa hiyo, validation inategemea hasa comparison ya kuona ya curves.

Jaribio la oedometer la udongo wa Kaolin limeigwaje?

Kwa udongo wa Kaolin, ilitumika elementi moja ya mchemraba yenye upande wa 50 mm. One-dimensional compression condition yenye lateral freedoms zilizozuiwa ilitumika na loading history ikaongezewa unloading-reloading cycle.

Imeelezwa wazi kwamba parameters zilibainishwa kwa trial and error. Simulation ilinasa kwa ujumla nonlinear relationship kati ya hydrostatic pressure na void ratio.

Hadi takribani 9 kPa ilionekana elastic slope, huku katika eneo zaidi ya 10 kPa kukawa na mpito kuelekea plastic response. Mpito ulitokea kwa ulaini badala ya kink kali.

Unloading na reloading paths hazikuingiliana kikamilifu. Katika pressure level ileile, reloading path ilitoa void ratio ya juu zaidi na mwelekeo huu ulionekana kuwa sawa na tabia katika data za majaribio.

Modeli ya jaribio la triaxial imeundwaje?

Katika majaribio ya triaxial, sampuli ya udongo ya silinda kwanza iliwekewa constant confining pressure kutoka pande zote, kisha axial displacement ikawekwa kwenye uso wa juu. Confining pressure katika majaribio inawakilisha takribani homogeneous stress state inayoundwa na maji yenye shinikizo.

Silinda ya kihesabu iligawanywa katika elementi 2.600 za mwisho za hexahedral za mstari. Kielelezo 5 kinaonyesha sampuli ikiwa na urefu wa takribani 100 mm na kipenyo cha 50 mm.

Katika hatua ya kwanza, confining pressure iliwekwa kwenye nyuso za juu na pembeni, na equilibrium displacement ya uso wa juu ikabainishwa kwa iteration. Katika hatua zinazofuata, axial displacement iliongezwa kwa mstari hadi kufikia %10 ya original height.

Equivalent stress q ilihesabiwa kwenye middle node ya upper element layer kwa ku-project integration-point values kwenda nodes.

Ni parameters gani zilizotumika katika majaribio ya triaxial?

Ilidhaniwa kwamba confining pressures tofauti zinawakilisha soil depths na loading histories tofauti. Kwa hiyo Young modulus, pressure threshold na hardening modulus ziliwekwa tofauti kwa kila pressure level.

Slope ya critical-state line ilibainishwa kuwa takribani −1,2 kwa kuchora confining pressure na estimated maximum q values katika p-q plane, na M = 1,2 ikachaguliwa. Tension threshold ilikadiriwa kwanza kuwa 50 kPa, lakini 52 kPa ikatumika ili kupata ulinganifu bora na experimental curves.

Comparisons za majaribio ya triaxial zinaonyesha nini?

Katika confining pressures za 100, 200 na 300 kPa, equivalent stress iliongezeka pamoja na axial strain. Kadiri confining pressure ilivyoongezeka, q value iliyopatikana katika axial deformation ileile pia iliongezeka.

Modeli ilinasa general trend ya experimental curves kwa pressure levels zote tatu. Hata hivyo, katika experimental data mpito kutoka elastic region kwenda plastic region ni laini zaidi. Curves za numerical model zinaonyesha mabadiliko ya slope yaliyo wazi zaidi, hasa katika transition region.

Kwa kuwa parameters ziliwekwa tofauti kwa kila confining pressure, hii inaonyesha kwamba modeli inaweza kuwakilisha majaribio, lakini haithibitishi kwamba prediction imefanywa katika wide pressure range kwa single fixed parameter set.

Kwa nini boundary-condition study ilifanywa?

Katika theoretical triaxial test evaluations, sampuli hudhaniwa kuwa chini ya homogeneous stress. Katika jaribio halisi, friction kati ya loading platens na sampuli inaweza kuzuia lateral movement kwenye nyuso za juu na chini na kusababisha barrel-shaped deformation.

Watafiti walilinganisha extreme conditions mbili kwa confining pressure ya 100 kPa:

• Homogeneous boundary condition: Lateral movement ya nyuso za juu na chini imeachwa huru.
• Non-homogeneous boundary condition: Lateral movement ya nyuso za juu na chini imefungwa kabisa ili kuwakilisha upper bound ya athari inayowezekana ya friction halisi.

Kwa sababu hakuna actual friction coefficient au contact law, complete lateral fixation ni conservative assumption inayokadiria friction zaidi ya uhalisia.

Athari ya boundary conditions ilikuwa nini?

q–axial strain curves za modeli mbili zilipatikana kuwa karibu sana. Kwa hiyo, kwa jaribio na parameters zilizochunguzwa, athari ya end-surface friction kwenye global stress-response curve ilichukuliwa kuwa ndogo kiasi cha kupuuzwa.

Hata hivyo, plastic deformation distributions katika takribani %10 axial strain hazifanani. Modeli yenye lateral movement huru inatoa axial plastic strain field iliyo sawia zaidi, huku katika modeli yenye lateral movement iliyofungwa kukiwa na tofauti wazi za rangi na thamani kati ya maeneo ya juu, kati na chini.

Matokeo haya yanaonyesha kwamba single global stress-strain curve huenda isionyeshe tofauti zote katika local deformation distribution.

Nguvu kuu za utafiti ni zipi?

• Modified Cam Clay yield surface imeformulishwa ndani ya multiplicative kinematics inayofaa finite deformation.
• Stresses, hardening force, plastic flow na hardening rate zimetolewa kutoka muundo mmoja wa thermodynamic.
• Total dissipation imeelezwa kwa uwazi.
• Tabia inayofanana na classical exponential hardening imetokea kiasili kutokana na uchaguzi wa energy na yield function.
• Exponential map na implicit Euler integration zimetumika kwa local plastic update.
• Milinganyo ya msingi ya finite-element implementation yenye consistent tangent operator imetolewa.
• Modeli imetathminiwa kwa data za theluji na Kaolin clay.
• Oedometer na triaxial loading types zimechunguzwa kwa pamoja.
• Confining pressures tatu tofauti zimetumika katika jaribio la triaxial.
• Athari za loading-platen boundary conditions kwenye global na local results zimechunguzwa kando.

Mapungufu makuu ya utafiti ni yapi?

• Utafiti ni preprint ambayo haijapitia peer review.
• Hakuna jaribio jipya la maabara lililofanywa; data za majaribio zilizochapishwa hapo awali zimetumika.
• Parameters za theluji na udongo zimebainishwa kwa kiasi kikubwa kwa curve fitting au trial and error.
• Hakuna optimization, uncertainty analysis au confidence interval iliyotolewa kwa parameter identification.
• Simulation-experiment agreement haijatathminiwa kwa quantitative error metrics.
• Oedometer samples zimeigizwa kwa finite element moja tu.
• Different E, p_c na H values zimetumika kwa kila confining pressure ya triaxial tests.
• Predictive ability ya single parameter set katika wide stress range haijaonyeshwa.
• Pore-water pressure na solid-fluid interaction hazijajumuishwa katika modeli.
• Drained na undrained behavior hazijaigwa tofauti.
• Degree of saturation au partial-saturation effect haipo.
• Soil anisotropy na rotational hardening hazijatathminiwa.
• Cyclic, dynamic na earthquake-type loadings hazijachunguzwa.
• Viscosity, loading-rate na creep effects hazipo.
• Softening, localization na shear-band formation hazijashughulikiwa.
• Kwa theluji, damage, fracture na healing mechanisms hazipo katika modeli hii.
• Elastic model haiwakilishi moja kwa moja ongezeko la stiffness kadiri compression inavyoongezeka.
• Hardening modulus H bado haijaunganishwa na standard experimental parameters kama virgin compression na swelling indices.
• Hakuna open-access link iliyotolewa kwa source code au finite-element implementation.
• Hakuna conflict-of-interest na data-availability statements.

Utafiti unaunga mkono nini?

• Unaunga mkono kwamba Modified Cam Clay model inaweza kuundwa katika thermodynamically consistent multiplicative-plasticity structure.
• Unaonyesha kwamba quadratic hardening energy pamoja na yield function iliyochaguliwa inaweza kuzalisha pressure hardening inayofanana na exponential.
• Unaonyesha kwamba modeli inaweza kuwakilisha general trends za data za oedometer za theluji na Kaolin clay zilizochunguzwa.
• Unaonyesha kwamba modeli inaweza kufuatilia kwa takribani triaxial-test curves katika confining pressures za 100-300 kPa kwa parameters zinazofaa.
• Unaunga mkono kwamba katika triaxial test iliyochunguzwa end-surface constraints zina athari ndogo kwenye global q-strain response.
• Unaonyesha kwamba hata boundary conditions zikitoa global response inayofanana, local plastic-deformation fields zinaweza kubadilika.

Utafiti hauthibitishi nini?

• Hauthibitishi kwamba modeli ni halali kwa kila aina ya clay, sand, fill au theluji.
• Hauelezi undrained behavior ya saturated soils au pore-water pressure.
• Hauonyeshi earthquake, liquefaction au cyclic-mobility behavior.
• Hauthibitishi performance ya foundation, slope, tunnel au dam katika real field scale.
• Hauthibitishi kwamba parameters zilizobainishwa ndizo single na unique parameter set za nyenzo.
• Hauonyeshi kwamba modeli itatoa prediction ya kuaminika kwenye stress paths zisizo na experimental data.
• Hauelezi fracture, shear, temperature effects au rate-dependent behavior ya theluji.
• Hauthibitishi kwamba thermodynamic consistency yenyewe inatoa higher experimental accuracy.

Mbinu na Matokeo ya Utafiti

Muundo wa utafiti

Utafiti unajumuisha utoaji wa milinganyo ya constitutive, uchanganuzi wa thermodynamic, implicit time integration, finite-element implementation na numerical comparison na data za majaribio zilizochapishwa.

Vigezo vikuu vya modeli ya constitutive

Milinganyo ya msingi ya modeli

• Multiplicative decomposition: F = F^eF^p
• Hardening force: β = −Hα
• Plastic flow: Ḟ^p = λ̇(∂f_y/∂Σ^I)F^p
• Hardening: α̇ = λ̇(p_t − p)/Σ_y
• Dissipation: D = Σ^I:L^p + βα̇
• Time integration: F^p_n+1 exponential map, α_n+1 implicit Euler

Mbinu ya numerical solution

• Linear momentum balance imediscretize kwa Galerkin finite-element method.
• Global nodal equilibrium imetatuliwa kwa Newton-Raphson method.
• Katika kila integration point imetumika local Newton-Raphson loop ya pili kwa F^p, α na λ.
• Plastic deformation gradient imesasishwa kwa matrix exponential function.
• Consistent material tangent imejumuishwa katika global stiffness matrix.

Muhtasari wa mifano ya uthibitishaji

Matokeo makuu

• Yield function iliyopendekezwa na energy potential ziliweza kutumika pamoja ndani ya maximum dissipation framework.
• Hardening variable iliongezeka monotonically chini ya compression.
• Compression inayoongezeka ilizalisha tabia ya hardening inayofanana na exponential.
• Volume reduction katika sampuli za theluji iliakisiwa kwa ujumla katika loading paths tofauti.
• Void-ratio–pressure curve na unloading-reloading difference katika Kaolin clay zilinaswa.
• Katika triaxial tests, q value iliongezeka kadiri confining pressure ilivyoongezeka.
• Experimental elastic-plastic transition ilibaki laini zaidi kuliko numerical results.
• End-surface constraints zilibadilisha global q–axial strain curve kwa kiwango kidogo.
• End-surface constraints zilibadilisha kwa kiwango kikubwa local axial plastic-strain distribution.

Hali ya uthibitishaji wa takwimu na kiasi

Katika utafiti, simulation na experimental curves zimelinganishwa kwa graphical comparison. Hakuna mean error, root mean square error, coefficient of determination, confidence interval au parameter uncertainty iliyowasilishwa. Haijachunguzwa kama parameter calibrations nyingi zinaweza kuzalisha experimental curve ileile.

Kwa hiyo, maneno kama “ulinganifu mzuri” au “wa kutosha” yanategemea hasa tathmini ya kuona na ya kihandisi ya watafiti. Matokeo si matokeo ya statistical hypothesis test.

Mahitaji ya kazi za baadaye

• Kuongeza pore-water pressure kama independent degree of freedom
• Kuendeleza coupled solid-fluid na undrained analyses
• Kuunganisha parameter H na virgin compression na swelling indices
• Kuiga elastic stiffness inayoongezeka pamoja na pressure
• Kuthibitisha multiple stress paths kwa single parameter set
• Kuongeza anisotropy na rotational hardening mechanisms
• Cyclic na dynamic loading tests
• Kuchunguza rate dependence na viscoplastic behavior
• Kuongeza damage, softening na deformation localization
• Kubainisha parameters kwa automatic inverse analysis
• Validation kwa independent datasets na field scale
• Kushiriki source code na sample finite-element input files

Dokezo la Chanzo na Mbinu

Jina kamili asili la utafiti: Thermodynamically consistent formulation of the Modified Cam Clay yield surfaces in multiplicative kinematics

Waandishi: Jakob Platen; Ahmad Moeineddin; Michael Kaliske.

Mpangilio wa waandishi: Mpangilio ulio hapo juu ndio mpangilio asili wa waandishi katika utafiti.

Taarifa ya mchango sawa: Hakuna taarifa ya co-first author au equal contribution katika utafiti.

Mwandishi wa mawasiliano: Michael Kaliske.

Barua pepe ya mwandishi wa mawasiliano: Katika toleo hili la preprint hakuna barua pepe ya mawasiliano iliyotolewa.

Taarifa za ORCID:

• Jakob Platen: 0000-0001-6705-6023
• Ahmad Moeineddin: 0000-0002-0427-4119
• Michael Kaliske: 0000-0002-3290-9740

Taasisi: Institute for Structural Analysis, Georg-Schumann-Str. 7, 01187 Dresden, Saxony, Ujerumani. Jina rasmi linalolingana ni TUD Dresden University of Technology / Technische Universität Dresden Institute for Structural Analysis.

DOI:10.2139/ssrn.7202517

Jukwaa rasmi la uchapishaji: SSRN.

Kiungo rasmi cha utafiti:Ukurasa wa utafiti wa SSRN

Mwaka wa uchapishaji: 2026.

Jarida: Kwa toleo hili hakuna jina la jarida mahususi lililopitia peer review au taarifa ya acceptance.

Hali ya uchapishaji: Hati ina kauli “Preprint submitted to Elsevier”. Kauli hii haimaanishi kwamba imekubaliwa na jarida la Elsevier au kwamba mchakato wa peer review umekamilika.

Aina ya chanzo: Makala ya utafiti ya preprint inayojumuisha utoaji wa thermodynamic constitutive model, finite-element implementation na numerical validation kwa data za majaribio zilizochapishwa.

Hali ya peer review: Utafiti haujapitia peer review.

Ufadhili: Utafiti uliungwa mkono na German Research Foundation Deutsche Forschungsgemeinschaft kupitia mradi namba KA 1163/46-2, kitambulisho cha mradi 414936990.

Michango ya waandishi:

• Jakob Platen: Conceptualization, methodology, software, validation, investigation, writing original draft na visualization.
• Ahmad Moeineddin: Conceptualization, methodology, validation, investigation, review na editing.
• Michael Kaliske: Supervision, project administration, funding acquisition, review na editing.

Mgongano wa maslahi: Katika toleo hili hakuna conflict-of-interest statement.

Upatikanaji wa data na code: Data za majaribio zilizotumika zimetolewa kutoka machapisho ya awali. Hakuna open-access link iliyotolewa kwa model implementation, source code, finite-element input files au calibration data.

Maelezo haya ya Kiswahili yameandaliwa kwa kuchunguza derivations za kinematic na thermodynamic za utafiti uliopakiwa, yield function, free-energy equations, plastic-flow na hardening law, finite-element implementation, mifano ya oedometer ya theluji na Kaolin clay, comparisons za triaxial tests, boundary-condition analysis, tables na figures. Hakuna matokeo mapya ya kisayansi kutoka nje ya utafiti yaliyoongezwa. External verification imetumika tu kuthibitisha DOI, rekodi ya SSRN na taarifa za utambulisho wa taasisi.

Matokeo ya utafiti yanaunga mkono uwezekano wa kutumia formulasheni ya Modified Cam Clay iliyo thermodynamically consistent na inayofaa finite deformations. Hata hivyo, saturated-soil behavior, pore-water pressure, cyclic loading, earthquake effects, rate dependence, anisotropy na field-scale applications bado hazijathibitishwa. Parameters zinahitaji kuunganishwa moja kwa moja na experimental indices na quantitative error analysis kufanywa kwa independent datasets.