Discrete differential geometry-based model for nonlinear analysis of axisymmetric shells

In this paper, we propose a novel one-dimensional (1D) discrete differential geometry (DDG)-based numerical method for geometrically nonlinear mechanics analysis (e.g., buckling and snapping) of axisymmetric shell structures. Our numerical model leverages differential geometry principles to accurately capture the complex nonlinear deformation patterns exhibited by axisymmetric shells. By discretizing the axisymmetric shell into interconnected 1D elements along the meridional direction, the in-plane stretching and out-of-bending potentials are formulated based on the geometric principles of 1D nodes and edges under the Kirchhoff-Love hypothesis, and elastic force vector and associated Hession matrix required by equations of motion are later derived based on symbolic calculation. Through extensive validation with available theoretical solutions and finite element method (FEM) simulations in literature, our model demonstrates high accuracy in predicting the nonlinear behavior of axisymmetric shells. Importantly, compared to the classical theoretical model and three-dimensional (3D) FEM simulation, our model is highly computationally efficient, making it suitable for large-scale real-time simulations of nonlinear problems of shell structures such as instability and snap-through phenomena. Moreover, our framework can easily incorporate complex loading conditions, e.g., boundary nonlinear contact and multi-physics actuation, which play an essential role in the use of engineering applications, such as soft robots and flexible devices. This study demonstrates that the simplicity and effectiveness of the 1D discrete differential geometry-based approach render it a powerful tool for engineers and researchers interested in nonlinear mechanics analysis of axisymmetric shells, with potential applications in various engineering fields.Comment: 36 pages, 11 figure

A numerical method for fluid-structure interactions of slender rods in turbulent flow

This thesis presents a numerical method for the simulation of fluid-structure interaction (FSI) problems on high-performance computers. The proposed method is specifically tailored to interactions between Newtonian fluids and a large number of slender viscoelastic structures, the latter being modeled as Cosserat rods. From a numerical point of view, such kind of FSI requires special techniques to reach numerical stability. When using a partitioned fluid-structure coupling approach this is usually achieved by an iterative procedure, which drastically increases the computational effort. In the present work, an alternative coupling approach is developed based on an immersed boundary method (IBM). It is unconditionally stable and exempt from any global iteration between the fluid part and the structure part. The proposed FSI solver is employed to simulate the flow over a dense layer of vegetation elements, usually designated as canopy flow. The abstracted canopy model used in the simulation consists of 800 strip-shaped blades, which is the largest canopy-resolving simulation of this type done so far. To gain a deeper understanding of the physics of aquatic canopy flows the simulation data obtained are analyzed, e.g., concerning the existence and shape of coherent structures

Development of A Kinetic Model For Loop-Free Colonoscopy Technology

The colonoscope is an important tool in diagnosis and management of diseases of the colon. One of the ongoing challenges with this device is that the colonoscope may form a loop together with the colon during the procedure. The result of the loop is that further insertion of the scope in the colon may not be possible. The loop may also cause risks of perforation of the colon and pain in the patient. There are currently several existing devices to overcome loop formation in colonoscopy, some of which have been introduced in clinical work. However, empirical assessment shows that these devices do not work very well. This is the motivation for the research presented in this thesis. In this thesis, a new paradigm of thinking, “doctor-assisted colonoscopy,” is proposed to overcome loop formation. In this new approach, the physician’s role is enhanced with new information that is acquired by sensors outside the human body and inferred from the mathematical model. It is referred to as a kinetic model due to the fact that this model describes the kinetic behaviour of the scope. This thesis is devoted to development of this kinetic model. In this study, the model of the colonoscope and the model of the colon are developed based on the Timoshenko beam theory, and parameters in both models are determined by the experiments. The following conclusions then are made: (1) self-locking of the colonoscope is the most basic cause for a loop to occur, while structural instability of the colonsocope is dependent on the self-locking; (2) both the scope and the colon can be well represented with the Timoshenko beam elements and the Linear Complementary Problem (LCP) formulation derived from Signorini’s law, and Coulom’s law for representation of interactions between the colon and scope is adequate; (3) there are effects from the location, looping, and tip deflection of the scope on flexural rigidity of the scope. Approximately, the flexural rigidity of the CF-Q160L colonoscope ranges from 300 to 650 N•cm2, and its accuracy is proven by a good agreement between the model predicted result and experimental result; (4) Rayleigh damping for the CF-Q160L colonoscope depends more on the mass matrix [M] of the colonoscope than the stiffness matrix [K], which is evident by the large coefficient value of “alpha” (0.3864) and the small coefficient value of “beta” (0.0164). The contributions of this thesis are: (1) the finding that the main cause of the loop is not structural instability of the colonoscope but rather self-locking of the colonoscope, which could lead to design of a “new-generation” colonoscope to avoid the loop; (2) a systematic evaluation of the existing colonoscopy technologies based on the well-proven Axiomatic Design Theory (ADT), which will serve as a guideline for the development of future new colonoscopes in future; (3) an approach to developing a kinetic model of the colonoscope useful to modeling similar objects such as a catheter guide-wire; (4) a novel ex-vivo colonoscopy test-bed with the kinetic and kinematic measurements useful for validation of new designs in colonoscopy technology and also useful for training physicians who perform the colonoscopy procedure; and (5) a new paradigm of thinking for colonoscopy called “doctor-assisted colonoscopy,” which has potential applications to other medical procedures such as catheter-based procedures

