Static and Dynamic Responses of Micro-Structured Beams

In this study, we developed a one-dimensional Timoshenko beam model, embedded in a 3D space for static and dynamic analyses of beam-like structures. These are grid cylinders, that is, micro-structured bodies, made of a periodic and specifically designed three-dimensional assembly of beams. Derivation is performed in the framework of the direct 1D approach, while the constitutive law is determined by a homogenization procedure based on an energy equivalence between a cell of the periodic model and a segment of the solid beam. Warping of the cross-section, caused by shear and torsion, is approximatively taken into account by the concept of a shear factor, namely, a corrective factor for the constitutive coefficients of the equivalent beam. The inertial properties of the Timoshenko model are analytically identified under the hypothesis, and the masses are lumped at the joints. Linear static and dynamic responses of some micro-structured beams, taken as case studies, are analyzed, and a comparison between the results given by the Timoshenko model and those obtained by Finite-Element analyses on 3D frames is made. In this framework, the effectiveness of the equivalent model and its limits of applicability are highlighted.


Introduction
Recent advances in 3D printing techniques and associated micro-fabrication technologies have made (periodic) micro-structured systems and architected the material object of study in several fields of engineering, ranging from aerospace to biomechanics [1]. Interest in research is particularly focused on some peculiar characteristics of such objects, far exceeding those of conventional bulk materials, like a very low lightweight-stiffness ratio, accompanied with significant mechanical properties, such as resilience and impact resistance, which usually allow for important energy absorption capability, see for example, [2,3], also through micro-scale buckling phenomena. A unit-cell of such micro-structured bodies, whose characteristic dimension is of the order of some millimeters, consists of a lattice of fibers, which interacts at the internal connections. The intrinsic extensional, shear, flexural, and torsional stiffnesses of the fibers, together with the topology of the lattice, is able to confer to the resulting structure, which is the above-mentioned properties at the macro-scale.
The mixed up-scaling technique, adopted in the present work, has also been successfully applied to model tower buildings and multi-storey frames as generalized (continuous) beams in [26][27][28][29][30][31][32][33][34][35], to describe their overall mechanical behavior in statics [28,30,34,35], dynamics [28,33,35], buckling [29,32], and aeroelasticity [26,27,31] with a lower computational effort, with respect to, for example, Finite-Element models of the correspondent fine structures. These beam-like structures have been macroscopically modeled as shear-shear-torsional beams, that is, where macro-bending has been neglected, or Timoshenko beams, by identifying the floors and the columns of the periodic frames, with the cross-sections and longitudinal fibers of the continuous beam, respectively. In a subset of them [26][27][28][29][30]32,33], identification has been carried out by assuming the floor could behave as a rigid body, both in its plane and out-of-plane; whereas in [25,34,35], this hypothesis has been partially removed and out-of-plane deformability of the floor has been taken into account on an energy average ground. In particular, in [26,27], an equivalent shear-shear-torsional beam model embedded in a 3D-space has been developed to describe the nonlinear aeroelastic behavior of tower buildings. In [30], a rational justification of the main assumptions, heuristically adopted in [27], is provided. In this framework, the limits of applicability of the model have been discussed in terms of nonlinear static response of some buildings. In [33], under a small displacements assumption, the same model was adopted to analyze linear free and forced dynamics of multi-storey buildings, for which both external and internal (proportional) damping forms have been considered, thus allowing a discussion of the damping features of the shear beam model. In [28], a Timoshenko equivalent beam model, embedded in a 3D-space, was developed to also account for macro-bending of multi-storey buildings in linear statics and dynamics. In [29,32], a generalization of this latter model was provided for buckling analyses, including for the geometric effects both in the equilibrium equations and in the constitutive law-in this context, the effect of uniform compression ( [29]), as well as that of non-uniform compression and the interactions with elastic soil ( [32]), are discussed. In [34], by limiting the analysis to planar frames, it has been shown that it is possible to consider the out-of plane true flexibility of the cross-section, even when rigid cross-section 1D (unwarpable) models, such as the shear and the Timoshenko ones, are heuristically adopted. Indeed, corrective factors have been analytically derived in [34] for the case of geometrically regular planar frames, to suitably reduce the shear and flexural stiffnesses of the macro-model by accounting on (energy) averages for the neglected warping. This analysis has been extended in [35], where an energy-based numerical algorithm, sided on a Finite-Element (FE) analysis of a single cell, has been again applied to planar frames, in order to detect the constitutive properties of the equivalent beam model. Finally, in [25], a review of these works was presented, together with a systematization of the homogenization procedure there adopted.
In this paper, the analysis so far carried out in [25,34,35] is extended to the spatial behavior of micro-structured beams, here referred to as grid beams, which are namely cylinders made of a cubic micro-structure [36], embedded in a 3D space, consisting of a pattern of three orders of orthogonal micro-beams periodically arranged. In particular, an equivalent 1D Timoshenko beam model is formulated through a direct approach. Moreover, the out-of-plane deformability of the cross-section, induced by shear and torsion, is approximatively taken into account by introducing corrective factors for the constitutive coefficients of the equivalent beam model. Attention is focused on statics and free dynamics of such structures. It is worth noticing that the main goal of the paper relies on the interpretation of the mechanical behaviour of such a complex micro-structured and three-dimensional object, and not on the algorithmic formulation. In particular, it is of interest to prove that the observation on a large scale of the mechanics of a periodic system can be pursued (with a good approximation) by using the classical tools of Continuum Mechanics.
The paper is organized as follows. In Section 2 the equations of the Timoshenko beam model, embedded in a 3D space, are recalled. In Section 3 the identification procedure is described. In Section 4, comparisons between results provided by Finite-Element models of grid beams and the equivalent Timoshenko model, for which exact solutions are viable, are developed. In Section 5, some conclusions and perspectives of this work are drawn.

