On the mechanics of fishscale structures

Biological and manmade structures often share the same specifications and design constraints: structural support, lightweight or protection against specific threats. In this context, the structure of fishscales, consisting of small rigid plates growing out of the skin of a majority of fish species, are characterized by a large variety of shape, size and properties in order to achieve particular functions. The present study introduces a basic two-dimensional micromechanical model that permits to establish a correlation between the flexural response of a scaled skin and the nature of its underlying structure, including both geometric and material aspects. The model is used to predict trends in the structure’s response and illustrates the fact that the scale design, arrangement and properties can be tailored to achieve a wide spectrum of response. In particular, fishscale structure possesses an inherent strain-stiffening response that can be suppressed or magnified by certain structural features. This particularity, shared by most biological materials, ensures that the structure provides both a structural and protective support for the animal.     

水凝胶多场耦合计算力学

作为一种具有多场耦合特性的智能柔体材料,水凝胶的制备技术、性能表征与结构应用得到迅速发展。本文在分析水凝胶本构理论和结构设计的基础上,提出了水凝胶多场耦合计算力学的基本方法和范式,包括微观粗粒化分子动力学模拟和宏观耦合有限元方法等,计算了化学-力学耦合作用下水凝胶材料与结构的变形和应力,给出了多个数值算例与结果比较。研究指出多场耦合计算力学将成为水凝胶材料和结构分析的主要手段,并推动水凝胶等这类智柔材料的性能设计与工程应用。     

非均质材料动力分析的广义多尺度有限元法

自然界和工程中的大部分材料都具有多尺度特征,当考察尺度小到一定程度后,都将表现出非均质性.针对非均质材料的动力问题,提出了一种广义多尺度有限元方法,其基本思想是利用静态凝聚法以及罚函数法构造能够反映单元内部材料非均质特性的多尺度位移基函数.与传统扩展多尺度有限元法中的基函数构造方式不同,广义多尺度有限元法的基函数无需通过在子网格域上多次求解椭圆问题得到,而可直接通过矩阵运算获得.其主要步骤如下:利用数值基函数将一个非均质单胞等效为一个宏观单元,进而形成整个结构的等效刚度矩阵,并得到宏观网格的节点位移,最后再次利用数值基函数得到微观尺度上的位移结果.该广义多尺度有限元法是扩展多尺度有限元法的一种新的拓展,可模拟具有更加复杂几何的非均质单胞的力学行为.通过数值算例,模拟了非均质材料的静力问题、广义特征值问题以及瞬态响应问题,计算结果表明:在边界条件一样的情况下,广义多尺度有限元法的计算结果与传统有限元的计算结果保持高度一致.与传统有限元相比,该方法在保证计算精度的同时极大地提高了计算效率.研究结果表明,广义多尺度有限元法能够很好地模拟非均质单胞的力学行为,具有良好的工程应用潜力.      

零泊松比手风琴蜂窝等效模量

柔性蒙皮是变形机翼和风力机叶片的关键组成部分。一维变形的柔性蒙皮不仅要求其支撑结构具有良好的面内变形和面外承载能力,还需要具有零泊松比特性。手风琴蜂窝具有零泊松比特性,可用作一维变形柔性蒙皮支撑。为全面分析其面内外弹性变形特性,综合考虑结构的内力弯矩、轴力和剪力,通过卡氏第二定理对其x向等效弹性模量和x-y面内等效剪切模量进行了推导;利用最小余能原理和最小势能原理确定了x-z面内的等效剪切模量;此外还推导了其y和z向的等效弹性模量以及y-z面内的等效剪切模量;然后通过ANSYS有限元仿真对等效模量理论公式进行了验证;最后将本文理论模型与现有模型进行了比较。结果表明,理论公式和有限元仿真吻合较好,在结构设计时采用较大的斜梁高度系数h和斜梁间距系数g,较小的厚度系数t以及较大的竖直梁厚度系数η,有望获得具有较小面内刚度和较大面外刚度的手风琴蜂窝结构。该结果可用于手风琴蜂窝面内外等效模量的快速预测,为一维变形柔性蒙皮的结构设计提供相应的参考。此外,本文理论模型相比传统模型更为精确且具有更加广泛的应用范围。     