A novel musculoskeletal joint modelling for orthopaedic applications

This thesis was submitted for the degree of Docter of Philosophy and awarded by Brunel University.The objective of the work carried out in this thesis was to develop analytical and computational tools to model and investigate musculoskeletal human joints. It was recognised that the FEA was used by many researchers in modelling human musculoskeletal motion, loading and stresses. However the continuum mechanics played only a minor role in determining the articular joint motion, and its value was questionable. This is firstly due to the computational cost and secondly due to its impracticality for this application. On the other hand, there isn’t any suitable software for precise articular joint motion analysis to deal with the local joint stresses or non standard joints. The main requirement in orthopaedics field is to develop a modeller software (and its associated theories) to model anatomic joint as it is, without any simplification with respect to joint surface morphology and material properties of surrounding tissues. So that the proposed modeller can be used for evaluating and diagnosing different joint abnormalities but furthermore form the basis for performing implant insertion and analysis of the artificial joints. The work which is presented in this thesis is a new frame work and has been developed for human anatomic joint analysis which describes the joint in terms of its surface geometry and surrounding musculoskeletal tissues. In achieving such a framework several contributions were made to the 6DOF linear and nonlinear joint modelling, the mathematical definition of joint stiffness, tissue path finding and wrapping and the contact with collision analysis. In 6DOF linear joint modelling, the contribution is the development of joint stiffness and damping matrices. This modelling approach is suitable for the linear range of tissue stiffness and damping properties. This is the first of its kind and it gives a firm analytical basis for investigating joints with surrounding tissue and the cartilage. The 6DOF nonlinear joint modelling is a new scheme which is described for modelling the motion of multi bodies joined by non-linear stiffness and contact elements. The proposed method requires no matrix assembly for the stiffness and damping elements or mass elements. The novelty in the nonlinear modelling, relates to the overall algorithmic approach and handling local non-linearity by procedural means. The mathematical definition of joint stiffness is also a new proposal which is based on the mathematical definition of stiffness between two bodies. Based on the joint stiffness matrix properties, number of joint stiffness invariants was obtained analytically such as the centre of stiffness, the principal translational stiffnesses, and the principal rotational stiffnesses. In corresponding to these principal stiffnesses, their principal axes have been also obtained. Altogether, a joint is assessed by six principal axes and six principal stiffnesses and its centre of stiffness. These formulations are new and show that a joint can be described in terms of inherent stiffness properties. It is expected that these will be better in characterising a joint in comparison to laxity based characterisation. The development of tissue path finding and wrapping algorithms are also introduced as new approaches. The musculoskeletal tissue wrapping involves calculating the shortest distance between two points on a meshed surface. A new heuristic algorithm was proposed. The heuristic is based on minimising the accumulative divergence from the straight line between two points on the surface and the direction of travel on the surface (i.e. bone). In contact and collision based development, the novel algorithm has been proposed that detects possible colliding points on the motion trajectory by redefining the distance as a two dimensional measure along the velocity approach vector and perpendicular to this vector. The perpendicular distance determines if there are potentially colliding points, and the distance along the velocity determines how close they are. The closest pair among the potentially colliding points gives the “time to collision”. The algorithm can eliminate the “fly pass” situation where very close points may not collide because of the direction of their relative velocity. All these developed algorithms and modelling theories, have been encompassed in the developed prototype software in order to simulate the anatomic joint articulations through modelling formulations developed. The software platform provides a capability for analysing joints as 6DOF joints based on anatomic joint surfaces. The software is highly interactive and driven by well structured database, designed to be highly flexible for the future developments. Particularly, two case studies are carried out in this thesis in order to generate results relating to all the proposed elements of the study. The results obtained from the case studies show good agreement with previously published results or model based results obtained from Lifemod software, whenever comparison was possible. In some cases the comparison was not possible because there were no equivalent results; the results were supported by other indicators. The modelling based results were also supported by experiments performed in the Brunel Orthopaedic Research and Learning Centre