Model
The grid beam consists of a grid cylinder, that is, in a 3D frame, made of a periodic and specifically designed three-dimensional assembly of equal micro-beams. It is made of micro-beams, here referred to as fibers, forming a periodic cubic micro-structure, of side h (Figure 1a), whose pattern is repeated n x = n times in x-direction. Two families of fibers are detected: (i) those parallel to the beam x-axis, referred to as longitudinal fibers, which are n y in y-direction and n z in z-direction; (ii) those lying in the n + 1 planes orthogonal to the axis-line, referred to as transverse fibers, parallel to the (y, z)-plane. These latter are arranged in a square orthogonal pattern that is reinforced through bracing diagonal elements, which guarantee the indeformability of the cylinder in the (y, z)-plane. It is assumed that the (x, y)and (x, z)-plane are of symmetry for the grid beam, and the origin of the coordinate system is at the lower end A (see Figure 1). All the fibers are modeled as Timoshenko micro-beams, whose cross-section properties are taken as equal in the two principal inertia planes, where square, circular, and ring cross-sections, or analogous, are considered. They have axial, shear, bending, and torsional stiffnesses EA, GA * , EI, and GJ, respectively. An internally unconstrained model of beam, embedded in a 3D space, see, for example, [37][38][39][40][41][42], is formulated (Figure 1b). The beam is made of material points, equipped with orientation, which accounts for the attitude of the cross-section. It is referred to as a material abscissa x ∈ [0, ], spanning the interval (A, B) of total length := n h. The following strain-displacement relationships hold: which link the longitudinal displacement u (x, t), the transverse displacements v (x, t), w (x, t), the in y-direction and z-direction, respectively, and the rotations of the cross-section θ ν (x, t), with respect to ν-axis (ν = x, y, z), to the elongation ε (x, t), the shear strains γ y (x, t), γ z (x, t), the torsional curvature κ x , and the flexural curvatures κ y (x, t), κ z (x, t); here, t is the time, and a prime denotes space differentiation.
The following geometric boundary conditions are considered in what follows (subscripts A, B denote evaluation at the ends x = 0 and x = , respectively):

