Planar auxeticity from various inclusions

Size: px
Start display at page:

Download "Planar auxeticity from various inclusions"

Transcription

1 Planar auxeticity from various inclusions Artur A. Poźniak 1 Krzysztof W. Wojciechowski 2 Joseph N. Grima 3 Luke Mizzi 4 1 Institute of Physics, Poznań Universty of Technology, Piotrowo 3, Poznań 2 Institute of Molecular Physics PAS, M. Smoluchowskiego 17, Poznań 3 Department of Chemistry, Faculty of Science, University of Malta, Msida, MSD 28 4 Metamaterials Unit, Faculty of Science, University of Malta, Msida, MSD th September 215 AUXETICS 215 Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

2 Motivation What is an auxetic? ν < (Negtative Poisson s Ratio) is for mechanics what negative refractive index is for optics n <. Auxetics are metamaterials. Beneficial features from NPR Resistance to shape change and indentation; crack resistance; better vibration absorption (including acoustic one); synclastic curvature, different dynamics. Applications of NPR materials Medicine (stents, bandages, implants), defence (energy absorption), furniture industry (better mattresses indentation), automotive industry and sports (safety belts). Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

3 Poisson s ratio (stretching along x axis) W H Formal definition ν xy = ǫ yy ǫ xx = s yyxx s xxxx W H Engineering definition ν = H H W W Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

4 Motivation rotating units Grima, J. N., Alderson, A., & Evans, K. E. (25). Auxetic behaviour from rotating rigid units. Physica Status Solidi (B), 242(3), Grima, J. N., Chetcuti, E., Manicaro, E., Attard, D., Camilleri, M., Gatt, R., & Evans, K. E. (211). On the auxetic properties of generic rotating rigid triangles. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 468(2139), Grima, J., Zammit, V., Gatt, R., Alderson, A., & Evans, K. (27). Auxetic behaviour from rotating semi-rigid units. Physica Status Solidi (B), 244(3), Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

5 Motivation rotating units Grima, J. N., Alderson, A., & Evans, K. E. (25). Auxetic behaviour from rotating rigid units. Physica Status Solidi (B), 242(3), Grima, J. N., Chetcuti, E., Manicaro, E., Attard, D., Camilleri, M., Gatt, R., & Evans, K. E. (211). On the auxetic properties of generic rotating rigid triangles. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 468(2139), Grima, J., Zammit, V., Gatt, R., Alderson, A., & Evans, K. (27). Auxetic behaviour from rotating semi-rigid units. Physica Status Solidi (B), 244(3), Our result from 21 Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

6 Motivation rotating units Grima, J. N., Alderson, A., & Evans, K. E. (25). Auxetic behaviour from rotating rigid units. Physica Status Solidi (B), 242(3), Grima, J. N., Chetcuti, E., Manicaro, E., Attard, D., Camilleri, M., Gatt, R., & Evans, K. E. (211). On the auxetic properties of generic rotating rigid triangles. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 468(2139), Grima, J., Zammit, V., Gatt, R., Alderson, A., & Evans, K. (27). Auxetic behaviour from rotating semi-rigid units. Physica Status Solidi (B), 244(3), Our result from 21 Taylor, M., Francesconi, L., Gerendás, M., Shanian, A., Carson, C., & Bertoldi, K. (214). Advanced Materials, 26(15), Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

7 Motivation rotating units Grima, J. N., Alderson, A., & Evans, K. E. (25). Auxetic behaviour from rotating rigid units. Physica Status Solidi (B), 242(3), Grima, J. N., Chetcuti, E., Manicaro, E., Attard, D., Camilleri, M., Gatt, R., & Evans, K. E. (211). On the auxetic properties of generic rotating rigid triangles. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 468(2139), Grima, J., Zammit, V., Gatt, R., Alderson, A., & Evans, K. (27). Auxetic behaviour from rotating semi-rigid units. Physica Status Solidi (B), 244(3), Our result from 21 Taylor, M., Francesconi, L., Gerendás, M., Shanian, A., Carson, C., & Bertoldi, K. (214). Advanced Materials, 26(15), Why not fill the gaps? Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

8 Structures a) Geometry: Periodically arranged ellipses Alternating axes (perpendicular) b) Repeating unit cell (RUC): 4 rotating units Periodic mesh c) Boundary conditions within reduced RUC All boundaries remain straight lines Less computational resources needed Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

9 Simulations Finite Elements 1 Abaqus/CAE via Python interface generates inp files (wc -l ellipses-uniform-mesh.py gives 64 lines), 2 Abaqus/STANDARD performs simulations generating odb files, 3 Abaqus/CAE extracts important information from each odb file. One run is enough to calculate ν as well as E in one direction: V y the displacement of (an arbitrary) top node E elas elastic deformation energy T denotes the plane-stress thickness. U x known deformation ν eff = V y/h U x /W. (1) E eff = 2 Eelas W T U 2 x H. (2) Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

