非圆截面弹性细杆的平衡稳定性与分岔

本文研究存在初始曲率或挠率的非圆截面弹性细杆的平衡及稳定性问题,在两端受力矩单儿作用的条件下,杆的平衡微分方程可转换为用欧拉角表述的一阶自治系统,并有可能利用相平面的奇点理论分析弹性细杆平衡状态的稳定性,文中对杆截面的对称性,以及杆的初始曲率和挠率对平衡状态性的影响进行了定性分析,导出了解析形式的稳定性判据,揭示了杆平衡状态的列态分岔现象。     

Mechanical analysis of phase transition experiments of the bacterial flagellar filament

Bacterial flagellar filaments can undergo a polymorphic phase transition in both vitro and vivo environments. Each bacterial flagellar filament has 12 different helical forms which are macroscopically represented by different pitch lengths and helix radii. For external mechanical force induced filament phase transitions, there is so far only one experiment performed by Hotani in 1982, who showed a very beautiful cyclic phase transition phenomenon in his experiment on isolated flagellar filaments. In the present paper, we give a detailed mechanical analysis on Hotani＇s experiments. Through theoretical computations, we obtained a phase transition rule based on the phase transition mechanism. The theoretical analysis provides a foundation facilitating the establishment of phase transition theory for bacterial flagellar filaments.     

Non-linear elastic analysis of thin rods subjected to bending with arbitrary kinetic conditions of their cross-sections

In this paper problems concerning the non-linear analysis of thin rods due to pure bending with constant initial curvatures and twist and with arbitrary kinetic conditions of their cross-sections are presented. Couples are not considered as being applied to the rods except at their ends. The solutions developed in this paper, which determine the curvature components and the twist of the rod after deformation, are exact in the form of elliptic integrals.     

Simulation of bacterial flagellar phase transition by non-convex and non-local continuum modeling

Bacterial flagellar filament can undergo a stress-induced polymorphic phase transition in both vitro and vivo environments. The filament has 12 different helical forms (phases) characterized by different pitch lengths and helix radii. When subjected to the frictional force of flowing fluid, the filament changes between a left-handed normal phase and a right-handed semi-coiled phase via phase nucleation and growth. This paper develops non-local finite element method (FEM) to simulate the phase transition under a displacement-controlled loading condition (controlled helix-twist). The FEM formulation is based on the Ginzburg-Landau theory using a one-dimensional non-convex and non-local continuum model. To describe the processes of the phase nucleation and growth, viscosity-type kinetics is also used. The non-local FEM simulation captures the main features of the phase transition: two-phase coexistence with an interface of finite thickness, phase nucleation and phase growth with interface propagation. The non-local FEM model provides a tool to study the effects of the interfacial energy/thickness and loading conditions on the phase transition.     

非圆截面弹性细杆的螺旋线平衡及稳定性

本文研究端部受力和力矩作用，且存在初曲率和初扭率的非圆截面弹性细杆的螺旋线平衡及其稳定性。描述弹性细杆平衡状态的Kirchhoff方程存在与杆的螺旋线平衡状态相对应的特解。直杆和圆环杆为螺旋线状态的两种特例。文中分析了螺旋线的几何特性与作用力和力矩之间的相互关系，并导出螺旋线平衡的一次近似解析形式稳定性判据。分析表明，松弛状态下弹性杆可处于螺旋线状态，直杆只有在轴向压力的作用下才能保持螺旋线平衡。无初曲率和初扭率弹性杆的螺旋线平衡稳定性必要条件是杆截面绕副法线轴的抗弯刚度大于或等于绕法线轴的抗弯刚度。此条件也适用于带初扭率的圆环杆及更普遍情形。无初曲率和初扭率的圆截面杆的螺旋线平衡恒稳定。     

A Two-Strand Ply Hanging Under Its Own Weight

We consider an idealised model for a plied structure such as may form when a straight rod or filament is subjected to a high twisting moment. Examples are found in textile yarns, interwound DNA molecules and bacterial macrofibres. Plied structures, generally composed of more than two strands, are also used widely in engineering (mooring ropes, cables in lift shafts), although here the strands are often not intrinsically straight. The ply is assumed to consist of two strands of thin circular elastic rod winding around each other while touching on a straight line of contact. Each strand is therefore constrained to lie on a cylinder. Using a variational approach we give an unconstrained Hamiltonian formulation for this problem. We also derive an exact expression for the contact force acting between the two strands. We study the symmetry-breaking effect of gravity on the ply configuration as well as on the contact force.     

