08.02.2013 Views

Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP

Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP

Journal of Mechanics of Materials and Structures vol. 5 (2010 ... - MSP

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

MATRIX OPERATOR METHOD FOR THERMOVISCOELASTIC ANALYSIS OF COMPOSITE STRUCTURES 847<br />

Even though, for comparison purposes, the parameters Ĉ1, Ĉ2, Êi, <strong>and</strong> ˆγ i could be chosen arbitrarily,<br />

we chose them via a least squares fitting <strong>of</strong> the functions aT (T ) <strong>and</strong> E(ξ) experimentally obtained for a<br />

modified asphalt binder (PG58-40). (It is not essential to have perfect fitting <strong>of</strong> experimental data with<br />

the new functions âT (T ) <strong>and</strong> Ê(ξ); we just need to use the same functions in solving the problem using<br />

the two approaches being compared.)<br />

The values <strong>of</strong> the parameters in<strong>vol</strong>ved in âT (T ) are<br />

Those appearing in Ê(ξ) are<br />

Ĉ1 = 41.156640, Ĉ2 = 2.276346 ◦ C −1 .<br />

(in MPa) Ê0 = 38.598647, Ê1 = 349.348243, Ê2 = 195.939985, Ê3 = 126.336051,<br />

(in sec −1 ) ˆγ 0 = 0, ˆγ 1 = −0.517286 · 10 9 , ˆγ 2 = −0.265315 · 10 8 , ˆγ 3 = −0.137463 · 10 7 .<br />

We next present the parameters <strong>of</strong> the composite cylinder, chosen according to an ABCD specimen<br />

described in detail in the next subsection. The specimen consisted <strong>of</strong> an elastic ring with thickness<br />

r1 − r0, surrounded by a viscoelastic binder with thickness r2 − r1 (compare Figure 3 on page 845).<br />

These geometric parameters are<br />

while the remaining parameters are<br />

r0 = 23.75 · 10 −3 m, r1 = 25.4 · 10 −3 m, r2 = 31.75 · 10 −3 m, (36)<br />

T0 = 18 ◦ C, C0 = −1 ◦ C/hour,<br />

α1 = 1.4 · 10 −6 ◦ C −1 ,<br />

α2 = 2.0 · 10 −4 ◦ C −1 ,<br />

E1 = 141 GPa,<br />

ν1 = 0.3,<br />

ν2 = 0.33.<br />

Figure 4 compares the analytical plane strain solution for σ bind<br />

θθ (r1) obtained with the use <strong>of</strong> the Laplace<br />

����������������������� �� � ��������������������<br />

���<br />

���<br />

���<br />

���<br />

���<br />

���<br />

���<br />

�������������������������<br />

����������������������<br />

��������������������������<br />

��������������������������������<br />

��� ��� ��� ��� � �� ��<br />

��������������<br />

Figure 4. Comparison <strong>of</strong> circumferential stress σ bind<br />

θθ (r1) in the composite cylinder problem<br />

obtained analytically (Laplace transform method) <strong>and</strong> numerically by the matrix<br />

representation <strong>of</strong> the relaxation operator ˜E <strong>of</strong> (22)–(24).<br />

(37)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!