Linear wave-induced dynamic structural stress analysis of a 2D semi-flexible closed fish cage

Closed fish cages in the sea are proposed as a new concept in marine aquaculture, replacing the conventional net cages in order to meet ecological challenges related to fish lice and escapes. A closed fish cage can be compared to a floating tank structure with an internal free surface. Several types of closed cages have been suggested, and they are categorised according to structural properties as flexible membrane structures (fabric), semi-flexible structures (glass fibre) and rigid structures (steel or concrete). To be able to develop safe and reliable structures, more knowledge is required on the seakeeping behaviour of closed cages in waves and the structural response to the wave loads. This paper builds on a theory presented in Strand and Faltinsen (2019) on the linear wave loads on a 2D closed flexible fish cage. A modelling error has been found in Strand and Faltinsen (2019), however, all the main conclusions are in hold. The error has been corrected in the model in the present paper. The present paper extends the model to include bending in the structural model to be able to handle semi-flexible structures where bending stiffness is significant. In this paper, the linear theory of a 2D semi-flexible closed fish cage in waves is developed and analysed to investigate the structural response of the semi-flexible closed cage in waves. We have compared a quasi-static analysis with a fully coupled hydroelastic analysis to investigate if it is a valid and conservative assumption to assume that the stresses in the structure can be assumed quasi-static. If a hydroelastic analysis is necessary or not, is dependent on the stiffness of the structure. We have investigated what happens with the stress in the curved beam part of the closed fish cage for increasing and decreasing stiffness relative to a reference composite structure. One stiffer and two softer cases have been analysed. One major concern for the structural stresses in a closed cage is the effect of sloshing. Sloshing is internal wave motion inside the cage and have multiple resonance periods. The results indicate that to use the quasi-static assumption in structural stress calculation is conservative within the given frequency range for all examined stiffnesses and frequencies, except the frequencies very close to the second sloshing frequency. Close to the second sloshing frequency for all the examined stiffnesses, a localised peak can be observed in the coupled hydroelastic results. The second sloshing frequency is a frequency connected to a symmetric sloshing mode. Rigid body motion is not affected at the symmetric sloshing frequency for an assumed rigid structure, and are therefore also not visible in the stress results from the quasi-static analysis. The structural stress in irregular sea was calculated. These results show no indication of increased stress close to the second sloshing frequency. However, this is not a surprising result since the stress peak is very localised in frequency, and the accumulated effect on the stress standard deviation is therefore small.     

Stability of hard plates on soft substrates and application to the design of bioinspired segmented armor

Flexible natural armors from fish, alligators or armadillo are attracting an increasing amount of attention from their unique and attractive combinations of hardness, flexibility and light weight. In particular, the extreme contrast of stiffness between hard plates and surrounding soft tissues give rise to unusual and attractive mechanisms, which now serve as model for the design of bio-inspired armors. Despite a growing interest in bio-inspired flexible protection, there is little guidelines as to the choice of materials, optimum thickness, size, shape and arrangement for the protective plates. In this work, we focus on a failure mode we recently observed on natural and bio-inspired scaled armors: the unstable tilting of individual scales subjected to off-centered point forces. We first present a series of experiments on this system, followed by a model based on contact mechanics and friction. We condense the result into a single stability diagram which capture the key parameters that govern the onset of plate tilting from a localized force. We found that the stability of individual plates is governed by the location of the point force on the plate, by the friction at the surface of the plate, by the size of the plate and by the stiffness of the substrate. We finally discuss how some of these parameters can be optimized at the design stage to produce bio-inspired protective systems with desired combination of surface hardness, stability and flexural compliance.     

重型车驾驶室结构特性仿真分析与研究

