Construction of periodic adapted orthonormal frames on closed space curves

The construction of continuous adapted orthonormal frames along C1 closed–loop spatial curves is addressed. Such frames are important in the design of periodic spatial rigid–body motions along smooth closed paths. The construction is illustrated through the simplest non–trivial context — namely, C1 closed loops defined by a single Pythagorean–hodograph (PH) quintic space curve of a prescribed total arc length. It is shown that such curves comprise a two–parameter family, dependent on two angular variables, and they degenerate to planar curves when these parameters differ by an integer multiple of π. The desired frame is constructed through a rotation applied to the normal–plane vectors of the Euler–Rodrigues frame, so as to interpolate a given initial/final frame orientation. A general solution for periodic adapted frames of minimal twist on C1 closed–loop PH curves is possible, although this incurs transcendental terms. However, the C1 closed–loop PH quintics admit particularly simple rational periodic adapted frames

Characterization of Spherical and Plane Curves Using Rotation Minimizing Frames

In this work, we study plane and spherical curves in Euclidean and Lorentz-Minkowski 3-spaces by employing rotation minimizing (RM) frames. By conveniently writing the curvature and torsion for a curve on a sphere, we show how to find the angle between the principal normal and an RM vector field for spherical curves. Later, we characterize plane and spherical curves as curves whose position vector lies, up to a translation, on a moving plane spanned by their unit tangent and an RM vector field. Finally, as an application, we characterize Bertrand curves as curves whose so-called natural mates are spherical.Comment: 8 pages. This version is an improvement of the previous one. In addition to a study of some properties of plane and spherical curves, it contains a characterization of Bertrand curves in terms of the so-called natural mate

Mapping rational rotation-minimizing frames from polynomial curves on to rational curves

Given a polynomial space curve r(ξ) that has a rational rotation–minimizing frame (an RRMF curve), a methodology is developed to construct families of rational space curves r˜(ξ) with the same rotation–minimizing frame as r(ξ) at corresponding points. The construction employs the dual form of a rational space curve, interpreted as the edge of regression of the envelope of a family of osculating planes, having normals in the direction u(ξ)=r′(ξ)×r″(ξ) and distances from the origin specified in terms of a rational function f(ξ) as f(ξ)/‖u(ξ)‖. An explicit characterization of the rational curves r˜(ξ) generated by a given RRMF curve r(ξ) in this manner is developed, and the problem of matching initial and final points and frames is shown to impose only linear conditions on the coefficients of f(ξ), obviating the non–linear equations (and existence questions) that arise in addressing this problem with the RRMF curve r(ξ). Criteria for identifying low–degree instances of the curves r˜(ξ) are identified, by a cancellation of factors common to their numerators and denominators, and the methodology is illustrated by a number of computed examples

기하학적으로 정밀한 비선형 구조물의 아이소-지오메트릭 형상 설계 민감도 해석

