We present a generalization of Signorini's method to the case of live loads which allows us to derive approximate solutions to some pure traction-value problems in finite elastostatics. The boundary-value problems and the corresponding compatibility conditions are formulated in order to determine the displacement of the system up to the second-order of approximation. In particular, we consider the case of homogeneous and isotropic elastic bodies and we solve the following two traction-value problems with live loads:(i) a sphere subjected to the action of a uniform pressure field;(ii)a hollow circular cylinder whose inner and outer surfaces are subjected to uniform pressures. Then, starting from these solutions, we suggest experiments to determine the second-order constitutive constants of the elastic body. Expressions of the second-order material constants in terms of displacements and Lame' coefficients are determined. / Master of Science

The singular solutions for linear elastostatics at corners are essential in both theory and computation. In this thesis, we seek new singular solutions for corners with the fixed (displacement), the free stress (traction) boundary conditions, and their mixed types, and to explore their corner singularity and provide the algorithms and error estimates in detail. The singular solutions of linear elastostatics are derived, and a number of new models of corner and crack singularity are proposed. Effective numerical methods, such as the collocation Trefftz methods (CTM), the method of fundamental solutions (MFS), the method of particular solutions (MPS) and their combinations: the so called combined method, are developed. Such solutions are useful to examine other numerical methods for singularity problems in linear elastostatics. This thesis consists of three parts, Part I: Basic approaches, Part II: Advanced topics, and Part III: Mixed types of displacement and traction conditions. Contents of Parts I and II have been published in [47,82]. In Part I, the collocation Trefftz methods are used to obtain highly accurate solutions, where the leading coefficient has 14 (or 13) significant digits by the computation with double precision. In part II, two more new models (symmetric and anti-symmetric) of interior crack singularities are proposed, for the corner and crack singularity problems, the combined methods by using many fundamental solutions, but by adding a few singular solutions are proposed. Such a kind of combined methods is significant for linear elastostatics with corners (i.e., the L-shaped domain), because the singular solutions can only be obtained by seeking the power £hk of r£hk numerically. Hence, only a few singular solutions used may greatly simplify the numerical algorithms; Part III is a continued study of Parts I and II, to explore mixed type of displacement and free traction boundary conditions. To our best knowledge, this is the first time to provide the particular solutions near the corner with mixed types of boundary conditions and to report their numerical computation with different boundary conditions on the same corner edge in linear elastostatics. This thesis explores corner singularity and its numerical methods, to form a systematic study of basic theory and advanced computation for linear elastostatics.

Realistic and real-time computational simulation of biological organs (e.g., human kidneys, human liver) is a necessity when one tries to build a quality surgical simulator that can simulate surgical procedures involving these organs. Currently deformable models, spring-mass models, or finite element models are widely used to achieve the realistic simulations and/or the real-time performance. It is widely agreed that continuum mechanics based numerical techniques are preferred over deformable models or spring-mass models, but those techniques are computationally expensive and hence the higher accuracy offered by those numerical techniques come at the expense of speed. Hence there is a need to study the speed of different numerical techniques, while keeping an eye on the accuracy offered by those numerical techniques. Such studies are available for the Finite Element Method (FEM) but rarely available for the Boundary Element Method (BEM). Hence the present work aims to conduct a study on the viability of BEM for the real-time simulation of biological organs, and the present study is justified by the fact that BEM is considered to be inherently efficient when compared to mesh based techniques like FEM. A significant portion of literature on the real-time simulation of biological organs suggests the use of BEM to achieve better simulations. When one talks about the simulation of biological organs, one needs to have the geometry of a biological organ in hand. Geometry of biological organs of interest is not readily available many a times, and hence there is a need to extract the three dimensional (3D) geometry of biological organs from a stack of two dimensional (2D) scanned images. Software packages that can readily reconstruct 3D geometry of biological organs from 2D images are expensive. Hence, a novel procedure that requires only a few free software packages to obtain the geometry of biological organs from 2D image sequences is presented. The geometry of a pig liver is extracted from CT scan images for illustration purpose. Next, the three dimensional geometry of human kidney (left and right kidneys of male, and left and right kidneys of female) is obtained from the Visible Human Dataset (VHD). The novel procedure presented in this work can be used to obtain patient specific organ geometry from patient specific images, without requiring any of the many commercial software packages that can readily do the job. To carry out studies on the speed and accuracy of BEM, a source code for BEM is needed. Since the BEM code for 3D elasticity is not readily available, a BEM code that can solve 3D linear elastostatic problems without accounting for body forces is developed from scratch. The code comes in three varieties: a MATLAB version, a Fortran version (sequential version), and a Fortran version (parallelized version). This is the first free and open source BEM code for 3D elasticity. The developed code is used to carry out studies on the viability of BEM for the real-time simulation of biological organs, and a few representative problems involving kidneys and liver are found to give accurate solutions. The present work demonstrates that it is possible to simulate linear elastostatic behaviour in real-time using BEM without resorting to any type of precomputations, on a computer cluster by fully parallelizing the simulations and by performing simulations on different number of processors and for different block sizes. Since it is possible to get a complete solution in real-time, there is no need to separately prove that every type of cutting, suturing etc. can be simulated in real-time. Future work could involve incorporating nonlinearities into the simulations. Finally, a BEM based simulator may be built, after taking into account details like rendering.