针对某重型车驾驶室结构特性分析与优化的问题进行研究,采用CATIA软件对驾驶室进行几何建模,利用HYPERMESH软件建立有限元模型,并进行自由模态仿真分析;利用试验模态和计算模态分析对比的方法,验证了有限元模型的有效性;结合此款驾驶室的承载特点,利用MSC.Patran/Nastran对驾驶室的弯曲和扭转工况进行刚度和强度仿真计算,利用Fatigue对扭转工况下驾驶室的疲劳强度进行仿真计算,给出了驾驶室的易疲劳破坏位置。研究结果表明:驾驶室固有频率为37.826Hz,避开了激振频率,避免了共振现象;在扭转、弯曲两种工况下,驾驶室产生的最大应力分别为108MPa、38.8MPa,均小于材料许用应力,弯曲、扭转刚度满足要求,驾驶室抵抗变形的能力较好,结构设计合理。     

关于极端力学

随着前沿科学和新技术不断发展,工程材料与结构的超常规尺度、密度、硬度、刚度等性能以及在超常规温度、速度、场强和恶劣天气等极端服役环境中的力学响应规律,需要力学提供更为有效的理论和方法. 本报告从极端力学的基本定义和科学内涵出发,结合重大工程问题和大科学问题,从极端性能、极端载荷、学科发展等三个方面系统介绍了极端力学的研究现状,并总结了极端力学的特点及其对力学理论、计算方法和实验技术的挑战,最后对极端力学未来的发展进行了展望.      

高强钢波纹夹芯结构的力学性能研究

本文利用ABAQUS仿真软件分析了高强钢波纹夹芯结构在三点弯曲以及面内压缩时的力学响应。将计算结果与实验值和理论值进行了对比,验证了数值模拟的准确性。分析了芯板材料属性（高强钢DP900和普通低碳钢DC01）以及波纹芯板排列方向(横向、纵向)对夹芯结构力学性能的影响。结果表明DP900高强钢波纹芯板夹芯结构抗弯强度是DC01低碳钢夹芯结构的2.39倍,抗压强度是低碳钢夹芯结构的1.40倍,纵向波纹夹芯结构比横向波纹夹芯结构弯曲刚度高11.63倍。     

ANCF索梁单元应变耦合问题与模型解耦

索粱结构在土木工程、航空航天等领域有着广泛的应用.在各类索梁动力学建模方法中,由于绝对节点坐标方法(absolute nodal coordinate formulation,ANCF)能够描述柔性体的大变形和大转动问题,因此非常适合大变形索梁结构的动力学建模.对绝对节点坐标索梁单元的应变进行分析可知,弯曲变形会引起单元内部轴向应变的不均匀分布,即单元轴向应变与弯曲应变相互耦合.这种应变耦合效应使单元产生伪应变能,导致单元刚度增大,造成单元失真.分析不同弯曲角下的单元应变及应变能可知,弯曲变形越大,单元失真越严重.通过构造等效一维杆单元重新描述轴向应变,实现了轴向应变与弯曲应变解耦.在此基础上推导广义弹性力,得到了绝对节点坐标索梁单元的应变解耦模型.对解耦前后的两种梁模型进行静力学和动力学仿真,结果表明;解耦模型消除了单元伪应变,相比原模型表现出更好的收敛性和曲率连续性,在相同单元数目下具有更高的精度.同时由于解耦模型降低了单元刚度,因此相比原模型,速度曲线中不再有高频振动.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

广义协调元在薄壳自由振动分析中的应用

本文采用２０个自由度的矩形平板型壳单元分析薄筒壳自由振动问题。在文献（１）的基础上，将广义协调元的应用范围扩大到了壳体分析中，目的是找出简单实用的方法分析板壳这类具有特殊性质的结构。该单元由平面应力单元和平板弯曲单元经简单叠加而成，是一种性能良好的单元。沿用常规作法非常便利，程序容易实现。经算例验证，该单元自由度少，精度较好，实用方便，适合于工程应用。     

