next up previous contents
Next: INPUT Up: Using the penalty method Previous: Using the penalty method

  
Axisymmetric solution, TUBE7C1

Mesh: Ten quadrilaterals--see appendix A.3.
Material properties: $\alpha=10^{-5}$ 1/K, $E=2\times 10^5$ MPa, $\nu=0.3$. 



$\dot\epsilon_c=a_1 +a_2\sigma_e+a_3\sigma_e^2+a_4\sigma_e^3~,~~~a_1=a_2=a_3=0~,~~~a_4=2\times 10^{-28}\,\mbox{1/h}$
Support: Plane strain conditions are enforced by setting v=0 for all the nodes 1 to 53, where v is the displacement component parallel to the axis of revolution.
Loading: The internal pressure p=100 MPa.
Execute from prompt:
>rmd2 tube7c1.I1
>rpd2 tube7c1.I2
>srh2 tube7c1.I3
>fefs tube7c1.I4
>hpp2 tube7c1.IP
>hpls tube7c1.IL
>str2 tube7c1.I5