Continuum Modeling of Lattice Structures
Continuum Modeling of Lattice Structures
Ahmed K Noor
George Washington University, NASA Langley Research Center, Hampton VA
Appl Mech Rev vol 41, no 7, July 1988 285 © Copyright 1988 American Society of Mechanical Engineers
Periodic structure approaches: Not efficient for complicated Combine with substitute
Substructuring, combined use configurations and/or for continuum approach
of finite elements, and transient response calculations
transfer matrix methods
Exact representation of Limited to beamlike lattices
stiffness of individual with simply supported ends or
members ring configuration with rigid
central mast
the lattice and either solving the resulting difference equations ing for large lattice structures. The present survey paper is an
directly, or using truncated Taylor series expansions to replace attempt to fill this void. Specifically, the objectives of this paper
the difference equations by differential equations (see, for exam- are:
ple, Refs 27, 92, and 96, and the three monographs, Refs 19, 28,
and 121). This approach works well for simple lattice configura- (a) to assess the effectiveness of the currently used ap-
tions, but becomes quite involved for lattices with complex proaches for continuum modeling;
geometry. (b) to present a simple and rational approach for develop-
The third group of methods, referred to as periodic structure ment of continuum models for large lattice structures:
approaches, are based on either: (1) the combined use of finite (c) to identify the future directions of research which have
elements and transfer matrix methods, which is efficient only high potential for realizing the advantages of continuum
for rotationally periodic (ie, cyclically symmetric) structures or modeling.
lattices with simple geometries (Refs 67, 123, and 124), or (2)
the exact representation of the stiffness of an individual mem-
The scope of the present survey includes thermoelastic stress
ber from which the analysis of beamlike lattices with simply
analysis, buckling, free vibration, and geometrically nonlinear
supported edges can be performed (Refs 4-6).
problems of large lattice structures. Beamlike and platelike
The fourth approach is based on using a smearing technique repetitive lattices with pin and rigid joints are considered.
to replace the actual lattice by a substitute continuum model Continuum modeling of lattices with flexible joints is also
that is equivalent to the original structure in some sense, such as discussed.
the constitutive relations, the strain energy, and/or kinetic
energy (see, for example, Refs 3, 22, 23, 29, 30, 33, 40, 41, 43,
44, 46, 47, 48, 49, 51, 55, 58-61, 69-73, 75, 77-81, 83, 85, 88, 2. WHY CONTINUUM MODELING?
95, 98,100, 104, 105,108-111,113, 122,127, and 128). The use
of continuum models to simulate the behavior of planar lattice Before an assessment is made of the different approaches for
beams dates back to the previous century (see Ref 116, p 483). developing continuum models, the following three advantages
It has gained popularity only in recent years and has been of using the continuum modeling approach for analyzing repe-
applied to a variety of other discrete systems and phenomena titive lattice structures are identified:
including solid and liquid crystals, dislocations and defects,
composite materials and biological systems (see Refs 21, 62, 87, (1) It offers a practical and efficient approach for analyzing
106, and 115). large lattice structures. This is particularly true for
The number of publications on continuum modeling of beamlike and platelike lattices, wherein a dimensionality
. repetitive lattice structures has been steadily increasing. Despite reduction can result in a substantial reduction in the
the availability of some survey papers (for example, Ref 3) and number of degrees of freedom.
monographs (Refs 60 and 107), there is a need to broaden (2) Provides a simple means of comparing structural, ther-
awareness among practicing engineers and research workers of mal, and dynamic characteristics of lattices with differ-
the recent developments in various aspects of continuum model- ent configurations and assessing the sensitivity of their
responses to variations in material and geometric proper- The last two characteristics are perhaps the most important
ties. in terms of recent developments, and are discussed in the next
(3) Provides an effective tool for parameter and system two sections.
identification, and for feedback control system design of
lattice structures. This is because in these situations only
the global behavior of the large lattice structure is of 3.1. Local deformations
interest, and the detailed deformation in each member is
The local deformations of two axially loaded planar trusses
not needed.
are shown in Fig 2. The first truss has double lacing and a
single-bay repeating cell. The second truss has single lacing and
It is interesting to note that in the precomputer and first- a double-bay repeating cell. The chord members of the first
generation computer eras an opposite path to continuum mod- truss remain straight as shown on the top sketch. On the other
eling was followed, namely, replacing small segments of the hand, the actual deformation of the single-laced truss has the
two- (and three-) dimensional continuum by discrete analogues zig-zag pattern shown on the top right sketch. On the average,
consisting of one-dimensional members (ie, repetitive lattices, however, the chord members remain straight. Early continuum
see, for example, Refs 50, 68, 74, 99, and 130). models averaged these deformations, thereby substantially over-
estimating the axial stiffness. Recent continuum models, for
lattices with more than one bay in their repeating cells, do
3, DEFINITION AND KEY ELEMENTS OF account for the local deformations (see Refs 77, 78 and 81).
A SUBSTITUTE CONTINUUM MODEL
A number of definitions have been given for the substitute
continuum model. Herein an effective continuum model is de- 3.2. Ordinary versus micropolar continua
fined as a continuum which has the following characteristics: A contrast between the ordinary and micropolar continua is
made in Fig 3. For an axially loaded, pin-jointed truss member
(1) the same amount of thermoelastic strain and kinetic the transverse motion is completely characterized by the joint
energies are stored in it as those of the original lattice displacements. The member rotation t|/ is related to the joint
structure when both are deformed identically; displacements w, and w-. Therefore, the appropriate continuum
(2) the temperature distribution, loading and boundary con- to use in modeling pin-jointed trusses is the ordinary continuum
ditions of the continuum simulate those of the original for which the displacement field completely characterizes the mo-
lattice structure; tion of the structure.
(3) for beamlike and platelike lattices the continuum models On the other hand, for a rigid-jointed member, the trans-
are one-dimensional beams and two-dimensional plates, verse motion is characterized by both the joint displacements w,
respectively (see Fig 1); and vv as well as the joint rotations #, and 0- which are
(4) local deformations are accounted for; independent degrees of freedom. Therefore, the appropriate
(5) lattices with pin joints are modeled as classical continua, continuum to use in modeling rigid-jointed flexural members is
and lattices with rigid (and/or flexible) joints are gener-
ally modeled as micropolar continua.
^m
Deformation
Beamlike lattices Continuum beam mode! predicted by classical
continuum (based on
uniform strain state
within repeating cell)
i L v
9p ~~^S$X j
Rotations
Continuum plate model V .i<w,».,
L 1 1 Bj,6- (joint, rotations)
2. expansion of each of the nodal displacements, strain Phase 2—Generation of thermoelastic stiffnesses and density
components, and temperature in a Taylor series; parameters of the beam and plate continua
3. generation of three transformation matrices which relate
nodal displacements, axial strains and temperatures of Step 4. For beamlike and platelike lattices reduction in
individual members of the repeating cell to the displace- dimensionality is achieved by introducing kinematic and tem-
ments, strain and temperature parameters of the con- perature assumptions regarding the variation of the displace-
tinuum. ment components, strain components, and temperature in the
plane of the cross section (for beamlike lattices) or in the
thickness direction (for platelike lattices).
A schematic representation of the sequence of steps involved
The vectors of nodal displacements and strain components
in the generation of the stiffness coefficients for simplified one-
and the temperature of member k are expressed in terms of the
and two-dimensional continuum models is given in Fig 4. The
corresponding continuum parameters by means of transforma-
application of the procedure to thermoelastic stress and vibra-
tion matrices as follows:
tion analyses of lattices with pin joints consists of the three
major phases described in the following three sections:
(6)
€
22
€
33
[M]C = - Z [ir]'[m<*>)[r<*>] (10)
l-}( (2) hers
members
Yl2
Yl3
1
E
( M c = o (aEALyk)[R^]'[T^]{T}c, (11)
members
[* ( A , ] = [/ 2 Im In (3)
where 0 is a characteristic geometric property of the repeating
cell (length of the repeating cell for beamlike lattices and
€ u , e22, and £33 are the axial strains in the coordinate direc-
planform of the repeating cell for platelike lattices).
tions; yu, yu, and y23 are the shearing strains; and (/, m, n)
are the direction cosines of the member.
For simplicity, in the present study the strain state is as-
sumed to be uniform within each repeating cell. Variation of the Phase 3—Generation of properties of simplified beam
strain state within the repeating cell can be accounted for by and plate continua
expanding {e} ( A ) in a Taylor series about the center of the
repeating cell. The number of terms in the Taylor series expan- Step 6. Simplified continuum models are obtained by parti-
sion is equal to the number of independent deformation modes tioning the displacement and strain parameters of the con-
of the repeating cell. tinuum as follows:
Step 3. The stiffness matrix and the thermal load vector of
the repeating cell are generated using the following equations:
(12)
w<.=
[K]= E (£4LW]'[*< (4)
members
{"}, (13)
{Pr = £ (aEAL) (A) [R ( *>]'r< (5)
members
where {« a } are the strain parameters to be retained in the
where E, A, L, and a are the elastic modulus, cross-sectional simplified model, {tc2} are the strain parameters associated
area, length, and coefficient of thermal expansion of member k, with the neglected stress resultants (or internal forces), {e c3 }
and superscript t denotes transposition. The thermoelastic stiff- are the strain parameters to be omitted, and { uA } and {uc2 }
nesses of the equivalent three-dimensional classical continuum are the displacement parameters associated with the retained
are obtained by dividing the right-hand sides of Eqs. (4) and (5) and neglected mass coefficients, respectively. The stiffness ma-
by the volume of the repeating cell. trix, matrix of density parameters, and thermal load vector of
the continuum are partitioned accordingly: matrices [C], [m], and {Cr} numerically using the FOR.
TRAN program listed in Ref 86 on an IBM PC.
Kn Kn K
n 3. Higher-order continuum models can be obtained by using
*1,= Kn K23 (14) more strain and displacement parameters than {e(1} and
symtn K)3 { u d } . The additional parameters may be needed for the
accurate prediction of the detailed displacement and stress
Mn Mn distribution in the lattice.
[M]c = (15)
symtn M22 4. Similar procedures for generating the characteristics of
the continuum model are presented in Refs 29, 30, 70,
P-n) and 71. In the first two cited references the use of
{Pr}c = - P-ri)- (16) transformation matrices in developing the continuum
PTJ characteristics is discussed. In Refs 70 and 71 the thermo-
elastic stiffness coefficients of the continuum model are
The effective thermoelastic coefficients of the simplified con- obtained by applying the tensor transformation relations
tinuum models are obtained by deleting the rows and columns to the unidirectional stiffness and thermal coefficients of
associated with {«,..,} and by expressing {ec2} in terms of each of the individual members constituting the repeating
{c rf }(using static condensation). The resulting constitutive ma- cell. The three-dimensional continuum coefficients gener-
trix and thermal load vector have the following forms: ated are equal to the stiffness matrix and thermal load
vector [K] and {PT} of Eqs (4) and (5), each divided by
[C]=[Ku]-[Kl2][K22]-l[K21], (17) the volume of the repeating cell.
{CT} = {PT1} -[K^iK^y'iP^}. (18)
The density parameters of the simplified continuum are
obtained by neglecting the terms associated with {u t . 2 }, that is, 6. APPLICATION OF CONTINUUM MODELS TO
STABILITY PROBLEMS
M = [Mn (19)
The dominant criterion in the design of lightweight lattice
Step 7. The thermoelastic strain energy density of the con- trasses is the onset of elastic instability. For finely divided
tinuum can be written in the following form: lattice trusses the failure mode is usually that of general (sys-
£/(.= U0- Ux-U2, (20) tem) instabilily rather than local snap-through of a single joint
or member buckling. The approach presented in the preceding
where U0 is the isothermal strain energy density and Ul and U2 sections can be extended to the general (global) stability analy-
are the contributions to the strain energy of the linear and sis of lattice structures. This is accomplished by equating the
quadratic terms in the temperature parameters. The expressions potential energy due to initial stresses of the original lattice to
for U0 and Ux are that of the continuum, from which expressions for the geomet-
^ = iKi}'[C]{ecl}, ric stiffness coefficients of the continuum are obtained in terms
(21) of the geometric and elastic properties of the original lattice
structure. The procedure is outlined subsequently.
The expression for U2 is not presented herein since U2 is
inconsequential in the development of the properties of the 1. A repeating cell is isolated from the lattice grid. The
continuum models. geometric stiffness matrix of a typical member k of the
Step 8. The kinetic energy density of the continuum is given repeating cell referred to the coordinates of the repeating
by element is given by [G<k)].
2. The nodal displacements of member k are expressed in
® = *{«c-l}'['»]{« r f }, (22) terms of the continuum displacement parameters and
their derivatives with respect to the spatial coordinates by
where a dot refers to derivative with respect to time.
means of the following transformation:
The foregoing approach has been applied in Ref. 86 to the
development of continuum models for beamlike and platelike { « r = [riA,]{A},, (23)
lattice trusses with pin joints.
where { A }c is the vector ofjiisplacement parameters and
its spatial derivatives and [Tlk)] is a transformation ma-
trix. Equation (23) is obtained by expanding the nodal
5. COMMENTS ON THE FOREGOING APPROACH FOR displacements of member A: in a Taylor series about the
DEVELOPING CONTINUUM MODELS center of the repeating cell.
The following comments regarding the procedure for devel- 3. The geometric stiffness matrix of the continuum is given
oping continuum models seem to be in order: by
associated with the retained and neglected geometric stiff- ! Coupling terms included
FIG. 7. Vierendeel girder and its repeating cell (see Ref. 80).
Cross
section Length Designation
area
Longerons A L The error reduces to less than 1.5% when the number of bays
i
increases to 30. On the other hand, neglecting the coupling
Diagonals Ad d
terms results in overestimating the stiffness of the structure,
Battens A
b b thereby increasing the vibration frequencies and reducing the
maximum displacements. This effect is very pronounced for the
vibration frequencies associated with the higher modes. For the
b=L=5m 10-bay truss, when the coupling terms are neglected, the maxi-
A = 2.359x10~4m2 mum errors in the lowest 10 frequencies are over 65%.
The second set of results is for the bifurcation buckling of
P = 1743.8 kg/m3
axially-loaded spatial vierendeel girder with rigid joints shown
E = 2.482x10 Pa (longerons and battens) in Fig 7. The continuum model for this structure is a micropolar
= 2.557 x10 11
Pa (diagonals) beam whose stiffness characteristics are given in Ref 80. Figure
8 shows the lowest two bifurcation buckling modes of a 10-bay
FIG 5. Orthogonal tetrahedral beamlike lattice. lattice. The predictions of the micropolar beam model are given
» " 3 x X ) = 1.330 Nnwtons (1.330) K>"3 x X 5 = 1.807 Newtons (1.823) > (1.248) Sy'r- >2 = 2.064(2.087)
,K4^
\-
rrr—
A ; = 8 0 x 10 m . E - 7 1 7 x 1 0 " Pa L De S 1 g
"i
Ab = 50x10"6m2 T
—
A d -. 10 x 1 0 - 6 m 2 r
. ... I_i = 6.93 m
2
A v , , 0 » , 0 - 6 m „ t . 1 V P d - P v
L =7 5 m = 2768 k g - m 3
i i -1
l_2 = 8.5 m
h = 7.5 m
to =(25.916) : (36.966)
8. OTHEM REPORTED APPLICATIONS OF 21 25.710 Hz 22 37.152 Hz
CONTINUUM MODELING
FIG. 12. Vibration mode shapes for simply supported hexagons'
Among the reported applications of continuum modeling grid shown in Fig 11. Numbers in parentheses are exact frequencies
not discussed in the preceding sections are the following. (Hz); other numbers are continuum plate frequencies.
8.1. Geometrically nonlinear problems of beamlike lattices 9.1. Needed extensions and new applications
with pinjoints
The needed extensions include incorporating effects of joint
In Refs 1 and 17 incremental stiffness matrix of the equiv- flexibility and damping into the continuum model.
alent continuum is used to account for the large-displacement The new applications include:
nonlinear effects in beamlike lattices. Applications of this ap-
proach are made to static, large deflection analysis, large-ampli- 1. random periodic structures (structures with random varia-
tude free vibration analysis, as well as buckling and postbuck- tions in material properties, geometric characteristics, or
ling analysis. periodicity),
2. nonlinear structural dynamics and wave propagation
problems,
8.2. Beamlike lattices with material damping
3. multidisciplinary optimization problems (eg, simulta-
Planar beamlike lattices with members having viscous damp- neous optimization of structures and control systems).
ing are modeled as continuum Timoshenko beams in Refs 2 and
112. The global damping characteristics are determined in terms
of the damping coefficients and dimensions of the truss mem- 9.2. Computational strategies and modeling techniques
bers. With regard to computational strategies and modeling tech-
niques, the following two areas appear to have high potential
for research:
8.3. Stress, free vibrations, and buckling problems of beamlike
lattices with rigid and iexible joints
(a) Hybrid modeling approaches include combination of
Micropolar continuum models are developed in Refs 11-14, either (1) continuum/repetitive structure approach for
53, 54, 79, and 80 for the static, free vibrations, and buckling handling lattice structures with arbitrary boundary con-
problems of planar- and spatial beamhke lattices with rigid ditions, or (2) continuum/discrete structure approach
joints. The study made in Ref 80 shows that ordinary shear- for specialized problems such as stress concentration.
flexible continuum beam models are found to be adequate for (b) Application of operator splitting/reduction technique for
predicting the global response characteristics of lattices which analyzing complicated continuum models. Operator split-
do not need to have rigid joints for their kinematic stability. An ting/reduction techniques are very effective techniques
exception to this is the case when very small diagonal members for generating the response of a complex structure (eg,
are used. The incorporation of joint flexibility in the continuum anisotropic continuum model) using large perturbations
model is discussed in Ref 97. from the response of a simpler structure (eg, correspond-
ing orthotropic model).
8.4. Stress and free vibration problems of single-layer grids Reduction techniques have been shown to be related to the
Finite-difference calculus has been extensively used in preconditioned conjugate gradient iterative method. The pre-
analyzing the stress, free vibrations, and buckling problems of conditioning matrix is taken to be the matrix of the simpler
planar and cylindrical gridworks (see, for example, Refs 31, 32, structure; and the preconditioned residuals provide sensitivity
35, 36, 65, 89, 90, 101, 117-120, 125, and 126). Ordinary information of the response with respect to the complicating
continuum plate models are developed in Refs 43 and 73 for factors (eg, anisotropy).
predicting the response of single-layered grids with rigid joints.
The transverse motion of the grids are predicted accurately by
these models. On the other hand, the accurate prediction of the 10. CONCLUDING REMARKS
in-plane motion of single-layer grids with rigid joints requires
The status and some recent developments in continuum
the use of micropolar continuum (see Refs 11-14, 34, 37-39,
modeling of large repetitive lattice structures are summarized.
53, 54, 56, 110, 118, 125, and 126).
Discussion focuses on a number of aspects including definition
of an effective substitute continuum, characterization of the
continuum model, and the different approaches for generating
8.5. Use of continuum models for system identification and
the constitutive matrix and the matrix of mass densities of the
control law design of lattice structures
continuum. Also, a simple approach is presented for generating
In view of the significant reduction in the number of struc- the continuum properties. The approach can be used in con-
tural parameters obtained by replacing the actual lattice struc- junction with computerized symbohc manipulation to generate
ture by a continuum model, a number of recent studies have analytic expressions for the continuum properties, or in a small
performed the control and/or system identification analysis to FORTRAN program to determine the numerical values of these
either the governing partial differential equations, or the sys- properties. Sample numerical results are presented to demon-
tems of ordinary differential equations resulting from discretiza- strate the effectiveness of the continuum modeling approach.
tions of the simple continuum model (see Refs 10, 16, 52, 63, The status of continuum modeling for repetitive lattice struc-
and 129). tures can be summarized in the following:
Simplified continuum models (associated with known en- 17. Berry, D T. and Yang, T Y, Simplified lattice beam elements tor
gineering theoiies such as Timoshenko type beam theory geometrically nonlinear static, dynamic and postbuckling anah si.,
AIAA 24(8), Aug 1986, 1346-1347.
and Reissner-Mindlin type plate theoiy) are adequate for 18. Bibliography on latticed structures, J Struc Div ASCE 98(ST7), Jul
many practical problems. Continuum models associated 1972, 1545-1566.
with higher-order theories may be required for predicting 19. Bleich, F, and Melan, E, Die gewohnlichen und partiellen Difjeren-
zengleichungen der Baustatik, Springer, Berlin, 1927.
localized effects. 20. Borrego, J, Space grid structures—skeletal frameworks and stressed-skm
2. For lattices with rigid joints, ordinary continuum models systems, MIT Press, Cambridge, MA, 1968.
are available for the accurate prediction of the global 21. Brulin, O. and Hsieh, R K T, Continuum models of discrete systems 4,
in Proceedings of the fourth international conference on continuum oioJeh
response of lattices with rigid joints, provided that the of discrete systems, Stockholm, Sweden, 29 J u n - 3 Jul, 1981, North-Hol-
response is not dominated by local member deformation land, Amsterdam, 1981.
(eg, transverse motion of single-layered grids). Micropolar 22. Chen, C C and Sun, C T, Transient analysis of large frame structures
continuum models exist for predicting the in-plane mo- by simple models, J Astronaut Sci 31(3), 1983, 359-379.
23. Collings, A G, Analysis of some double layer grids, ./ Struc Din ASCE
tions of lattices with rigid joints. 102(ST8), Aug 1976, 1591-1608.
24. Dean, D L, On the statics of latticed shells, IABSE Publ, Zurich, V)f>S.
vol 25, pp 6 5 - 8 1 .
Needed extensions of continuum modeling include incorpo- 25. Dean, D L, and Ugarte, C P, Field solutions for two-dimensional
ration of joint flexibihty and damping. Also, application of frameworks, Int J Mech Sci 10, 1968, 315-339.
continuum modeling to multidisciplinary optimization prob- 26. Dean, D L, and letter, F R, Analysis for truss buckling, ./ Struc Die
ASCE 98(ST8), Aug 1972, 1893-1897.
lems need more development. These extensions are recom- 27. Dean, D L, and Avent, R R, State-of-the-art of discrete field analysis of
mended as future directions for research. New computational space structures, in Proceedings of the second international conference un
and modeling strategies that have high potential for use in space structures, W J Supple, Ed, Univ of Surrey, Guildford, UK.. Sep
1975, pp 7-16.
conjunction with continuum modeling include hybrid con- 28. Dean, D L, Discrete field analysis of structural systems, Springer, New
tinuum/repetitive structure; hybrid continuum/discrete struc- York, 1976.
tures approaches; and the application of operator splitting/ 29. Dow, J O, Su, Z W , Feng, C C, and Bodley, C, Equivalent continuum
reduction technique for analyzing complicated continuum representation of structures composed of repeated elements, AIAA J
23(10), Oct 1985, 1564-1569.
models. 30. Dow, J O and Huyer, S A, An equivalent continuum analysis procedure
for space station lattice structures, Part I—A collection of technical
papers, A I A A / A S M E / A S C E / A H S 28th structures, structural dy-
namics and materials conference, Apr 1987, Monterey, CA, pp 110-122.
ACKNOWLEDGMENT
31. Ellington, J P, and McCallion, H, Moments and deflections of a
simply-supported beam grillage, Aeronaut Q 8(part 4), Nov 1V57.
The present work is supported by a NASA grant no NAG1- 360-368.
740. 32. Ellington, J P, and McCallion, H, The free vibrations of grillages. ./
Appl Mech 26, Trans ASME Ser E 81, Dec 1959, 603-607.
33. Flower, W R and Schmidt, L C, Analysis of space truss as equivalent
REFERENCES plate, J Struc Div ASCE 97(ST12), Dec 1971, 2777-2789.
34. Frackiewicz, II, Buckling of lattice plates. Arch Mech Stosowanej 23i4).
1. Abrate, S, and Sim, C T, Dynamic analysis of geometrically nonlinear 1971,495-515.
truss structures, Comput Struct 17(4), 1983, 491-497. 35. GangaRao, H V S, and Smith, J C, Dynamic field analysis of torsion-
2. Abrate, S, and Sun, C T, Continuum modeling of damping in large less grids, J Eng Mech Div ASCE 98(EM3), Jun 1972, 679-693.
space structures, in Proceedings of the second international conference on 36. GangaRao, H V S, Vibration analysis of grid-works, Shock Vib Dig
recent advances in structural dynamics, M Petyt and H F Wolfe, Eds, vol 8(9), Sep 1976, 25-30.
II, Inst of Sound and Vib Res, Univ of Southampton, Apr 1984, pp 37. Gutkowski, W, Statics and stability of prismatic frame-lattice shells.
877-885. Bull Acad Pol Sci Ser Sci Tech IX(5), 1961, 297-303.
3. Abrate, S, Continuum modeling of lattice structures, Shock Vib Dig 38. Gutkowski, W, Cylindrical grid shells, Bull Acad Pol Sci Ser Sci Tech
17(1), Jan 1985, 15-21. XIII(3), 1965, 133-141.
4. Anderson, M S, Buckling of periodic lattice structures, AIAA J 19(6), 39. Gutkowski, W, Mechanical problems of elastic lattice structures. Pro-
Jun 1981, 782-788. gress in aerospace sciences, Vol. 15, D Kuchemann et al, Eds. Per-
5. Anderson, M S, Vibration of prestressed periodic lattice structures, gamon, New York, 1977, pp 181-216.
AIAA J 20(4), Apr 1982, 551-555. 40. Haug, E, Dowlatyari, P, and Dubois, J, Equivalent continuum analysis
6. Anderson, M S, and Williams, F W, Natural vibration and buckling of of large space trusses, in Spacecraft structures, Proc. of an international
general periodic lattice structures, AIAA J 24(1), Jan 1986, 163-169. conference held at Centre National d'Etudes Spatiales, Toulouse.
7. Anderson, M S, and Williams, F S, BUNVIS-RG: Exact frame buck- France, Dec 1985, European Space Agency, ESA SP-238, Apr 1986, pp
ling and vibration program, with repetitive geometry and substructur- 57-62.
ing, J Spacecraft Rockets 24(4), Jul/Aug 1987, 353-361. 41. Hefzy, M S, and Nayfeh, A H, Shear deformation plate continua of
8. Artola, M, and Duvaut, G, Elasticite homogeneisation d'une plaque large double layered space structures, Int J Solids Struc 22(12), l u 86.
renforcee, C R Seances Acad Sci 284A(12), 21 Mar 1977, 707-710. 1455-1469.
9. Babuska, I, Homogenization and its application. Mathematical and 42. Heki, K and Fujitani, Y, The stress analysis of grids under the action
computational problems, in Numerical solution of partial differential of bending and shear, in Space structures, R M Davies, Ed, Blacks ell,
equations—III, SYNSPADE 1975, B Hubbard, Ed. Academic, New Oxford and Edinburgh, 1967, pp 33-43.
York, 1976, pp 89-116. 43. Heki, D, On the effective rigidities of lattice plates, in Recent resean lies
10. Banks, FI T, and Crowley, J M, Parameter identification in continuum of structural mechanics, H Tanaka and S Kawamata, Eds, Uno Shoten,
models, NASA CR-172132, May 1983. Tokyo, 1968, pp. 31-46.
11. Bazant, Z P, Micropolar medium as a model for buckling of grid 44. Heki, K, and Saka, T, Stress analysis of lattice plates as anisotropic
frameworks, Developments in mechanics, proc. of the 12th midwestern continuum plates, in Proceedings of the 1971 Pacific symposium, Pan II
mechanics conference, University of Notre Dtmie, Notre Dame, IN, on tension structures and space frames, Tokyo and Kyoto, Paper no 7-7,
16-18 Aug 1971, vol 6, pp 587-593. 1972, Architectural Inst of Japan, Tokyo, pp 663-674.
12. Bazant, Z P, and Christensen, M, Analogy between micropolar con- 45. Heki, K, The effect of shear deformation on double layer lattice plates
tinuum and grid frameworks under initial stress, Int J Solids Struc, 8, and shells, in Second international conference on space structures, Dcpt
Mar 1972, 327-346. of Civil Engineering, University of Surrey, Guildford, UK, Sep 1975,
13. Bazant, Z P, and Christensen M, Long-wave extensional buckling of pp. 189-198.
large regular frames, J Struc DivASCE 98(ST10), Oct 1972, 2269-2289. 46. Heki, K, The effective strength of double-layer grids in continuum
14. Bazant, Z P, and Christensen, M, Continuum solutions for long-wave treatment, in IASS world congress on shell and spatial structures, Madrid.
extensional buckling of regular frames, Int J Eng Sci 11, 1973, Spain, Sep 24-28, 1979, vol 2, pp 4.123-4.137.
1255-1263. 47. Heki, K, Examples of effective rigidities of lattice plates, appendix of
15. Benoussan, A, Lions, J L, and Papanicolaou, G, Asymptotic analysis for Analysis, design and realization of space frames, Bull IASS 84/ <V5.
periodic structures, North-Holland, Amsterdam, 1978. 1984, 96-106.
16. Berry, D T, Yang, T Y, and Skelton, R E, Dynamics and control of 48. Heki, K, The continuum treatment of discrete systems and its applica-
lattice beams using simplified finite element models, J Guidance Con- tion to the analysis of lattice structures, Bull IASS, XXVI-3(89), 19x5.
trol Dyn 8(5), Sep-Oct 1985, 612-619. 19-26.
49. Heki, K, The continuum analogy of lattice structures from the view- a/ialysis and optimization issues for large space structures, Larry D
point of difference equations, in Shells, membranes and space frames, Pinson, Anthony K Amos, and V B Venkayya, compilers, NASA
Proc. IASS symposium on membrane structures and space frames, CP-2258, 1983, pp. 5-32.
Osaka, Japan, K. Heki, Eds, Elsevier 1986, vol 3, pp 65-72. 84. Noor, A K, and Peters, J M, Recent advances in reduction methods for
50. Hrennikoff, A, Solution of problems of elasticity by the framework instability analysis of structures, Comput Struct 16(1-4), Jan 1983,
method, J Appl Mech, Dec 1941, A-169-A-175. 67-80.
51. Hussey, M J L, Tarzi, A I, and Theron, W F D, Simply supported 85. Noor, A K, and Russell, W C, Anisotropic continuum models for
rectangular double-layered grids, ./ Struc Din ASCII 97(ST3), Mar beamiike lattice trusses, Comput Methods Appl Mech Eng 57(3), Sep
1971, 753-764. 1986, 257-277.
52. Juang, J N, and Sun, C T, System identification of large flexible 86. Noor, A K, and Mikulas, M M, Jr, Continuum modeling of large lattice
structures by using simple continuum models, J Astronaut Set 31(1), structures — status and projections, NASA TP-2767, Feb 1988.
1983, 77-98. 87. Proven, J W, Ed, Continuum models of discrete systems, Proceedings of
53. Kanatani, K -I, A micropolar continuum model for vibrating grid the second international symposium on continuum models of discrete
frameworks, hit J Eng Set 17, 1979, 409-418. systems, Mont Gabriel, Quebec, Canada, June 26-July 2, 1977 Univ of
54. Kim, K S and Piziali, R L, Continuum models of materials with beam Waterloo Press, Waterloo ON, Canada, 1978.
microstructure, hit J Solids Struct 23(11), 1987, 1563-1578. 88. Pshenitchnov, G I, Analysis of cylindrical gridwork shells, Izd Akad
55. Kleiber, M, The approximate methods in the theory of elastic lattice- Nauk SSSR, Moscow, 1961 (in Russian).
type shells. Arch Mech 25(2), 1973, 195-211. 89. Renton, J D, A finite difference analysis of the flexural-torsional behav-
56. Kleiber, M, Statics of elastic lattice-type shells. Arch Mech 25(2), 1973. ior of grillages, hit J Mech Sci 6, 1964, 209-224.
179-194. 90. Renton, J D , On the gridwork analogy of plates, J Mech Phvs Solids
57. Kollar, L, and Hcgedus, I, Solution of double-layer space trusses of 13, 1965, 413-420.
general triangular grid by the equivalent continuum method. Acta Tech 91. Renton, J D, On the analysis of triangular mesh grillages, hit J Solids
Acad Set Hung 74(3-4), 1973, 363-381. Struct 2, 1966, 307-318.
58. Kollar, L, Continuum method of analysis for double-layer space trusses 92. Renton, J D, The related behaviour of plane grids, space grids and
with upper and lower chord planes of different rigidities, Acta Tech plates. Space structures, R M Davies, Ed, Blackwell, Oxford and Edin-
Acad Sci Hung 76(1-2), 1974, 53-63. burgh, 1967, pp 19-32.
59. Kollar, L, Analysis of double-layer space trusses with diagonally square 93. Renton, J D, Behavior of Howe, Pratt, and Warren trusses, J Struct
mesh by the continuum method, Acta Tech Acad Sci Hung 76(3-4), Div ASCE 95 (ST2), Feb 1969, 183-202.
1974, 273-292. 94. Renton, J D, General properties of space grids, hit J Mech Sci 12,
60. Kollar, L, and Hegedus, I, Analysis and design of space frames by the 1970, 801-810.
continuum method, Elsevier, Amsterdam, 1985. 95. Renton, J D, The formal derivation of simple analogies for space
61. Kollar, L, Continuum method of analysis for double-layer space trusses frames, in IASS symposium on tension structures and space frames,
of "hexagonal over triangular" mesh. Acta Tech Acad Sci Hung K6I1-2), Kyoto, Japan, 1971. Y. Yokoo, et al, Eds, Architectural Institute of
1978, 55-77. Japan, Tokyo.
62. Kroner, E, Anthony, K-H, and Leipholz, H H E, Continuum models of 96. Renton, J D, The beamlike behavior of space trusses, AIAA J 22(2),
discrete systems (CMDS3), in Proceedings of the third international Feb 1984, 273-280.
symposium on continuum models of discrete systems, Frcudenstadt, FRG, 97. Russell, W C, Continuum modeling theories for repetitive lattice space
Jun 1979, Univ of Waterloo Press, Waterloo ON, Canada, 1980. structures, MS thesis, George Washington Univ, Oct 1985.
63. Lamberson, S E, and Yang, T Y, Integrated design of space structures 98. Saka, T, and Heki, K, On the effective rigidities of b-ariike plane
using lattice plate finite elements, ./ Guidance Control Dyn 9(4), Jul-Aug trusses, in Mem Fac Eng Osaka City Univ 22 Dec 1981, 167-173.
1986, 478-484. 99. Schnobrich, W C, and Melin, J W, Model analogue for cylindrical
64. Latticed structures: State-of-the-art report, J Struc Die ASCE shells, in Proceedings of the world conference on shell structures, Nat Res
102(ST11), Nov 1976, 2197-2230. Council publ no 1187, San Francisco, 1962.
65. Lederer, F, Framework systems and gridwork. Space structures, R M 100. Shilin, D, and Hengxi, X, The solutions and their comparisons of the
Davies, Ed, Blackwell, Oxford and Edinburgh, 1967, pp 154-159. sandwich plate analogy method in two kinds of simply supported
66. Lutoborski, A, Homogenization of linear elastic shells, J Elast 15(1), two-way orthogonal-diagonal lattice grids, in Shells, membranes and
Mar 1985, 69-87. space frames, Proc. IASS symposium, Osaka, Japan, K Heki, Ed,
67. McDaniel, T J, and Chang, K J, Dynamics of rotationally periodic Elsevier, Amsterdam, 1986, pp 33-48.
large space structures, J Sound and Vib 68(3), Feb 1980, 351-368. 101. Smith, J C, and Dean, D L, Field analysis of torsionless bridge grids, J
68. McHenry, D, A lattice analogy for the solution of plane stress prob- Struct Div ASCE 99(ST9), Sep 1973, 1837-1849.
lems, J hist Civil Eng 21, 1943, 59-82. 102. Soare, M, Application des equations aux differences fmies au calcul des
69. Mitchell, L H, A shell analogy for framed domes, Sm note 208, reseaux spatiaux planaire carres, Rev Roum Sci Tech Mec Appl 14(3),
Aeronaut Res Labs, Melbourne, Australia. 1969, 595-628.
70. Nayfeh, A H, and Hefzy, M S, Continuum modeling of three dimen- 103. Soare, M, The cylindrical bending of double-layer diagonal square
sional truss-like space structures, AIAA J 16(8), Aug 1978, 779-787. mesh grids, Acta Tech Acad Sci Hung 69(3-4), 1970, 433-466.
71. Nayfeh, A H, and Flefzy, M S, Continuum modeling of the mechanical 104. Soare, M, Application of the equivalent continuum method to the
and thermal behavior of discrete large structures, AIAA J 19(6), Jun analysis of double-layer parallel square mesh grids, Bui Stiint Inst
1981, 766-773. Construct Bucuresti 14(1-2), 1971, 251-269.
72. Nayfeh, A H, and Hefzy, M S, Effective constitutive relations for large 105. Soare, M, Statica si dynamica retelelor spatiale planar patrate simple
repetitive frame-like structures, hit J Solids Struct 18, 1982, 975-987. prin aplicarea metodei mediului continuu echivalent, Stud Cercetari
73. Nemeth, M P, Continuum models for repetitive lattice structures with Mech Apl 31(3), 1972, 673-702.
rigid joints, MS thesis, George Washington Univ, Sep 1979. 106. Spencer, A J M, Ed, Continuum models of discrete systems. Proceedings
74. Newmark, N M, Numerical methods of analysis in bars, plates and of the fifth international symposium, Nottingham, UK, Jul 1985, Else-
elastic bodies, in Numerical methods of analysis in engineering, L E vier-North-Holland, Amsterdam, 1985.
Grinter, Ed, Macmillan, New York, 1949. 107. Sumec, J, Regularne mriezkove dosky a skrupiny (Regular lattice plates
75. Noor, A K, Thermal stress analysis of double-layered grids, J Struc Die and shells), 1st ed, VEDA Pub, Bratislava, Czechoslovakia, 1984.
ASCE 104(ST2), Feb 1978, 251-262. 108. Sumec, J, and Zingali, A E, A smdy of the influence of initial shape
76. Noor, A K, and Andersen, C M, Computerized symbolic manipulation imperfections on the stability of lattice shells by direct and shell
in structural mechanics — Progress and potential, Comput Struct analogy methods, Space Struct 2(4), 1986/87, 223-230.
10(1/2), Apr 1979, 95-118. 109. Sun, C T, and Yang, T Y, A continuum approach toward dynamics of
77. Noor, A K, Anderson, M S, and Greene, W H, Continuum models for gridworks, ,/ Appl Mech, Trans ASME Ser E 40(1), Mar 1973, 186-192.
beam- and platelike lattice structures, AIAA J 16(12), Dec 1978, 110. Sun, C T, and Yang, T Y, A couple-stress theory for gridwork-rein-
1219-1228. forced media, J Elast 5(1), Mar 1975, 45-58.
78. Noor, A K, and Andersen, C M, Analysis of beam-like lattice trusses, HI. Sun, C T, Kim, B J, and Bogdanoff, J A, On the derivation of
Comput Methods Appl Mech Eng 20(1), Oct 1979, 53-70. equivalent simple models for beam- and platelike structures in dynamic
79. Noor, A K, and Nemeth, M P, Micropolar beam models for lattice- analysis, in A collection of technical papers, part 2, A I A A / A S M E /
grids with rigid joints, Comput Methods Appl Mech Eng 21(2), Feb A S C E / A H S 22nd structures, structural dynamics and materials con-
1980, 249-263. ference, and AIAA dynamics specialists conference, 1981, pp 523-532.
80. Noor, A K, and Nemeth, M P, Analysis of spatial beamlike lattices 112. Sun, C T, and Juang, J N, Modeling global structural damping in
with rigid joints, Comput Methods Appl Mech Eng 24(1), Oct 1980, trusses using simple continuum models, in A collection of technical
35-59. papers part 1, A I A A / A S M E / A S C E / A H S 24th structures, structural
81. Noor, A K, and Weisstein, L S, Stability of beamlike lattice trusses, dynamics and materials conference. Lake Tahoe, NV, May 1983, pp
Comput Methods Appl Mech Eng 25(2), 1981, 179-193. 722-729.
82. Noor, A K, On making large nonlinear problems small, Comput Meth- 113. Sun, C T, and Kim, B J, Continuum modeling of periodic truss
ods Appl Mech Eng 34, 1982, 955-985. structures, in Damage mechanics and continuum modeling, N Stubbs
83. Noor, A K, Assessment of current state-of-the-art in modeling tech- and D Krajcinoric, Ed, ASCE, New York, 1985, pp 57-71.
niques and analysis methods for large space structures, in Modeling, 114. Suzuki, E, Kitamura, H, and Yamada, M, The analysis of the space
truss plate by dift'erence equations, Space structllres, R M Davies, Ed, rotationally periodic structures, Int J IVumer Methods Eng 23(4).
Blackwell, Oxford and Edinburgh, 1967, pp 136-144. 1986, 609-622.
115. Symposium on continuum models of discrete systems, Jodlowy Dwor, 124. Williams, F W, Exact eigenvalue calculations for structures with
Jun 1975, Arch Mech Stos. 28(3), 1976, pp 255-581. tionally periodic SUbstructures, Int J IVllmer Methods Eng 23(4) Apri!
116. Timoshenko, S P, and Gere, J M, Theory of elastic stability, 2nd ed, 1986, 695- 706.
McGraw-Hill, New York, 1961. 125. Wozniak, Cz, Load-carrying structures of the dense lattice type.
117. Wah, T, Natural frequencies of uniform grillages, J Appl Mech, Trans plane problem, Arch Mech Stosowanej 5(18),1966,581-597.
ASME, 30, Dec 1963, 571-578. 126, Wozniak, Cz, Bending and stability problems of plates with latllcc
118. Wah, T, Analysis of laterally loaded gridworks, J Eng Mech Diu structure, Arch JlIech Stosowanej 6(18),1966781-796.
ASCE, 9O(EM2), Apr 1964, 83-106. 127. Wright, D T, Membrane forces and buckling in reticulated shells,
119. Wah, T, The buckling of gridworks, J Mech Phys Solids 13, 1965, 1-16. Stract Diu ASCE 91(STl), Feb 1965, 173-20l,
120. Wah, T, Vibration of cylindrical gridwork shells, AJAA J 3(8), Aug 128. Wright, D T, A continunm analysis for double-layer space frame
1965,1467-1475. Publ Int Assoc Bridge and Struct Eng, Zurich, vol 26, 1966, pp 593 -6l{).
121. Wah, T, and Calcote, L R, Strllctural ana(vsis by .finite difference 129. Yang, T Y, and Lamberson, S E, Optimization using lattice
caleuills, Van Nostrand Reinhold, New York, 1970. clements for feedback control of space structures, in A
122. Weisstein, L S, Introduction and survey on continuum models' [or technical papers, pari 2, AIAA/ASME/ASCE/AHS structures,
repetitive lattice structures, in Proceedings of the workshop on applica- tural dynamics and materials conference, Orlando FL, Apr 1985.
tions of distributed system theory to the control of large space structures, 743- 750,
G Rodriguez, Ed, JPL publ 83-46, California Inst Tech, Jul 1983, pp 130. Yettram, A L, and Husain, H M, Grid-framework method for
63-70, t1exure, .I Eng Mech Diu, Proc ASCE, Jun 1965 91( Elvf3)
123. Williams, F W, An algorithm for exact eigenvalue calculations for
Continuum modeling facilitates the design of feedback control systems in lattice structures by focusing on the global behavior of the lattice rather than the detailed deformation in each member. This simplification allows for a more efficient analysis and design of control systems, as it compresses complex lattice behavior into a manageable continuum representation suitable for parameter and system identification .
Recent developments in continuum models address the challenge of axial stiffness overestimation by incorporating the local deformations of the lattice structures. Early models that averaged deformations tended to overestimate axial stiffness. The inclusion of local deformations allows these models to more accurately reflect the behavior of each repeating cell, particularly in lattices with more than one bay, thereby providing a realistic representation of stiffness .
The continuum modeling approach does not predict local member instability or snap-through buckling because it assumes that individual members of the lattice remain straight and stable during buckling. This limitation arises from the model's focus on global stability analysis, which does not account for specific joint or localized deformation behaviors that lead to local instability .
In using transformation matrices to develop the characteristics of continuum models, it is essential to consider how these matrices relate nodal displacements of lattice members to continuum displacement parameters. These transformations allow for the translation of discrete lattice behavior into continuous forms appropriate for analysis. For more complex configurations, employing numerical methods to generate transformation matrices ensures accurate reflection of the lattice structure's mechanical behavior .
The concept of an equivalent continuum is crucial for the static and dynamic analysis of space truss structures as it enables the approximation of the structures' behavior using continuum mechanics principles. By treating lattice space trusses as equivalent continua, one can simplify the system into beamlike or platelike models that capture the essential mechanical behaviors, such as vibrations and buckling, allowing for systematic analysis and result predictions .
In continuum models for beamlike and platelike lattice structures, the kinetic energy density is expressed in terms of the velocity of the continuum model's points, represented by a derivative with respect to time. The continuum approach simplifies the complex variable behavior of lattice structures into a framework where their kinetic energy can be systematically analyzed .
Higher-order continuum models enhance the prediction of displacement and stress distributions in lattice structures by incorporating additional strain and displacement parameters beyond the basic set ({e(1} and {u d}). These added parameters allow for a more detailed and accurate representation of the behavior captured in the lattice, addressing complexities that cannot be sufficiently captured by lower-order models .
Substitute continuum models differ from traditional continuum models by accounting for local deformations and appropriately modeling joint behavior based on the type of joints present. Traditional models tended to average deformations, which could lead to overestimations of axial stiffness, whereas recent models factor in the complexities of local deformations for more accuracy . Furthermore, pin-jointed lattices are modeled as classical continua, while those with rigid or flexible joints are modeled as micropolar continua .
In the stability analysis of lattice structures using continuum models, the geometric stiffness matrix is crucial as it provides the geometric stiffness coefficients of the continuum. This matrix is obtained by equating the potential energy due to initial stresses of the original lattice to that of the continuum. The matrix is computed by expanding the nodal displacements of the lattice members in a Taylor series about the center of their repeating cell and combining transformations for each member .
Numerical results have shown the effectiveness of continuum models in simulating lattice structures by accurately capturing static, free vibration, and buckling responses. For instance, Timoshenko-type beam models applied to lattice structures like orthogonal tetrahedral trusses demonstrate accurate vibrational characteristics and stability predictions matching physical structures, thus validating the continuum approach as a reliable simulation tool .