•
Clamped at A and free at B: Simply-supported and torsionally restrained in A, B: The equilibrium is ruled by the following equations: where N (x, t), T y (x, t) and T z (x, t) are the axial and shear internal forces, respectively, M x (x, t) is the torsional moment, M y (x, t) and M z (x, t) are the bending moments, p ν (x, t), c ν (x, t) are the external forces and couple per unit-length (ν = x, y, z). Moreover, m and I ν are the mass and second inertia moments per unit-length, evaluated with respect to the centroid (ν = x, y, z), and a dot denotes time-differentiation. When no external load and no lumped masses are present at the ends of the beam, the mechanical boundary conditions read: • Clamped at A and free at B: Simply supported and torsionally restrained in A, B: It is assumed that the beam material is hyperelastic, so that the existence of a strain energy function, quadratic in the strains, is postulated, namely: where the elastic constants c ii (i = 1, . . . , 6) must be identified according to the micro-structure of the grid beam. The constitutive law follows from the Green law as: Then, equilibrium (2) is written in terms of displacements, and by using Equations (4) and (1), the elasto-dynamic problem is obtained: Equation (5)-a is a system of partial differential equations, sided by geometric and mechanical boundary conditions at A and B, Equation (5)-b,c, and by initial conditions, prescribing that the system is initially at rest (not reported here). Moreover, where M is the mass matrix, K i (i = 0, 1, 2) are symmetric or skew-symmetric stiffness matrices, and u, p are displacement and load vectors, respectively. Boundary operators, involved in Equation (5)-b,c, are defined as: • Clamped at A and free at B: • Simply-supported and torsionally restrained in A, B: It is worth noticing that in the static case, that is, when inertial effects are disregarded, the elastic problem (5) consists in a system of ordinary differential equations and boundary conditions in the unknown displacement field. They can be easily solved in closed-form due to the uncoupled constitutive law, entailing that the equations of the equivalent model are the classical (partially uncoupled) equations of the Timoshenko beam model [43]. Moreover, the free dynamics problem associated with Equation (5) also admits a closed-form solution, since the mass matrix is diagonal; see [43][44][45].

Identification Algorithm for the Elastic Constants
The identification of the elastic constants follows the lines of [25], according to which a suitable energy equivalence between the (fine) 3D grid beam model and the (coarse) 1D Timoshenko beam model is enforced (time is understood in what follows). The analysis here presented extends to the 3D case, the results of [25], where a planar grid beam was considered.
The first step of the identification consists in defining a discrete map, which links the displacements of the fine model to that of the coarse one, at the selected abscissas x i := i h (i = 0, 1, . . . , n), h being the period of the grid beam. At each x i , the fine model possesses a cross-section D i occupying a domain of the (y, z)-plane, made of the nodal points of the (y, z)-fibers; two adjacent cross-sections D i , D i+1 bound a cell. In order to define the displacement map, a kinematic constraint is enforced-the sampled cross-section D i remains plane, while the displacements inside the cells are unconstrained. The rigid displacement of D i (the tilde denoting, from now on, a quantity relevant to the 3D body), can be written as: The second step requires enforcement of the equivalence of the energies stored by a 3D cell and a segment of equal length of the 1D beam, when the two models undergo the same displacements at the ends. This (so-called) cell analysis leads to the evaluation of the averaged energy density of the 3D cell, namelyφ :=Ũ/h, withŨ being the elastic energy of the 3D cell, as a function of the configuration variables of the 1D beam (according to Equation (9)). The energy density of the cell reads: The last step of the identification consists in a comparison between the energy densitiesφ and φ, in Equations (3) and (10), respectively. To this end, the strain-displacement relationships (1) must be integrated in the (x i , x i+1 ) interval, thus giving a suitable strain test field. By taking the (simplest) constant strain field, the end displacements (rigid motions removed) read: Once Equation (11) is substituted in Equation (10), by the equality φ =φ ∀ ε, γ y , γ z , κ x , κ y , κ z , the elastic constants C := diag (c ii ), i = 1, . . . , 6 follow. It is important to remark that the cell's boundary displacements, defined in Equation (11), can be read as the result of a superposition of six independent deformation modes, namely: • An extensional mode (EX), in which D i+1 translates axially; • Two shear modes (SH y , SH z ), in which D i+1 translates transversely along the yand z-directions, respectively; • two flexural modes (FL y , FL z ), in which D i+1 rotates around the yand zaxes, respectively, and translates transversely along the zand yaxes; • A torsional mode (TO), in which D i+1 rotates around the xaxis.
Moreover, it is worth noticing that the evaluation of the cell's elastic energy is the most difficult step in the identification procedure and, in general, it cannot be pursued in closed form except for few cases, where for example, the fibers are arranged in a simple and regular pattern and/or the cross-section is assumed to be rigid (see [28,32,34]). On the contrary, a numerical identification procedure, such as that based on a Finite Element analysis, is always viable, as done in [35].