10 Meshes: uniform vs. non-uniform The discretization (meshing) is an essential and most challenging step on FE analysis. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

11 Meshes: uniform vs. non-uniform The discretization (meshing) is an essential and most challenging step on FE analysis. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

12 Mesh & convergence Figure: Nonuniform exemplary mesh. Each edge was assigned a bias seed in order to obtain possibly dense mesh in the narrowest (bottleneck) regions. In this case L gap/d min = 8 was chosen in order to make he figure readable. In the simulations this ratio was 1. Fig. on the right shows that the convergence tests were performed up to the value d min = 1 4 Figure: Convergence of Poisson s ratio as a function of d min at fixed d max =.1. Here ν matr = ν incl =.99, E matr = 1. The geometric parameters are: a =.9, b =.89. The convergence displays various characteristics depending on the value of E incl. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

13 Mesh & convergence Figure: Nonuniform exemplary mesh. Each edge was assigned a bias seed in order to obtain possibly dense mesh in the narrowest (bottleneck) regions. In this case L gap/d min = 8 was chosen in order to make he figure readable. In the simulations this ratio was 1. Fig. on the right shows that the convergence tests were performed up to the value d min = 1 4 Figure: Convergence of Poisson s ratio as a function of d min at fixed d max =.1. Here ν matr = ν incl =.99, E matr = 1. The geometric parameters are: a =.9, b =.89. The convergence displays various characteristics depending on the value of E incl. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

14 Convergence E incl /E matr Figure: Convergence of ν eff and E eff /E matr as obtained within non-uniform meshing for a =.9, b =.89, ν incl = ν matr =.99 and increasing mesh density. The upper, monotonic curves refer to the right (Young s modulus, red) axis, while the non monotonic curves refer to the left (Poisson s ratio, blue) axis. The numbers in the legend stand for the L gap/d min ratio. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23 Eeff/Ematr

15 Detailed convergence ν / n l E incl = 1 4 uniform mesh Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

16 Detailed convergence ν E incl = 1 4 non-uniform mesh uniform mesh 1/ n l d max PR PR dmin dmin PR dmin PR dmin Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

17 Detailed convergence ν E incl = 1 4 non-uniform mesh uniform mesh 1/ n l d max 1 PR PR dmin dmin PR dmin PR dmin Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

18 Detailed convergence E /n l Eeff/Ematr E incl = 1 4 non-uniform mesh uniform mesh d max Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

19 Required time vs. nuber of FE layers UNIFORM N-UNIFORM, d max =.2 N-UNIFORM, d max =.1 N-UNIFORM, d max = t [s] L gap /d min or L gap /d u Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

20 Memory consumed vs. nuber of FE layers 1 9 UNIFORM N-UNIFORM, d max =.2 N-UNIFORM, d max =.1 N-UNIFORM, d max =.5 memory [MB] L gap /d min or L gap /d u Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

21 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

22 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

23 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

24 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

25 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

26 Results Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

27 E incl limits, anisotropy L gap=.1 L gap=.2 L gap=.5 L gap=.1 L gap= a/ b Void limit, E incl c) The dependence of Poisson s ratio as a function of a/b for various values of L gap. Here both ν matr and ν incl are equal. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

28 E incl limits, anisotropy L gap=.1 L gap=.2 L gap=.5 L gap=.1 L gap= a/ b Void limit, E incl L gap =.1 L gap =.2 L gap =.5 L gap =.1 L gap = a/b Rigid limit, E incl c) The dependence of Poisson s ratio as a function of a/b for various values of L gap. Here both ν matr and ν incl are equal. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

29 Rigid inclusions, f (L gap ) a/ b=.5 a/ b=.1 a/ b=.5 a/ b=.1 a/ b=.2 a/ b=.3 a/ b= L gap Figure: The dependence of Poisson s ratio as a function of L gap for various values of a/b within the rigid inclusion limit. Here both ν matr and ν incl are equal. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

30 Figure: The dependence of effective Poisson s ratio as a function of the matrix Poisson s ratio in the limit of infinitely hard inclusions (E incl ). The value of ν incl in this limit has no impact on the effective properties since the inclusions do not deform. The inset figure helps in estimating of the a/b ratio, for which ν eff changes sign. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23 Rigid inclusions, f (ν matr ) a =.9,b =.89 a =.14,b =.84 a =.19,b =.79 a =.24,b =.74 a =.29,b =.69 a =.34,b =.64 a =.39,b =.59 a =.44,b =.54 a =.49,b =.49 E incl νmatr = log(a/ b) ν matr

