Researcher profile

Sundararajan Natarajan

Sundararajan Natarajan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2023arXiv

On the fractional transversely isotropic functionally graded nature of soft biological tissues

This paper focuses on the origin of the poroelastic anisotropic behaviour of the meniscal tissue and its spatially varying properties. We present confined compression creep test results on samples extracted from three parts of the tissue (Central body, Anterior horn and Posterior horn) in three orientations (Circumferential, Radial and Vertical). We show that a poroelastic model in which the fluid flow evolution is ruled by non-integer order operators (fractional Darcy's law) provides accurate agreement with the experimental creep data. The model is validated against two additional sets of experimental data: stress relaxation and fluid loss during the consolidation process measured as weight reduction. Results show that the meniscus can be considered as a transversely isotropic poroelastic material. This behaviour is due to the fluid flow rate being about three times higher in the circumferential direction than in the radial and vertical directions in the body region of the meniscus. In the anterior horn, the elastic properties are transversely isotropic, with the aggregate modulus higher in the radial direction than in the circumferential and vertical directions. The 3D fractional poroelastic model is implemented in finite element software and quantities such as flux of interstitial fluid during the consolidation process, a non-trivial experimental measure, are determined.

preprint2020arXiv

Adaptive phase field modelling of crack propagation in orthotropic functionally graded materials

In this work, we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials (FGMs). A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement. The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency. The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations. The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle. The study is then extended to the analysis of orthotropic FGMs. It is observed that, if the gradation in fracture properties is neglected, the material gradient plays a secondary role, with the fracture behaviour being dominated by the orthotropy of the material. However, when the toughness increases along the crack propagation path, a substantial gain in fracture resistance is observed.

preprint2020arXiv

Continuum modelling of stress diffusion interactions in an elastoplastic medium in the presence of geometric discontinuity

Chemo-mechanical coupled systems have been a subject of interest for many decades now. Previous attempts to solve such models have mainly focused on elastic materials without taking into account the plastic deformation beyond yield, thus causing inaccuracies in failure calculations. This paper aims to study the effect of stress-diffusion interactions in an elastoplastic material using a coupled chemo-mechanical system. The induced stress is dependent on the local concentration in a one way coupled system, and vice versa in a two way coupled system. The time-dependent transient coupled system is solved using a finite element formulation in an open-source finite element solver FEniCS. This paper attempts to computationally study the interaction of deformation and diffusion and its effect on the localization of plastic strain. We investigate the role of geometric discontinuities in scenarios involving diffusing species, namely, a plate with a notch/hole/void and particle with a void/hole/core. We also study the effect of stress concentrations and plastic yielding on the diffusion-deformation. The developed code can be from https://github.com/mrupeshkumar/Elastoplastic-stress-diffusion-coupling

preprint2013arXiv

Finite element computations on quadtree meshes: strain smoothing and semi-analytical formulation

This short communication discusses two alternate techniques to treat hanging nodes in a quadtree mesh. Both the techniques share similarities, in that, they require only boundary information. Moreover, they do not require an explicit form of the shape functions, unlike the conventional approaches, for example, as in the work of Gupta \cite{gupta1978} or Tabarraei and Sukumar \cite{tabarraeisukumar2005}. Hence, no special numerical integration technique is required. One of the techniques relies on the strain projection procedure, whilst the other is based on the scaled boundary finite element method. Numerical examples are presented to demonstrate the accuracy and the convergence properties of the two techniques.

preprint2013arXiv

Smoothed Finite Element and Genetic Algorithm based optimization for Shape Adaptive Composite Marine Propellers

An optimization scheme using the Cell-based Smoothed Finite Element Method (CS-FEM) combined with a Genetic Algorithm (GA) framework is proposed in this paper to design shape adaptive laminated composite marine propellers. The proposed scheme utilise the bend-twist coupling characteristics of the composites to achieve the required performance. An iterative procedure to evaluate the unloaded shape of the propeller blade is proposed, confirming the manufacturing requirements at the initial stage. The optimization algorithm and codes developed in this work were implemented under a variety of parameter settings and compared against the requirement to achieve an ideally passive pitch varying propeller. Recommendations for the required thickness of the propeller blade to achieve optimum bend-twist coupling performance without resulting large rake deformations are also presented.

preprint2013arXiv

Vibration of functionally graded material plates with cutouts & cracks in thermal environment

In this paper, the effect of a centrally located cutout (circular and elliptical) and cracks emanating from the cutout on the free flexural vibration behaviour of functionally graded material plates in thermal environment is studied. The discontinuity surface is represented independent of the mesh by exploiting the partition of unity method framework. A Heaviside function is used to capture the jump in the displacement across the discontinuity surface and asymptotic branch functions are used to capture the singularity around the crack tip. An enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction. The effective properties of the functionally graded material are estimated using the Mori- Tanaka homogenization scheme and the plate kinematics is based on the first order shear deformation theory. The influence of the plate geometry, the geometry of the cutout, the crack length, the thermal gradient and the boundary conditions on the free flexural vibration is numerically studied.

preprint2012arXiv

Efficient recovery-based error estimation for the smoothed finite element method for smooth and singular linear elasticity

An error control technique aimed to assess the quality of smoothed finite element approximations is presented in this paper. Finite element techniques based on strain smoothing appeared in 2007 were shown to provide significant advantages compared to conventional finite element approximations. In particular, a widely cited strength of such methods is improved accuracy for the same computational cost. Yet, few attempts have been made to directly assess the quality of the results obtained during the simulation by evaluating an estimate of the discretization error. Here we propose a recovery type error estimator based on an enhanced recovery technique. The salient features of the recovery are: enforcement of local equilibrium and, for singular problems a "smooth+singular" decomposition of the recovered stress. We evaluate the proposed estimator on a number of test cases from linear elastic structural mechanics and obtain precise error estimations whose effectivities, both at local and global levels, are improved compared to recovery procedures not implementing these features.

preprint2011arXiv

Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method (XFEM) allow discontinuities to be simulated independently of the mesh [1]. This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains [2] to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem, show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code.

preprint2011arXiv

On the approximation in the smoothed finite element method (SFEM)

This letter aims at resolving the issues raised in the recent short communication [1] and answered by [2] by proposing a systematic approximation scheme based on non-mapped shape functions, which both allows to fully exploit the unique advantages of the smoothed finite element method (SFEM) [3, 4, 5, 6, 7, 8, 9] and resolve the existence, linearity and positivity deficiencies pointed out in [1]. We show that Wachspress interpolants [10] computed in the physical coordinate system are very well suited to the SFEM, especially when elements are heavily distorted (obtuse interior angles). The proposed approximation leads to results which are almost identical to those of the SFEM initially proposed in [3]. These results that the proposed approximation scheme forms a strong and rigorous basis for construction of smoothed finite element methods.