弹性杆基因模型的力学问题

概述弹性杆静力学与刚体动力学之间相似性的Kirchhoff理论．讨论其在分子生物学的弹性杆基因模型中的应用，以及与分析力学和运动稳定性理论有关的若干问题．     

Kirchhoff弹性杆分析动力学的准坐标表达

作为DNA的力学模型,依据Kirchhoff动力学比拟思想建立的弹性细杆的分析力学方法已从静力学深入到动力学。由于静力学平衡微分方程与刚体动力学相当,因此,弹性细杆动力学的分析力学方程必是以弧坐标和时间为双自变量的偏微分方程。以横截面的形心速度以及弯扭度和角速度沿主轴的分量为准速度,定义了准坐标,导出了准坐标的微分和变分运算的交换关系。从Hamilton原理出发,利用准坐标的微分和变分运算的交换关系,导出了Kirchhoff弹性杆动力学准坐标下的Boltzmann-Hamel方程,并由此导出Lanrange方程。指出了Boltzmann-Hamel方程显式即为弹性杆动力学的Kirchhoff方程。定义关于弧坐标和时间的正则变量和Hamilton函数,导出Boltzmann-Hamel方程的正则形式。本文结果是以弹性杆静力学和刚性杆动力学为其特例。作为例子,建立了垂挂的在重力作用下作平面运动的弹性细杆的动力学微分方程以说明本文方法的应用。     

Influence of transverse curvature on the stability of pre-stressed helical structures

This study is concerned with the stability characteristics of helix shaped structures made of anisotropic, pre-stressed, thin flanges arranged in such a way as to enable and develop multi-stability. Previous research on similar structures assumed the structural response of the flanges to be one-dimensional due to the narrow width of the pre-stressed members in comparison to their length. In this work, a refined two-dimensional model of the flanges is employed to model the influence of transverse curvature as well as the membrane strain energy associated with the non-zero Gaussian curvature deformations. While longitudinal curvature changes and twist are inherent to the geometry of the helices; the transverse curvature results from a consideration of boundary effects and the minimisation of the (expensive) membrane elastic energy. A qualitative study of the changes in transverse curvature reveals ways of simplifying the two-dimensional model into a simpler, closed form, one-dimensional version applicable to helices with relatively narrow flanges. Correlation is found between experimental results, finite element modelling and analytical predictions for the two models.     

Meshless methods for the inverse problem related to the determination of elastoplastic properties from the torsional experiment

The problem of determining the elastoplastic properties of a prismatic bar from the given experimental relation between the torsional moment M and the angle of twist per unit length of the rod’s length θ is investigated as an inverse problem. The proposed method to solve the inverse problem is based on the solution of some sequences of the direct problem by applying the Levenberg-Marquardt iteration method. In the direct problem, these properties are known, and the torsional moment is calculated as a function of the angle of twist from the solution of a non-linear boundary value problem. This non-linear problem results from the Saint-Venant displacement assumption, the Ramberg–Osgood constitutive equation, and the deformation theory of plasticity for the stress–strain relation. To solve the direct problem in each iteration step, the Kansa method is used for the circular cross section of the rod, or the method of fundamental solutions (MFS) and the method of particular solutions (MPS) are used for the prismatic cross section of the rod. The non-linear torsion problem in the plastic region is solved using the Picard iteration.     

Helical and Localised Buckling in Twisted Rods: A Unified Analysis of the Symmetric Case

We review the geometric rod theory for the case of a naturallystraight, linearly elastic, inextensible, circular rod suffering bendingand torsion but no shear. Our primary focus is on the post-bucklingbehaviour of such rods when subjected to end moment and tension.Although this is a classic problem with an extensive literature, datingback to Kirchhoff, the usual approach tends to neglect the physicalinterpretation of solutions (i.e., rod configurations) to the modelsproposed. Here, we explicitly compute geometrical properties of buckledrods. In a unified approach, making use of Kirchhoff's dynamic analogy,both the classical helical and the more recently investigated localisedbuckling are considered. Special attention is given to a consistenttreatment of concepts of link, twist and writhe.     

On the Bending and Twisting of Rods with Misfit

We derive a one-dimensional variational problem representing the elastic energy of a rod with misfit, starting from a nonlinear, three-dimensional elastic energy with nontrivial preferred strain. Our approach to dimension reduction is to find a Gamma-limit as the thickness of the rod tends to 0. The limiting energy is a quadratic function of the rates at which the rod bends and twists, and we give explicit expressions for the preferred curvature and twist in the special case of isotropic elastic moduli.     

