[MAP Logo]

Materials Algorithms Project
Program Library



Program MAP_COMP_COAXCY

  1. Provenance of code.
  2. Purpose of code.
  3. Specification.
  4. Description of subroutine's operation.
  5. References.
  6. Parameter descriptions.
  7. Error indicators.
  8. Accuracy estimate.
  9. Any additional information.
  10. Example of code
  11. Auxiliary subroutines required.
  12. Keywords.
  13. Download source code.
  14. Links.

Provenance of Source Code

C.M. Warwick and T.W. Clyne,
Composites and Coatings Group,
Department of Materials Science and Metallurgy,
University of Cambridge,
Cambridge, U.K.

E-mail
T.W. Clyne: twc10@cus.cam.ac.uk

Released: January 1998.
Added to MAP: September 1999.

Top | Next

Purpose

This program calculates the radial distribution of axial, radial and hoop stresses in a series of coaxial cylinders (the inner cylinder being solid), due to a temperature change from the initial (stress-free) state and/or an externally applied axial and/or radial stress.

Top | Next | Prev

Specification

Language:N/A
Product form:Executable files for use on almost any Apple Macintosh or PC.

Complete Program.

Top | Next | Prev

Description

This program calculates the radial distribution of axial, radial and hoop stresses in a series of coaxial cylinders, the inner cylinder being solid, due to a temperature from the initial (stress-free) state change and/or an externally applied axial and/or radial stress. The calculations are carried out using a series of analytical solutions, which are presented in reference 1.

The number of component cylinders is user-defined. The outer surface of the outermost cylinder is free (unconstrained). Input data are required in the form of the elastic constants, thermal expansivities and thicknesses of each cylinder, and the externally applied stress and temperature change. These are screen input in response to a series of prompts. There is an option to use a default material system or to use material property data for any particular cylinder from a library. This library contains data for SiC, Mg, Mg-10wt%Li, Al, Y, Ti, graphite, Ni, W, glass, epoxy and rubber.

Three output files are produced which contain the following output data:

*.rad - radial position (microns) and radial stress (MPa) (sigmar)
*.hoo - radial position (microns) and hoop stress (MPa) (sigmah)
*.axi - radial position (microns) and axial stress (MPa) (sigmaz)
(* refers to is a user-supplied name for the output files.)

Downloading and running the program
Compiled versions of the programs have been produced as stand-alone applications. They are run simply by double-clicking on the icon concerned. They should run on virtually any Apple Macintosh or PC. Data input is via the screen or by selection of default cases which are offered. Data output is to files which are named by the user. These are produced as files for the plotting application "Kaleidagraph", but they can be read as text files from many other plotting or spreadsheet applications. These output files are normally created within the currently-active folder. The program quits after each complete set of computations. For further use, it is necessary to double-click on the icon again.

The executable files for downloading have been compressed using STUFFIT EXPANDER on the Macintosh and WINZIP or PKUNZIP on the PC. These decoders can be downloaded from the following websites:

STUFFIT EXPANDER at http://www.aladdinsys.com/expander/
WINZIP at http://www.winzip.com/

The run time depends on the number of cylinders specified (and on machine speed), but even on a relatively slow machine and with the maximum number of cylinders (10) specified, the program only takes a second or two to run.

Top | Next | Prev

References

  1. C.M. Warwick and T.W. Clyne, 1991, J. Mat. Sci., 26, 3817-3827.

Top | Next | Prev

Parameters

Input parameters

The data library contains material parameters for SiC, Mg, Mg-10wt%Li, Al, Y, Ti, graphite, Ni, W, glass, epoxy and rubber.

Output parameters

r - real
The radial position (microns).

sigmar - real
Radial stress (MPa).

sigmah - real
Hoop stress (MPa).

sigmarz - real
Axial stress (MPa).

The output data is written to the output files:
*.rad - radial position (microns) and radial stress (MPa)
*.hoo - radial position (microns) and hoop stress (MPa)
*.axi - radial position (microns) and axial stress (MPa)

Top | Next | Prev

Error Indicators

None.

Top | Next | Prev

Accuracy

No information supplied.

Top | Next | Prev

Further Comments

Further information about this program can be obtained from the Composites and Coating Group website at http://www.msm.cam.ac.uk /mmc/mmc.html and is one of a series of programs produced by Bill Clyne and co-workers in the Materials Science Department at Cambridge.

It should be noted:
Nothing is expected of anyone downloading a program and there is no obligation to use results obtained from it in any particular manner. Equally, there is no liability on the part of the supplier and no guarantee that the programs do not incorporate errors or invalid data. In general, however, the offer is aimed at researchers and is designed to stimulate collaborative work. Anyone downloading a program is therefore invited to give their address to the supplier and is also welcome to enter into contact if they wish to explore any details. In the event that results obtained using any of the programs are published in any form, it would be appreciated if their source could be acknowledged.