31 Figure: Effective Poisson s ratio, ν eff, (blue colour) and effective Young s modulus, E eff /E matr, (red colour) as functions of E incl /E matr for strongly auxetic matrices and inclusions. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23 More and more auxetic constituents ν matr =.9,ν incl =.9 ν matr =.99,ν incl =.99 ν matr =.999,ν incl =.999 ν matr =.9999,ν incl = E incl /E matr Eeff/Ematr

32 Large deformations % 5% 1% ε 15% 2% Poisson s ratio Effective Poisson s ratio, ν eff, and effective Young s modulus, E eff /E matr as functions of nominal strain for ν matr = and void inclusions (the material of inclusions is entirely removed). The numbers in the legend correspond to the following geometry parameter pairs: 1) a =.9,.89, 2) a =.14,.84, 3) a =.19,.79, 4) a =.24,.74, 5) a =.29,.69, 6) a =.34,.64, 7) a =.39,.59, 8) a =.44,.54, 9) a =.49,.49. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

33 Large deformations % 5% 1% ε 15% 2% Poisson s ratio Eeff/Ematr % 5% 1% ε 15% 2% Young s modulus Effective Poisson s ratio, ν eff, and effective Young s modulus, E eff /E matr as functions of nominal strain for ν matr = and void inclusions (the material of inclusions is entirely removed). The numbers in the legend correspond to the following geometry parameter pairs: 1) a =.9,.89, 2) a =.14,.84, 3) a =.19,.79, 4) a =.24,.74, 5) a =.29,.69, 6) a =.34,.64, 7) a =.39,.59, 8) a =.44,.54, 9) a =.49,.49. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

34 Star-like inclusions, parameters a and b a =.5, b =.945 Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

35 Star-like inclusions, parameters a and b a =.5, b =.945 Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

36 Conclusions For strongly anisotropic inclusions the effective Poisson s ratio, which is close to for very low E incl /E matr, grows and after reaching a maximum (of negative value) decays to the matrix values close to E incl /E matr 1 and then grows again to the large E incl /E matr limit. void homogenous rigid Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

37 Conclusions For strongly anisotropic inclusions the effective Poisson s ratio, which is close to for very low E incl /E matr, grows and after reaching a maximum (of negative value) decays to the matrix values close to E incl /E matr 1 and then grows again to the large E incl /E matr limit. In the vicinity of E incl /E matr for which the effective PR shows its maximum, for the most auxetic inclusions the effective Young s modulus is essentially larger than for less auxetic inclusions. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

38 Conclusions For strongly anisotropic inclusions the effective Poisson s ratio, which is close to for very low E incl /E matr, grows and after reaching a maximum (of negative value) decays to the matrix values close to E incl /E matr 1 and then grows again to the large E incl /E matr limit. In the vicinity of E incl /E matr for which the effective PR shows its maximum, for the most auxetic inclusions the effective Young s modulus is essentially larger than for less auxetic inclusions. For highly anisotropic inclusions of very large Young s modulus, the effective PR of the composite can be negative for nonauxetic matrix and inclusions. These is a very simple example of an auxetic structure being not only entirely continuous, i.e. without empty regions (cuts, holes) introduced into the system, but with E incl E matr, i.e. inclusions (much) harder than the matrix Lgap =.1 Lgap =.2 Lgap =.5 Lgap =.1 Lgap = a/ b Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

39 The end Thank you for your attention This work was supported by the (Polish) National Centre for Science under the grant NCN 212/5/N/ST5/1476. Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

40 Poźniak, Wojciechowski, Grima, Mizzi Eeffν/ eff Ematr Eeffν/ eff Ematr Eeffν/ eff Ematr Eeffν/ eff Ematr νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = Planar auxeticity from various inclusions νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = νincl =.99 νincl =.49 νincl =. νincl = Sep / 23

41 ν ν ν a=.9 b=.89, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 a=.24 b=.74, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 a=.39 b=.59, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 ν ν ν a=.14 b=.84, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 a=.29 b=.69, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 a=.44 b=.55, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 ν ν a=.19 b=.79, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 a=.34 b=.64, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl = νmat =.99 νmat =.49 νmat =. νmat =.49 νmat =.99 νmat =.49 νmat =. νmat =.49 νmat =.99 νmat =.49 νmat =. νmat = νmat =.99 νmat =.49 νmat =. νmat =.49 νmat =.99 νmat =.49 νmat =. νmat =.49 νmat =.99 νmat =.49 νmat =. νmat = νmat =.99 νmat =.49 νmat =. νmat = νmat =.99 νmat =.49 νmat =. νmat =.49 ν a=.49 b=.49, gwiazdki νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νincl =.999 νincl =.99 νincl =.9 νincl =.49 νincl =. νincl =.49 νmat =.99 νmat =.49 νmat =. νmat =.49 Poźniak, Wojciechowski, Grima, Mizzi Planar auxeticity from various inclusions 15 Sep / 23