Eulerian on Lagrangian Cloth Simulation

This thesis introduces a novel Eulerian-on-Lagrangian (EoL) approach for simulating cloth. This approach allows for the simulation of traditionally difficult cloth scenarios, such as draping and sliding cloth over sharp features like the edge of a table. A traditional Lagrangian approach models a cloth as a series of connected nodes. These nodes are free to move in 3d space, but have difficulty with sliding over hard edges. The cloth cannot always bend smoothly around these edges, as motion can only occur at existing nodes. An EoL approach adds additional flexibility to a Lagrangian approach by constructing special Eulerian on Lagrangian nodes (EoL Nodes), where cloth material can pass through a fixed point. On contact with the edge of a box, EoL nodes are introduced directly on the edge. These nodes allow the cloth to bend exactly at the edge, and pass smoothly over the area while sliding. Using this ‘Eulerian-on-Lagrangian’ discretization, a set of rules for introducing and constraining EoL Nodes, and an adaptive remesher, This simulator allows cloth to move in a sliding motion over sharp edges. The current implementation is limited to cloth collision with static boxes, but the method presented can be expanded to include contact with more complicated meshes and dynamic rigid bodies

Macro-mechanical modelling and simulation of textile fabric and clothing with S-FEM

Tese de Doutoramento Programa Doutoral em Engenharia TêxtilEsta tese propõe um método de elementos finitos, designado por S-FEM (Smoothed Finite Element Method), para modelação e análise mecânica de estruturas têxteis planas. Neste enquadramento teórico, supõe-se que a estrutura têxtil não-tecida é um material isotrópico elástico, enquanto a estrutura têxtil tecida é um material elástico com anisotropia ortotrópica, para os quais as leis constitutivas utilizam propriedades mecânicas de baixa pressão (low stress) com base na Medição Objetiva de Tecidos (FOM - Fabric Objective Measurement). As formulações de elementos finitos de baixa ordem baseadas em deslocamento quando aplicadas a elementos finitos de placas (plate/shell) quadriláteras de 4 nós, incluindo campos de tensão de cisalhamento transversal, baseiam-se nas contribuições de Raymond Mindlin e por Eric Reissner, no que agora se designa teoria de deformação por cisalhamento de primeira ordem (first-order shear deformation, do inglês, ou FSDT de forma abreviada), ou simplesmente teoria de Mindlin-Reissner, e nas abordagens MITC (Mixed Interpolation of Tensorial Components), são nesta tese combinadas com a técnica de suavização do/da gradiente/tensão nos termos dos modelos S-FEM por forma a mitigar problemas como são o caso da distorção de elementos finitos, da granularidade grosseira da malha, bem como dos bem conhecidos fenómenos de bloqueio. As malhas de quadriláteros são utilizadas nesta tese devido à sua capacidade de representar geometrias complexas de tecidos em resultado de deformações mecânicas como são os casos da recuperação face à pressão planar, flexão, deformação, vibração, drapejamento, etc. Refira-se que foi desenvolvido e implementado em Matlab um software para os novos modelos de elementos finitos, em grande medida devido à inexistência de modelos S-FEM em softwares de análise de elementos finitos (finite element analysis ou FEA), lacuna esta que ocorre quer em softwares comerciais, quer não comerciais, e até em códigos abertos. Exemplos numéricos para as aplicações básicas de engenharia no que respeita à modelação mecânica de folhas de tecido fino e de folhas de tecido de espessura média em estudos de casos típicos, como é o caso da recuperação face a pressão planar, flexão, deformação e comportamento livre de vibrações, indicam que os elementos finitos (plate/shell) desenvolvidos com a técnica de suavização de tensão e MITC acabam por aliviar os efeitos de distorção dos elementos, a granularidade grosseira da malha e efeito de bloqueio na modelação e análise mecânica de tecidos muito finos e até mesmo de tecidos de espessura média. Os modelos de elementos finitos de placas (plate/shell) desenvolvidos durante o trajeto desta tese, bem como as suas propriedades mecânicas de baixa tensão em termos de FOM, são, portanto, bem adaptados à modelação e análise numérica de deformação macro-mecânica de folhas de tecido muito fino e de folhas de tecido de espessura média, incluindo ao mesmo tempo análise de deformação mecânica simples e complexa.An S-FEM (Smoothed Finite Element Method) for mechanical analysis and modelling of the textile fabrics is proposed. In this theoretical framework, one assumes that the non-woven fabric is an elastic isotropic material, while the woven fabric is an elastic with orthotropic anisotropy for which the constitutive laws formulated are using low-stress mechanical properties based on FOM (Fabric Objective Measurement). The displacement-based low-order finite element formulations for four-node quadrilateral plate/shell finite element, including assumed transverse shear strain fields, are based on the contributions of Raymond Mindlin and by Eric Reissner as FSDT (first-order shear deformation theory and so-called the Mindlin-Reissner theory) together with MITC (Mixed Interpolation of Tensorial Components) approaches, which are combined with the gradient/strain smoothing technique in terms of S-FEM models contributed by G. R. Liu et al. in order to mitigate problems as element distortion, mesh coarseness as well as the well-known locking phenomena. Quadrilateral meshes are used due to ability to represent complicated geometries of complex mechanical deformation of the fabric such as plane stress recovery, bending, buckling, vibration, draping behavior, etc. The finite element computer codes were developed in MATLAB for the new formulated plate/shell finite element models due to the lack of FEM (Finite Element Method) packages for S-FEM models in both commercial and non-commercial FEA (Finite Element Analysis) computer applications, and even from open-source platforms. Numerical examples for the basic engineering applications of mechanical modelling of thin to moderately thick fabric sheet in the typical case studies such as in-plane stress recovery, bending, buckling and free-vibration behavior, indicate that the developed plate/shell finite elements with assumed strain smoothing technique and MITC, do alleviate element distortion, mesh coarseness, and locking effect even for mechanical analysis and modelling very thin to moderately thick fabric. The developed plate/shell finite element models and low-stress mechanic properties in terms of FOM are, therefore, well adapted for numerical analysis and modelling of macro-mechanical deformation of the thin to moderately thick fabric sheet including both simple and complex mechanical deformation analysis.EMECW L12 MOBILITY GRANT AWARD CONTRACT BTG_559 Grant agreement n 2009/1661-001 001EC

