ElAM: A computer program for the analysis and representation of anisotropic elastic properties

详细信息    查看全文

• 作者：Arnaud Marmier ; Zoe A.D. Lethbridge ; Richard I. Walton ; Christopher W. Smith ; Stephen C. Parker ; Kenneth E. Evans
• 关键词：Elastic properties ; Anisotropy ; Visualisation ; Database
• 刊名：Computer Physics Communications
• 出版年：2010
• 出版时间：December 2010
• 年：2010
• 卷：181
• 期：12
• 页码：2102-2115
• 全文大小：485 K

文摘

The continuum theory of elasticity has been used for more than a century and has applications in many fields of science and engineering. It is very robust, well understood and mathematically elegant. In the isotropic case elastic properties are easily represented, but for non-isotropic materials, even in the simple cubic symmetry, it can be difficult to visualise how properties such as Young's modulus or Poisson's ratio vary with stress/strain orientation. The ElAM (Elastic Anisotropy Measures) code carries out the required tensorial operations (inversion, rotation, diagonalisation) and creates 3D models of an elastic property's anisotropy. It can also produce 2D cuts in any given plane, compute averages following diverse schemes and query a database of elastic constants to support meta-analyses.

Program summary

Program title: ElAM1.0

Catalogue identifier: AEHB_v1_0

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEHB_v1_0.html

Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland

Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html

No. of lines in distributed program, including test data, etc.: 43 848

No. of bytes in distributed program, including test data, etc.: 2 498 882

Distribution format: tar.gz

Programming language: Fortran90

Computer: Any

Operating system: Linux, Windows (XP, Vista)

RAM: Depends chiefly on the size of the arrays representing elastic properties in 3D

Classification: 7.7

Nature of problem: Representation of elastic moduli and ratios, and of wave velocities, in 3D; automatic discovery of unusual elastic properties.

Solution method: Stiffness matrix (6×6) inversion and conversion to compliance tensor (3×3×3×3), tensor rotation, dynamic matrix diagonalisation, simple optimisation, postscript and VRML output preparation.

Running time: Dependent on angular accuracy and size of elastic constant database (from a few seconds to a few hours). The tests provided take from a few seconds for test0 to approximately 1 hour for test4.