The ATILA element library isdescribed on the following pages. The information in brackets following theelement name refers to the list of Section I.B and specifies one or several references concerning the element design.
The SPRI02AM element is a two-node spring element used tomodel the memory effect of a shape memory material.
Ux, Uy, and Uz (3 translations).
MATERIALS entry parameters
These correspond to a shape memory alloy with a memorybehavior (Section III.D).
GEOMETRY entry parameters
S
where S is the cross section area of the rod.
GEOMETRY
1 = 0.0005 * Section 1
2 = 0.0015 * Section 2
…
ELEMENTS
SPRI02AM SMA1 2 = 4 5 * The rod
SPRI02AM SMA1 1 = 5 1 / 5 6 / 5 8 * The “umbrella”
QUAD08E STEEL
1 3 8 10 2 6 7 9
…
END
4 2 6 * Same Oy displacement as that of node 6
5 2 6 * Idem
The SPRI02AS element is a two-node spring element used tomodel the superelastic effect in a shape memory material.
Ux, Uy, and Uz (3 translations).
MATERIALS entry parameters
These correspond to a shape memory alloy with a superelasticbehavior (Section III.D).
GEOMETRY entry parameters
S
where S is the cross section area of the rod.
GEOMETRY
1 = 0.0005 * Section 1
2 = 0.0015 * Section 2
…
ELEMENTS
SPRI02AS SMA2 2 = 4 5 * The rod
SPRI02AS SMA2 1 = 5 1 / 5 6 / 5 8 * The “umbrella”
QUAD08E STEEL
1 3 8 10 2 6 7 9
…
END
4 2 6 * Same Oy displacement as that of node 6
5 2 6 * Idem
The HEXA20C element is a twenty-node isoparametrichexahedron used to model multi-layer composite materials. The formulation ofthis element relies on the Halpin and Tsai model. The element must be a partof a thick plate or of a thick cylindrical shell having a circularcross-section. It must be simply curved with a constant thickness and aconstant curvature in the element Oxy plane. No curvature or thicknessvariation is allowed in the element Oyz plane.
Ux, Uy, and Uz, (3 translations).
(fiber skew angle PHI1 shown here)
MATERIALS entry parameters
Fiber and resin characteristics. See the detaileddescription of the MATERIALS entry, (Section III.D).
GEOMETRY entry parameters
NLAY PHI1 PHI2 ... PHIN
where NLAY is the number of layers, PH1 to PHIN are, foreach layer, the values of the fiber skew angle referenced to the element Ozaxis. All layers have the same thickness, deduced from element size and NLAY.
GEOMETRY
2 = 4 30. -30. 30. -30.
…
…
ELEMENTS
HEXA20C PLC3 2
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
These eight-node or six-node curved shell elements are usedto model multi-layer composite materials. The formulation of this elementrelies on the Halpin and Tsai model.
Ux, Uy, Uz, qx, and qy (3 translations, 2 rotations).
(Fiber skew angle PHI1 shown here)
MATERIALS entry parameters
Fiber and matrix characteristics. See the detaileddescription of the MATERIALS entry, (Section III.D).
GEOMETRY entry parameters
NLAY PHI1 PHI2 ... PHIN T
where NLAY is the number of layers, PH1 to PHIN are, foreach layer, the values of the fiber skew angle referenced to the element Oxaxis, T is the shell thickness.
Degrees-of-freedom qx and qy areexpressed in local axis. For each node, if is the vector normal to the element,
,
and
are the directions of the global axis, qx is the rotation around the direction
and qy the rotation around thedirection
. If
and
are parallel, qx is the rotationaround the direction
and qy the rotation around the direction
.
GEOMETRY
2 = 4 30. -30. 30. -30. 0.01
…
…
ELEMENTS
SHEL08C PLC3 2
6 1 8 3 4 7 2 5
The QUAD08C element is an eight-node isoparametric elementused to model composite materials in axisymmetrical analysis. The formulationof this element relies on the Halpin and Tsai model. Ox is the axis ofsymmetry. The shell can be doubly curved.
Ux, and Uy (2 translations).
MATERIALS entry parameters
Fiber and matrix characteristics. See the detaileddescription of the MATERIALS entry, (Section III.D).
GEOMETRY entry parameters
NLAY PHI1 PHI2 ... PHIN
where NLAY is the number of layers, PH1 to PHIN are, foreach layer, the values of the fiber skew angle referenced to the element Ozaxis.
Local axes (directions 1-2, 1-3, and global Oz) of all 2Delements must form either a uniquely direct or a uniquely inverted system.ATILA does not detect inconsistent axes.
GEOMETRY
2 = 4 30. -30. 30. -30.
…
…
ELEMENTS
QUAD08C PLC3 2
6 1 8 3 4 7 2 5
These isoparametric elements are used to model isotropicelastic materials, with or without material losses.
Ux, Uy, and Uz, (3 translations).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
Not required.
Local axes (directions 1-2, 1-3, and 1-5 for HEXA20E andPYRA13E, directions 1-2, 1-3 and 1-4 for PRIS15E and TETR10E) should form adirect system. ATILA detects inverted elements and renumbers themautomatically.
ELEMENTS
HEXA20E AU4G
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
These flat plate elements are used to model isotropicelastic materials with or without material losses. They rely on the classicalLove-Kirchhoff hypotheses and must be flat.
Uz, qx, and qy (1 translation, 2 rotations).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
T
where T is the plate thickness
GEOMETRY = 10 = 0.065 * Thickness
…
…
ELEMENTS
PLAT08E AU4G 10
6 1 8 3 4 7 2 5
The FACE08E element is an eight-node quadrilateral thinshell element, used for homogeneous isotropic elastic materials, with orwithout material losses. It is plane and enables the facet modeling ofshells. It relies on the classical Love-Kirchhoff hypotheses and must be flat (facet).
Ux, Uy, Ux, qx, and qy (3 translations, 2 rotations).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
T
where T is the shell thickness.
GEOMETRY = 10 = 0.065 * Thickness
…
…
ELEMENTS
FACE08E AU4G 10
6 1 8 3 4 7 2 5
The SPRI02E element is a two-node linear spring element,which is only able to transmit to its two limiting nodes a restoring force dueto a change of its length. The spring constitutive material may or may notinclude material losses. Its classical use is the modeling of prestress rods.
Ux, Uy, and Uz (3 translations).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
S
where S is the cross section area of the rod.
GEOMETRY
1 = 0.0005 * Section 1
2 = 0.0015 * Section 2
…
ELEMENTS
SPRI02E STEEL 2 = 4 5 * The rod
SPRI02E STEEL 1 = 5 1 / 5 6 / 5 8 * The “umbrella”
QUAD08E STEEL
1 3 8 10 2 6 7 9
…
END
4 2 6 * Same Oy displacement as that of node 6
5 2 6 * Idem
These isotropic elements are used to model isotropic elasticmaterials with or without material losses. They can be used for plane-stress,plane-strain or axisymmetrical analyses. For axisymmetrical problems, theglobal Ox axis is the axis of symmetry.
Ux, and Uy (2 translations).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
T
where T is the thickness of the element. Not required andignored for plane-strain or axisymmetrical analyses.
Local axes (directions 1-2, 1-3, and global Oz) of all 2Delements must form either a uniquely direct or a uniquely inverted system.ATILA does not detect inconsistent axes.
CLASS PLSTRESS
GEOMETRY = 10 = 0.065 * Thickness
…
…
ELEMENTS
QUAD08E AU4G 10
6 1 8 3 4 7 2 5
The SHEL03E element is a three-node axisymmetrical thinshell element, with or without material losses, the formulation of which takesaccount of its double curvature. It relies on the classical Love-Kirchhoffhypotheses. The global Ox axis is the symmetry axis.
Ux, Uy, and qz (2 translations, 1 rotation).
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicelastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
T R
where T is the thickness of the element and R its radius ofcurvature.
GEOMETRY
4 = 0.065 0.7 * Thickness + radius
5 = 0.065 0.9 * Thickness + radius
…
…
ELEMENTS
SHEL03E AU4G 5
25 23 22
SHEL03E AU4G 4
16 19 25
These isoparametric elements are used to modelelectrostrictive materials.
Ux, Uy, Uz , and F (3 translations, 1electrical potential).
MATERIALS entry parameters
The elastic, electrostrictive, dielectric tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
1 = 0. 0. 0. *
…
ELEMENTS
HEXA20ES PMN 1
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
These isoparametric elements are designed to model anyelectrostrictive material. They can be used for plane-strain or axisymmetricalanalyses.
Ux, Uy, and
MATERIALS entry parameters
The elastic, electrostrictive, dielectric tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to an isotropic elastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
1 = 0. 0. 0.
…
ELEMENTS
QUAD08P PMN 1
6 1 8 3 4 7 2 5
These isoparametric elements are used to model homogeneousfluid media, with or without losses.
The pressure P.
MATERIALS entry parameters
These correspond, depending upon the case, to a fluid, withor without material losses (Section III.D).
GEOMETRY entry parameters
Not required.
Local axes (directions 1-2, 1-3, and 1-5 for HEXA20F andPYRA13F, directions 1-2, 1-3 and 1-4 for PRIS15F and TETR10F) should form adirect system. ATILA detects inverted elements and renumbers them automatically.
ELEMENTS
HEXA20F WATER
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
These isoparametric elements are used to model homogeneousfluid media, with or without material losses. These elements can be used forplane-strain or axisymmetrical analyses.
The pressure P.
MATERIALS entry parameters
These correspond, depending upon the case, to a fluid, withor without material losses (Section III.D).
GEOMETRY entry parameters
Not required.
Local axes (directions 1-2, 1-3, and global Oz) of all 2Delements must form either a uniquely direct or a uniquely inverted system.ATILA does not detect inconsistent axes.
ELEMENTS
QUAD08F WATER
6 1 8 3 4 7 2 5
These isoparametric elements are used to model isotropicmagnetic media.
The magnetic potential f.
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicmagnetic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
Not required.
ELEMENTS
HEXA20G VACUUM
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
These isoparametric elements are designed to model anyisotropic magnetic media. These elements can be used for plane-strain oraxisymmetrical analyses.
The magnetic potential f.
MATERIALS entry parameters
These correspond, depending upon the case, to an isotropicmagnetic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
Not required.
ELEMENTS
QUAD08G VACUUM
6 1 8 3 4 7 2 5
These isoparametric elements are used to ensure the matchingbetween solid and fluid meshes along their interface. An element includes 8solid nodes and 8 fluid nodes that have the same coordinates as the solidnodes. It has no thickness.
Ux, Uy, and Uz for each solid node (3 translations), and P (the pressure)for each fluid node.
MATERIALS entry parameters
Not required.
GEOMETRY entry parameters
Not required.
ELEMENTS
QUAD16I
9 14 11 16 1 6 3 8 12 10 15 13 4 2 7 5
The LINE06I element is a six-node isoparametric element usedfor plane-strain (two-dimensional) or axisymmetrical analyses. It ensures thematching between solid and fluid meshes along their interface. It includes 3solid nodes and 3 fluid nodes that have the same coordinates as the solidnodes. It has no thickness.
Ux, and Uy for each solid node (2 translations), and P (the pressure) for each fluid node.
MATERIALS entry parameters
Not required.
GEOMETRY entry parameters
Not required.
ELEMENTS
QUAD08E AU4G
1 3 6 8 2 4 5 7
LINE06I
9 14 3 8 12 5
QUAD08F WATER
9 11 14 16 10 12 13 15
These isoparametric elements are used to modelmagnetostrictive materials, with or without material losses.
Ux, Uy, Uz and f (3 translations, 1 magneticpotential).
MATERIALS entry parameters
The elastic, piezomagnetic, magnetic tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to a magnetostrictive material, with or without material losses (Section III.D).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
1 = 0. 90. 0. * polarization along global Oz
…
ELEMENTS
HEXA20M TERFENOL 1
1 3 6 8 13 15 18 20 2 4 5 7 9 10 11 12 14 16 17 19
The QUAD16M element is a sixteen-node isoparametric elementused to model magnetostrictive thin films, with or without material losses.Because the thickness is small compared to other dimensions, a HEXA20M elementis not suitable. This element is made of two node layers having the samedisplacements but different magnetic potential degrees-of-freedom.
Ux, Uy, Uz and f (3 translations, 1 magneticpotential).
MATERIALS entry parameters
The elastic, piezomagnetic, magnetic tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to a magnetostrictive material, with or without material losses (Section III.D).
GEOMETRY entry parameters
ALPHA BETA GAMMA T
These are the three Euler angles that define the naturalaxes of the material with respect to the global coordinate system, followed bythe element thickness (see the detailed description of the GEOMETRYPOLARIZATION CARTESIAN entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
6 = 0. 0. 0. 0.0001 * polarization along global Ox
…
ELEMENTS
QUAD16M TERFENOL 6
9 14 11 16 1 6 3 8 12 10 15 13 4 2 7 5
...
END
14 4 6 * Same magnetic potential
12 4 4 * on the magnetostrictive film edges
9 4 1
These isoparametric elements are designed to model anymagnetostrictive material, with or without material losses. These elements canbe used for plane-strain or axisymmetrical analyses.
Ux, Uy, and
MATERIALS entry parameters
The elastic, piezomagnetic, magnetic tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to a magnetostrictive material, with or without material losses (Section III.D).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
1 = 90. 0. 0. * polarization along global Oy
…
ELEMENTS
QUAD08M TERFENOL 1
6 1 8 3 4 7 2 5
These isoparametric elements are used to model piezoelectricmaterials with or without material losses.
Ux, Uy, Uz , and F (3 translations, 1electrical potential).
MATERIALS entry parameters
The elastic, piezoelectric, dielectric tensors, and the density(see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to an isotropic elastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
The TRIL08P element is an eight-node quadrilateral used tomodel piezoelectric trilaminars. These trilaminars are composed of a metalliccore sandwiched between two plates of a piezoelectric material. The elementformulation relies, for the mechanical part, on the classical Love-Kirchhoffhypotheses. The polarization and the electrical field are takenperpendicular to the plane and uniform on the piezoelectric plates.
Uz, qx, qy, and
MATERIALS entry parameters
Young’s modulus, Poisson’s ratio, and density for the metalliclayer, elastic, piezoelectric, dielectric tensors, and density for thepiezoelectric part. See the detailed description of the materials entry,Section III.D.
GEOMETRY entry parameters
TM TP
where TM is the thickness of the metallic layer and TP thepiezoelectric plate thickness.
GEOMETRY = 4 = 0.065 0.003 * Thicknesses
…
…
ELEMENTS
TRIL08P AG5X31 4
6 1 8 3 4 7 2 5
These isoparametric elements are designed to model anypiezoelectric material, with or without material losses. These elements can beused for plane-strain or axisymmetrical analyses.
Ux, Uy, and
MATERIALS entry parameters
The elastic, piezoelectric, dielectric tensors, and thedensity (see the detailed description of the MATERIALS entry in Chapter III). These correspond, depending upon the case, to an isotropic elastic material, with or without material losses (Section III.D).
GEOMETRY entry parameters
These are the coordinates of the O'x'y'z' system origin andthe Euler angles that define the natural axes of the material with respect tothe global coordinate system (see the detailed description of the GEOMETRYPOLARIZATION entry, Section III.D).
GEOMETRY POLARIZATION CARTESIAN
1 = 90. 0. 0. * polarization along global Oy
…
ELEMENTS
QUAD08P X51 1
6 1 8 3 4 7 2 5
These isoparametric elements are used to prescribe amonopolar or dipolar radiation condition. They are attached to the outsidesurface of a three-dimensional fluid mesh.
The pressure P.
MATERIALS entry parameters
These correspond to a fluid, withoutmaterial losses (Section III.D).
GEOMETRY entry parameters
R
where R is the element radius of curvature.
GEOMETRY = 3 = 7.5 * Radius
…
…
ELEMENTS
QUAD08R WATER 3
6 1 8 3 4 7 2 5
The LINE03R element is a three-node isoparametric linearelement used to prescribe a monopolar or dipolar radiation condition. It mustbe circular and attached on the outer surface of a two-dimensional fluiddomain. It can be used for plane-strain and axisymmetrical analyses.
The pressure P.
MATERIALS entry parameters
These correspond to a fluid, withoutmaterial losses (Section III.D).
GEOMETRY entry parameters
R
where R is the element radius of curvature.
GEOMETRY = 4 = 0.75 * Radius
…
…
ELEMENTS
LINE03R WATER 4
16 18 17
These isoparametric elements are used to couple an
Ux, Uy, and Uz (3 translations).
MATERIALS entry parameters
Not required for BEM coupling. Not used but required forlocal impedance.
GEOMETRY entry parameters
Not required for BEM coupling.
N D CN0R CN0I CN1R CN1I ... CN(N-1)R CN(N-1)I&
CD0R CD0I CD1R CD1I ... CD(D-1)R CD(D-1)I
where N and D are the number of terms in the numerator and denominator (polynomial orders minus one), CNiR + j CNiIis the complex constant of the ith power of pulsation
The local impedance Z(w) is :
CHIEF
ELEMENTS
QUAD08Z
6 1 8 3 4 7 2 5
GEOMETRY
2 = 2 2 & * This geometry to set Z(
0. 0. 0. 100. & * that is, the radiation impedance ofa pulsating sphere
1. 0. 0. 6.7156E-05 * of radius = 0.1 in water.
ELEMENTS
QUAD08Z DUMMY 2
6 1 8 3 4 7 2 5
The LINE03Z element is a three-node isoparametric elementused either to couple an ATILA finite-element modelwith a CHIEF-type or EQI-type boundary element model, or to load a solid structure with a frequency-dependentlocal mechanical impedance. In the first use (BEM coupling), it is superimposedon the radiating surfaces of the mesh of the solid structure. These elementsmust be ordered in the file in the same way as the mutual impedance matrixdescribed in the JOB.ZRAD file. In the second use (local impedance), the localmechanical impedance is provided by means of a fluid materials constants and ageometry number corresponding to the coefficients of a ratio of two frequencydependent polynomials.
Ux, and Uy (2 translations).
MATERIALS entry parameters
Not required for BEM coupling. Not used but required forlocal impedance.
GEOMETRY entry parameters
Not required for BEM coupling.
N D CN0R CN0I CN1R CN1I ... CN(N-1)R CN(N-1)I&
CD0R CD0I CD1R CD1I ... CD(D-1)R CD(D-1)I
where N and D are the number of terms in the numerator anddenominator (polynomial orders minus one), CNiR + j CNiIis the complex constant of the ith power of pulsation
The local impedance Z(w) is :
CHIEF
ELEMENTS
QUAD08E AU4G
9 11 4 6 10 7 8 5
LINE03Z
4 6 5
GEOMETRY
2 = 2 2 & * This geometry to set Z(
0. 0. 0. 100. & * that is, the radiation impedance ofa pulsating sphere
1. 0. 0. 6.7156E-05 * of radius = 0.1 in water.
ELEMENTS
LINE03Z DUMMY 2
4 6 5
![]() |
![]() |
![]() |