Poisson s ratio of rectangular anti-chiral lattices with disorder

Poisson s ratio of rectangular anti-chiral lattices with disorder Poisson s ratio of rectangular anti-chiral lattices with disorder Artur A. Poźniak 1 Krzysztof W. Wojciechowski 2 1 Institute of Physics, Poznań Universty of Technology, ul. Nieszawska 13A, 60-695 Poznań,

More information

ANOMALOUS DEFORMATION OF CONSTRAINED AUXETIC SQUARE

ANOMALOUS DEFORMATION OF CONSTRAINED AUXETIC SQUARE Anomalous Rev.Adv.Mater.Sci. deformation 23(2010) of constrained 169-174 auxetic square 169 ANOMALOUS DEFORMATION OF CONSTRAINED AUXETIC SQUARE Artur A. PoŸniak 1, Henryk Kamiñski 2, Przemys³aw Kêdziora

More information

548 Advances of Computational Mechanics in Australia

548 Advances of Computational Mechanics in Australia Applied Mechanics and Materials Online: 2016-07-25 ISSN: 1662-7482, Vol. 846, pp 547-552 doi:10.4028/www.scientific.net/amm.846.547 2016 Trans Tech Publications, Switzerland Geometric bounds for buckling-induced

More information

Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements

Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements Acta Mech. Sin. (2016) 32(2):295 300 DOI 10.1007/s10409-015-0534-2 RESEARCH PAPER Analysis of the geometrical dependence of auxetic behavior in reentrant structures by finite elements V. H. Carneiro 1

More information

ENGN 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 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 information

Finite Element Method in Geotechnical Engineering

Finite 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 information

Extreme Cosserat elastic cube structure with large magnitude of negative Poisson s ratio

Extreme Cosserat elastic cube structure with large magnitude of negative Poisson s ratio Extreme Cosserat elastic cube structure with large magnitude of negative Poisson s ratio Carlos Andrade Chan Soo Ha +,androderics.lakes? December 7, 2018 Department of Engineering Physics, Department of

More information

Computer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature

Computer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature Computer simulation of the Poisson's ratio of soft polydisperse discs at zero temperature JAKUB NAROJCZYK KRZYSZTOF W. WOJCIECHOWSKI Institute of Molecular Physics Polish Academy of Sciences ul. M. Smoluchowskiego

More information

Three-dimensional stiff cellular structures with negative Poisson's ratio

Three-dimensional stiff cellular structures with negative Poisson's ratio ! Three-dimensional stiff cellular structures with negative Poisson's ratio Dong Li a, Jie Ma a, Liang Dong b and Roderic S. Lakes c a College of Sciences, Northeastern University, Shenyang 110819, PR

More information

Minimization of Poisson s ratio in anti-tetra-chiral two-phase structure

Minimization of Poisson s ratio in anti-tetra-chiral two-phase structure IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Minimization of Poisson s ratio in anti-tetra-chiral two-phase structure To cite this article: E Idczak and T Strek 2017 IOP Conf.

More information

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.

Example 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60. 162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides

More information

FINITE ELEMENT APPROACH AND MATHEMATICAL FORMULATION OF VISCOELASTIC AUXETIC HONEYCOMB STRUCTURES FOR IMPACT MITIGATION

FINITE ELEMENT APPROACH AND MATHEMATICAL FORMULATION OF VISCOELASTIC AUXETIC HONEYCOMB STRUCTURES FOR IMPACT MITIGATION Journal of Engineering Science and Technology Vol. 1, No. (017) 471-490 School of Engineering, Taylor s University FINITE ELEMENT APPROACH AND MATHEMATICAL FORMULATION OF VISCOELASTIC AUXETIC HONEYCOMB

More information

Mechanical Properties of Materials

Mechanical Properties of Materials Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of

More information

Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure

Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure Vibrations in Physical Systems Vol. 27 (2016) Abstract Computational Modelling of Vibrations Transmission Loss of Auxetic Lattice Structure Eligiusz IDCZAK Institute of Applied Mechanics, Poznan University

More information

arxiv: v1 [physics.comp-ph] 12 Apr 2015

arxiv: v1 [physics.comp-ph] 12 Apr 2015 arxiv:1504.03021v1 [physics.comp-ph] 12 Apr 2015 Elastic properties of the degenerate f.c.c. crystal of polydisperse soft dimers at zero temperature J. W. Narojczyk 1 and K. W. Wojciechowski 2 Institute

More information

FINITE ELEMENT ANALYSIS OF EFFECTIVE MECHANICAL PROPERTIES, VIBRATION AND ACOUSTIC PERFORMANCE OF AUXETIC CHIRAL CORE SANDWICH STRUCTURES

