Assignment 2

CH5130: Rheology of Complex Materials

Due on Friday (14/2/2025) just before the class

1. For a certain material we measure G’ and G” as

where ω is in units of rad/s.

(a) Sketch logG (ω) and logG(ω) versus log ω at frequencies between 0.001 to 100 rad/s.
(b) What is G(t) for this material?
(c) What is η0?
(d) What is the shear stress as a function of time after this material is subjected to a sudden
step strain of γ0 = 5 at time t = 0?

2. If G(t) = G0tp, p>-1, then find J(t) using the convolution theorem.

Useful information:

3. Maxwell mode fitting:

Fit a multimode Maxwell to the following data. Show the attempt of fitting for all
modes. Plot experimental modulus versus modulus obtained from fitting. Also, plot
the error as a function of number of modes. Once you have identified Gi and τi, write
expression for G(t) and plot log G versus log t.

Based on the last digit of your roll number, use the appropriate dataset from the table
given below:

Last digit of your Roll number Dataset to fit

0,7 A
1,9 B
2,6 C
3,5 D
4,8 E

1
Dataset A:
ω G' ω G"
0.041217 0.018802 0.009869 0.16679
0.058956 0.037321 0.020781 0.339406
0.084331 0.074077 0.028681 0.499669
0.11858 0.140875 0.042315 0.690454
0.172519 0.285691 0.060377 0.974344
0.246733 0.555102 0.087636 1.435084
0.352982 1.12556 0.122959 2.024824
0.487987 2.095027 0.178445 2.919391
0.733432 3.660739 0.250408 4.207869
1.031134 6.814865 0.363349 5.938918
1.472413 10.69878 0.509643 8.029667
2.06653 15.75303 0.739162 10.63073
2.998129 21.76453 1.036605 14.06998
4.273869 27.02511 1.477693 17.47075
6.19283 31.48282 4.264567 20.05068
8.823834 36.67011 6.171662 19.69501
12.77975 40.07179 8.782731 19.34264
18.20636 45.68962 12.50042 19.40596
25.48904 47.82876 18.09337 19.47252
36.30103 52.25743 18.09337 19.47252
52.57558 57.10505 26.18057 18.7236

Dataset B:
ω G' ω G"
0.099389 1.210891 0.099157 5.458416
0.12764 2.256436 0.127122 5.981188
0.161231 3.269307 0.162938 6.438614
0.20364 4.249505 0.208868 6.928713
0.261441 5.19703 0.263292 7.320792
0.330029 6.013861 0.331898 7.712871
0.409771 6.79802 0.418197 7.974257
0.5432 7.484158 0.535559 8.170297
0.685112 8.039604 0.87776 8.366337
1.089016 8.921782 1.087008 8.366337
1.395088 9.215842 1.390552 8.235644
1.786794 9.444554 1.749651 8.072277
2.250662 9.607921 2.274611 7.843564
2.881966 9.771287 2.86139 7.614851
3.629362 9.869307 3.59954 7.386139
4.570088 9.934653 4.528602 7.190099
5.850079 10 5.791312 6.961386
7.48857 10.06535 7.407712 6.79802

2
 9.429598 10.13069 9.319685 6.60198
11.87116 10.13069 11.92218 6.471287
15.19769 10.22871 15.00587 6.405941
19.45003 10.22871 19.20038 6.340594
24.49146 10.29406 24.17184 6.340594
30.83961 10.35941 38.95354 6.405941
39.47289 10.39208 38.95354 6.405941
50.52846 10.45743 49.0609 6.536634
82.82316 10.68614 62.84962 6.830693
102.6007 10.78416 81.76007 6.79802

