SummaryThe Burgers model is a simple mathematical model which illustrates the nature of real mechanical behaviour of amorphous polymers and other materials with both elastic and viscous nature. The model consists of a Maxwell element joined in series with a Voigt (or Kelvin) element. Maxwell's element can be visualized as a Hookean spring and a Newtonian dashpot in series, thus having viscoelastic properties by combining one elastic parameter (E_{1}) and one viscous parameter (η_{2}). Similarly, the Voigt element obtains viscoelasticity as a result of placing a Hookean spring and a Newtonian dashpot in parallel. Using the Burgers model the program simulates creep. At time t = 0 a fixed tension (tensile stress) is applied to a polymer sample which can be represented by a Burgers model with known parameters. Each of the three elements in series will elongate in a different way. The stress and parameter input values used in the program listing below are adopted from problem 3 of the preliminary examination example Polymeren II. For a better understanding of the calculations that are being executed, consult the solution help at page 4 of this document. You can download Polymeren II from here. 
Program Listing 
Elucidation and Help 
PROGRAM:CREEP  Title 
Normal:Float  > MODE > Normal, Enter 
Func:Connected  
Sequential  Or Simul 
GridOn:AxesOn:LabelOff  > 2nd, FORMAT > .. 
PlotsOff  > 2nd, STAT PLOT > PlotsOff 
FnOff  > VARS > YVARS > On/Off... > FnOff 
ClrHome  > PRGM > I/O > ClrHome 
Disp "E.G.: 1.5E6", "1.0E8 9.0E10", "4.0E7 1.9E9",""  > PRGM > I/O > Disp " : ALPHA, " keys E : 2nd, EE keys 
Input: 

Input "SIGMA? ",S  sigma, σ : tensile stress (Pa) for example 1.5E6 Pa (1.5∙10^{6} Pa) 
Input of the characteristic parameters: 

Disp "E1,ETA2,E3,ETA3?"  
Input A  E_{1} : Young's modulus (Pa) of the Hookean spring 1 e.g. 1.0E8 Pa 
Input B  η_{2} : viscosity coefficient (Pa s) of the Newtonian dashpot 2 e.g. 9.0E10 Pa s 
Input C  E_{3} : Young's modulus (Pa) of the Hookean spring 3 in the Voigt element 3 e.g. 4.0E7 Pa 
Input D  η_{3} : viscosity coefficient (Pa s) of the Newtonian dashpot 3 in the Voigt element 3 e.g. 1.9E9 Pa s 
0 → Xmin : 800 → Xmax  Plot the curves from t = 0 to t = 800 s "→" Symbolizes the STO> key > VARS > Window... > Xmin 
100 → Xscl  
0 → Ymin:.07 → Ymax  Display elongations between 0 and 7 % 
.01 → Yscl  
1 → Xres  
"S/A" → Y5  > VARS > YVARS > Function... > Y5 
"SX/B" → Y6  
"S*(1–e^(–XC/D))/C" → Y7  
"Y5+Y6+Y7" → Y8  
GraphStyle(5,1)  > PRGM > CTL > GraphStyle( 
GraphStyle(6,1)  
GraphStyle(7,1)  
GraphStyle(8,1)  
0 → TblStart  > VARS > Table > TblStart 
50 → ΔTbl  > VARS > Table > ΔTbl (DeltaTbl) 
DispGraph  > PRGM > I/O > DispGraph 
After running the program with the given example input values, the graph is displayed in a window with time values along the horizontal xaxis ranging from 0 to 800 s. The elongation (or strain) values along the yaxis range from ε = 0 to ε = 0,07 (7 %). To retrieve elongation and time values from the four displayed lines, use the TRACE key and the CALC function (2nd, CALC keys). Time is expressed in seconds. Using the TABLE function is another option. Just press the 2nd, TABLE keys. The right arrow key enables viewing the elongations of the Voigt element in the fourth column of the table and the elongations of the Burgers model as a whole in the fifth column. You may wish to call upon the TBLSET function to adjust the table settings. Press CLEAR to return to the normal screen. After running the program with other input values, the presented graph may need some changes for the better. For instance, if the total strain at 800 s exceeds 7 %, use the WINDOW key to adjust Ymax in the WINDOWmenu, then press the GRAPH key to draw the new diagram. Now, if the default Yscl value of 1 % isn't satisfying any more, press the WINDOW key again, and so on. 
Hein ten Horn
Top  Programs Texas Instruments  Home