Top | Next | Prev

Example

1. Program text

Complete program.

2. Program data

Stress distribution in a set of Coaxial Cylinders
  CM Warwick & TW Clyne, Cambridge University, 1991
ref: J Mat Sci, vol.26 (1991) p.3817-3827
(Unlimited distribution version - please acknowledge when publishing)

Stresses from differential thermal contraction and/or from axial & radial applied stress
System is transversely isotropic

Enter number of data points required per cylinder (2 to 30) 2
Default system (graphite core (r=15 |m), SiC sheath r=75 |m))? (Y/N)  n
Enter number of cylinders (2 to 10)  2

Cylinder number 1
Enter outer radius in |m 50
Choose material from library? (Y/N)  y

Materials available:
A - SiC          B - Mg    C - Mg/10wt%Li   D - Al         E - Y              F - Ti
G - graphite     H - Ni   I - W            J - glass      K - epoxy          L - rubber

Enter a letter (A-L) G
Axial Youngs modulus (GPa) = 2.3E+0002
Transverse Youngs modulus (GPa) = 2.2E+0001
Poisson ratio (nu12) = 4.2E-0001
Poisson ratio (nu13) = 3.0E-0001
Axial thermal expansivity (microstrain/K) = -1.5E+0000
Transverse thermal expansivity (microstrain/K) = 2.7E+0001

Cylinder number 2
Enter outer radius in |m 80
Choose material from library? (Y/N)  y

Materials available:
A - SiC          B - Mg    C - Mg/10wt%Li   D - Al         E - Y              F - Ti
G - graphite     H - Ni   I - W            J - glass      K - epoxy          L - rubber

Enter a letter (A-L) A
Axial Youngs modulus (GPa) = 4.0E+0002
Transverse Youngs modulus (GPa) = 4.0E+0002
Poisson ratio (nu12) = 1.9E-0001
Poisson ratio (nu13) = 1.9E-0001
Axial thermal expansivity (microstrain/K) = 4.3E+0000
Transverse thermal expansivity (microstrain/K) = 4.3E+0000

Enter temperature change in K (-ive for cooling) 20
Enter applied axial stress in MPa 50
Enter applied radial stress in MPa 30
Computation complete.  Output is stored in 3 files.
Enter prefix for these date files coaxcy

3. Program results

coaxcy.axi
          r(microns)         sigmaz(MPa)
 0.00000000000E+0000 4.26580302358E+0001
 1.25000000000E+0001 4.26580302358E+0001
 2.50000000000E+0001 4.26580302358E+0001
 3.75000000000E+0001 4.26580302358E+0001
 5.00000000000E+0001 4.26580302358E+0001
 5.00000000000E+0001 5.47063908745E+0001
 5.75000000000E+0001 5.47063908745E+0001
 6.50000000000E+0001 5.47063908745E+0001
 7.25000000000E+0001 5.47063908745E+0001
 8.00000000000E+0001 5.47063908745E+0001

coaxcy.hoo
          r(microns)         sigmah(MPa)
 0.00000000000E+0000-5.41722689735E+0000
 1.25000000000E+0001-5.41722689735E+0000
 2.50000000000E+0001-5.41722689735E+0000
 3.75000000000E+0001-5.41722689735E+0000
 5.00000000000E+0001-5.41722689735E+0000
 5.00000000000E+0001 1.10823928048E+0002
 5.75000000000E+0001 9.66508571708E+0001
 6.50000000000E+0001 8.70942248194E+0001
 7.25000000000E+0001 8.03469070426E+0001
 8.00000000000E+0001 7.54067011505E+0001

coaxcy.rad
          r(microns)         sigmar(MPa)
 0.00000000000E+0000-5.41722689735E+0000
 1.25000000000E+0001-5.41722689735E+0000
 2.50000000000E+0001-5.41722689735E+0000
 3.75000000000E+0001-5.41722689735E+0000
 5.00000000000E+0001-5.41722689735E+0000
 5.00000000000E+0001-5.41722689735E+0000
 5.75000000000E+0001 8.75584397970E+0000
 6.50000000000E+0001 1.83124763311E+0001
 7.25000000000E+0001 2.50597941079E+0001
 8.00000000000E+0001 3.00000000000E+0001

Top | Next | Prev

Auxiliary Routines

None.

Top | Next | Prev

Keywords

composite, laminate, coating, cylinder, axial, radial, hoop, stress, temperature

Top | Next | Prev

Download

Download MAP files
Download executable file for the PC
Download executable file for the Macintosh

Top | Prev


MAP originated from a joint project of the National Physical Laboratory and the University of Cambridge.

MAP Website administration / map@msm.cam.ac.uk

Top | Index | MAP Homepage Valid HTML 3.2!