FINITE ELEMENT ANALYSIS OF EFFECTIVE MECHANICAL PROPERTIES, VIBRATION AND ACOUSTIC PERFORMANCE OF AUXETIC CHIRAL CORE SANDWICH STRUCTURES Clemson University TigerPrints All Theses Theses 8-2013 FINITE ELEMENT ANALYSIS OF EFFECTIVE MECHANICAL PROPERTIES, VIBRATION AND ACOUSTIC PERFORMANCE OF AUXETIC CHIRAL CORE SANDWICH STRUCTURES Hrishikesh

More information

Computational Analysis for Composites

Computational 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 information

p s s On the auxetic properties of rotating rhombi and parallelograms: A preliminary investigation solidi physica status

p s s On the auxetic properties of rotating rhombi and parallelograms: A preliminary investigation solidi physica status phys. stat. sol. () 45, No. 3, 5 59 (8) / DOI./pss.77775 On the auxetic properties of rotating rhomi and parallelograms: A preliminary investigation physica Joseph N. Grima *, Pierre-Sandre Farrugia, Ruen

More information

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design Fracture mechanics fundamentals Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design Failure modes Failure can occur in a number of modes: - plastic deformation

More information

Forced Response of Plate with Viscoelastic Auxetic Dampers

Forced Response of Plate with Viscoelastic Auxetic Dampers Vibrations in Physical Systems 2018, 29, 2018003 (1 of 9) Abstract Forced Response of Plate with Viscoelastic Auxetic Dampers Tomasz STREK Poznan University of Technology, Institute of Applied Mechanics

More information

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST

More information

The in-plane linear elastic constants and out-ofplane bending of 3-coordinated ligament and cylinder-ligament honeycombs.

The in-plane linear elastic constants and out-ofplane bending of 3-coordinated ligament and cylinder-ligament honeycombs. University of Bolton UBIR: University of Bolton Institutional Repository IMRI: Journal Articles (Peer-Reviewed) Institute for Materials Research and Innovation 200 The in-plane linear elastic constants

More information

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution

More information

Basic Equations of Elasticity

Basic Equations of Elasticity A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 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 information

UWE has obtained warranties from all depositors as to their title in the material deposited and as to their right to deposit such material.

UWE has obtained warranties from all depositors as to their title in the material deposited and as to their right to deposit such material. Mohanraj, H., Filho Ribeiro, S. L. M., Panzera, T. H., Scarpa, F., Farrow, I. R., Jones, R., Davies-Smith, A., Remillat, C. D. L., Walters, P. and Peng, H.-X. (2016) Hybrid auxetic foam and perforated

More information

Towards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr.

Towards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr. Towards Efficient Finite Element Model Review Dr. Richard Witasse, Plaxis bv (based on the original presentation of Dr. Brinkgreve) Journée Technique du CFMS, 16 Mars 2011, Paris 1/32 Topics FEA in geotechnical

More information

ICM11. Homogenization of periodic auxetic materials

ICM11. Homogenization of periodic auxetic materials Available online at www.sciencedirect.com Procedia Engineering 10 (2011) 1847 1852 ICM11 Homogenization of periodic auxetic materials J. Dirrenberger a,, S. Forest a, D. Jeulin a,b, C. Colin a a Centre

More information

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Elasto Plastic 3D Finite Element Contact Analysis of a Hole Containing a Circular Insert in a Fatigue Test Coupon

Elasto Plastic 3D Finite Element Contact Analysis of a Hole Containing a Circular Insert in a Fatigue Test Coupon Elasto Plastic 3D Finite Element Contact Analysis of a Hole Containing a Circular Insert in a Fatigue Test Coupon Witold Waldman Aerospace Division Defence Science and Technology Group DST-Group-TR-3140

More information

DYNAMIC MODELING OF SPOT WELDS USING THIN LAYER INTERFACE THEORY

DYNAMIC MODELING OF SPOT WELDS USING THIN LAYER INTERFACE THEORY Hamid Ahmadian: 1 DYNAMIC MODELING OF SPOT WELDS USING THIN LAYER INTERFACE THEORY H. Ahmadian 1, H. Jalali 1, JE Mottershead 2, MI Friswell 3 1 Iran University of Science and Technology, Faculty of Mechanical

More information

Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Cellular 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 information

BECAS - an Open-Source Cross Section Analysis Tool

BECAS - an Open-Source Cross Section Analysis Tool BECAS - an Open-Source Cross Section Analysis Tool José P. Blasques and Robert D. Bitsche Presented at DTU Wind Energy stand at the EWEA 2012 conference, Copenhagen, 16.4.2012 BECAS-DTUWind@dtu.dk Motivation

More information

Elasticity in two dimensions 1

Elasticity in two dimensions 1 Elasticity in two dimensions 1 Elasticity in two dimensions Chapters 3 and 4 of Mechanics of the Cell, as well as its Appendix D, contain selected results for the elastic behavior of materials in two and