剪力墙的面外刚度对建筑结构计算结果的影响

剪力在钢筋混凝土建筑结构中有着广泛的应用,在计算时一般都忽略剪力面外刚度的贡献。本文利用一种同时具有面内刚度和面外刚度的剪力墙模型,研究了剪力墙的面外刚度对整体结构的控制位移和固有频率的影响。结果表明,剪力面外刚度的贡献与剪力墙之间的耦合程度以及是要用刚性楼板假定有关,在大型复杂的纵横剪力墙结构中,剪力墙面外刚度的贡献不可忽略。     

双曲率壳的力学研究进展

双曲率壳是飞机、汽车以及船舶上常见的薄壁结构，其中性面可看作是一条动曲线沿着另一条曲线扫掠所形成的曲面．双曲率壳的非线性理论不断更新推动着双曲率壳力学行为的研究．随着工程实际应用的不断改进，如功能梯度材料(FGM)，加筋壳，弹性地基模型等的引入，双曲率壳在强度、变形和稳定性等方面的研究得到了进一步促进．本文首先回顾了双曲率壳结构非线性力学基本理论发展过程，主要阐述了经典的二维板壳理论，如Donnell 薄壳理论，一阶剪切变形壳理论，高阶剪切变形壳理论，和三维板壳理论的理论体系及基本公式，并对几种理论之间的联系和应用进行了总结和讨论，简述了近几十年来国内外学者在双曲率壳非线性弯曲、稳定性和振动等方面的最新研究成果，最后对双曲率壳体研究目前的局限性和未来的研究方向进行了探讨．     

基于Ω形蜂窝的夹芯式大变形柔性蒙皮结构设计与性能分析

变体飞机根据不同飞行状态改变机翼形状以获得最佳的气动效能，从而提高在多工况下执行多目标任务的能力。表面结构的柔性蒙皮设计是实现变体结构的关键，需要设计具有足够协同变形能力的同时，还能够承受表面的气动载荷、保持光滑气动外形性能的柔性蒙皮结构。针对这一需要，本文提出一种基于Ω形蜂窝的夹芯式大变形柔性蒙皮结构（Ω形柔性蒙皮）。利用参数优化技术以最小化构型变形方向的刚度为目标确定Ω构型的拓扑结构；推导出Ω形蜂窝柔性蒙皮结构的面内等效弹性性能预测公式，并通过数值仿真分析了该结构的面内变形性能和面外承载能力。通过与蛇形蜂窝和正余弦蜂窝两种夹芯式柔性蒙皮结构的面内变形能力的对比分析，验证了该柔性蒙皮的良好性能。结果表明，所提出的Ω形柔性蒙皮结构具有良好的面内可恢复的大变形能力和面外承载能力。     

Modeling fracture in large scale shell structures

We consider a problem of modeling fracture and failure preceded by large scale yielding of ductile shells from the point of view of large-scale structural analysis. We place a special emphasis on the computational efficiency of the constitutive formulation. In this context, we seek the formulation embedded in the shell mechanics framework, which is both theoretically sound and easily implementable into a large-scale explicit dynamic finite element code without precluding vectorization or parallelization. This is achieved through the elasto-plastic damage constitutive model for finite-element analysis of plates and shells. The proposed damage model is purely phenomenological with a scalar damage parameter, which has no physical interpretation, except that it represents on a global scale the micromechanical changes the material undergoes during the process of necking and fracture. The localization leading to softening and fracture is represented by the damage calibration function with exponential damage growth after the onset of necking. The proposed phenomenological damage model uses a general plasticity and shell mechanics frameworks which makes it general and easily implementable into existing finite element codes. The proposed formulation has been implemented into the explicit dynamic finite element software code EPSA (Atkatsh et al., 1980, Atkatsh et al., 1983).     

参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

AN INTERVAL FINITE ELEMENT METHOD BASED ON THE NEUMANN SERIES EXPANSION 1)

