Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
Techniques for Identification of Bending and Extensional Elastic Stiffness Matrices on Thin Composite Material Plates Based on Virtual Field Method (VFM): Theoretical and Numerical Aspects
Written By
Fabiano Bianchini Batista and Éder Lima de Albuquerque
Submitted: 27 October 2010Published: 23 August 2011
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
Federal University of São João Del ReiUniversity of Brasília, Brazil
Éder Lima de Albuquerque
Federal University of São João Del ReiUniversity of Brasília, Brazil
*Address all correspondence to:
1. Introduction
Nowadays, the demand and the necessity for the use of materials with specificcharacteristics have increased in many engineering fields. Due to this necessity of making new materials, composite materials have been an alternative, or maybe the unique option, to attempt a large number of design requirements suchas high strength-to-weight ratio,high resistance to mechanical shocks, chemical attacks, corrosion, and fatigue, that cannot be obtained only from the commonly used structural materials (metals, ceramics, polymers and wood). Because of this, their applications are present in the main industries such as aerospace, automotive, marine, and sportive
Thanks to their flexibility characteristicsthere are many combinations and arrangements and, consequently, constitutive properties, that are possible to be achieved. This particularity represents one of the main advantages of these materials. Nevertheless, some factors related to arrangements as the number of layers and the orientation of the fibers can introduce abehavior called anisotropy that, in the most of the cases, is not required. The anisotropy makes the structural analysis more complex due to increasing the number of independent variables, as for example, the number of bending and extensional elastic stiffness constants.
Recently, wide part of works presented by scientific literature whose goal is to identify constitutive parameters of materials (these being composites or not) is based in the called “inverse problems”. Experimental data such as geometry, resultant forces and strain (or stress) fields are used as input data, and the unknown variables are the required constitutive parameters. In general, the solution is basically associated to two methods: iterative (or also called indirect methods) and non-iterative (or also called of direct methods). The first one is related to optimization problems where the design variables are constitutive parameters and the objective function represents, in general, a residue (or error) between experimental and numerical (generally obtained by finite elements) data. For the numerical simulations, it is considered structures that have the same geometrical characteristics and boundary conditions of the real ones. For each step, required parameters are checked out, and, the optimum represents the iteration whose residue (objective function) has its lowest value. Unlike indirect methods, the direct methods are ones where the required parameters are computed from the solution of constitutive equation(s) that are functions of these parameters.
In general, according to the type of experimental test, it is possible to separate the methods of elastic property identification in two categories: static (destructives and non-destructives) and dynamic methods (non-destructive), as shown in Fig.1. A large number of identification techniques that use data from these categories of tests have been proposed, especially ones dedicated to composite materials. It is possible to say that these techniques identify effective properties of the entire material. The way as each formulation is built, and, the adopted procedures and devices are the main differences among the many proposed methodologies.
Static tests with monotonic load are experimental tests that were more commonly used in the last years, and maybe the simplestones, for this material property identification. Despite the simplicity of these tests, some aspects render them less attractive than dynamic tests, such as the fact of requiring a number of samples with fiber orientations according to standard norms, for example, American Society for Testing and Materials (ASTM), which, in the most of cases, are not in accordance to the real characteristics of the required material. Furthermore, some variables difficult of controlling during the tests can contribute to worsen the experimental results, e.g., the presence of non-uniform stress fields near the ends of the sample from the clamped boundary conditions. For these reasons, dynamic tests have been considered an interesting alternative. In general, they are tests that combine experimental data with numerical methods, and allow the identification of elastic constants from only one unique sample or even from composite material part. Sample is usually thin plate (that reflect Kirchhoff’s hypotheses), cylindrical shell, or beam. In many cases, input data of the numerical methods are natural frequencies and/or mode shapes.
Figure 1.
Methodologies more used to identify elastic properties of materials.
The majority of works identifies only the bending stiffness matrix or directly the engineering elastic constants. However, the extensional elastic stiffness matrix is also needed to model composite materials under multi-axial loads. In general, these stiffness matrices are independent. The extensional stiffness matrix relates the in-plane resultant forces to the midplane strains, and, the bending stiffness matrix relates the resultant moments to the plate curvatures. In a laminate composite, if only the stacking sequence of layers is changed, the bending matrix is changed but the extensional matrix remains the same. In other words, different laminates can have different bending stiffness matrices and the same extensional stiffness matrix. It is not possible to obtain the extensional matrix from the bending matrix. It will be possible only if the stacking sequence of layers and their thickness are known and, also, if the material is the same for all lamina.
Sometimes, it is more convenient to use effective laminate engineering constants rather than the laminate stiffness. These effective laminated engineering constants may be easily obtained from the extensional elastic constants. However, due to difficulties on experimental in-plane modal analysis, such as the necessity of using specific devices to measure in-plane displacements and to excite high frequencies, the identification of extensional elastic stiffness constants using modal testing is less attractive. The main challenge to perform in-plane vibration testing is the excitation and measurement ofonly in-plane and not out-of-plane vibration modes. Today there are some new techniques that are suitable for this kind of problems, for example, the excitation by piezoelectric (PZT) and measurements by digital image correlation.
In this chapter, a review about the VFM applied to compute bending elastic stiffness constants proposed by Grédiac& Paris, 1996 is presented. Furthermore, a formulation based on the VFM is proposed in order to identify the extensional elastic stiffness matrix ofKirchhoff’s thin plates. The linear system of equations that provides the required elastic constants is obtained from differential equations that govern the forced vibration of anisotropic, symmetric and non-damped plates under in-plane loads. The common procedures to find the weak form (or integral form) of these equations are applied here. The correct choice of weighting functions (which are the virtual fields) and mode shapes representsa key characteristic to the accuracy of the results. Numerical simulations using anisotropic, orthotropic, quasi-isotropic plates are carried out to demonstratethe accuracy of the methodology.
The VFM has been developed for extracting constitutive parameters from full-field measurements and it is associated to problems of identification of parameters from constitutive equations. Two cases are clearly distinguished: constitutive equations depending linearly on the constitutive parameters and non-linear constitutive equations. The type of constitutive equations is chosen a priori for its relevancy and objective of it is to determine the parameters which govern the constitutive equations. The main difficulty comes from the fact that the measured displacement or strain components are generally not directly related to the unknown parameters (Grédiac et al., 2006), and no closed-form solution for the displacement, strain and stress fields is available.
Mathematically, the VFM is based on the principle of virtual work and can be written as:
-∫Vσ:ε*dV+∫SfT.u*dS+∫Vf.u*dV=∫VργudVE1
where V is volume of the solid, σ is the actual stress tensor, ε* is the virtual strain tensor, T is the distribution vector of loading tractions acting on the boundary, Sf is the part of the solid boundary where the tractions are applied, u* is the virtual displacement vector, f is the distribution of volume forces acting on V, ρ is the density and γ the acceleration. Eq. (1) is verified for any kinematically admissible virtual field (u*,ε*). Kinematically admissible means that u* must be continuous across the whole volume and it must be equal to the prescribed displacement on the boundary Sf where displacements are prescribed. Let’s introduce the constitutive equations in the general case as:
σ=g(ε)E2
whereg is a function of the actual strain. Thus, when constitutive equations are introduced and volume forces are disregarded, Eq. (1) can be rewritten as:
−∫Vg(ε):ε*dV+∫SfT.u*dS=∫Vργu*dVE3
It is possible to see in Eq. (3) that each virtual field originates a new equation involving the constitutive parameters. The VFM relies on this important property. It is a method based on setting virtual fields that provide a set of equations. This set of equations is used to extract the required unknown constitutive parameters. The correct choice of the virtual fields that combine to actual fields in Eq. (3) is the key issue of the method. Their number and their type depend on the nature of g in Eq. (3).
2.2. Review of the identification method of bending elastic stiffness matrix
The method proposed by Grédiac&Paris,1996, consists of obtaining elastic constants based on the partial differential equation that governs the transversal vibration of an anisotropic thin plate (Kirchhoff’s plate). This equation is given by:
where Dij are thin plate bending stiffness constants; ρis the mass density of the material; h is the plate thickness; x and y are coordinates of the plate; t is time; and w(x,y,t) is the deflection function that represents the transversal displacement of a point of the plate at an instant t. Eq. (4) doesn’t state the global equilibrium of the plate since the excitation force and damping are not considered. However, for many composite materials, as for example, aeronautic carbon epoxy tested in this work, the damping is low enough to disregard its contribution in the formulation. Besides, if the input data refer to resonant response of the plate, the work provided by the excitation is balanced by internal dissipation of the plate. A detailed discussion about when excitation and damping should be considered in Eq. (4) can be found in Giraudeau&Pierron, 2006.
After some mathematical manipulations in Eq. (4), Grédiac& Paris,1996 obtained a linear system in which the unknown variables are the elastic constants. Briefly, the sequence of operations is as follows: (a) multiply both sides of Eq. (4) by an arbitrary weighting function; (b) integrate twice by parts along the plate domain; (c) eliminate the boundary integrals by applying the free-edge boundary conditions; (d) decompose the displacement function w(x, y, t) as a product of the deflection amplitude Φ and sin(ѡt), where w is the out-of-plane natural frequency of a particular mode shape of the plate; and (e) choose appropriate weighting functions and mode shapes to build the matrix of the linear system. At this point, as Grédiac& Paris, 1996 explain, the choice of mode shapes associated with the weighting function is extremely important for the accuracy of this method. Three particular modes are strongly dependent on the required coefficients Dij: a twisting mode that strongly depends on terms D66,D16, and D26; a bending mode along direction 1 that strongly depends on terms D11, D12, and D16, and a bending mode along direction 2 that strongly depends on D22, D12, and D26. These modes present smooth curvatures and are generally among the first modes, with lower frequencies. If these modes are not found, it is recommended to use modes that have similar shapes to them. Furthermore, they are modes that can be approximated by quadratic functions with constant curvatures: x2, y2, and xy. For this reason, these quadratic functions were the weighting functions chosen by Grédiac& Paris, 1996. Thus, using these previous quadratic-weighting functions, the following simplified system of equations can be obtained:
where, considering L as the number of modes used in the linear system of equations, K is a 3L x 6 matrix, D is a 6 x 1 matrix, and C is a 3L x 1 matrix. As can be seen, Eq. (7) is an overdetermined system of equations. The solution can be found by least squares:
D=(KTK)−1(KTC)E8
2.3. Identification method of the extensional elastic stiffness matrix
In the general case of composite laminates, each lamina is assumed to have orthotropic material properties. After the assembly, the behavior can be anisotropic due to the interaction of different laminas. Considering a plate under plane state of stress and using Hooke’s generalized law, stresses can be integrated over its thickness yielding the following force-deformation equations:
{N⋯M}=[A⋮B⋯⋮⋯B⋮D]{ε⋯κ}E9
where N and M are vectors that contain normal forces and resultant moments, respectively, A is the extensional elastic stiffness matrix, B is the coupling elastic stiffness matrix (B is a null matrix in the case of a symmetric laminate), D is the bending elastic stiffness matrix, ε and κ are vectors that contain middle plane linear strains and rotations, respectively. Considering a symmetrical (B = [0]) and fully anisotropic laminate under free-edge in-plane vibration (the plate is not under bending) and using the equilibrium relations, the following equations can be written:
E10
E11
where Aij are the elements of matrix A (i,j = 1,2,6); ρ is the mass density; h is the plate thickness; x and y are the coordinates in the plate plane; t is the time, and u(x,y,t) and v(x,y,t) are functions that represent the displacements along x and y direction, respectively, of a point with coordinates (x,y) of the plate at an instant t. Multiplying Eq. (10) by a weighting function W(x, y) and integrating along the domain of the plate, we can obtain:
E12
whereΩ is the plate domain. Using chain rule, i.e.,
ddx(f.g)=fdgdx+gdfdxE13
wheref and g are any two continuous functions, and Green’s theorem, i.e.,
∫Ω∂f∂xdΩ=∫ΓfnxdxE14
,
where Γ is the boundary domain, and nx is the component of the normal unity vector in directions x, the left hand side terms of Eq. (12) can be written as:
E15
whereny is the component of the normal unity vector in direction y. From the constitutive equation, one has:
whereNx and Ny are axial forces per unit length along directions x and y, respectively, Nxy is the shear force per unit length along plane xy, εx and εy are the middle surface axial strain along directions x and y, respectively, and γxy is the shear angular strain along plane xy. Rewriting Eq. (16), one obtains:
Multiplying Eq. (17) by a weighting function W and by nx, which is the x component of unit normal vector n, and integrating on the boundary Γ, one obtains:
Multiplying Eq. (19) by a weighting function W and by ny, which is the y component of the unit normal vector n, and integrating on the boundary Γ, yields:
Now, if free-edge boundary conditions are considered, boundary integrals of Eq. (22) vanish. Considering that the plate is vibrating, functions u and v can be written as:
u(x,y,t)=U(x,y)sin(ϖt)E23
,
v(x,y,t)=V(x,y)sin(ϖt)E24
,
whereϖ is the in-plane natural frequency associated to any in-plane mode, and U(x,y) and V(x,y) are the displacement amplitudes along directions x and y, respectively, of a point with coordinate (x,y). In this sense, the amplitude is only a function of coordinates x and y. Eliminating the boundary integrals of Eq. (22), substituting Eqs. (23) and (24) into Eq. (22), and considering any mode j, one obtains:
Eqs. (25) and (26) are theoretically valid for isotropic, orthotropic, or anisotropic plates, provided that the laminate is symmetrical. As it can be seen, the function W is arbitrary since itself and its first order derivative is continuous in the domain Ω. Amplitudes U(x,y) and V(x,y), and frequencies ϖj are obtained from dynamic tests. Dimensions of the plate and parameters ρ and h can also be easily measured on the sample plate. Thus, the next steps are the choice of a suitable numerical method to compute the derivatives and integrals in Eqs. (25) and (26). Furthermore, suitable mode shapes and weighting functions should be chosen. In this work, finite differences and Gauss-Legendre numerical integration scheme are used to compute these derivatives and integrals, respectively. For numerical reasons, modes with several sign changes in the mode shape are avoided because their numerical derivatives and integrals are more sensitive to errors (Grédiac& Paris, 1996). Generally, first modes present more smooth curvatures and are, at the same time, easier to be obtained experimentally. It is worth noting that in-plane modal analysis presents much higher frequencies than transverse modal analysis (bending modes). This is because stiffness along direction x and y is much higher than stiffness along the transversal direction of the plate. For numerical reasons, smooth mode shapes associated with smooth weighting functions are preferred. For all these reasons, and in order to simplify Eqs. (25) and (26), the following group of weighting functions are proposed:
W(x, y) = x2, which applied to Eqs. (25) and (26) provides the following integral equations:
Eq. (39) can also be rewritten in a compact form as:
KA=CE40
,
where, considering L modes, K is a 6L x 6 matrix in Eq. (40), A is a 6 x 1 matrix, and C is a 6L x 1 matrix in Eq. (40). Eq. (40) is an over determined system. Thus, the solution can be found by least squares:
A=(KTK)−1(KTC)E41
,
from where the extensional elastic constants Aij are computed.
A commercial finite element code (ANSYS 11.0) was used to give particular mode shapes and their corresponding natural frequencies from both in-plane and out-of-plane numerical modal analysis. Element SHELL99 was used and plates under free-edge boundary conditions were considered.
To exemplify the method proposed by Grédiac and Paris (1996), it was used an anisotropic plate with dimensions 0.450 x 0.350 x 0.0021 m and density 1500 kg/m3. It was used a laminate with 8 plies, [0 45 90 135]S, and the following engineering elastic constants by ply: E1=120 GPa (Young’s module along the principal direction 1), E2=10 GPa (Young’s module along the principal direction 2), G12=4.9 GPa (shear module along the plane 1-2), and ν12=0.3 (Poisson’s ratio along the plane 1-2). A mesh of 651 nodes was used. Fig. 2 shows the three modes used to identify the required properties. As can be seen, depending on type of material anisotropy, it is not possible to find all three modes necessary to apply the method. In this case, it is necessary to find the more approximated ones.
Table 1 shows the bending elastic stiffness constants computed using the engineering constants and the classical theory of laminates, and it also shows the errors computed after applying the identification method. As can be observed, the technique is able to find very satisfactory results when it is used the correct modes. The problem of this technique is the high sensitivity to noise presence because of second-order derivatives. More results and comments about this method can be found in Grédiac& Paris, 1996.
In order to verify the accuracy of the extensional elastic stiffness identification method, it was used six graphite/polymer symmetric laminated plates (Table 2): a fully anisotropic with all Aij different from zero (i,j = 1, 2, and 6); a cross-ply orthotropic with A11 = A22, A16 = A26 = 0, and A11 – A12 ≠ 2A66; a 0° unidirectional orthotropic; a 30° unidirectional orthotropic (generallyorthotropic); a +30°/-30 angle-ply orthotropic; and a quasi-isotropic with A11 = A22, A16 = A26 = 0, and A11 – A12 = 2A66. These laminates have 8 plies with the following engineering elastic constants by ply: E1=155 GPa, E2=12.10 GPa, G12=4.4 GPa, and ν12=0.248. Dimensions considered were 0.450 x 0.350 x 0.003 m (x, y, and z plate coordinate axis, respectively), for the rectangular plate, 0.350 x 0.350 x 0.003 m, for the square plate, and, the density material was 1500 kg/m3. The plates were modeled using a mesh with 651 nodes, for the rectangular plates, and with 441 nodes, for the square plate. The extensional elastic stiffness constants are shown in Table 2. The terms “Aniso”, “Ortho” and “Quasi-iso” are simplifications of “Anisotropic”, “Orthotropic” and “Quasi-isotropic”, respectively. Only the first fifteen mode shapes were analyzed.
As can be seen in Table 2, the constants A12 and A66 for the orthotropic laminates I and IIare equals. They also are equals to laminates with the same characteristics but with 90° unidirectional fibers. 0° and 90° unidirectional laminated plates are only different in relation to A11 and A22 terms. These are inverted: A11 term of the 0° laminate is equal to A22 term of the 90° laminate, and vice-verse. For the generally orthotropic laminate III and orthotropic laminate IV the difference are only the A16 and A26 constants: they are nulls for the laminate IV and non-nulls for the laminate III. For laminate III all extensional elastic constants are non-nulls, similar to fully anisotropic laminates, what, consequently, originates to full extensional elastic stiffness matrix.
Bending elastic constants
N x mm
Errors (%)
D11
64363.9
0.02
D22
24155.8
0.04
D12
8875.1
0.02
D66
10032.7
1.22
D16
6019.6
0.63
D26
6019.6
0.64
Table 1.
Bending elastic stiffness constants of the tested anisotropic plate and the computed errors.
The key point of this technique of identification is related with the correct choice of mode shapes together to the weighting functions (virtual fields). The correct mode shapes are called here by “suitable modes” and the correct combination between these modes and the weighting functions are called by “suitable combinations”. The identification of the suitable modes is not difficult, as it will be shown in the next topics. But, the suitable combinations are more difficult because they depend on the type and the geometry of the material. Fortunately, there are some aspects that help finding the best choice. Unlike the bending stiffness identification method originally proposed, for this method there are a lot of modes and suitable combinations that give satisfactory results.
Aij 108 [N/m]
Aniso [90 0 0 45]S
Ortho I [90 0 90 0]S
Ortho II [0]4S
Ortho III [30]4S
Ortho IV [30 -30]2S
Quasi-iso [90 45 0 -45]S
A11
2.7865
2.5186
4.6724
2.7840
2.7840
1.9776
A12
0.3610
0.0905
0.0905
0.9020
0.9020
0.6315
A16
0.2692
0
0
1.4012
0
0
A22
1.7096
2.5186
0.3648
0.6301
0.6301
1.9776
A26
0.2692
0
0
0.4641
0
0
A66
0.4025
0.1320
0.1320
0.9436
0.9436
0.6730
Table 2.
Laminated plates used for the verification of the method
[90 0 0 45]S anisotropic plate
Table 3 shows some errors computed for the anisotropic plate, rectangular and square. It was considered only the first fifteen in-plane modes shapes. Anisotropic plates, in general, give very satisfactory results using the combinations among suitable modes. This factor can be justified by the fact of these combinations be hardly involved with all required extensional elastic constants Aij’s.
The numerical contribution of each mode to the computation of a specific constant cannot be jeopardized by numerical contribution of another mode during the solution of the system given by Eq. (40).The suitable modes are those that when associated with weighting functions do not null or give very low values for integrals of the right (K matrix) and/or left (C matrix) hand sides of Eq. (40). The suitable combinations are one composed by suitable modes and that give more accurate results. In the majority of the cases, combinations using a higher number of suitable modes can be suitable combinations. According to Table 3 is possible to see that using combinations with only two suitable modes very satisfactory results can be obtained. Satisfactory results would also be obtained using combinations with any modes since the number of suitable modes among all used modes is higher than non-suitable modes. But the accuracy of these results cannot be guaranteed for all combinations.
[0 90 0 90]Sorthotropic plate (ortho I)
In general, for the orthotropic and isotropic materials is more difficult to find the suitable combinations when it is compared to fully anisotropic materials. It is necessary to take care to correctly identifying the combinations that give the best results. In these types of materials not all combinations are among suitable modes that can be considered as being suitable combinations. According to values found to terms of the K and C matrices, Eq. (40), and using combinations among suitable modes, it is possible to see the following types of systems:
Errors computed for the anisotropic plate using some suitable combinations.
whereI1, I2, I3, I4, I5, I6, and I7 are the integral values of Eq. (40). These integrals can be negatives or positives depending on strain direction and reference coordinate axis. As can be observed, for these two types of systems, Eq. (42) and Eq. (43), all elastic constants are involved, and, however, combinations associated to only one type of system can be suitable combinations and sufficient to give correct results. For this plate were found the following suitable modes: 2, 3, 6, 8, 11, and 14.
For the square plate: for this plate, the suitable modes are: 2, 3, 6, 7, 11, and 12. The systems of equations are full, even though of the additional terms, that are null in rectangular plate, to be low in this square plate. It is observed in this plate that for each suitable mode there is another identical but out-of-phase at 90°: modes 2 and 3, 6 and 7, and, 11 and 12.
Table 4 shows errors computed to some suitable combinations for these orthotropic plates. As can be seen, very satisfactory results can be obtained using correct combinations of modes. For this type of orthotropy, it can be more difficult to compute an accurate value for constant A12. The identification of the suitable combinations is not so clear. As this technique of identification is associated to the solution of an equation system having different modal contributions, it is difficult to identify which modes compose a correct combination.
Ortho I rectangular plate – suitable modes: 2, 3, 6, 8, 11, and 14
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-3-6-8-11-14
0.14
3.80
0.12
1.22
-2.7 x 104
0.8 x 104
2-3-6-8
0.11
2.90
0.24
0.74
-2.5 x 104
0.6 x 104
3-6-8-11
0.17
1.58
0.31
1.67
-1.4 x 104
0.5 x 104
2-8-11
0.17
0.29
0.23
1.46
-1.5 x 104
0.6 x 104
6-8-11
0.18
0.89
0.34
1.91
-4.8 x 104
1.0 x 104
3-8
0.02
5.75
0.91
0.13
-1.4 x 104
8.4 x 103
8-11
0.25
0.42
0.40
2.17
-2.6 x 104
1.1 x 104
Ortho I square plate – suitable modes: 2, 3, 6, 7, 11, and 12
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-3-6-7-11-12
0.09
5.81
0.05
1.89
2.7 x 104
-1.6 x 104
2-3-6-11
0.10
5.09
0.06
1.92
-4.4 x 104
-8.7 x 104
3-6-7-12
0.38
2.33
0.08
2.92
6.4 x 104
-1.7 x 104
2-6-11
0.13
5.47
0.11
2.51
-1.0 x 105
-1.8 x 104
2-7-12
0.22
5.37
0.04
1.95
1.1 x 105
0.5 x 105
2-11
0.02
9.04
0.18
0.56
-1.3 x 104
-1.9 x 104
6-11
0.18
5.63
0.30
3.46
-1.1 x 105
-2.7 x 104
Table 4.
Errors computed for the ortho I plate using some suitable combinations.
[0]4Sorthotropic plate (ortho II)
For the rectangular plate: for this plate, modes 3, 10, 11, and 14 originate systems of type 1, Eq. (42), and modes 2, 5, 9, and 15 originate systems of type 2, Eq. (43).
For the square plate: for this plate, modes 4, 10, 12, and 14 originate systems of type 1, Eq. (51), and modes 2, 5, 9, and 13 originate systems of the type 2, Eq. (43).
Table 5 shows the errors computed to some suitable combinations for these orthotropic plates. Using correct combinations very satisfactory results can be obtained. Similar to ortho I plate, for this type of orthotropy, constant A12is more difficult of being accurately computed.
Ortho II rectangular plate – suitable modes: 2, 3, 5, 9, 10, 11, 14, and 15
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-3-5
2.77
6.06
0.08
0.38
2.4 x 104
0.5 x 104
11-14-15
1.14
2.15
0.92
1.15
2.3 x 105
0.1 x 105
2-3
2.36
0.57
0.19
0.72
1.7 x 104
0.1 x 104
2-10
2.10
0.82
0.60
2.42
-2.0 x 105
-0.4 x 105
2-14
2.04
5.22
0.37
1.12
-3.9 x 104
-1.2 x 104
11-14
1.64
2.14
0.87
0.98
-1.8 x 105
-0.1 x 104
Ortho II square plate – suitable modes: 2, 4, 5, 9, 10, 12, 13, and 14
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-4-5
4.13
8.63
0.39
0.45
-3.5 x 104
-0.1 x 104
2-4-14
2.43
0.78
0.72
2.89
4.2 x 104
1.1 x 104
2-4
3.35
0.10
0.21
1.25
-3.2 x 104
-0.1 x 104
2-14
1.63
1.59
1.30
4.72
1.1 x 105
0.3 x 105
5-14
2.67
13.54
2.21
2.82
5.6 x 104
-1.0 x 104
9-14
2.12
62.05
0.85
3.89
1.1 x 105
0.2 x 105
Table 5.
Errors computed for the ortho II plate using some suitable combinations.
[30]4Sorthotropic plate (ortho III)
For these plates, rectangular and square the computed matrices K and C, Eq.(40), are full matrices, similar ones of the anisotropic plates. Thus, the majority of combinations among suitable modes are suitable combinations. Combinations that are not suitable present high errors for all constants, what, consequently, make them easy to be identified. Table 6 shows the errors computed to some suitable combinations. Using correct combinations very satisfactory results can be obtained.
Ortho III rectangular plate – suitable modes: 2, 3, 6, 8, 10, 12, and 13
Suitable combinations
Errors (%)
A11
A12
A22
A66
A16
A26
3-6-12
0.25
0.51
0.10
1.05
0.69
0.99
8-10-12
4.26
5.42
3.67
3.43
4.60
5.62
2-3
0.85
1.72
1.22
1.39
0.10
0.96
3-6
0.63
0.32
0.04
0.22
0.53
0.14
6-12
0.71
1.06
0.05
1.39
1.19
1.30
10-13
4.42
2.05
0.71
4.07
5.21
0.63
Ortho III square plate – suitable modes: 2, 4, 6, 8, 9, 12, and 13
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-4-6
2.62
3.94
3.18
0.59
2.02
3.23
2-4-8
1.20
2.40
1.64
1.27
0.11
1.19
2-4
2.35
1.11
0.02
5.46
4.26
2.90
2-6
0.76
2.76
3.52
0.01
0.53
2.90
2-13
1.75
2.75
1.51
1.92
0.33
1.25
12-13
3.08
2.41
2.63
0.53
1.85
2.54
Table 6.
Errors computed for the ortho III plate using some suitable combinations.
[30 -30 30 -30]Sorthotropic plate (ortho IV)
For the rectangular plate: it was found the following suitable modes: 2, 4, 5, 8, 11, 12 and 14. For this plate, the majority of combinations among the suitable modes are suitable combinations.
For the square plate: it was found the following suitable modes:1, 4, 6, 9, 11, 12, and 14. Similar to rectangular plate, here the most of combinations among the suitable modes are suitable combinations.
Table 7 shows errors computed to some suitable combinations. Using correct combinations, very satisfactory results can be obtained.
Ortho IV rectangular plate – suitable modes: 2, 4, 5, 8, 11, 12, and 14
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
2-4-5
1.52
2.15
1.28
1.54
-1.7 x 104
-0.6 x 104
2-4-8
1.50
2.09
0.83
1.14
1.4 x 104
0.9 x 104
5-11-12
0.27
0.09
0.04
1.11
-1.2 x 104
1.9 x 104
2-4
1.18
1.83
1.02
2.87
-1.2 x 104
-0.4 x 104
2-5
0.49
1.01
1.01
1.41
-1.8 x 104
-0.7 x 104
5-11
0.90
0.37
0.04
1.05
3.1 x 104
0.2 x 104
Ortho IV square plate – suitable modes: 1, 4, 6, 9, 11, 12, and 14
Suitable combinations
Errors (%)
Differences (N/m)
A11
A12
A22
A66
A16
A26
1-6-12
0.16
0.30
0.13
3.07
1.7 x 104
2.5 x 104
6-9-12
0.18
0.35
0.30
2,68
1.3 x 104
4.8 x 104
9-12-14
0.18
1.87
0.40
0.68
-6.2 x 104
8.1 x 104
1-12
0.59
0.26
0.74
2.79
2.0 x 104
2.9 x 104
4-11
0.94
1.55
2.77
0.23
2.2 x 104
0.8 x 104
12-14
0.26
2.17
0.19
0.90
-9.7 x 104
9.4 x 104
Table 7.
Errors computed for the ortho IV plate using some suitable combinations.
[90 45 0 -45]Squasi-isotropic plate
For the rectangular plate: it was found the following suitable modes: 1, 4, 6, 8, 11, 12, and 14. For this plate the most of combinations among the suitable modes are suitable combinations.
For the square plate: it was found the following suitable modes: 2, 3, 7, 8, 10, and 11. Similar to rectangular plate, here the majority of combinations among the suitable modes are suitable combinations.
Table 8 shows errors computed to some suitable combinations. Using correct combinations, very satisfactory results can be obtained.
The identification of elastic properties using VFM has shown to be a very efficient technique since the correct combinations among mode shapes and weighing functions are used. This factor is the key point to find the correct results. However, the identification of these suitable combinations is not so simple in some situations, mainly to the extensional elastic stiffness identification method. Fortunately, there are some characteristics that can help to find such combinations, as it was shown here. A great advantage of this method is related to the large number of possibilities to make combinations able to give very satisfactory results.
References
1.AlfanoM. .PagnottaL.2007A non-destructive technique for the elastic characterization of thin isotropic plates.NDT Internatinal, 40112120
2.AraújoA. L.MotaSoares. C. M.MoreiraFreitas. M. J.PedersenP.HerskovitsJ.2000Combined numerical-experimental model for the identification of mechanical properties of laminated structures.Composite Structures, 50363372
3.AvrilS.HuntleyJ. M.PierronF.SteeleD. D.2008D heterogeneous stiffness identification using MRI and the virtual fields method. IEEE Transactions on Medical Imaging, 484479494
4.AvrilS. .PierronF.2007General framework for the identification of elastic constitutive parameters from full-field measurements.International Journal of Solids and Structures, 4449785002
5.AyorindeE. O.GibsonR. F.1995Improved method for in-situ elastic constants of isotropic and orthotropic composite materials using plate modal data with trimodal and hexamodal Rayleigh formulations.Journal of Vibration and Acoustics, 117180186
6.AyorindeE. O.YuL.1999On the use of diagonal modes in the elastic identification of thin plates.Journal of Vibration and Acoustics, 1213340
7.AyorindeE. O.YuL.2005On the elastic characterization of composite plates with vibration data.Journal of Sound and Vibration, 283243262
8.Bledzki, A.K.; Kessler, A.;Rikards, R. &Chate, A. (1999). Determination of elastic constants of glass/epoxy unidirectional laminates by the vibration testing of plates. Composites Science and Technology, 59:2015-2024
9.BrunoL.FeliceG.PagnottaL.PoggialiniA. .StiglianoG.2008Elastic characterization of orthotropic plates of any shape via static testing.International Journal of Solids and Structures, 45908920
10.ChalalH.AvrilS.PierronF. .MeraghniF.2006Experimental identification of a nonlinear model for composites using the grid technique coupled to the virtual fields method. Composite Part A: Applied Science and Manufacturing, 372315325
11.ChakrabortyS. .MukhopadhyayM.2000Estimation of in-plane elastic parameters and stiffener geometry of stiffener plates.Journal of Sound and Vibration, 231199124
12.CugnoniJ.GmürT. .SchorderetA.2007Inverse method based on modal analysis for charactering the constitutive properties of thick composite plates. Computers and Structures, 8513101320
13.DeobaldL. R.GibsonR. F.1988Determination of elastic constants of orthotropic plates by modal analysistechnique.Journal of Sound and Vibration, 1242269283
14.DiveyevB. .Butiter(2008I.2008Identifying the elastic moduli of composite plates by using high-order theories. 1: theoretical approach. Mechanics of Composite Materials, 4412536
15.DiveyevB.ButiterShcherbinaI. N.2008Identifying the elastic moduli of composite plates by using high-order theories. 2: theoretical-experimental approach. Mechanics of Composite Materials, 442139144
16.GibsonR. F.2000Modal vibration response measurements for characterization of composite materials and structures.Composite Science and Technology, 6027692780
17.GiraudeauA.GuoB. .PierronF.2006Stiffness and damping identification from full-field measurements on vibration plates.Experimental Mechanics, 466777787
18.GiraudeauA. .PierronF.2003Simultaneous identification of stiffness and damping properties of isotropic materials from forced vibration plates. ComptesrendusMécanique, 3314259264
19.GiraudeauA. .PierronF.2005Identification of stiffness and damping properties of thin isotropic vibrating plates using the Virtual Fields Method.Theory and simulations.Journal of Sound and Vibration, 2843-5757781
20.GrédiacM.1996On the direct determination of invariant parameters governing the bending of anisotropic plates.International Journal of Solids and Structures, 333969033982
21.GrédiacM.2004The use of full-field measurement methods in composite material characterization: interest and limitations. Composites Part A, 35751761
22.GrédiacM.AuslenderF. .PierronF.2001Applying the virtual fields method to determine the through-thickness moduli of thick composites with a nonlinear shear response. Composites Part A: Applied Science and Manufacturing, 321217131725
23.GrédiacM.FounierN.ParisP. A. .SurrelY.1999aDirect identification of elastic constants of anisotropic plates by modal analysis: experiments and results. Journal of Sound and Vibration, 210645659
24.Grédiac, M.;Founier, N.;Surrel, Y. &Pierron, F. (1999b).Direct measurement of invariant parameters of composite plates.Journal of Composite Materials, Vol. 33, pp. 1939-1965
25.GrédiacM.ParisP. A.1996Direct identification of elastic constants of anisotropic plates by modal analysis: theoretical and numerical aspects. Journal of Sound and Vibration, 195401415
26.GrédiacM. .PierronF.1998A T-shaped specimen for the direct characterization of orthotropic materials. International Journal for Numerical Methods in Engineering, 41293309
27.GrédiacM. .PierronF.2006Applying the Virtual Fields Method to the identification of elasto-plastic constitutive parameters.International Journal of Plasticity, 22602627
28.GrédiacM.PierronF.AvrilS.ToussaintE.2006The virtual fields method for extracting constitutive parameters from full-field measurements: a review. Strain: an International Journal for Experimental Mechanics, 42233253
29.GrédiacM.ToussaintE. .PierronF.2002aSpecial virtual fields for direct determination of material parameters with the virtual fields method. 1-Principle and definition.International Journal of Solids and Structures, 391026912705
30.GrédiacM.ToussaintE. .PierronF.2002bSpecial virtual fields for direct determination of material parameters with the virtual fields method. 2-Application to in-plane properties.International Journal of Solids and Structures, 391027072730
31.Grédiac, M.; Toussaint, E. &Pierron, F. (2003). Special virtual fields for direct determination of material parameters with the virtual fields method. 2-Application to the bending rigidities of anisotropic plates.International Journal of Solids and Structures, Vol. 40, No. 10, pp. 2401-2419
32.HwangS.ChangC.2000Determination of elastic constants of materials by vibration testing.Composite Structures, 49193190
33.LaiT. C.LauT. C.1993Determination of elastic constants of a generally orthotropic plate by modal analysis.International Journal of Analytical and Experimental Modal Analysis, 81533
34.LauwagieT.SolH.HeylenW. .RoebbenG.2004Determination of in-plane elastic properties of different layers of laminated plates by means of vibration testing and model updating.Journal of Sound and Vibration, 274529546
35.LauwagieT.SolH.RoebbenG.HeylenW.ShiY.Van der BiestO.2003Mixed numerical-experimental identification of elastic properties of orthotropic metal plates. NDT e International, 36487495
36.LeeC. R. .KamT. Y.2006Identification of mechanical properties of elastically restrained laminated composite plates using vibration data. Journal of Sound and Vibration, 2959991016
37.PagnottaL. .StiglianoG.2008Elastic characterization of isotropic plates of any shape via dynamic tests: theoretical aspects and numerical simulations. Mechanics Research Communications, 35351360
38.PedersenP. .FrederiksenP. S.1992Identification of orthotropic material moduli by a combined experimental /numerical method. Measurement, 10113118
39.PierronF. .GrédiacM.2000Identification of the through-thickness moduli of thick composite from whole-field measurements using the Iosipescu fixture.Composites Part A: Applied Science and Manufacturing, 314309318
40.PierronF.ZhavoronokS. .GrédiacM.2000Identification of the through-thickness properties of thick laminated tubes using the virtual fields method. Internatinal Journal of Solids and Structures, 373244374453
41.PierronF.VertG.BurgueteR.AvrilS.RotinatR. .WisnomM.2007Identification of the orthotropic elastic stiffness of composites with the virtual fields method: sensitivity study and experimental validation. Strain: an International Journal for Experimental Mechanics, 433250259
42.ReverdyF. .AudoinB.2001Ultrasonic measurement of elastic constants of anisotropic materials with laser source and laser receiver focused on the same interface. Journal of Applied Physics, 90948294835
43.RikardsR. .ChateA.1998Identification of elastic properties of composites by method of planning of experiments.Composite Structures, 42257263
44.RikardsR.ChateA. .GailisG.2001Identification of elastic properties of laminates based on experiment design.International Journal of Solids and Structures, 3850975115
45.RikardsR.ChateA.SteinchenW.KesslerA.BledzkiA. K.1999Method for identification of elastic properties of laminates based on experiment design.Composites: Part B, 30279289
46.ToussaintE.GrédiacM. .PierronF.2006The virtual fields method with piecewise virtual fields. International Journal of Mechanical Sciences, 48256264
Written By
Fabiano Bianchini Batista and Éder Lima de Albuquerque
Submitted: 27 October 2010Published: 23 August 2011