학위논문 (박사)-- 서울대학교 대학원 : 공과대학 조선해양공학과, 2019. 2. 조선호.In this thesis, a continuum-based analytical adjoint configuration design sensitivity analysis (DSA) method is developed for gradient-based optimal design of curved built-up structures undergoing finite deformations. First, we investigate basic invariance property of linearized strain measures of a planar Timoshenko beam model which is combined with the selective reduced integration and B-bar projection method to alleviate shear and membrane locking. For a nonlinear structural analysis, geometrically exact beam and shell structural models are basically employed. A planar Kirchhoff beam problem is solved using the rotation-free discretization capability of isogeometric analysis (IGA) due to higher order continuity of NURBS basis function whose superior per-DOF(degree-of-freedom) accuracy over the conventional finite element analysis using Hermite basis function is verified. Various inter-patch continuity conditions including rotation continuity are enforced using Lagrage multiplier and penalty methods. This formulation is combined with a phenomenological constitutive model of shape memory polymer (SMP), and shape programming and recovery processes of SMP structures are simulated. Furthermore, for shear-deformable structures, a multiplicative update of finite rotations by an exponential map of a skew-symmetric matrix is employed. A procedure of explicit parameterization of local orthonormal frames in a spatial curve is presented using the smallest rotation method within the IGA framework. In the configuration DSA, the material derivative is applied to a variational equation, and an orientation design variation of curved structure is identified as a change of embedded local orthonormal frames. In a shell model, we use a regularized variational equation with a drilling rotational DOF. The material derivative of the orthogonal transformation matrix can be evaluated at final equilibrium configuration, which enables to compute design sensitivity using the tangent stiffness at the equilibrium without further iterations. A design optimization method for a constrained structure in a curved domain is also developed, which focuses on a lattice structure design on a specified surface. We define a lattice structure and its design variables on a rectangular plane, and utilize a concept of free-form deformation and a global curve interpolation to obtain an analytical expression for the control net of the structure on curved surface. The material derivative of the analytical expression eventually leads to precise design velocity field. Using this method, the number of design variables is reduced and design parameterization becomes more straightforward. In demonstrative examples, we verify the developed analytical adjoint DSA method in beam and shell structural problems undergoing finite deformations with various kinematic and force boundary conditions. The method is also applied to practical optimal design problems of curved built-up structures. For example, we extremize auxeticity of lattice structures, and experimentally verify nearly constant negative Poisson's ratio during large tensile and compressive deformations by using the 3-D printing and optical deformation measurement technologies. Also, we architect phononic band gap structures having significantly large band gap for mitigating noise in low audible frequency ranges.본 연구에서는 대변형을 고려한 휘어진 조립 구조물의 연속체 기반 해석적 애조인 형상 설계 민감도 해석 기법을 개발하였다. 평면 Timoshenko 빔의 선형화된 변형률의 invariance 특성을 고찰하였고 invariant 정식화를 선택적 축소적분(selective reduced integration) 기법 및 B-bar projection 기법과 결합하여 shear 및 membrane 잠김 현상을 해소하였다. 비선형 구조 모델로서 기하학적으로 정밀한 빔 및 쉘 모델을 활용하였다. 평면 Kirchhoff 빔 모델을 NURBS 기저함수의 고차 연속성에 따른 아이소-지오메트릭 해석 기반 rotation-free 이산화를 활용하여 다루었으며, 기존의 Hermite 기저함수 기반의 유한요소법에 비해 자유도당 해의 정확도가 높음을 검증하였다. 라그랑지 승수법 및 벌칙 기법을 도입하여 회전의 연속성을 포함한 다양한 다중패치간 연속 조건을 고려하였다. 이러한 기법을 현상학적 (phenomenological) 형상기억폴리머 (SMP) 재료 구성방정식과 결합하여 형상의 프로그래밍과 회복 과정을 시뮬레이션하였다. 전단변형을 겪는 (shear-deformable) 구조 모델에 대하여 대회전의 갱신을 교대 행렬의 exponential map에 의한 곱의 형태로 수행하였다. 공간상의 곡선 모델에서 최소회전 (smallest rotation) 기법을 통해 국소 정규직교좌표계의 명시적 매개화를 수행하였다. 형상 설계 민감도 해석을 위하여 전미분을 변분 방정식에 적용하였으며 휘어진 구조물의 배향 설계 변화는 국소 정규직교좌표계의 회전에 의하여 기술된다. 최종 변형 형상에서 직교 변환 행렬의 전미분을 계산함으로써 대회전 문제에서 추가적인 반복 계산없이 변형 해석에서의 접선강성행렬에 의해 해석적 설계 민감도를 계산할 수 있다. 쉘 구조물의 경우 면내 회전 자유도 및 안정화된 변분 방정식을 활용하여 보강재(stiffener)의 모델링을 용이하게 하였다. 또한 본 연구에서는 휘어진 영역에 구속되어있는 구조물에 대한 설계 속도장 계산 및 최적 설계기법을 제안하며 특히 곡면에 구속된 빔 구조물의 설계를 집중적으로 다룬다. 자유형상변형(Free-form deformation)기법과 전역 곡선 보간기법을 활용하여 직사각 평면에서 형상 및 설계 변수를 정의하고 곡면상의 곡선 형상을 나타내는 조정점 위치를 해석적으로 표현할 수 있으며 이의 전미분을 통해 정확한 설계속도장을 계산한다. 이를 통해 설계 변수의 개수를 줄일 수 있고 설계의 매개화가 간편해진다. 개발된 방법론은 다양한 하중 및 운동학적 경계조건을 갖는 빔과 쉘의 대변형 문제를 통해 검증되며 여러가지 휘어진 조립 구조물의 최적 설계에 적용된다. 대표적으로, 전단 강성 및 충격 흡수 특성과 같은 기계적 물성치의 개선을 위해 활용되는 오그제틱 (auxetic) 특성이 극대화된 격자 구조를 설계하며 인장 및 압축 대변형 모두에서 일정한 음의 포아송비를 나타냄을 3차원 프린팅과 광학적 변형 측정 기술을 이용하여 실험적으로 검증한다. 또한 우리는 소음의 저감을 위해 활용되는 가청 저주파수 영역대에서의 밴드갭이 극대화된 격자 구조를 제시한다.Abstract 1. Introduction 2. Isogeometric analysis of geometrically exact nonlinear structures 3. Isogeometric confinguration DSA of geometrically exact nonlinear structures 4. Numerical examples 5. Conclusions and future works A. Supplements to the geometrically exact Kirchhoff beam model B. Supplements to the geometrically exact shear-deformable beam model C. Supplements to the geometrically exact shear-deformable shell model D. Supplements to the invariant formulations E. Supplements to the geometric constraints in design optimization F. Supplements to the design of auxetic structures 초록Docto