Analytical identification of elastic constants of the grid beam
The illustrated algorithm is applied to the grid beam. A cell made of n y × n z longitudinal fibers, comprised between two adjacent cross-sections D i , D i+1 , is considered. The fibers belonging to the cross-sections are taken of half of their stiffnesses, and only their nodal points are assumed to remain on the same plane, as displayed in Figure 2. However, the same fibers are free to warp around the plane-this deformation is here referred to as micro-warping.
The previously introduced deformation modes, shown in Figure 2, are then assigned to the cell under the following assumptions: Hypothesis 1. In the shear modes SH y , SH z all nodes rotate of the unknown angles ϕ y and ϕ z , respectively (see [34]); Hypothesis 2. In the flexural modes FL y , FL z the effect of micro-warping is neglected, since the associated elastic energy is small with respect to the extension of the longitudinal fibers (see [34]); Hypothesis 3. In the torsional mode TO the following energy contributions, related to the fibers belonging to the cross-section plane, are neglected: • That due to their bending in the cross-section plane; • That due to their torsion.
It is worth noticing that Hypothesis 3 has been numerically validated in what follows on case studies, as it will be shown ahead.
The energyŨ EX , associated to the EX-mode reads: in which the only contribution is due to the extension of the longitudinal fibers. The energiesŨ SH y ,Ũ SH z associated to the correspondent shear modes can be determined via Clapeyron's theorem, that is, as with the deformation work spent by the end-forces of each fiber in the correspondent (assigned) nodal displacements. They read (see [34] for details): where the two contributions on the right-end sides of Equation (13) are related to bending of the longitudinal fibers and micro-warping of the transverse fibers, respectively. Moreover, α := 12EI/ GA * h 2 is a nondimensional stiffness ratio. By requiring ∂Ũ SHy ∂ϕ y = 0, ∂Ũ SHz ∂ϕ z = 0 the unknowns ϕ y and ϕ z are eliminated, and Equation (13) reads: where: are nondimensional quantities, here referred to as shear factors, which, as shown in [34], account for the micro-warping of the fibers (i.e., a warping which leaves the joints aligned) of the fine model in the same way they account for warping of the solid cross-section in the Timoshenko beam model. The energiesŨ FL y ,Ũ FL z associated to the correspondent flexural modes are: where the two contributions on the right-end sides of Equation (16) accounts for bending and extension of the longitudinal fibers. Here, according to that discussed in [34], micro-warping of yand z-fibers is neglected. The energyŨ TO , associated with the torsional mode TO, is the sum of the torsional energies: (i) of its longitudinal fibers, which, according to the Hypothesis 3, experience the same twist angle; and (ii) of the two orthogonal orders of planar frames parallel to y and z. Concerning this latter contribution, it is found by evaluating the torsional energy of a single planar frame, undergoing a shear mode, namely γ yj = −κ x z j (y -parallel) or γ zk = κ x y k (z-parallel), analogously to what was done to determine Equation (14). Accordingly,Ũ TO reads: where the terms in the square brackets are the torsional stiffnesses of the longitudinal fibers, and of the frames parallel to the (x, z)and (y, z)-planes, respectively. Moreover, sinceŨ =Ũ EX +Ũ SH y +Ũ SH z +Ũ FL y +Ũ FL z +Ũ TO , Equation (10) suppliesφ; then, the equality φ =φ ∀ ε, γ y , γ z , κ x , κ y , κ z is enforced and the Green law supplies the following elastic constants, Equation (4), namely: c 11 = n y n z EA, , c 55 = n y EI n z + EA As a final remark, it is important to observe that, since bracing elements are assumed to be hinged, they do not energetically contribute toŨ. Accordingly, their stiffness does not enter in the definition of the elastic constants (18).

Identification of the Inertial Properties
The inertial properties of the grid beam are identified in this section for dynamical analysis. A cell is taken, made of: (i) transverse fibers, belonging to a generic cross-section D i , and (ii) half of the x-fibers below it and half above it. The mass of the fibers and that of the bracings is assumed to be lumped at joints. The kinetic energy of the cell,T , is then determined and equated to the kinetic energy T of a segment of beam of same length h, subjected to a rigid field of velocity: where M jk are the lumped masses at the (j, k)-joint, of coordinates y jk , z jk , andu jk ,v jk ,ẇ jk the correspondent translational velocities. By using the displacement map (9), it follows that:

Numerical Results
Numerical results are referred to as two grid beams, taken as case studies, made of a thermoplastic polymer, having elastic modulus E = 2180 N/mm 2 , Poisson coefficient ν = 0 and mass density ρ = 1040 kg/m 3 . The layout of the longitudinal fibers for the two case studies is shown in Figure 3, namely: in case study I (Figure 3a) a squared (or "compact") layout is considered, while in case study II (Figure 3b) a rectangular (or "thin") layout is taken. All the fibers and bracing elements have equal square cross-sections, of dimensions 1.5 mm × 1.5 mm. The elasto-geometric characteristics of the cross-section are reported in Table 1 (A is the area, A * the shear area, I the inertia moment, J the torsional inertia).  The considered boundary conditions are: (i) clamped at the ground, case studies denoted with "-a", or (ii) simply supported and torsionally restrained at the ends, case studies denoted with "-b". The number of cells is taken as n = 40 or n = 100 and, in all the cases, h = 10 mm. The static loads consist in horizontal point forces, namely, F = 0.5 N (placed as in Figure 3), applied at each cross-section, except for the free end cross-section of the clamped case, where F = 0.25 N.
Both static response, under lateral loads, and free dynamics characteristics, accounting for self-mass, are of interest. The analytical (exact) solutions relevant to the equivalent beam model are compared with the numerical results obtained by FE models, in which the grid beam is modeled as a 3D frame. Moreover: (i) a Timoshenko beam (exact) finite-element has been adopted to model the fibers of the grid beam; (ii) each fiber, comprised between two nodes of the frame, has not been meshed in static analyses, while in dynamics, it has been divided into three finite-elements.
The elasto-geometric characteristics of the equivalent beam model are reported in Table 2 for each case study. Moreover, the distributed load and couples assigned to the coarse model are p y = F/h = 0.05 N/mm and c x = (F/h) 3h = 1.5 N. In the following Figures, the green dots represent the FE solution, and the continuous blue curves are relevant to the solution of the Timoshenko beam model, respectively, while the displacement unit is mm.

Statics
The static analysis here developed leads to a comparison of the fine and coarse models in terms of: (i) lateral displacement v (x) and torsional angle θ x (x), and (ii) vertical displacements of the grid's cross-sections at selected abscissa x =x and z = 0, namely,ũ (x, y, 0). The aim of this analysis is to check the effectiveness of the Timoshenko beam model in reproducing the static behavior of the grid beam.
A very good accordance is found in cases I-a, irrespectively from the number of cells, that is, n = 40 ( Figure 4) and n = 100 ( Figure 5). In the first case, the deformed shape of the equivalent model resembles that of a shear-deformable beam, while in the second case, it resembles that of a slender (Euler) beam. In the former case, the very good agreement between the coarse and fine models is also shown in Figure 6 by comparing their deformed shapes (displacements are amplified by a factor equal to 10 in the Figure). Here, only the joints of the fine model (green dots) are represented, while the blue lines are related to the current configuration of the equivalent model, whose displacements are evaluated by using Equation (9). In Figure 6a, the lateral view at z = 0 is displayed and some of the cross-sections of the Timoshenko beam, at sampled abscissas x i , are shown to compare the two models. Figure 6b,c show two top views of the same deformed shape, at x = /2 and x = , respectively, and only the edges of the cross-section of the coarse model are represented for comparison. The order of magnitude of the percentage error, evaluated at x = , is about: 4% for v (x) and 5% for θ x (x), when n = 40; 0.5% for v (x) and 7% for θ x (x), when n = 100. In all the cases, the vertical displacements, given by the FE analysis, at different levels, are very well-described by the Timoshenko beam, by remembering the fact that the coarse model is a rigid cross-section model and, therefore, it furnishes the mean-value of the local rotation of the deformable sections of the grid beam (see Figures 4b  and 5b), which, indeed, also undergo macro-warping (i.e., misalignment of the joints).  The same very good accordance is also found in case I-b, that is, for different boundary conditions, as apparent in Figure 7. Here, the order of magnitude of the percentage error at x = /2 is about 2% for both v (x) and θ x (x). In contrast, some important quantitative differences occur for the torsional displacement in case II-a, when n = 40 (see Figure 8): indeed, here the error at x = is about 0.5% for v (x) and 32% for θ x (x). This occurrence is attributed to the fact that, in case study II, the thinner cross-sections of the grid beam are much more prone to warping, with respect to the (ticker) ones of case I, where the fibers lying in a cross-section plane constitute a more stiff lattice. Therefore, significant macro-warping is experienced from the grid beam in case study II, which cannot be well-captured by a rigid cross-section equivalent model; on the contrary, to account for this effect, more refined beam theories are needed to be able to correctly describe the deformation modes of the cross-sections, as with, for example, the Generalized Beam Theory (see, e.g., [46][47][48]). As proof of what conjectured, a numerical evaluation of the cell's energy, sided on a Finite-Element analysis, has been carried out in cases study I and II, by extending to the 3D case what was done in [35] for a planar cell, namely: (i) the nodal displacements at the boundary joints, compatible with the deformation modes, have been assigned; (ii) the reactive forces at the same joints have been evaluated; and (iii) the elastic energy has been computed via the Clapeyron Theorem. In particular, the cell's energy has been here evaluated by assigning (i) the TO mode, and (ii) an axially-unrestrained version of the TO mode (au-TO), where the kinematic constraint enforcing planarity of the cross-sections D i , is removed and warping is allowed. It was found that: • case I: c 44 = 1.90 × 10 6 N × mm 2 (TO mode) and c 44 = 1.76 × 10 6 N × mm 2 (au-TO mode); • case II: c 44 = 4.27 × 10 5 N × mm 2 (TO mode) and c 44 = 2.36 × 10 5 N × mm 2 (au-TO mode).
It can be seen that, in case I, the two torsional stiffnesses are almost the same, while in case II, the c 44 given by the axially unrestrained FEM model is about half of the restrained one, thus confirming the importance of the macro-warping effect. Finally, the numerical analysis of the cell, also corroborates the assumptions here made on the analytical identification of the torsional stiffness: indeed, the percentage error made in the analytical evaluation of the (constrained model) coefficients (remember Table 2), with respect to the numerical ones, is about 2% in both the case studies.