Simulating Humans: Computer Graphics, Animation, and Control

People are all around us. They inhabit our home, workplace, entertainment, and environment. Their presence and actions are noted or ignored, enjoyed or disdained, analyzed or prescribed. The very ubiquitousness of other people in our lives poses a tantalizing challenge to the computational modeler: people are at once the most common object of interest and yet the most structurally complex. Their everyday movements are amazingly uid yet demanding to reproduce, with actions driven not just mechanically by muscles and bones but also cognitively by beliefs and intentions. Our motor systems manage to learn how to make us move without leaving us the burden or pleasure of knowing how we did it. Likewise we learn how to describe the actions and behaviors of others without consciously struggling with the processes of perception, recognition, and language

A survey on personal computer applications in industrial design process

Thesis (Master)--Izmir Institute of Technology, Industrial Design, Izmir, 1999Includes bibliographical references (leaves: 157-162)Text in English, Abstract: Turkish and Englishxii, 194 leavesIn this thesis, computer aided design systems are studied from the industrial designer's point of view. The study includes industrial design processes, computer aided design systems and the integration aspects.The technical issues are priorly studied, including current hardware and software technologies. The pure technical concepts are tried to be supported with real-world examples and graphics. Several important design software are examined, whether by personal practice or by literature research, depending on the availability of the software.Finally, the thesis include a case study, a 17" LCD computer monitor designed with a set of graphic programs including two-dimensional and three-dimensional packages.Keywords: Computers, industrial design methods, design software, computer aided design