A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites
|
|
- Aldous Heath
- 5 years ago
- Views:
Transcription
1 Materials Science and Engineering A 412 (2005) A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites Harald Berger a,, Sreedhar Kari a, Ulrich Gabbert a, Reinaldo Rodríguez-Ramos b, Julián Bravo-Castillero b,raúl Guinovart-Díaz b a Institute of Mechanics, Otto-von-Guericke-University of Magdeburg, Universitaetsplatz 2, D Magdeburg, Germany b Facultad de Matemática y Computación, Universidad de La Habana, San Lázaro y L, CP Vedado, Habana 4, Cuba Received in revised form 3 June 2005 Abstract This work deals with the modelling of periodic composites made of piezoceramic (lead zirconate-titanat) fibres embedded in a soft nonpiezoelectric matrix (polymer). The goal is to predict thective coefficients of such periodic transversely isotropic piezoelectric fibre composites by use of a representative volume element or unit cell. The solution is based on a numerical approach using finite element method (FEM). The necessary basic equations for the piezoelectric material are introduced and the special concept for definition of generalised periodic boundary conditions for the unit cell is explained. For a composite with square arrangements of cylindrical fibres the algorithm is demonstrated and the extension to other fibre arrangements is shown. For different fibre volume fractions the results are compared with analytical solutions Elsevier B.V. All rights reserved. Keywords: Homogenisation; Representative volume element; Piezoelectricity; Finite element method 1. Introduction Piezoelectric materials have the property of converting electrical energy into mechanical energy and vice versa. This reciprocity in the energy conversion makes piezoelectric ceramics such as lead zirconate-titanat (PZT) very attractive materials towards sensors and actuators applications. But bulk piezoelectric materials have several drawbacks, hence composite materials are often a better technological solution in the case of a lot of applications such as ultrasonic transducers, medical imaging, sensors, actuators, and damping. For the last 20 years, composite piezoelectric materials have been developed by combining piezoceramics with passive non-piezoelectric polymers. Superior properties have been achieved by these composites by taking advantage of most profitable properties of each constituents. Recently, due to the miniaturisation of piezoelectric composites and the use of PZT fibres instead of piezoelectric bars, homogenisation techniques are necessary to describe the behaviour of piezoelectric composites. Corresponding author. Tel.: ; fax: address: berger@mb.uni-magdeburg.de (H. Berger). Even if analytical and semi analytical models have been developed to homogenise piezoelectric composites, they are often reduced to specific cases. Numerical methods, such as the finite element method, seem to be a well-suited approach to describe the behaviour of these materials, because there are no restrictions to the geometry, the material properties, the number of phases in the piezoelectric composite, and the size. The prediction of the mechanical and electrical properties of piezoelectric fibre composites became an active research area in recent years. Micro mechanical methods provide an overall behaviour of piezoelectric fibre composites from known properties of their constituents (fibre and matrix) through an analysis of a periodic representative volume element (RVE) or a unit cell model. In the macro mechanical approach, on the other hand, the heterogeneous structure of the composite can then be replaced by a homogeneous medium with anisotropic properties. A number of methods have been developed to predict and to simulate the coupled piezoelectric and mechanical behaviour of composites. Basic analytical approaches have been reported [1,2], which are not capable of predicting the response to general loading. This restriction can be overcome by employing periodic micro field approaches (commonly referred to as unit cell models) where the fields are typically solved numerically with high resolution, e.g. by the finite element method [3]. Sometimes /$ see front matter 2005 Elsevier B.V. All rights reserved. doi: /j.msea
2 54 H. Berger et al. / Materials Science and Engineering A 412 (2005) special finite element formulations were introduced to solve the problem [4]. In such models, the representative unit cell and the boundary conditions are designed to capture a few special load cases, which are connected to specific deformation patterns [5,6]. But many of these approaches allow the prediction of only a few key material parameters and not the full set due to restrictions in load-boundary combinations for the numerical unit cell model. The aim of this paper is to predict the full set of material moduli, i.e. to determine the complete tensors associated with the overall elastic, dielectric, and piezoelectric behaviour. This means that the linear response to any mechanical and electrical load, or any combination of both, will be determined. In the present paper, the FEM based on a micro mechanical analysis is applied to uniaxial periodic piezoelectric fibre composites subjected to different loading conditions with different boundary T T 22 T T 23 T 31 = T 12 D 1 D 2 D 3 conditions to predict thective coefficients of transversely isotropic piezoelectric fibre composites (1 3 periodic). The advantage of this technique lies in using a standard FE code in connection with auxiliary routines written in FORTRAN for modelling of unit cells with the appropriate boundary and loading conditions and its possibility to extend it to complex composites with arbitrary types of inclusions and non-uniform distributions. The algorithm is demonstrated for a composite where the embedded fibres has a square arrangement. For verification of results various fibre volume fractions and different fibre arrangements are investigated and a comparison is done with values calculated by an analytical solution based on the asymptotic homogenisation method (AHM) reported in [7,8]. 2. Piezoelectricity and piezoelectric composites Coupled piezoelectric problems are those in which an electric potential gradient causes deformation (converse piezoelectric effect), while mechanical strains cause an electric potential gradient in the material (direct piezoelectric effect). The coupling between mechanical and electric fields is characterised by piezoelectric coefficients. Those materials respond linearly to changes in the electric field, the electrical displacements, or mechanical stresses and strains. These assumptions are compatible with the piezoelectric ceramics, polymers, and composites in current use [9]. Therefore, the behaviour of the piezoelectric medium is described by the following piezoelectric constitutive equations, which correlate stresses T, strains S, electric fields E, and electrical displacements D as follows: [ ] [ ][ ] T C e T S =, (1) D e ε E where C is the elasticity matrix, ε the permittivity matrix, and e is the piezoelectric strain coupling matrix. For a transversely isotropic piezoelectric solid, the stiffness matrix, the piezoelectric matrix, and the dielectric matrix simplify so that there remain independent coefficients. In the case of aligned fibres made of a transversely isotropic piezoelectric solid (PZT), embedded in an isotropic polymer matrix, the resulting composite is a transversely isotropic piezoelectric material too. Consequently, the constitutive Eq. (1) can be written as S S S S 23 S 31. (2) S 12 Ē 1 Ē C44 eff C66 eff ε eff ε eff ε eff In this matrix, the general variables of the coupled electromechanical problem were replaced by the appropriate values for the homogenised structure. So, Cij eff, eeff ij, εeff ij denote thective material coefficients and S ij, Ē i, T ij, D i denote average values. These relations represent the basis for the further considerations based on the unit cell. Ē 3 3. Representative volume element (unit cell) In general, the object under consideration is regarded as a large-scale/macroscopic structure. The common approach to model the macroscopic properties of 3D piezoelectric fibre composites is to create a representative volume element (RVE) or a unit cell that captures the major features of the underlying microstructure. The mechanical and physical properties of the constituent material is always regarded as a smallscale/microstructure. One of the most powerful tools to speed up the modelling process, both the composite discretisation and the computer simulation of composites in real conditions, is the homogenisation method. The main idea of the method is to find a globally homogeneous medium equivalent to the original composite, where the strain energy stored in both systems is approximately the same. In this paper, we limit ourselves to a quasistatic analysis of periodic (1 3) structures with perfectly bonded continuous fibres, which are aligned and poled along the x 3 axis as shown in Fig. 1. With help of symmetry, such a regular piezoelectric fibre composite may be analysed by using a representative volume element or unit cell. A unit cell is the smallest
3 H. Berger et al. / Materials Science and Engineering A 412 (2005) Fig. 1. Schematic diagram of a periodic 1 3 composite (a) and unit cell (b), picked from the original composite. part that contains sufficient information on the above-mentioned geometrical and material parameters at the microscopic level to deduce thective properties of the composite. Fig. 1 shows the unit cell that is picked from the periodic piezoelectric fibre composite. It has infinite lengths in all three directions. For demonstration of the algorithm we consider a composite with a square fibre arrangement. It is assumed that the material properties are the same in the first two directions (i.e. along x 1 and x 2 axis). All the fibres are assumed to be straight and poled in the third direction (i.e. along axis x 3 ). Fig. 1(b) shows the schematic diagram of the unit cell picked from the considered composite. 4. Numerical solution using finite element method 4.1. Periodic boundary conditions to RVE Composite materials can be represented as a periodical array of RVEs. Therefore, periodic boundary conditions must be applied to the RVE models. This implies that each RVE in the composite has the same deformation mode and there is no separation or overlap between the neighbouring RVEs. These periodic boundary conditions on the boundary of RVE are given in [10] u i = S ij x j + v i. (3) In the above Eq. (3) S ij are the average strains and v i is the periodic part of the displacement components u i on the boundary surfaces (local fluctuation), which is generally unknown and is dependent on the applied global loads. The indices i and j denote the global three-dimensional coordinate directions in the range from 1 to 3. A more explicit form of periodic boundary conditions, suitable for square RVE models can be derived from the above general expression. For the RVE as shown in Fig. 1(b), the displacements on a pair of opposite boundary surfaces (with their normal along the x j axis) are ui K+ = S ij xj K+ + v K+ i, u K i = S ij xj K + v K i (4) where index K + means along the positive x j direction and K means along the negative x j direction on the corresponding surfaces A /A +, B /B +, and C /C + (see Fig. 1(b)). The local fluctuations v K+ i and v K i around the average macroscopic value are identical on two opposing faces due to periodic conditions of RVE. So, the difference between the above two equations is the applied macroscopic strain condition and electric field condition, respectively u K+ i u K i = S ij (xj K+ xj K ), Φ K+ Φ K = Ē i (x K+ j x K j ). (5) It is assumed that the average mechanical and electrical properties of a RVE are equal to the average properties of the particular composite. The average stresses, strains, electric fields, and electrical displacements in a RVE are defined by S ij = 1 S ij dv, T ij = 1 T ij dv, V V V V Ē i = 1 E i dv, D i = 1 D i dv (6) V V V V where V is the volume of the periodic representative volume element Finite element modelling All finite element calculations were made with FE package ANSYS. For modelling of the RVE three-dimensional multifield eight node brick elements with displacement degrees of freedom (DOF) and additional electric potential (voltage) degree of freedom were used. These allows for fully coupled electromechanical analyses. To obtain the homogenised effective properties we apply the macroscopic boundary conditions (Eq. (5)) to the RVE by coupling opposite nodes on opposite boundaries. In order to apply these periodic boundary conditions in the FE analysis, the mesh on the opposite boundary surfaces must be the same. For each pair of displacement components at two corresponding nodes with identical in-plane coordinates on the two opposite boundary surfaces a constraint condition (periodic boundary condition Eq. (5)) is imposed. The prescription of these constraint conditions to all opposite nodes at opposite boundary surfaces directly in ANSYS is a very difficult task because of too many number of nodes. Consequently, we developed a FORTRAN program, which generates all required constraint conditions automatically. First the finite element mesh is created from the ANSYS preprocessor. Then based on the generated nodal coordinates the appropriate nodal pairs are selected by the FORTRAN program and a partial ANSYS input file is created containing the constraint conditions. Using this input file ANSYS continues with assigning the constraint equations and finally with solving the
4 56 H. Berger et al. / Materials Science and Engineering A 412 (2005) such a way that, except the strain in the x 3 direction ( S ), all other mechanical strains and gradients of electric potential (Ē i ) become zero. This can be achieved by constraining the normal displacements at all surfaces to zero except of surface C + (see Fig. 1(b)). At surface C + the periodic boundary condition corresponding to surface C must be applied. Due to applied zero displacements to surface C in x 3 direction (u C 3 = 0) the periodic boundary condition in this direction according to Eq. (5) simplifies to u C+ 3 = S (x3 C+ x3 C ). (7) Fig. 2. Periodic boundary conditions for a pair of nodes on opposite surfaces A and A +. problem. As an example Fig. 2 shows the constraint equations for a pair of nodes on opposite surfaces A and A +. For demonstration of calculation of effective coefficients we consider a piezoceramic (PZT-5) fibre embedded in a soft nonpiezoelectric material (polymer) with a square arrangement like shown in Fig. 1. We are assuming that the fibres and the matrix are ideally bonded and that the fibres are straight and parallel to the x 3 axis. The fibre cross section is circular and the unit cell is having a square cross section. The piezoelectric fibres are uniformly poled along the x 3 direction. To find thective coefficients special load cases with different boundary conditions must be constructed in such a way that for a particular load case only one value in the strain/electric field vector (see Eq. (2)) is non-zero and all others become zero. Then from one row in Eq. (2) the corresponding effective coefficient can be determined using the calculated average non-zero value in the strain/electric field vector and the calculated average values in the stress/electrical displacement vector. In the next chapters, the special models for calculating the different effective coefficients are explained in detail. But due to limitation of this paper only for one load case the FE model and deformed contour plots are shown. The used material properties of polymer and PZT-5 are listed in Table 1, where elastic properties, piezoelectric constants and permittivities are given in N/m 2, C/m 2, and F/m, respectively. 5. Numerical calculation of the different effective coefficients 5.1. Calculation of and Ceff For the calculation of thective coefficients and Ceff the boundary conditions have to be applied to the RVE in Because now this equation is independent of u C 3, instead of using constraint equations, an arbitrary constant prescribed displacement can be applied on surface C + to produce a strain in x 3 direction. To ensure that no charge transfer is between opposite surfaces the voltage degree of freedom on all surfaces is set to zero. For the calculation of the total average values S, T, and T according to Eq. (6) the integral was replaced by a sum over averaged element values multiplied by the respective element volume. Using these total average values the coefficients C eff and can be calculated from the matrix Eq. (2). Due to zero strains and electric fields, except S the first row becomes T = S. Then C44 eff can be calculated as the ratio of T / S. Similarly, can be evaluated as the ratio of T / S from the third row of matrix Eq. (2) Calculation of and Ceff 12 For the calculation of thective coefficients C eff and Ceff 12 we have similar conditions like for the calculation of thective coefficients and Ceff. But now a prescribed displacement in x 1 direction must be applied on surface A + and the normal displacements on all other surfaces must be set to zero. Also the electric potential DOF at all surfaces is set to zero. For the calculation of these coefficients the in-plane behaviour is relevant only. Using the total average strain S and stress T from first row of matrix Eq. (2) we get T = C eff S. From this relation C eff can be calculated as the ratio of T / S. Similarly, C12 eff can be found as the ratio of T 22 / S using the second row of matrix Eq. (2) Calculation of 66 and Ceff 44 To evaluate thective coefficient C66 eff the in-plane shear strain S 12 may have a non-zero value in strain/electric field vector Table 1 Material properties of the composite constituents fibre (PZT-5) and matrix (polymer) C (10 10 ) C 12 (10 10 ) C (10 10 ) C (10 10 ) C 44 (10 10 ) C 66 (10 10 ) e e e ε (10 9 ) ε (10 9 ) PZT Polymer
5 H. Berger et al. / Materials Science and Engineering A 412 (2005) Fig. 3. Finite element mesh with shear loads and out-of-plane supports. of Eq. (2) only. For the evaluation of C44 eff the out-of-plane shear strain S 23 or S 31 may have a non-zero value only. This can be achieved by applying appropriate shear forces to produce a pure shear stress state in the desired directions, which consequently, results in a pure shear strain state. Beside the applied shear forces the displacement boundary conditions must be chosen in the right way to get these states. That means that a symmetric deformation mode respective to the diagonal line of the shear plane must be ensured; see also Fig. 2, which visualizes such a pure shear stress state in a plane view. To receive such pure shear stress state a support normal to the diagonal at one edge of the finite element model was applied. To avoid rigid body movement the opposite edge was fixed in the shear plane. For calculation of C66 eff the x 1 x 2 plane was used as shear plane and for calculation of C44 eff the x 2 x 3 plane was chosen. Fig. 3 shows the finite element model with shear forces and the above-mentioned supports for the out-of-plane case. Furthermore, normal displacements at parallel surfaces to the shear plane must be set to zero to ensure zero strains perpendicular to the shear plane. Also the electric potential DOF on all surfaces must be set to zero. Using the calculated non-zero average strain and stress values from the first model, the sixth row in the matrix Eq. (2) becomes T 12 = C66 eff S 12, and consequently, C66 eff can be evaluated as the ratio of T 12 / S 12. With the second model C44 eff can be found as the ratio of T 23 / S 23. Figs. 4 and 5 show the calculated strains and stresses for this coefficient in deformed shape. Fig. 4. Strain distribution S 23. the electric potential difference in x 3 direction the fibre will try to expand and stresses as well as electric fields are produced in x 3 direction. Consequently, the total average stress and electrical displacement values as well as the total average electric field Ē 3 can be calculated. Then the third row of the matrix Eq. (2) reduces to T = Ē3, so that = T /Ē 3. Similarly, is the ratio of T /Ē 3. For the evaluation of ε eff the last row of matrix Eq. (2) is used which becomes D 3 = ε eff Ē3. From this equation ε eff can be calculated by εeff = D 3 /Ē Calculation of Thective coefficient can be calculated as the ratio of D 2 / S 23 from the seventh row of the matrix Eq. (2). In order to evaluate this coefficient, the boundary conditions should be applied in such a way that except of the strain S 23 in the strain/electric field vector in Eq. (2), all other strains and electric 5.4. Calculation of, eeff, and εeff In order to evaluate thective coefficients, eeff,and εeff, all surfaces must be constrained to have zero normal displacements. A zero electric potential is applied to the surface C and a non-zero electric potential is applied to the surface C +. The electric potential at all other surfaces is set to zero. Since all surfaces are constrained to have zero displacements and because of Fig. 5. Stress distribution T 23.
6 58 H. Berger et al. / Materials Science and Engineering A 412 (2005) Fig. 6. Hexagonal fibre arrangement and unit cell. Fig. 9. Variation of. Fig. 7. Random fibre arrangement and unit cell. fields are made to zero. Similarly, to the calculation of C44 eff outof-plane shear forces are applied in the x 2 x 3 plane. The electric potential in all directions is made to zero by constraining the voltage degree of freedom to zero. Consequently, in the matrix Eq. (2) the eighth row becomes D 2 = S 23. From this equation can be calculated by the ratio eeff = D 2 / S 23. Fig. 10. Variation of Calculation of ε eff ε eff can be evaluated as the ratio of D 2 /Ē 2. In order to solve these equations, all mechanical strains ( S ij ) should be made equal to zero. This can be achieved by constraining the displacements at all surfaces to zero. The electric potential difference should be applied in x 1 direction. Then the eighth row of the matrix Eq. (2) gives the ε eff coefficient. 6. Discussion of results The introduced method was applied to the calculation of piezoelectric fibre composites with square fibre arrangement shown in Fig. 1 and with hexagonal fibre arrangement shown in Fig. 6. Hexagonal arrangement is characterised by a constant distance between one fibre and all surrounding six fibres. So a repeatable part and consequently a unit cell can be found like shown in Fig. 6. A more complex composite with random fibre distribution and fibre diameter is shown in Fig. 7. The unit cell is here an appropriate repeatable part of the composite. Results Fig.. Variation of. are not presented in this paper. But the algorithm is in the same manner applicable like for the regular fibre arrangements. All effective coefficients have been calculated for 6 discrete fibre volume fractions. Fig. 8 demonstrates the used fractions for the square cell. Since only the volume fibre fraction has an influence on the results the size of the RVE was chosen with unit length in all directions. From these discrete fractions graphs were interpolated. For selected coefficients the behaviour related to the volume fraction is shown in the diagrams Figs as comparison between asymptotic homogenisation Fig. 8. Investigated ratios of fibre volume fraction shown for square arrangement.
7 H. Berger et al. / Materials Science and Engineering A 412 (2005) Fig. 12. Variation of 44. Fig. 16. Variation of ε eff. Fig.. Variation of 66. method (AHM) [6,7] and the introduced numerical method (FEM). FEM-Squ denotes the results for square arrangement and FEM-Hex for hexagonal arrangement. The coefficients, and εeff are not plotted because no difference is identifiable between the curves. The results show in general a good coincidence between calculation by AHM and FEM. Thective coefficients, Ceff, Fig. 14. Variation of.,, eeff, eeff, εeff, and εeff calculated by FEM show a very good agreement with the coefficients obtained by AHM in all range of the volume fraction. Qualitatively, the behaviour of the coefficients C eff, Ceff 12, Ceff 66, and Ceff 44 is the same although for a greater value of the volume fraction (0.4 and higher) the curves are not so close. For reliability the results were checked with the universal relation of Schulgasser []. It was satisfied with a maximum error less than 6%. Furthermore, the results were tested against the Hashin/Shtrikman lower and upper bounds [12]. It could be stated that the calculated coefficients lie inside these bounds over the full range of volume fraction. 7. Conclusions A numerical approach for predicting the homogenised properties of piezoelectric fibre composites has been presented. The numerical approach is based on the finite element method. Longitudinal and transversal elastic and piezoelectric effective coefficients have been calculated with the introduced finite element numerical approach and compared with analytical solutions based on the asymptotic homogenisation method. This permits us to estimate the range of validity of our approach. It can be stated that this comprehensive numerical homogenisation technique provides reliable results. In comparison to many other approaches our technique can be easily extended to determination of homogenized material data for more complex composites with arbitrary inclusions and arbitrary distribution (not only uniaxial fibres, also angle plied fibres, ellipsoids, spheres, etc.). The developed tool is based on using a standard FE package like ANSYS by interfacing with corresponding FORTRAN routines. It reduces the manual work and time and can be used as a template to determinective coefficients. Acknowledgement This work has been supported by DFG Graduiertenkolleg 828 Micro Macro Interactions in Structured Media and Particle Systems. References Fig.. Variation of. [1] M.L. Dunn, M. Taya, Int. J. Sol. Struct. 30 (1993) [2] P. Bisegna, R. Luciano, J. Mech. Phys. Solids 44 (1996)
8 60 H. Berger et al. / Materials Science and Engineering A 412 (2005) [3] P. Gaudenzi, Comput. Struct. 65 (1997) [4] M.-L. Feng, C.-C. Wu, Composites Sci. Technol. 61 (2001) [5] J.R. Brockenbrough, S. Suresh, Scr. Metall. Mater. 24 (1990) [6] S.R. Gunawardena, S. Jansson, Acta Metall. Mater. 41 (1993) [7] R. Guinovart-Díaz, J. Bravo-Castillero, R. Rodríguez-Ramos, R. Martínez-Rosado, F. Serranía, M. Navarrete, Mater. Lett. 56 (2002) [8] H. Berger, U. Gabbert, H. Köppe, R. Rodriguez-Ramos, J. Bravo- Castillero, R. Guinovart-Diaz, J.A. Otero, G.A. Maugin, J. Comput. Mech. (2003) [9] E.C.N. Silva, J.S.O. Fonseca, N. Kikuchi, Comput. Methods Appl. Eng. 9 (1998) [10] P. Suquet, in: E. Sanchez-Palencia, A. Zaoui (Eds.), Homogenisation Techniques for Composite Media, Springer-Verlag, Berlin, 1987, pp [] K. Schulgasser, J. Mech. Phys. Solids 40 (1992) [12] Z. Hashin, J. Appl. Mech. 50 (1983)
Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites
Computational Materials Science 39 (2007) 198 204 www.elsevier.com/locate/commatsci Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites S. Kari
More informationPLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by:[otto-von-guericke-universitaet] On: 0 September 007 Access Details: [subscription number 77067] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More informationFINITE ELEMENT ANALYSIS APPLIED TO EVALUATION OF EFFECTIVE MATERIAL COEFFICIENTS FOR PIEZOELECTRIC FIBER COMPOSITES
29 Brazilian Symposium on Aerospace Eng. & Applications Copyright 29 by AAB September 14-16, 29, S. J. Campos, SP, Brazil FINITE ELEMENT ANALYSIS APPLIED TO EALUATION OF EFFECTIE MATERIAL COEFFICIENTS
More informationInternational Journal of Engineering Science
International Journal of Engineering Science 49 (2011) 322 332 Contents lists available at ScienceDirect International Journal of Engineering Science ournal homepage: www.elsevier.com/locate/iengsci On
More informationInternational Journal of Engineering Science
International Journal of Engineering Science 63 (2013) 1 9 Contents lists available at SciVerse ScienceDirect International Journal of Engineering Science journal homepage: www.elsevier.com/locate/ijengsci
More informationRVE Analysis of Light Weight Carbon Nanotubes Embedded Piezoelectric Fibre Composites
Nanoscience and Nanotechnology 2016, 6(1): 11-16 DOI: 10.5923/j.nn.20160601.03 RVE Analysis of Light Weight Carbon Nanotubes Embedded Piezoelectric Fibre Composites V. K. Srivastava 1,*, H. Berger 2, U.
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationEVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD
U.P.B. Sci. Bull., Series D, Vol. 79, Iss. 3, 2017 ISSN 1454-2358 EVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD Mohamed said BOUTAANI 1, Salah
More informationModule 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method
Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without
More informationVariational bounds for anisotropic elastic multiphase composites with different shapes of inclusions
Arch Appl Mech DOI 1.17/s419-8-246-1 ORIGINAL R. Rodríguez-Ramos R. Guinovart-Díaz J. Bravo-Castillero F. J. Sabina H. Berger S. Kari U. Gabbert Variational bounds for anisotropic elastic multiphase composites
More informationA multiscale-multiphysics strategy for numerical modeling of thin piezoelectric sheets
A multiscale-multiphysics strategy for numerical modeling of thin piezoelectric sheets Claudio Maruccio 1*, Laura De Lorenzis 2, Luana Persano 3 and Dario Pisignano 3,4 1 Department of Innovation Engineering,
More informationModule 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains
Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model
More informationELECTROMECHANICAL RESPONSE OF PIEZOELECTRIC FOAMS
18 TH INTRNATIONAL CONFRNC ON COMPOSIT MATRIALS LCTROMCHANICAL RSPONS OF PIZOLCTRIC FOAMS K.S. Challagulla 1 *, T.A. Venkatesh 1 School of ngineering, Laurentian University, Sudbury, Canada, Department
More informationA coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites.
A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites. P.F.Giddings, C.R.Bowen, H.A.Kim University of Bath, UK Dept. Mech. ng, University
More informationEnhancement of magnetoelectric coupling in multiferroic composites via FEM simulation
Enhancement of magnetoelectric coupling in multiferroic composites via FEM simulation *Artjom Avakian 1), Andreas Ricoeur 2) 1), 2) Institute of Mechanics, University of Kassel, Kassel 34125, Germany 1)
More informationMoving screw dislocations in piezoelectric bimaterials
phys stat sol (b) 38 No 1 10 16 (003) / DOI 10100/pssb00301805 Moving screw dislocations in piezoelectric bimaterials Xiang-Fa Wu *1 Yuris A Dzenis 1 and Wen-Sheng Zou 1 Department of Engineering Mechanics
More informationDepartment of Mechanical Engineering, National Institute of Technology Rourkela, Rourkela , India
Composites, Article ID 629175, 12 pages http://dx.doi.org/10.1155/2014/629175 Research Article Effect of Fiber Geometry and Representative Volume Element on Elastic and Thermal Properties of Unidirectional
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f
MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown
More informationNumerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions
TECHNISCHE MECHANIK, 33, 2, (2013), 97 103 submitted: December 11, 2012 Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions R. Glüge, M. Weber
More informationMicromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading
Indian Journal of Engineering & Materials Sciences Vol. 15, October 2008, pp. 382-390 Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading K Sivaji Babu a *, K Mohana
More informationModeling and analysis of the electromechanical behavior of surface-bonded piezoelectric actuators using finite element method
Modeling and analysis of the electromechanical behavior of surface-bonded piezoelectric actuators using finite element method Huangchao Yu and Xiaodong Wang Abstract Piezoelectric actuators have been widely
More informationSimulation of Piezoelectric Induced Lamb Waves in Plates
Simulation of Piezoelectric Induced Lamb Waves in Plates C. WILLBERG, J. M. VIVAR-PEREZ, Z. AHMAD and U. GABBERT ABSTRACT The use of Lamb waves for damage detection and non-destructive evaluation have
More informationEffective properties evaluation for smart composite materials
Universidade de ão Paulo Biblioteca Digital da Produção Intelectual - BDPI Departamento de Engenharia Aeronáutica - EEC/AA Artigos e Materiais de Revistas Científicas - EEC/AA 2012 Effective properties
More informationBENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION
BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION Ahmed Elgamal and Jinchi Lu October 07 Introduction In this study: I) The response
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationValidation of High Displacement Piezoelectric Actuator Finite Element Models
Validation of High Displacement Piezoelectric Actuator Finite Element Models Barmac Taleghani * Army Research Laboratory Vehicle Technology Directorate NASA Langley Research Center Hampton, VA ABSTRACT
More informationThermal deformation compensation of a composite beam using piezoelectric actuators
INSTITUTE OF PHYSICS PUBLISHING Smart Mater. Struct. 13 (24) 3 37 SMART MATERIALS AND STRUCTURES PII: S964-1726(4)7973-8 Thermal deformation compensation of a composite beam using piezoelectric actuators
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationMechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach
Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Sandhyarani Biswas Department of Mechanical Engineering, N.I.T Rourkela, INDIA Abstract: Now-a-days,
More informationEffect of fibre shape on transverse thermal conductivity of unidirectional composites
Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and
More informationExploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description. Igor Dobovšek
Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description Igor Dobovšek University of Ljubljana Faculty of Mathematics and Physics Institute of Mathematics Physics and Mechanics
More informationThickness Optimization of a Piezoelectric Converter for Energy Harvesting
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Thickness Optimization of a Piezoelectric Converter for Energy Harvesting M. Guizzetti* 1, V. Ferrari 1, D. Marioli 1 and T. Zawada 2 1 Dept.
More informationASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS
ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS Dieter H. Pahr and Helmut J. Böhm Institute of Lightweight Design and Structural
More informationApparent thermal conductivity of periodic two-dimensional composites
Computational Materials Science 25 (2002) 329 338 www.elsevier.com/locate/commatsci Apparent thermal conductivity of periodic two-dimensional composites M. Jiang a, I. Jasiuk b, *, M. Ostoja-Starzewski
More informationA FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE
A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,
More informationSupplementary figures
Supplementary figures Supplementary figure 1: The fang model (longitudinal section) used for the stiffness analysis. The base of the fang is considered fixed (i.e. zero displacement constraint) and the
More informationPiezo materials. Actuators Sensors Generators Transducers. Piezoelectric materials may be used to produce e.g.: Piezo materials Ver1404
Noliac Group develops and manufactures piezoelectric materials based on modified lead zirconate titanate (PZT) of high quality and tailored for custom specifications. Piezoelectric materials may be used
More informationAuthors: Correspondence: ABSTRACT:
Multi-Scale Modeling of the Impact and Failure of Fiber Reinforced Polymer Structures using DIGIMAT to LS-DYNA Interface Authors: L. Adam, R. Assaker & R. Ramaya e-xstream engineering S.A. 7, Rue du Bosquet.
More informationDEVELOPMENT OF A MICROMECHANICAL MODEL IN INTERACTION WITH PARAMETERS RELATED TO THE MICROSTRUCTURE OF CARBON/EPOXY COMPOSITES
ECCM16-16 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, 22-26 June 2014 DEVELOPMENT OF A MICROMECHANICAL MODEL IN INTERACTION WITH PARAMETERS RELATED TO THE MICROSTRUCTURE OF CARBON/EPOXY
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationEffect of Crystal Orientation on Piezoelectric Response of Single-Crystal Piezoelectric Layers
Copyright American Scientific Publishers All rights reserved Printed in the United States of America Journal of Computational and Theoretical Nanoscience Vol. 9, 667, Effect of Crystal Orientation on Piezoelectric
More informationANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM BEHAVIOR UNDER THE EFFECT OF EXTERNAL SOLICITATIONS
Third International Conference on Energy, Materials, Applied Energetics and Pollution. ICEMAEP016, October 30-31, 016, Constantine,Algeria. ANALYSIS AND NUMERICAL MODELLING OF CERAMIC PIEZOELECTRIC BEAM
More informationA TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH
9th HSTAM International Congress on Mechanics Limassol, Cyprus, - July, A TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH Cristian Dascalu, Bertrand François, Laboratoire Sols
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationA multiscale approach for composite materials as multifield continua
Materials Science Forum Vols. 539-543 (27) pp. 2551-2556 online at http://www.scientific.net (27) Trans Tech Publications, Switzerland A multiscale approach for composite materials as multifield continua
More informationMicromechanical modeling and simulation of piezoceramic materials
Micromechanical modeling and simulation of piezoceramic materials B. Delibas 1, A. Arockia Rajan 1 & W. Seemann 2 1 Workgroup for Machine Dynamics, Technical University of Kaiserslautern, Germany 2 Institut
More informationThermal load-induced notch stress intensity factors derived from averaged strain energy density
Available online at www.sciencedirect.com Draft ScienceDirect Draft Draft Structural Integrity Procedia 00 (2016) 000 000 www.elsevier.com/locate/procedia 21st European Conference on Fracture, ECF21, 20-24
More informationShock enhancement of cellular structures under impact loading: Part II Analysis
Shock enhancement of cellular structures under impact loading: Part II Analysis S. Pattofatto, 1 I. Elnasri, 1 H. Zhao, 1* H. Tsitsiris 1 F. Hild 1 and Y. Girard 2 1 Laboratoire de Mécanique et Technologie
More informationAddition 1. Shear Stack Piezoelectric Elements and Shear Effect Basics
120 Addition 1 Shear Stack Piezoelectric Elements and Shear Effect Basics Introduction The STM scanner built up in this work is a Besocke type scanner (see room temperature STM instrumental chapter). The
More informationPIEZOELECTRIC TECHNOLOGY PRIMER
PIEZOELECTRIC TECHNOLOGY PRIMER James R. Phillips Sr. Member of Technical Staff CTS Wireless Components 4800 Alameda Blvd. N.E. Albuquerque, New Mexico 87113 Piezoelectricity The piezoelectric effect is
More informationPROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS
PROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS Abstract Matthew T. Tilbrook, Robert J. Moon and Mark Hoffman School of Materials Science and Engineering University of New South Wales,
More informationAnalytical homogenization method for periodic composite materials
Analytical homogenization method for periodic composite materials Ying Chen and Christopher A. Schuh* Department of Materials Science and Engineering, Massachusetts Institute of Technology, 77 Massachusetts
More informationHomogenization methods for multi-phase elastic composites: Comparisons and benchmarks
TECHNISCHE MECHANIK, 3, 4, (), 374 386 submitted: October 6, 9 Homogenization methods for multi-phase elastic composites: Comparisons and benchmarks B. Klusemann, B. Svendsen Usually homogenization methods
More informationFinite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements
Proceedings o the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. Finite Element Modeling o Residual Thermal Stresses in Fiber-Reinorced Composites Using Dierent Representative
More informationAccepted Manuscript. Numerical simulation of the guided Lamb wave propagation in particle reinforced
Accepted Manuscript Numerical simulation of the guided Lamb wave propagation in particle reinforced composites Ralf Weber, Seyed Mohammad Hossein Hosseini, Ulrich Gabbert PII: S0263-8223(12)00180-8 DOI:
More informationTHE ANALYSIS OF FEM RESULTS CONVERGENCE IN MODELLING PIEZOELECTRIC ACTIVE SHELL STRUCTURES
Dragan Marinković Manfred Zehn Zoran Marinković ISSN 1333-1124 eissn 1849-1391 THE ANALYSIS OF FEM RESULTS CONVERGENCE IN MODELLING PIEZOELECTRIC ACTIVE SHELL STRUCTURES Summary UDC 624.04+531:519.6 The
More informationPiezoelectric Bimorph Response with Imperfect Bonding Conditions
Copyright c 28 ICCES ICCES, vol.6, no.3, pp.5-56 Piezoelectric Bimorph Response with Imperfect Bonding Conditions Milazzo A., Alaimo A. and Benedetti I. Summary The effect of the finite stiffness bonding
More informationA Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core
Commun. Theor. Phys. 56 774 778 Vol. 56, No. 4, October 5, A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core JIANG
More informationGeneralized Plane Piezoelectric Problem: Application to Heterostructure Nanowires
Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires H. T. Mengistu and A. Garcia-Cristobal Institute of Material Science, University of Valencia, P.O.Box 22085, E-46071, Valencia,
More informationMICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING
MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING N. Krishna Vihari 1, P. Phani Prasanthi 1, V. Bala Krishna Murthy 2* and A. Srihari Prasad 3 1 Mech. Engg. Dept., P. V. P. Siddhartha
More informationNUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING
NUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING Evangelos V. Liarakos Postdoctoral researcher School of Architecture, Technical
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
More informationTesting and analysis of high frequency electroelastic characteristics of piezoelectric transformers
Arch. Mech., 59, 2, pp. 119 131, Warszawa 2007 Testing and analysis of high frequency electroelastic characteristics of piezoelectric transformers F. NARITA, Y. SHINDO, F. SAITO, M. MIKAMI Department of
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationAnalysis of high loss viscoelastic composites
Analysis of high loss viscoelastic composites by C. P. Chen, Ph.D. and R. S. Lakes, Ph.D. Department of Engineering Physics Engineering Mechanics Program; Biomedical Engineering Department Materials Science
More informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationStructural Health Monitoring Using Smart Piezoelectric Material
Structural Health Monitoring Using Smart Piezoelectric Material Kevin K Tseng and Liangsheng Wang Department of Civil and Environmental Engineering, Vanderbilt University Nashville, TN 37235, USA Abstract
More informationCellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).
1 Cellular solid structures with unbounded thermal expansion Roderic Lakes Journal of Materials Science Letters, 15, 475-477 (1996). Abstract Material microstructures are presented which can exhibit coefficients
More informationFinite element analysis of longitudinal debonding between fibre and matrix interface
Indian Journal of Engineering & Materials Sciences Vol. 11, February 2004, pp. 43-48 Finite element analysis of longitudinal debonding between fibre and matrix interface K Aslantaş & S Taşgetiren Department
More informationPiezoelectric composite lattices with high sensitivity
Piezoelectric composite lattices with high sensitivity Roderic Lakes Department of Engineering Physics, Engineering Mechanics Program, Department of Materials Science, Rheology Research Center, University
More informationDae Won Ji and Sang-Joo Kim
Journal of the Korean Ceramic Society Vol. 55, No. 1, pp. 7~73, 2018. https://doi.org/10.4191/kcers.2018.55.1.03 Communication Construction and Application of Experimental Formula for Nonlinear Behavior
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More informationEstimating the overall properties of heterogeneous Cosserat materials
Modelling Simul. Mater. Sci. Eng. 7 (1999) 829 840. Printed in the UK PII: S0965-0393(99)03041-7 Estimating the overall properties of heterogeneous Cosserat materials S Forest, R Dendievel and G R Canova
More information3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials
3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials Abstract: Talapady S. Bhat and T. A. Venkatesh Department of Material Science and Engineering Stony Brook University,
More informationA direct evaluation of the Fabric Tensor in anisotropic porous media
A direct evaluation of the Fabric Tensor in anisotropic porous media Maria Cristina Pernice 1, Luciano Nunziante 1, Massimiliano Fraldi 1,2 1 Department of Structural Engineering, University of Naples
More informationGeneric Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic
More informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
More informationMechanics of Irregular Honeycomb Structures
Mechanics of Irregular Honeycomb Structures S. Adhikari 1, T. Mukhopadhyay 1 Chair of Aerospace Engineering, College of Engineering, Swansea University, Singleton Park, Swansea SA2 8PP, UK Sixth International
More informationDegenerated shell element for geometrically nonlinear analysis of thin-walled piezoelectric active structures
IOP PUBLISHING Smart Mater. Struct. 17 (2008) 015030 (10pp) SMART MATERIALS AND STRUCTURES doi:10.1088/0964-1726/17/01/015030 Degenerated shell element for geometrically nonlinear analysis of thin-walled
More informationIntelligent Materials and their applications. Autors: Jana Pintea, Ionut Balan INCDIE ICPE-CA Bucharest, ROMANIA
Intelligent Materials and their applications Autors: Jana Pintea, Ionut Balan INCDIE ICPE-CA Bucharest, ROMANIA What are intelligent materials? The intelligence of a material is a characteristic difficult
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationP program, version 2. Theoretical manual & User guide
Documentation for H E L P program, version 2 Theoretical manual & User guide Jan Vorel 1 May 22, 2015 1 Czech Technical University in Prague, Faculty of Civil Engineering, Department of Mechanics Contents
More informationStress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear
Xi an 2-25 th August 217 Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension compression or shear Heyin Qi 1 Mingming Chen 2 Yonghong Duan 3 Daxu
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationNanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation. Mark Hardiman
Nanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation Mark Hardiman Materials and Surface Science Institute (MSSI), Department of Mechanical and Aeronautical
More informationVariational Asymptotic Method for Unit Cell Homogenization
Variational Asymptotic Method for Unit Cell Homogenization Wenbin Yu 1 and Tian Tang 2 1 Utah State University, Logan, Utah 84322-4130 Wenbin.Yu@usu.edu 2 Utah State University, Logan, Utah 84322-4130
More informationDesign of in-plane piezoelectric sensors for static response by simultaneously optimizing the host structure and the electrode profile
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Design of in-plane piezoelectric sensors for static response by simultaneously optimizing the
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011
Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.
More informationZero Energy Modes in One Dimension: An Introduction to Hourglass Modes
Zero Energy Modes in One Dimension: An Introduction to Hourglass Modes David J. Benson March 9, 2003 Reduced integration does a lot of good things for an element: it reduces the computational cost, it
More informationTensile behaviour of anti-symmetric CFRP composite
Available online at www.sciencedirect.com Procedia Engineering 1 (211) 1865 187 ICM11 Tensile behaviour of anti-symmetric CFRP composite K. J. Wong a,b, *, X. J. Gong a, S. Aivazzadeh a, M. N. Tamin b
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationTransduction Based on Changes in the Energy Stored in an Electrical Field
Lecture 7-1 Transduction Based on Changes in the Energy Stored in an Electrical Field - Electrostriction The electrostrictive effect is a quadratic dependence of strain or stress on the polarization P
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationInfluence of the cover plate thickness on the Lamb wave propagation in honeycomb sandwich panels
CEA Aeronaut J (2013) 4:69 76 DOI 10.1007/s13272-012-0054-8 ORIGINAL PAPER Influence of the cover plate thickness on the Lamb wave propagation in honeycomb sandwich panels. M. H. Hosseini U. Gabbert R.
More informationInterfacial effects in electromagnetic coupling within piezoelectric phononic crystals
Acta Mech Sin (29) 25:95 99 DOI 1.17/s149-8-21-y RESEARCH PAPER Interfacial effects in electromagnetic coupling within pieoelectric phononic crystals F. J. Sabina A. B. Movchan Received: 14 July 28 / Accepted:
More informationStudies on effective fiber and matrix poling characteristics of 1-3 piezoelectric composites
INTERNATIONAL JOURNAL OF STRUCTURAL CHANGES IN SOLIDS Mechanics and Applications Volume 3, Number 3, 211, pp. 1-19 Studies on eective iber and matrix poling characteristics o 1-3 piezoelectric composites
More information1 Introduction. Abstract
Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical
More information