Dynamics
The free dynamics of the grid beam is now analyzed. By letting p = 0 and u =û (s) e λt in Equation (5), (λ,û (s)) being the eigenpairs of the system, a (spatial) boundary eigenvalue problem is obtained. Since the problem is conservative, the eigenvalues are purely imaginary, namely, λ j = ±iω j , where ω j is the angular frequency of the j-th mode. As already discussed, the boundary eigenvalue problem admits an exact solution, which, in what follows is compared with FE results, to check the effectiveness of the homogenized beam model in capturing frequencies and modal shapes of the real grid beam.
Comparison is relevant to the first six modes of the fine and coarse models and it is made for all the case studies, already discussed in statics. The results of modal analysis are presented in Table 3, where the frequencies, ω EQ and ω FE , relevant to the equivalent beam model and to FE solutions, respectively, are reported; moreover, % = 100 ω FE − ω EQ /ω FE is the percentage error. It is seen that the equivalent beam model furnishes a very accurate approximation of the first six frequencies of the fine model, particularly in the case studies I-a (n = 40, 100) and I-b (n = 40), where the maximum error is about 5%. Otherwise, in case study II-a, the percentage error increases, where just for the torsional modes (up to about 20%), this behavior is in agreement with what was already discussed in the static analysis.
Finally, a very good accordance is confirmed also by analyzing the the modal shapes of the first six modes in case studies I-a (n = 40, 100), and I-b (n = 40), reported in Figures 9-11, respectively.

Conclusions and Perspectives
Static and dynamic analyses of grid beams have been addressed in the present paper. To this end, a Timoshenko beam model was formulated in the framework of direct 1D modeling. Its constitutive law was determined via a homogenization procedure, based on the energy equivalence between the cell of the periodic grid beam and a segment of the equivalent beam, under the same point-wise displacement fields. The presented model also incorporates corrective factors, analytically determined, which accounts, on average, for the micro-warping effect. Moreover, inertial properties of the coarse model have been determined, under the assumption of lumped masses at the joints.
The effectiveness and limits of applicability of the Timoshenko beam model have been discussed by numerical analyses, carried out on sample grid beams, taken as case studies. Numerical results, concerning the exact solution of the coarse model in terms of static response and modal properties, have been compared with Finite-Element analyses.
The following conclusions have been drawn. 1.
An excellent accuracy in terms of translational displacements is given by the Timoshenko model, both in statics and in free dynamics (modal shapes), independently from the number of cells.

2.
An excellent agreement between coarse and fine models is also detected, both in statics and in free dynamics (modal shapes), for torsional displacements, where when the aspect ratio of the grid beam cross-section is close to 1, this circumstance reduces the effect of macro-warping. In contrast, when the cross-section is thin, significant errors due to macro-warping appear.

3.
Very good agreement is also detected in terms of natural frequencies. In this case, the macro-warping effect also entails worsening of the accuracy in the torsional modes.
The main perspective of this work relies, therefore, on the development of direct 1D models of beams, and is able to take into account the effects of macro-warping and to investigate the apparent similarity existing between solid and grid beams, when the aspect ratio is far from 1.