IMT Institutional Repository: No conditions. Results ordered -Date Deposited. 2024-03-28T12:24:10ZEPrintshttp://eprints.imtlucca.it/images/logowhite.pnghttp://eprints.imtlucca.it/2017-09-26T07:57:29Z2017-09-26T07:57:29Zhttp://eprints.imtlucca.it/id/eprint/3799This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/37992017-09-26T07:57:29ZA nonlocal adaptive discrete empirical interpolation method combined with modified hp-refinement for order reduction of molecular dynamics systemsModel order reduction is an emerging technique to tackle the computational complexities of molecular
dynamics (MD) simulations. Different strategies are required to adequately obtain the reduced solutions
of different classes of molecular dynamics systems. In this work, a proper orthogonal decomposition
(POD) is combined with the discrete empirical interpolation method (DEIM) to study atomic systems.
Due to the limitations of the DEIM in capturing the nonlocal response of the nonlinear force field of
MD systems, a nonlocal adaptive discrete empirical interpolation method (ADEIM) is proposed.
Furthermore, a modified hp-refinement algorithm is introduced to extend the application of the PODDEIM
approach to order reduction of multi-dimensional MD systems. In the DEIM, the distance between
atoms and hence the reduced internal force vector is estimated based on a local interpolation of the state
variables. The internal forces of a multi-dimensional MD system depend on the distance between the
atoms, represented in space by more than one coordinate. Therefore, the ADEIM approach seeks to obtain
a nonlocal interpolation of the state variables to accurately predict the distance between the interpolated
atoms and hence the reduced force vector. Simulation of MD systems with frequently changing neighbour
atoms leads to change in the system dynamics, which further leads to change of properties of the
snapshots. Therefore, the temporal domain is adaptively subdivided into smaller sub-domains using
the adopted hp-refinement procedure. The reduced system parameters are effectively derived over the
sub-domains. Considering the computational cost, a modified hp-refinement algorithm is developed in
this study, which is further coupled with the POD-ADEIM approach to obtain the reduced-order solution
of the MD systems. The results of the proposed approach demonstrate the efficiency and accuracy of the
reduced solutions.Saheed Olalekan Ojosaheed.ojo@imtlucca.itPattabhi R. Budarapupattabhi.budarapu@imtlucca.itMarco Paggimarco.paggi@imtlucca.it2017-01-30T11:29:16Z2017-09-21T14:32:15Zhttp://eprints.imtlucca.it/id/eprint/3648This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/36482017-01-30T11:29:16ZConcurrently coupled solid shell based adaptive multiscale method for fractureA solid shell-based adaptive atomistic–continuum numerical method is herein proposed to simulate complex crack growth patterns in thin-walled structures. A hybrid solid shell formulation relying on the combined use of the enhanced assumed strain (EAS) and the assumed natural strain (ANS) methods has been considered to efficiently model the material in thin structures at the continuum level. The phantom node method (PNM) is employed to model the discontinuities in the bulk. The discontinuous solid shell element is then concurrently coupled with a molecular statics model placed around the crack tip. The coupling between the coarse scale and the fine scale is realized through the use of ghost atoms, whose positions are interpolated from the coarse scale solution and enforced as boundary conditions to the fine scale model. In the proposed numerical scheme, the fine scale region is adaptively enlarged as the crack propagates and the region behind the crack tip is adaptively coarsened in order to reduce the computation costs. An energy criterion is used to detect the crack tip location. All the atomistic simulations are carried out using the LAMMPS software. A computational framework has been developed in MATLAB to trigger LAMMPS through system command. This allows a two way interaction between the coarse and fine scales in MATLAB platform, where the boundary conditions to the fine region are extracted from the coarse scale, and the crack tip location from the atomistic model is transferred back to the continuum scale. The developed framework has been applied to study crack growth in the energy minimization problems. Inspired by the influence of fracture on current–voltage characteristics of thin Silicon photovoltaic cells, the cubic diamond lattice structure of Silicon is used to model the material in the fine scale region, whilst the Tersoff potential function is employed to model the atom–atom interactions. The versatility and robustness of the proposed methodology is demonstrated by means of several fracture applications.Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itJosé ReinosoMarco Paggimarco.paggi@imtlucca.it2017-01-30T10:33:52Z2017-01-30T11:33:18Zhttp://eprints.imtlucca.it/id/eprint/3647This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/36472017-01-30T10:33:52ZLattice orientation and crack size effect on the mechanical properties of GrapheneThe effect of lattice orientation and crack length on the mechanical properties of Graphene are studied based on molecular dynamics simulations. Bond breaking and crack initiation in an initial edge crack model with 13 different crack lengths, in 10 different lattice orientations of Graphene are examined. In all the lattice orientations, three recurrent fracture patterns are reported. The influence of the lattice orientation and crack length on yield stress and yield strain of Graphene is also investigated. The arm-chair fracture pattern is observed to possess the lowest yield properties. A sudden decrease in yield stress and yield strain can be noticed for crack sizes <10 nm. However, for larger crack sizes, a linear decrease in yield stress is observed, whereas a constant yield strain of ≈≈0.05 is noticed. Therefore, the yield strain of ≈≈0.05 can be considered as a critical strain value below which Graphene does not show failure. This information can be utilized as a lower bound for the design of nano-devices for various strain sensor applications. Furthermore, the yield data will be useful while developing the Graphene coating on Silicon surface in order to enhance the mechanical and electrical characteristics of solar cells and to arrest the growth of micro-cracks in Silicon cells.Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itB. JavvajiV. K. SutrakarD. Roy MahapatraMarco Paggimarco.paggi@imtlucca.itG. ZiT. Rabczuk2016-12-27T09:06:58Z2016-12-27T09:06:58Zhttp://eprints.imtlucca.it/id/eprint/3617This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/36172016-12-27T09:06:58ZMechanical properties of Graphene: Molecular dynamics simulations correlated to continuum based scaling lawsIn this paper, the combined effect of domain size, lattice orientation and crack length on the mechanical properties of Graphene, namely the yield strength and strain, are studied extensively based on molecular dynamics simulations. Numerical predictions are compared with the continuum-based laws of size effect and multifractal scaling. The yield strength is found to vary with the specimen size as ≈L^{−1/3}, which is in agreement with the multifractal scaling law, and with the inverse square of the initial crack length as ≈a^{-1/2}, according to the Griffith’s energy criterion for fracture.Brahmanandam JavvajiPattabhi R. Budarapupattabhi.budarapu@imtlucca.itV. K. SutrakarD. Roy MahapatraMarco Paggimarco.paggi@imtlucca.itGoangseup ZiTimon Rabczuk2016-01-20T15:45:46Z2016-04-06T09:41:29Zhttp://eprints.imtlucca.it/id/eprint/3027This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/30272016-01-20T15:45:46ZCrack propagation in grapheneThe crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated.Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itBrahmanandam JavvajiV. K. SutrakarD. Roy MahapatraGoangseup ZiTimon Rabczuk2015-03-27T11:14:59Z2015-03-27T12:54:14Zhttp://eprints.imtlucca.it/id/eprint/2651This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26512015-03-27T11:14:59ZAero-Elastic Analysis of Stiffened Composite Wing StructureThe interaction of the Elastic, Inertia and Aerodynamic forces is a dynamic phenomena resulting in flutter. Dynamic Aero elasticity is critical for a high-speed subsonic class of AerialVehicles. In the presentwork, dynamic Aero-elastic analysis of the Unmanned Aerial Vehicles (UAVs) has been studied for a stiffened composite structure. The flutter speed and the corresponding flutter frequencies are computed using the Velocity-Damping (V -g) method. The V -g method is employed to estimate the flutter speed and flutter frequencies, for a high-speed subsonic composite wing structure, designed and developed at ADE. Also the improvement of the flutter frequencies over the existing metallic wing structure is discussed. The work can be further extended to develop an optimized composite structure with higher margins in flutter speeds.
Keywords: Aeroelasticity, Flutter, Velocity-damping methodPattabhi R. Budarapupattabhi.budarapu@imtlucca.itB. RammohanS. Vijay KumarD. Satish Babu R. Raghuatnhan2015-03-27T10:54:35Z2016-03-18T10:43:50Zhttp://eprints.imtlucca.it/id/eprint/2650This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26502015-03-27T10:54:35ZAn adaptive multiscale method for quasi-static crack growthThis paper proposes an adaptive atomistic- continuum numerical method for quasi-static crack growth. The phantom node method is used to model the crack in the continuum region and a molecular statics model is used near the crack tip. To ensure self-consistency in the bulk, a virtual atom cluster is used to model the material of the coarse scale. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively enlarged as the crack propagates and the region behind the crack tip is adaptively coarsened. An energy criterion is used to detect the crack tip location. The triangular lattice in the fine scale region corresponds to the lattice structure of the (111) plane of an FCC crystal. The Lennard-Jones potential is used to model the atom–atom interactions. The method is implemented in two dimensions. The results are compared to pure atomistic simulations; they show excellent agreement.Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itRobert GracieStéphane P.A. BordasTimon Rabczuk2015-03-27T10:51:27Z2015-03-27T12:54:14Zhttp://eprints.imtlucca.it/id/eprint/2649This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26492015-03-27T10:51:27ZEfficient coarse graining in multiscale modeling of fracture Abstract We propose a coarse-graining technique to reduce a given atomistic model into an equivalent coarse grained continuum model. The developed technique is tailored for problems involving complex crack patterns in 2D and 3D including crack branching and coalescence. Atoms on the crack surface are separated from the atoms not on the crack surface by employing the centro symmetry parameter. A rectangular grid is superimposed on the atomistic model. Atoms on the crack surface in each cell are used to estimate the equivalent coarse-scale crack surface of that particular cell. The crack path in the coarse model is produced by joining the approximated crack paths in each cell. The developed technique serves as a sound basis to study the crack propagation in multiscale methods for fracture. Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itRobert GracieShih-Wei YangXiaoying ZhuangTimon Rabczuk2015-03-27T10:40:25Z2015-03-27T12:54:13Zhttp://eprints.imtlucca.it/id/eprint/2648This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26482015-03-27T10:40:25ZA meshless adaptive multiscale method for fracture Abstract The paper presents a multiscale method for crack propagation. The coarse region is modelled by the differential reproducing kernel particle method. Fracture in the coarse scale region is modelled with the Phantom node method. A molecular statics approach is employed in the fine scale where crack propagation is modelled naturally by breaking of bonds. The triangular lattice corresponds to the lattice structure of the (111)plane of an {FCC} crystal in the fine scale region. The Lennard–Jones potential is used to model the atom–atom interactions. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively refined and coarsened as the crack propagates. The centro symmetry parameter is used to detect the crack tip location. The method is implemented in two dimensions. The results are compared to pure atomistic simulations and show excellent agreement. Shih-Wei YangPattabhi R. Budarapupattabhi.budarapu@imtlucca.itD. Roy MahapatraStéphane P.A. BordasGoangseup ZiTimon Rabczuk2015-03-27T10:38:14Z2016-04-06T08:49:47Zhttp://eprints.imtlucca.it/id/eprint/2647This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26472015-03-27T10:38:14ZVibration analysis of multi-walled carbon nanotubes embedded in elastic mediumWe propose a method to estimate the natural frequencies of the multi-walled carbon nanotubes (MWCNTs) embedded in an elastic medium. Each of the nested tubes is treated as an individual bar interacting with the adjacent nanotubes through the inter-tube Van der Waals forces. The effect of the elastic medium is introduced through an elastic model. The mathematical model is finally reduced to an eigen value problem and the eigen value problem is solved to arrive at the inter-tube resonances of the MWCNTs. Variation of the natural frequencies with different parameters are studied. The estimated results from the present method are compared with the literature and results are observed to be in close agreement.Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itSudhirSastry YbBrahmanandam JavvajiD. Roy Mahapatra2015-03-27T10:35:00Z2015-03-27T12:54:13Zhttp://eprints.imtlucca.it/id/eprint/2646This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26462015-03-27T10:35:00ZParametric studies on buckling of thin walled channel beams Abstract The lateral buckling analysis of cold-formed thin walled beams subjected to pure bending moments has been performed. The critical buckling loads are estimated based an optimization criteria. The estimated critical buckling stresses are compared with the published results, they show excellent agreement. The effect of the beam length, radius and thickness of the flanges and the length of the extended open flanges, on the critical buckling stresses have been studied for several combinations of the geometric parameters of the beam. Among the three beams, the critical buckling moments for the beam with the extended open flanges are found to be the maximum. However, considering the material and manufacturing costs, beams with rounded cross section are efficient in resisting the buckling loads. Y.B. SudhirSastryY. KrishnaPattabhi R. Budarapupattabhi.budarapu@imtlucca.it2015-03-27T10:32:33Z2015-03-27T12:54:13Zhttp://eprints.imtlucca.it/id/eprint/2645This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26452015-03-27T10:32:33ZBuckling analysis of thin wall stiffened composite panels Abstract We present the pre and post buckling analysis of stiffened composite panels based on the finite element models. Individual buckling studies are conducted on the stiffened composite panel madeup of woven fabric CFC/epoxy, E-glass/epoxy and the Kevlar/epoxy composites. Straight, T shaped and I shaped stiffeners are considered to stiffen the panel. The panel is fabricated with 8 layers and the stiffeners are madeup of 16 layers, of equal thickness arranged in different orientations. The panel is subjected to a uniform axial compression load of 10 kN. The distribution of the buckling stresses and the buckling loads with different panel and stiffener combinations are estimated for three different layup sequences. The variation of the buckling stresses and the buckling loads from the numerical model are compared with the experiment. The results are in excellent agreement. Y.B. SudhirSastryPattabhi R. Budarapupattabhi.budarapu@imtlucca.itN. MadhaviY. Krishna2015-03-27T10:27:08Z2015-03-27T12:54:14Zhttp://eprints.imtlucca.it/id/eprint/2644This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26442015-03-27T10:27:08ZStudies on ballistic impact of the composite panels Abstract The ballistic impact of the composite materials is studied using the numerical models. Individual impact studies are conducted on the composite plate made-up of woven fabric CFRP, E-glass/epoxy and the Kevlar/epoxy composites. The plate is fabricated with 8 layers of equal thickness arranged in different orientations. A spherical steel projectile is considered for the high velocity impact. The projectile is placed very close to the plate, at the center and impacted with a velocity of 100 m/s. The displacement and the stress distribution in each layer are studied for the layup sequence +45/−45/+45/−45/−45/+45/−45/+45. The variation of the kinetic energy, the increase in the internal energy of the laminate and the decrease in velocity of the projectile with time are also studied. Based on the results, the best layup sequence for the ballistic impact of each material is suggested. Y.B. Sudhir SastryPattabhi R. Budarapupattabhi.budarapu@imtlucca.itY. KrishnaS. Devaraj2015-03-27T10:18:27Z2015-03-27T12:54:13Zhttp://eprints.imtlucca.it/id/eprint/2643This item is in the repository with the URL: http://eprints.imtlucca.it/id/eprint/26432015-03-27T10:18:27ZDirectionality of sound radiation from rectangular panels Abstract In this paper, the directionality of sound radiated from a rectangular panel, attached with masses/springs, set in a baffle, is studied. The attachment of masses/springs is done based on the receptance method. The receptance method is used to generate new mode shapes and natural frequencies of the coupled system, in terms of the old mode shapes and natural frequencies. The Rayleigh integral is then used to compute the sound field. The point mass/spring locations are arbitrary, but chosen with the objective of attaining a unique directionality. The excitation frequency to a large degree decides the sound field variations. However, the size of the masses and the locations of the masses/springs do influence the new mode shapes and hence the sound field. The problem is more complex when the number of masses/springs are increased and/or their values are made different. The technique of receptance method is demonstrated through a steel plate with attached point masses in the first example. In the second and third examples, the present method is applied to estimate the sound field from a composite panel with attached springs and masses, respectively. The layup sequence of the composite panel considered in the examples corresponds to the multifunctional structure battery material system, used in the micro air vehicle (MAV) (Thomas and Qidwai, 2005). The demonstrated receptance method does give a reasonable estimate of the new modes. Pattabhi R. Budarapupattabhi.budarapu@imtlucca.itT.S.S. NarayanaB. RammohanTimon Rabczuk