More information

Plane and axisymmetric models in Mentat & MARC. Tutorial with some Background

Plane and axisymmetric models in Mentat & MARC. Tutorial with some Background Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,

More information

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of

More information

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925

More information

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2 Samantha Ramirez, MSE Stress The intensity of the internal force acting on a specific plane (area) passing through a point. Δ ΔA Δ z Δ 1 2 ΔA Δ x Δ y ΔA is an infinitesimal size area with a uniform force

More information

Low Porosity Metallic Periodic Structures with Negative Poisson's Ratio

Low Porosity Metallic Periodic Structures with Negative Poisson's Ratio Low Porosity Metallic Periodic Structures with Negative Poisson's Ratio The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters. Citation

More information

VIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES. 1. Introduction

VIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES. 1. Introduction ARCHIVES OF ACOUSTICS 31, 4 (Supplement), 53 58 (2006) VIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES J. CIEŚLIK, W. BOCHNIAK AGH University of Science and Technology Department of Robotics and Mechatronics

More information

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering

Alternative 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 information

Structural behaviour of traditional mortise-and-tenon timber joints

Structural behaviour of traditional mortise-and-tenon timber joints Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University

More information

COMPRESSION AND BENDING STIFFNESS OF FIBER-REINFORCED ELASTOMERIC BEARINGS. Abstract. Introduction

COMPRESSION AND BENDING STIFFNESS OF FIBER-REINFORCED ELASTOMERIC BEARINGS. Abstract. Introduction COMPRESSION AND BENDING STIFFNESS OF FIBER-REINFORCED ELASTOMERIC BEARINGS Hsiang-Chuan Tsai, National Taiwan University of Science and Technology, Taipei, Taiwan James M. Kelly, University of California,

More information

International Journal of Scientific & Engineering Research Volume 8, Issue 12, December-2017 ISSN

International Journal of Scientific & Engineering Research Volume 8, Issue 12, December-2017 ISSN 362 Study of the Effects of Random Inclusions in Composite Construction with Isotropic Negative Poisson s Ratio [Kundan Singh Deshwal1, a and Dr. Satyendra Singh2, b] 1 Assistant Professor, Department

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

Macroscopic theory Rock as 'elastic continuum'

Macroscopic theory Rock as 'elastic continuum' Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave

More information

Lecture 8. Stress Strain in Multi-dimension

Lecture 8. Stress Strain in Multi-dimension Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]

More information

Module I: Two-dimensional linear elasticity. application notes and tutorial. Problems

Module I: Two-dimensional linear elasticity. application notes and tutorial. Problems Module I: Two-dimensional linear elasticity application notes and tutorial Problems 53 selected excerpts from Read Me file for: ElemFin 1.1.1 Yannick CALLAUD in Symantec C++. 1 place of Falleron, 44300

More information

Prediction of Elastic Constants on 3D Four-directional Braided

Prediction 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 information

Constitutive models. Constitutive model: determines P in terms of deformation

Constitutive models. Constitutive model: determines P in terms of deformation Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function

More information

Crack Tip Plastic Zone under Mode I Loading and the Non-singular T zz -stress

Crack Tip Plastic Zone under Mode I Loading and the Non-singular T zz -stress Crack Tip Plastic Zone under Mode Loading and the Non-singular T -stress Yu.G. Matvienko Mechanical Engineering Research nstitute of the Russian Academy of Sciences Email: ygmatvienko@gmail.com Abstract:

More information

HPLP300 - Plate with Young modulus function of the temperature

HPLP300 - Plate with Young modulus function of the temperature Titre : HPLP300 - Plaque aec module d Young fonction de l[...] Date : 04/08/2011 Page : 1/6 HPLP300 - Plate with Young modulus function of the temperature Summarized: This thermoelastic test makes it possible

More information

International Journal of Advanced Engineering Technology E-ISSN

International Journal of Advanced Engineering Technology E-ISSN Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,

More information

IDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE

IDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE IDENTIFICATION OF THE ELASTIC PROPERTIES ON COMPOSITE MATERIALS AS A FUNCTION OF TEMPERATURE Hugo Sol, hugos@vub.ac.be Massimo Bottiglieri, Massimo.Bottiglieri@vub.ac.be Department Mechanics of Materials

More information

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric

More information

GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE

GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE GLOBAL AND LOCAL LINEAR BUCKLING BEHAVIOR OF A CHIRAL CELLULAR STRUCTURE Alessandro Spadoni, Massimo Ruzzene School of Aerospace Engineering Georgia Institute of Technology Atlanta, GA 30332 Fabrizio Scarpa

More information

Exercise: concepts from chapter 5