Patient-specific technology for in vivo assessment of 3-D spinal motion

One of the most common musculoskeletal problems affecting people is neck and low back pain. Traditional clinical diagnostic techniques such as fluoroscopic imaging or CT scans are limited due to their static and/or planar measurements which may not be able to capture all neurological pathologies. More advanced diagnostics have proven successful in assessing 3-D patient-specific spinal kinematics by combining a patient-specific 3-D spine model (CT or MRI) with bi-planar fluoroscopic imaging; however, custom, not clinically available advanced imaging equipment as well as an increase in radiation exposure is required to acquire a complete patient-specific spinal kinematic description. Hence, the purpose of this research was to develop a clinically viable bi-planar fluoroscopic imaging technique which acquires a complete patient-specific kinematic description of the spine with reduced radiation exposure. Development of the proposed technique required evaluating the accuracy of 3-D kinematic interpolation techniques in reconstructing spinal kinematic data in order to reduce radiation exposure from bi-planar fluoroscopic diagnostic techniques. Several interpolation and sampling algorithms were evaluated in reconstructing cadaveric lumbar (L2-S1) flexion-extension motion data; ultimately, a new interpolation algorithm was proposed. Similarly, the success of the interpolation algorithm was evaluated in reconstructing spine-specific kinematic parameters. Next, the interpolation algorithm was combined with a CT-based bi-planar fluoroscopic method. Accuracy of the proposed diagnostic technique was evaluated against previously validated work on an ex vivo optoelectronic 3-D kinematic assessment technique. Bi-planar fluoroscopic images were acquired during both flexion-extension and lateral bending motions of cadaveric cervical (C4-T1) and lumbar (L2-S1) spine. Registration of the bi-planar fluoroscopic images to the CT-based 3-D model was optimized using a gradient derived similarity function. Additionally, a stochastic approach, covariance matrix adaptive evolution strategy, was used as the optimizing function. The newly developed interpolation algorithm was used to reduce the sample size of the bi-planar fluoroscopic images which reduces radiation exposure. Experimental results illustrate the potential success of the technique, but ultimately improvements in registration and validation methods are needed before becoming clinically viable

Numerical Methods in Shape Spaces and Optimal Branching Patterns

The contribution of this thesis is twofold. The main part deals with numerical methods in the context of shape space analysis, where the shape space at hand is considered as a Riemannian manifold. In detail, we apply and extend the time-discrete geodesic calculus (established by Rumpf and Wirth [WBRS11, RW15]) to the space of discrete shells, i.e. triangular meshes with fixed connectivity. The essential building block is a variational time-discretization of geodesic curves, which is based on a local approximation of the squared Riemannian distance on the manifold. On physical shape spaces this approximation can be derived e.g. from a dissimilarity measure. The dissimilarity measure between two shell surfaces can naturally be defined as an elastic deformation energy capturing both membrane and bending distortions. Combined with a non-conforming discretization of a physically sound thin shell model the time-discrete geodesic calculus applied to the space of discrete shells is shown to be suitable to solve important problems in computer graphics and animation. To extend the existing calculus, we introduce a generalized spline functional based on the covariant derivative along a curve in shape space whose minimizers can be considered as Riemannian splines. We establish a corresponding time-discrete functional that fits perfectly into the framework of Rumpf and Wirth, and prove this discretization to be consistent. Several numerical simulations reveal that the optimization of the spline functional—subject to appropriate constraints—can be used to solve the multiple interpolation problem in shape space, e.g. to realize keyframe animation. Based on the spline functional, we further develop a simple regression model which generalizes linear regression to nonlinear shape spaces. Numerical examples based on real data from anatomy and botany show the capability of the model. Finally, we apply the statistical analysis of elastic shape spaces presented by Rumpf and Wirth [RW09, RW11] to the space of discrete shells. To this end, we compute a Fréchet mean within a class of shapes bearing highly nonlinear variations and perform a principal component analysis with respect to the metric induced by the Hessian of an elastic shell energy. The last part of this thesis deals with the optimization of microstructures arising e.g. at austenite-martensite interfaces in shape memory alloys. For a corresponding scalar problem, Kohn and Müller [KM92, KM94] proved existence of a minimizer and provided a lower and an upper bound for the optimal energy. To establish the upper bound, they studied a particular branching pattern generated by mixing two different martensite phases. We perform a finite element simulation based on subdivision surfaces that suggests a topologically different class of branching patterns to represent an optimal microstructure. Based on these observations we derive a novel, low dimensional family of patterns and show—numerically and analytically—that our new branching pattern results in a significantly better upper energy bound