Modeling of properties and motions of an inhomogeneous one-dimensional continuum of a complicated Cosserat-type microstructure

We consider an approach to modeling the properties of the one-dimensional Cosserat continuum [1] by using the mechanical modeling method proposed by Il’yushin in [2] and applied in [3]. In this method, elements (blocks, cells) of special form are used to develop a discrete model of the structure so that the average properties of the model reproduced the properties of the continuum under study. The rigged rod model, which is an elastic structure in the form of a thin rod with massive inclusions (pulleys) fixed by elastic hinges on its elastic line and connected by elastic belt transmissions, is taken to be the original discrete model of the Cosserat continuum. The complete system of equations describing the mechanical properties and the dynamical equilibrium of the rigged rod in arbitrary plane motions is derived. These equations are averaged in the case of a sufficiently smooth variation in the parameters of motion along the rod (the long-wave approximation). It was found that the average equations exactly coincide with the equations for the one-dimensional Cosserat medium [1] and, in some specific cases, with the classical equations of motion of an elastic rod [4–6]. We study the plane motions of the one-dimensional continuum model thus constructed. The equations characterizing the continuum properties and motions are linearized by using several assumptions that the kinematic parameters are small. We solve the problem of natural vibrations with homogeneous boundary conditions and establish that each value of the parameter distinguishing the natural vibration modes is associated with exactly two distinct vibration mode shapes (in the same mode), each of which has its own frequency value.     

A rod model for three dimensional deformations of single-walled carbon nanotubes

A continuum model for single-walled carbon nanotubes (SWCNT) is presented which is based on an extension to the special Cosserat theory of rods (Kumar and Mukherjee, 2011). The model allows deformation of a nanotube’s lateral surface in a one dimensional framework and hence is an efficient substitute to the commonly used two dimensional shell models for nanotubes. The model predicts a new coupling mode in chiral nanotubes – coupling between twist and cross-sectional shrinkage implying that the three deformation modes (extension, twist and cross-sectional shrinkage) are all coupled to each other. Atomistic simulations based on the density functional based tight binding method (DFTB) are performed on a (9, 6) SWCNT and the simulation data is used to estimate material parameters of this rod model. A peculiar behavior of the nanotube is observed when it is axially stretched – induced rotation of each cross-section is equal in magnitude but opposite to that of its two neighboring cross-sections. This is shown to be an effect of relative shift/inner-displacement between the two SWCNT sub-lattices.     

Non-linear dynamic intertwining of rods with self-contact

Twisted marine cables on the sea floor can form highly contorted three-dimensional loops that resemble tangles. Such tangles or ‘hockles’ are topologically equivalent to the plectomenes that form in supercoiled DNA molecules. The dynamic evolution of these intertwined loops is studied herein using a computational rod model that explicitly accounts for dynamic self-contact. Numerical solutions are presented for an illustrative example of a long rod subjected to increasing twist at one end. The solutions reveal the dynamic evolution of the rod from an initially straight state, through a buckled state in the approximate form of a helix, through the dynamic collapse of this helix into a near-planar loop with one site of self-contact, and the subsequent intertwining of this loop with multiple sites of self-contact. This evolution is controlled by the dynamic conversion of torsional strain energy to bending strain energy or, alternatively, by the dynamic conversion of twist (Tw) to writhe (Wr).     

DYNAMICAL EXPLANATION ON TOROFLUX'S PHENOMENA

本文分析了流动环特殊现象的力学原理。流动环是由细长钢带缠绕形成的环状螺旋玩具,具有瞬间几何形态突变和绕圆柱物体运动时迅速旋转两种独特现象。文中利用封闭曲杆的拓扑学规律确定其曲率和扭率,用于计算和比较流动环两种平衡状态的弹性变形势能。证实流动环的形态突变源于最小弹性势能原理。以组成流动环的单个圆环为分析对象,且考虑邻近圆环的牵拉效应。分析了流动环下落时相对圆柱体的滚动过程,以解释旋转现象的产生原因。导出了流动环下落旋转的角加速度公式。理论计算结果与实验数据较好吻合。     

Mechanics and chemical thermodynamics of phase transition in temperature-sensitive hydrogels