Exercise: concepts from chapter 5 Reading: Fundamentals of Structural Geology, Ch 5 1) Study the oöids depicted in Figure 1a and 1b. Figure 1a Figure 1b Figure 1. Nearly undeformed (1a) and significantly deformed (1b) oöids with spherulitic

More information

Modelling Anisotropic, Hyperelastic Materials in ABAQUS

Modelling Anisotropic, Hyperelastic Materials in ABAQUS Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,

More information

1 332 Laboratories 1. 2 Computational Exercises 1 FEA of a Cantilever Beam... 1 Experimental Laboratory: Tensile Testing of Materials...

1 332 Laboratories 1. 2 Computational Exercises 1 FEA of a Cantilever Beam... 1 Experimental Laboratory: Tensile Testing of Materials... 1 332 Laboratories Contents 1 332 Laboratories 1 2 Computational Exercises 1 FEA of a Cantilever Beam.......................................... 1 Experimental Laboratory: Tensile Testing of Materials..........................

More information

Advanced Structural Analysis EGF Section Properties and Bending

Advanced Structural Analysis EGF Section Properties and Bending Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear

More information

Finite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure

Finite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure Vibrations in Physical Systems Vol.26 (2014) Finite Element Analysis of Dynamic Properties of Thermally Optimal Two-phase Composite Structure Abstract Maria NIENARTOWICZ Institute of Applied Mechanics,

More information

Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions

Numerical 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 information

TC211 Workshop CALIBRATION OF RIGID INCLUSION PARAMETERS BASED ON. Jérôme Racinais. September 15, 2015 PRESSUMETER TEST RESULTS

TC211 Workshop CALIBRATION OF RIGID INCLUSION PARAMETERS BASED ON. Jérôme Racinais. September 15, 2015 PRESSUMETER TEST RESULTS Jérôme Racinais September 15, 215 TC211 Workshop CALIBRATION OF RIGID INCLUSION PARAMETERS BASED ON PRESSUMETER TEST RESULTS Table of contents 1. Reminder about pressuremeter tests 2. General behaviour

More information

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.

More information

BENCHMARK 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 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 information

ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A BODY STRUCTURE

ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A BODY STRUCTURE The 3rd International Conference on Computational Mechanics and Virtual Engineering COMEC 2009 29 30 OCTOBER 2009, Brasov, Romania ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

More information

Homogenization on Multi-Materials Elements: Application to Printed Circuit Boards and Warpage Analysis

Homogenization on Multi-Materials Elements: Application to Printed Circuit Boards and Warpage Analysis Homogenization on Multi-Materials Elements: Application to Printed Circuit Boards and Warpage Analysis Manuel Araújo 1,a, J. L. Alves 1, Paulo Silva 1 and Pedro Delgado 2 1 MEMS, University of Minho, Dep.

More information

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference

More information

Stress analysis of a stepped bar

Stress analysis of a stepped bar Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.

More information

ME 475 Modal Analysis of a Tapered Beam

ME 475 Modal Analysis of a Tapered Beam ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory

More information

Nomenclature. Length of the panel between the supports. Width of the panel between the supports/ width of the beam

Nomenclature. Length of the panel between the supports. Width of the panel between the supports/ width of the beam omenclature a b c f h Length of the panel between the supports Width of the panel between the supports/ width of the beam Sandwich beam/ panel core thickness Thickness of the panel face sheet Sandwich

More information

A 3 D finite element model for sound transmission through a double plate system with isotropic elastic porous materials

A 3 D finite element model for sound transmission through a double plate system with isotropic elastic porous materials Acoustics and Vibrations Group Université de Sherbrooke, QC CANADA Département génie mécanique Université de Sherbrooke Sherbrooke, QC CANADA Tel.: (819) 821-7157 Fax: (819) 821-7163 A 3 D finite element

More information

DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD

DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,

More information

ENGINEERING TRIPOS PART IIA 3C7: EXPERIMENTAL STRESS ANALYSIS

ENGINEERING TRIPOS PART IIA 3C7: EXPERIMENTAL STRESS ANALYSIS ENGINEERING TRIPOS PART IIA 3C7: EXPERIMENTAL STRESS ANALYSIS Experiment takes place in BNB-06 (follow downward stairs opposite Baker Building reception). OBJECTIVES To develop an appreciation of two different

More information

Assessment of boundaryelement method for modelling the structural response of a pavement

Assessment of boundaryelement method for modelling the structural response of a pavement Assessment of boundaryelement method for modelling the structural response of a pavement (PhD Student, Lecturer DEC of FCTUC) (Full Professor - IST) Pavements: Materials, design and performance a, LNEC,

More information

Keywords: Re-entrant; cellular solid; Negative Poisson s ratio; Voronoi; Energy absorption; Compression