实际工程问题中通常存在大量的不确定参数, 区间有限元方法是一种结合有限元数值计算工具对结构进行不确定性分析的区间方法. 区间有限元的目的是获得在含有区间不确定性参数条件下的结构响应上下边界, 其关键问题在于区间平衡方程组的求解, 而这属于一类往往很难求解的NP-hard问题. 本文归纳了一类工程实际中常见的结构不确定性问题, 即可线性分解式区间有限元问题, 并针对此提出一种基于Neumann级数的区间有限元方法. 在区间有限元分析中, 当区间不确定参数表示为一组独立区间变量线性叠加时, 若结构的刚度矩阵也可表示为这些独立区间变量的线性叠加形式, 则称此类区间有限元问题为可线性分解式区间有限元问题. 对于此类问题, 采用Neumann级数对其刚度矩阵的逆矩阵进行表示, 可获得结构响应关于区间变量的显式表达式, 从而可高效求解结构响应的上下边界. 最后通过两个算例验证了本文所提方法的有效性.     

On the influence of interfacial properties to the bending rigidity of layered structures

Layered structures are ubiquitous, from one-atom thick layers in two-dimensional materials, to nanoscale lipid bi-layers, and to micro and millimeter thick layers in composites. The mechanical behavior of layered structures heavily depends on the interfacial properties and is of great interest in engineering practice. In this work, we give an analytical solution of the bending rigidity of bilayered structures as a function of the interfacial shear strength. Our results show that while the critical bending stiffness when the interface starts to slide plastically is proportional to the interfacial shear strength, there is a strong nonlinearity between the rigidity and the applied bending after interfacial plastic shearing. We further give semi-analytical solutions to the bending of bilayers when both interfacial shearing and pre-existing crack are present in the interface of rectangular and circular bilayers. The analytical solutions are validated by using finite element simulations. Our analysis suggests that interfacial shearing resistance, interfacial stiffness and preexisting cracks dramatically influence the bending rigidity of bilayers. The results can be utilized to understand the significant stiffness difference in typical biostructures and novel materials, and may also be used for non-destructive detection of interfacial crack in composites when stiffness can be probed through vibration techniques.     

特定方向"零膨胀"的最小柔顺性结构优化设计

工程中很多承载结构必须面对苛刻的温度变化工作环境,如卫星天线、太空照相机和电子器件等。剧烈的温度变化引起较大的热变形,造成仪器信号失真,精度下降;同时温度应力也会造成结构破坏甚至失效,因此零膨胀材料的研制备受关注。近年来国内外很多学者对此进行了研究,设计出具有特定等效膨胀系数的微结构,但考虑到制备工艺的限制,这类具有复杂微结构的材料制备起来比较困难,成本较高;同时这类材料一般不具备足够的刚度,难以满足承载性能的要求。本文基于结构优化设计技术,采用拓扑优化方法直接设计出具备较高的承载性能和特定方向变形较少受热载荷影响的结构。本文提出采用多目标优化的方法设计圆环结构,使其具有较高的刚度和在热载荷下圆环内表面具有较好的热几何稳定性。由于用单相材料无法同时满足高刚度和低热膨胀的要求,因此假设结构由两种不同的材料构成,用连续体拓扑优化的方法设计三相材料(两种实体材料MAT-I、MAT-II和空材料)在设计域上的最优分布,使结构满足设计要求。由对称性,设计域取为圆环的一个扇面,将设计域离散成有限元网格,每个单元具有两个设计变量:实体材料的体分比和MAT-I在实体材料中所占的体分比,采用伴随法进行灵敏度分析,用GCMMA方法求解此问题,采用体积守恒的Heaviside密度过滤函数保证获得清晰的最优拓扑构型以及避免棋盘格式的出现。通过两个数值算例,表明使用本文提出的多目标优化模型能够得到特定方向"零膨胀"同时具有一定刚度的结构设计,且这种宏观结构尺度上的两种材料组成的拓扑构型相对易于制造。