# Applied stress in a creep analysis **URL:** \ **Category:** FLAC\ **Created:** [October 31, 2024, 7:13pm UTC](https://forum.itascainternational.com/t/applied-stress-in-a-creep-analysis/1781 "2024-10-31T19:13:23Z")\ **Posts on this page:** 3\ **Page:** 1
**Author:** ![agapanto](https://avatars.discourse-cdn.com/v4/letter/a/f4b2a3/32.png) [@agapanto](https://forum.itascainternational.com/u/agapanto)\ **Post date:** [October 31, 2024, 7:13pm UTC](https://forum.itascainternational.com/t/applied-stress-in-a-creep-analysis/1781/1 "2024-10-31T19:13:23Z")
Hi all, I’m performing a creep analysis in which I’d like to apply a series of nstress on a wall. Each applied nstress has to be applied on the wall for only 1 second. **Premise** : I’m not performing a dynamic analysis because my idea is to model this applied stress pattern over a very large time window (days, months or years) that can not be easily modeled in a dynamic configuration. I thus set the creep analysis in this way: > _set cr 0_ > _set fobl 10_ > _set fobu 1E9_ > _set lat 100_ > _set lmul 2_ > _set max 1E8_ > _set min 1E3_ > _set um 0.5_ > _set crdt 0.1_ > _set large_ This is a tipical “time series” of applied nstress: > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 1 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 2 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 3 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 4 > apply nstress -1694.29 var 0.0 -1694.29 from 258, 204 to 257, 234 > solve age 5 > apply nstress -15302.50 var 0.0 -15302.50 from 258, 204 to 257, 234 > solve age 6 > apply nstress -19457.20 var 0.0 -19457.20 from 258, 204 to 257, 234 > solve age 7 > apply nstress -14706.00 var 0.0 -14706.00 from 258, 204 to 257, 234 > solve age 8 > apply nstress -6932.76 var 0.0 -6932.76 from 258, 204 to 257, 234 > solve age 9 > apply nstress -1327.86 var 0.0 -1327.86 from 258, 204 to 257, 234 > solve age 10 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 11 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 12 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > solve age 13 As you can see, it is a kind of sinusoidal “time series” (of 13 seconds) in which the applied stress starts from 0 Pa, reach a maximum value and then go back to 0 Pa. If I plot the xdisplacement induced on a specific zone on which I apply the nstress it seems that the displacements are not affected by the increasing and decreasing stress trend. ![Cattura(3)](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/2X/6/657054e7841dd10fd7c245786e74236e37cbe375.png) I would expect, instead, that displacements would also increase and then decrease following the applied stress trend. It appears, in essence, that the applied stresses accumulate with each other instead of following an up-and-down trend. Otherwise, if I set up a “classic” mechanical analysis (i.e., no creep) in which I define the number of solving steps for each applied nstress the results look different: > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -1694.29 var 0.0 -1694.29 from 258, 204 to 257, 234 > step 100 > apply nstress -15302.50 var 0.0 -15302.50 from 258, 204 to 257, 234 > step 100 > apply nstress -19457.20 var 0.0 -19457.20 from 258, 204 to 257, 234 > step 100 > apply nstress -14706.00 var 0.0 -14706.00 from 258, 204 to 257, 234 > step 100 > apply nstress -6932.76 var 0.0 -6932.76 from 258, 204 to 257, 234 > step 100 > apply nstress -1327.86 var 0.0 -1327.86 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 > apply nstress -0.00 var 0.0 -0.00 from 258, 204 to 257, 234 > step 100 ![Cattura2](https://canada1.discourse-cdn.com/flex031/uploads/itasca/original/2X/f/f0fe53effc6631bcf92036024b8d50b9af7377f5.png) The induced displacements agree with the applied stress trend, first increasing and then decreasing. Therefore, it seems that in this case the applied stresses do not accumulate but are overwritten with each new stress application. Can anyone explain to me why this difference? Where am I going wrong with the creep analysis? Thank you in advance ---
**Author:** ![markus](https://yyz1.discourse-cdn.com/flex031/user_avatar/forum.itascainternational.com/markus/32/288_2.png) [@markus](https://forum.itascainternational.com/u/markus)\ **Post date:** [November 1, 2024, 3:14pm UTC](https://forum.itascainternational.com/t/applied-stress-in-a-creep-analysis/1781/2 "2024-11-01T15:14:06Z")
Is that a one-element model or something more complicated? Anyways I would suspect that your creep timestepping is wrong. Without your total script I cannot recreate it, but my guess would be that the creep timestep is so large, that you only do a few steps for each iteration, thereby not transfering the load changes. Remember that the creep in Itasca-codes is simply a stress correction at each mechanical step and you - in prinicple - have to make sure yourself, that each creep timestep is close to mechanical equilibrium. You can check this in principle, by e.g using a very small fixed timestep (e.g. 1e-4), so that it would take at least 1e4 steps for one “second” and then see, if the displacement behavior is more like the one you expect. ---
**Author:** ![agapanto](https://avatars.discourse-cdn.com/v4/letter/a/f4b2a3/32.png) [@agapanto](https://forum.itascainternational.com/u/agapanto)\ **Post date:** [November 1, 2024, 7:08pm UTC](https://forum.itascainternational.com/t/applied-stress-in-a-creep-analysis/1781/3 "2024-11-01T19:08:26Z")
Thank you very much Markus for your reply! I am new to FLAC and therefore I’m not an expert… I’m trying to reproduce the impact of sea-waves on a vertical rocky cliff. The code that I’m trying to develop (1wave\_cr.txt) restores a .sav file (VT\_j\_th\_12m\_C\_divided.sav) on which I carried out a thermal analysis. I assign a new fictitious cvisc model to the domain “locking” the viscosity-related parameters (vis, k\_vis, k\_shear) by assigning them very high values. As I told you, my intention is not to analyse the viscous behaviour of the model (which is why I set the viscosity-related parameters very high), but I simply want to be able to control the real time of the modelling (i.e. not the mechanical one) without having to perform a dynamic analysis. I set up the creep analysis configuration and I call a .txt file (apply\_1wave\_cr.txt) in which I stored the applied pressure “time series”. History 250 is related to the xdisplacement that I already shown to you. By changing the crdt to 0.0001 as you suggested seems to make no difference ☹ I’ll leave you a wetransfer link from where you can download the above files in case you have a chance to help me out (I’m working with FLAC v. 8.00.463): > **[agapanto\_FLAC](https://wetransfer.com/downloads/08209503314172cc4588b2022c9a4aca20241101184501/d9f5f0?t_exp=1730745901&t_lsid=89817555-6ff4-4e94-999e-609d8ee2b2a5&t_network=link&t_s=download_link&t_ts=1730486701)** > > 3 files sent via WeTransfer, the simplest way to send your files around the world Thank you!!