Keywords: Re-entrant; cellular solid; Negative Poisson s ratio; Voronoi; Energy absorption; Compression Negative Poisson's ratio in 2D Voronoi cellular solids by biaxial compression: a numerical study Preprint: Li, D. Dong, L., Yin, J. and Lakes, R. S., "Negative Poisson's ratio in 2D Voronoi cellular solids

More information

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS Djoković, J. M., et.al.: Influence of Temperature on Behavior of the Interfacial THERMAL SCIENCE: Year 010, Vol. 14, Suppl., pp. S59-S68 S59 INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK

More information

When materials are stretched in a particular direction, they

When materials are stretched in a particular direction, they pubs.acs.org/nanolett Intrinsic Negative Poisson s Ratio for Single-Layer Graphene Jin-Wu Jiang,*, Tienchong Chang, Xingming Guo, and Harold S. Park Shanghai Institute of Applied Mathematics and Mechanics,

More information

MECHANICS OF CHIRAL, ANTI-CHIRAL, AND HIERARCHICAL HONEYCOMBS

MECHANICS OF CHIRAL, ANTI-CHIRAL, AND HIERARCHICAL HONEYCOMBS MECHANICS OF CHIRAL, ANTI-CHIRAL, AND HIERARCHICAL HONEYCOMBS A Dissertation Presented By DAVOOD MOUSANEZHAD VIYAND to The Department of Mechanical and Industrial Engineering in partial fulfillment of

More information

rsc.li/soft-matter-journal As featured in: Highlighting research from the University of Wisconsin Madison, and Technion.

rsc.li/soft-matter-journal As featured in: Highlighting research from the University of Wisconsin Madison, and Technion. Highlighting research from the University of Wisconsin Madison, and Technion. As featured in: Auxetic multiphase soft composite material design through instabilities with application for acoustic metamaterials

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft 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 information

INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES

INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES INFLUENCE OF THE LOCATION AND CRACK ANGLE ON THE MAGNITUDE OF STRESS INTENSITY FACTORS MODE I AND II UNDER UNIAXIAL TENSION STRESSES Najah Rustum Mohsin Southern Technical University, Technical Institute-Nasiriya,

More information

[5] Stress and Strain

[5] Stress and Strain [5] Stress and Strain Page 1 of 34 [5] Stress and Strain [5.1] Internal Stress of Solids [5.2] Design of Simple Connections (will not be covered in class) [5.3] Deformation and Strain [5.4] Hooke s Law

More information

Thermal stresses in hybrid materials with auxetic inclusions

Thermal stresses in hybrid materials with auxetic inclusions Thermal stresses in hybrid materials with auxetic inclusions Elena Pasternak a,, Igor Shufrin a, b, 1 and Arcady V. Dyskin b, 3 a School of Mechanical and Chemical Engineering, The University of Western

More information

ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH

ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN

More information

Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost

Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost Comparison of Ply-wise Stress-Strain results for graphite/epoxy laminated plate subjected to in-plane normal loads using CLT and ANSYS ACP PrepPost 1 Mihir A. Mehta, 2 Satyen D. Ramani 1 PG Student, Department

More information

Module 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy

Module 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy The Lecture Contains: Torsion of Rods Torsional Energy This lecture is adopted from the following book 1. Theory of Elasticity, 3 rd edition by Landau and Lifshitz. Course of Theoretical Physics, vol-7

More information

3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials

3-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 information

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT

UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT 2004 AIMETA International Tribology Conference, September 14-17, 2004, Rome, Italy UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT Yuri KLIGERMAN( ), Yuri Kadin( ), Izhak ETSION( ) Faculty of

More information

Fabrication and Testing of 2-Dimensional Z-Expandable Auxetic Textile. Structures for Impact Protective Clothing Applications

Fabrication and Testing of 2-Dimensional Z-Expandable Auxetic Textile. Structures for Impact Protective Clothing Applications 2D Z-EXPANDABLE AUXETIC TEXTILES FOR IMPACT PROTECTION 1 Fabrication and Testing of 2-Dimensional Z-Expandable Auxetic Textile Structures for Impact Protective Clothing Applications ASTM Student Project

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

ABHELSINKI UNIVERSITY OF TECHNOLOGY

ABHELSINKI UNIVERSITY OF TECHNOLOGY ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural

More information

Analytical Modeling of the Stress-Strain Distribution in a Multilayer Structure with Applied Bending

Analytical Modeling of the Stress-Strain Distribution in a Multilayer Structure with Applied Bending Analytical Modeling of the Stress-Strain Distribution in a Multilayer Structure with Applied Bending Ana Neves Vieira da Silva Prof. Doutora Susana Isabel Pinheiro Cardoso de Freitas Dep. Physics, IST,

More information

3D Elasticity Theory

3D Elasticity Theory 3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.

More information

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002 student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity

More information

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Module 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 information

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

Module 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 information