This paper uses the thermodynamic data of aqueous solutions of uncrosslinked poly(N-isopropylacrylamide) (PNIPAM) to study the phase transition of PNIPAM hydrogels. At a low temperature, uncrosslinked PNIPAM can be dissolved in water and form a homogenous liquid solution. When the temperature is increased, the solution separates into two liquid phases with different concentrations of the polymer. Covalently crosslinked PNIPAM, however, does not dissolve in water, but can imbibe water and form a hydrogel. When the temperature is changed, the hydrogel undergoes a phase transition: the amount of water in the hydrogel in equilibrium changes with temperature discontinuously. While the aqueous solution is a liquid and cannot sustain any nonhydrostatic stress in equilibrium, the hydrogel is a solid and can sustain nonhydrostatic stress in equilibrium. The nonhydrostatic stress can markedly affect various aspects of the phase transition in the hydrogel. We adopt the Flory-Rehner model, and show that the interaction parameter as a function of temperature and concentration obtained from the PNIPAM-water solution can be used to analyze diverse phenomena associated with the phase transition of the PNIPAM hydrogel. We analyze free swelling, uniaxially and biaxially constrained swelling of a hydrogel, swelling of a core-shell structure, and coexistent phases in a rod. The analysis is related to available experimental observations. Also outlined is a general theory of coexistent phases undergoing inhomogeneous deformation.     

Using filament stretching rheometry to predict strand formation and “processability” in adhesives and other non-Newtonian fluids

The spinning of polymeric fibers, the processing of numerous foodstuffs and the peel and tack characteristics of adhesives are all associated with the formation, stability and, ultimately, the longevity of thin fluid `strands'. This tendency to form strands is usually described in terms of the tackiness of the fluid or by heuristic concepts such as `stringiness' (Lakrout et al. J Adhesion 1999). The dynamics of such processes are complicated due to spatially and temporally non-homogeneous growth of extensional stresses, the action of capillary forces and the evaporation of volatile solvents. We describe the development and application of a simple instrument referred to as a microfilament rheometer (MFR) that can be used to readily differentiate between the dynamical response of different pressure-sensitive adhesive fluid formulations. The device relies on a quantitative observation of the rate of extensional thinning or `necking' of a thin viscous or viscoelastic fluid filament in which the solvent is free to evaporate across the free surface. This high-resolution measurement of the radial profile provides a direct indication of the ultimate time to break up of the fluid filament. This critical time is a sensitive function of the rheological properties of the fluid and the mass transfer characteristics of the solvent, and can be conveniently reported in terms of a new dimensionless quantity we refer to as a processability parameter P. We demonstrate the usefulness of this technique by presenting our results in the form of a case study in which we measure the visco-elasto-capillary thinning of slender liquid filaments for a number of different commercial polymer/solvent formulations and relate this to the reported processing performance of the materials. We also compare the MFR observations with the prediction of a simple 1D theory derived from the governing equations that model the capillary thinning of an adhesive filament. Received: 22 December 1999/Accepted: 4 January 2000     

Analysis of coupled thermo-hydro-mechanical behavior of unsaturated soils based on theory of mixtures I

This paper analyzes a coupled thermo-hydro-mechanical behavior of unsaturated soils based on the theory of mixtures. Unsaturated soil is considered as a mixture composed of soil skeleton, liquid water, vapor, dry air, and dissolved air. In addition to the mass and momentum conservation equations of each component and the energy conservation equation of the mixture, the system is closed using other 37 constitutive （or restriction） equations. As the change in water chemical potential is identical to the change in vapor chemical potential, a thermodynamic restriction relationship for the phase transition between pore water and pore vapor is formulated, in which the impact of the change in gas pressure on the phase transition is taken into account. Six final govern- ing equations are given in incremental form in terms of six primary variables, i.e., three displacement components of soil skeleton, water pressure, gas pressure, and temperature. The processes involved in the coupled model include thermal expansions of soil skeleton and soil particle, Soret effect, phase transition between water and vapor, air dissolution in pore water, and deformation of soil skeleton.     

螺旋细杆的均匀扭转振动

讨论螺旋细杆的特殊形式扭转振动，即均匀扭转振动．以非圆截面杆和有原始曲率的圆截面杆为研究对象．杆作均匀扭转振动时各截面有相同的扭角变化规律，且杆中心线的几何形状不受振动过程的影响．研究表明，扭振来源于杆截面的非对称性及杆的原始曲率．杆的扭振规律与单摆运动相似，其动力学方程存在精确解．圆环杆的均匀扭振为螺旋杆的倾角为零时的特例．