Dataset C:
ω G' ω G"
0.012241 2.265375 0.01 26.33666
0.015672 3.483823 0.012803 34.09654
0.019619 5.474163 0.015672 42.28324
0.024561 8.239253 0.024561 65.02556
0.031446 11.13599 0.024561 65.02556
0.039366 16.40413 0.038492 97.87104
0.049281 23.14653 0.063096 147.3073
0.061694 31.28435 0.098883 207.8533
0.077233 43.20301 0.151526 306.1832
0.123789 77.24146 0.242866 432.0301
0.198409 135.1578 0.380617 596.6243
0.310945 226.5375 0.610051 772.4146
0.48731 379.6988 1 937.4813
0.763707 609.6025 1.532378 1162.574
1.224066 897.989 3.074715 1380.978
1.918346 1213.702 3.849159 1441.712
3.074715 1605.49 4.818666 1505.116
4.818666 1948.586 7.55177 1640.413
7.55177 2469.013 12.10394 1787.872
11.83507 2996.646 18.96919 1907.102
19.40014 3483.823 30.40373 2034.282
30.40373 4138.304 47.64845 1990.973
47.64845 4811.085 76.37071 2034.282
95.60658 5714.909 122.4066 1990.973
153.2378 6364.119 191.8346 1948.586
240.1529 7087.079 300.6414 1866.5
307.4715 7398.759 471.1625 1948.586
376.3655 7241.242

3
Dataset D:
ω G' ω G"
0.202299 0.001217 0.202299 0.381821
0.23393 0.001656 0.237353 0.431604
0.278481 0.002184 0.278481 0.535278
0.372373 0.004172 0.322023 0.586712
0.512601 0.007725 0.44329 0.848334
0.695459 0.014303 0.601423 1.153203
0.943547 0.026486 0.804198 1.567787
1.280134 0.049043 1.123242 2.130794
1.736791 0.088048 1.480291 2.808939
2.35635 0.158076 2.008349 3.818399
3.196921 0.275161 2.764651 5.190127
4.337345 0.47897 3.696778 6.841266
5.884588 0.808364 5.015514 9.299839
8.100603 1.364154 6.904253 12.64073
14.91082 3.433756 9.367179 17.18348
20.2299 5.618823 12.89467 21.95645
27.44643 8.643221 17.24222 28.94147
36.70023 12.89215 23.73528 38.1412
50.52078 19.22607 31.73783 51.85321
69.54586 28.67184 43.05954 66.26268
92.99386 40.20312 58.41999 84.6764
128.0134 59.95492 80.41976 108.1965
171.1743 79.02843 107.534 138.2767
264.675 117.7632 145.8941 171.3245
315.0814 141.5919 200.8349 212.25
433.7345 192.4576 272.4779 271.2321
369.6778 346.6047

Dataset E:
ω G' ω G"
0.195497 0.01 0.195497 0.491219
0.236275 0.013622 0.232857 0.591306
0.3027 0.020997 0.294005 0.757172
0.371211 0.032364 0.371211 0.969566
0.46869 0.049887 0.46869 1.167116
0.591769 0.0818 0.591769 1.494502
0.747167 0.12225 0.943373 2.45054
1.173871 0.290461 1.503887 3.777298
1.844263 0.62901 2.397435 5.645181

4
 2.940048 1.320704 3.712105 8.974704
4.755708 2.527462 5.917686 13.00451
7.471671 4.836863 9.433732 18.27027
11.91103 7.931017 15.03887 25.66823
23.28566 16.65239 23.6275 34.96425
37.66599 26.47395 37.12105 46.17741
59.17686 39.56539 59.17686 60.98667
94.33732 60.98667 94.33732 78.09396
148.2129 83.07361 146.0686 103.139
232.8566 120.3752 187.1337 109.7156
298.3208 128.0509

Suggestions for Fitting Maxwell Model to viscoelastic data:

1. Initial guess value of 𝜏 can be taken as:

𝜏max= 1/𝜔_𝑚ax
𝜏min= 1/𝜔_𝑚in

2. Prepare a regression model (the most popular being least squares regression) using inbuilt
functions in MATLAB (eg: lsqcurvefit) or EXCEL (eg: solver function, you can google about
how to add solver function in excel if it is not already present). You can define an objective
function as sum of all relative squared error.
3. If you’re getting negative values of gi’s you’re probably overfitting, reduce N.
4.. On increasing number of relaxation modes check the change in relative squared error,
increasing N without much change in it is not useful.