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
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)