两种五零能模式材料单胞构型及其性能分析

五零能模式材料是一种新型的人工超材料,虽属于弹性材料,但组成其单胞的特殊构型使其宏观静态表现为仅能承载一种受力状态,动态表现为仅能传播一种弹性波。本文首先构造了两种五零能模式材料的单胞构型,其具有不同的弹性特性,其中一种材料可传播弹性膨胀波,另一种可传播弹性剪切波。然后分别采用代表体元法和均匀化法分析这两种单胞的等效弹性模量。五零能模式材料的分析分为两步更直观,开始从单胞桁架模型入手,检验单胞构型是否满足五零能模式的定义,然后分析单胞实体模型,考察单胞构型的结构参数与其等效弹性模量的关系。研究表明对于这种低密度弹性材料的分析,代表体元法更适合。     

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

周期性多孔材料等效剪切模量与尺寸效应研究

提出了以圆筒扭转力学模型为基础, 预测周期性多孔材料等效剪切模量及其研究尺寸效 应的一种简单有效计算方法. 以方形孔和圆形孔两种典型多孔材料为例进行了数值计算求解; 同时, 建立了几何参数随体胞尺寸缩放因子$n$的解析关系, 证明了两种构型的周期性多孔材 料的等效剪切模量均随尺寸缩放因子$n$的增大而减小. 当$n \to \infty $时, 即体胞尺寸相对整体结构无限小时, 多孔材料的等效剪切模量趋近收敛于一个恒定值; 当体胞的材料体分比增大时, 多孔材料等效剪切模量也随之增大. 此外, 依据周 期性多孔材料的结构对称特性, 使用体胞子结构有限元计算模型进行等效剪切模量及其尺寸 效应的预测, 极大地提高了计算效率.     

高承载梯度分层点阵结构的拓扑优化设计方法

随着增材制造技术的迅速发展, 点阵结构由于其高比强度、高比刚度等优异力学性能受到广泛关注, 但其单胞分布设计大多基于均布式假设, 导致其承载能力相对较差. 基于拓扑优化技术提出了一种梯度分层的点阵结构设计方法. 首先, 基于水平集函数建立点阵单胞几何构型的显式描述模型, 引入形状插值技术实现点阵单胞的梯度构型生成; 其次, 构建基于Kriging的梯度点阵单胞宏观等效力学属性预测模型, 建立宏观有限单元密度与微观点阵单胞等效力学属性的内在联系; 然后, 以点阵结构刚度最大为优化目标, 结构材料用量和力学控制方程为约束条件, 构建点阵结构的梯度分层拓扑优化模型, 并采用OC算法进行数值求解. 算例结果表明, 所提方法可实现点阵结构的最优梯度分层设计, 充分提高了点阵结构的承载性能, 同时可保证不同梯度点阵单胞之间的几何连续性. 最后, 开展梯度分层点阵结构与传统均匀点阵结构和线性梯度点阵结构的准静态压缩仿真分析, 仿真结果表明, 与传统均匀点阵结构和线性梯度点阵结构相比, 梯度分层点阵结构的承载能力明显提高. 研究结果可为高承载点阵结构设计提供理论参考.      

PREDICTIONS OF EFFECTIVE OUT-PLANE SHEAR MODULUS AND SIZE EFFECT OF HEXAGONAL HONEYCOMB

给出了预测六边形蜂窝材料等效剪切模量及其尺寸效应的圆筒扭转力学模型和扭转能量法,建立了等效面外剪切模量G₁₃相对于材料体分比ν、周向单胞数n、圆筒半径r和单胞层数参数m变化的解析表达式;同时将扭转能量法、有限元数值模拟计算和G-A经典细观力学方法进行了比较,从理论上揭示并验证了尺寸效应的存在性. 结果表明,当蜂窝体胞尺寸相对结构尺寸无穷小时,预测结果趋近于细观力学方法的结果. 此外,利用周期性蜂窝材料的结构对称特性,使用体胞子结构有限元计算模型进行等效面外剪切模量及其尺寸效应的预测,在不影响计算结果的前提下极大地提高了计算效率.     

一种基于平均场理论的格栅材料均匀化方法

将颗粒材料中发展的一种基于平均场理论的解析均匀化方法应用于二维周期格栅材料;依据尺度分离原理和统计均匀表征元概念构建了格栅材料的两尺度均匀化模型,包括细观杆件单元的本构关系、细观位移-宏观应变关系式以及应力的细观力学表达式;推导了两种二维周期格栅材料等效弹性参数包括弹性模量、泊松比和剪切弹性模量的细观力学表达式。结果表明:等边三角形结构等效为各向同性连续体时,弹性参数表达式与文献中其他方法所得结果一致;正方形结构均匀化为正交各向异性连续体时,主平面内弹性模量等于杆件单元轴向刚度,泊松比和剪切弹性模量分别由杆件单元的泊松比和剪切刚度决定,符合正方形格栅材料的力学特性;对于非主平面内的正方形本构矩阵,选取坐标轴与材料主轴夹角为45°的方向为例进行推导,本文方法与坐标变换方法所得结果一致。以上结果均验证了本文所发展方法的有效性。     

二维类桁架材料结构弹塑性分析

类桁架材料所构成结构的弹塑性行为的精确建模分析保证非常耗时, 为了 在保证精度的前提下提高此类问题的求解效率, 本文利用类桁架材料基本构件长细比 较大的特点，将材料单胞简化为桁架模型. 考虑到微单胞空间分布的周期性，基于 数值均匀化理论提出了类桁架材料结构的宏微观两 级弹塑性求解格式. 原问题转化为宏观上一个非线性弹性连续体计算问题和微观上多个小规 模桁架系统的弹塑性计算问题. 两个数值算例分别考虑了简单加载，非单调加载，规则宏观 结构和具有非完整单胞的较复杂宏观结构等问题. 与实际结构计算结果在精度和时间等方面 的比较验证了求解格式的有效性. 最后还探讨了算法的适用范围.     

利用松质骨理想化模型对骨小梁刚度的研究

利用结构张量的概念建立了考虑松质骨结构各向异性的理想模型，将杆型与板型两种传统模型统一起来．采用均匀化理论，通过有限单元的数值计算建立了松质骨弹性常数与结构张量、固体体积比之间的数值关系，并预测了骨小梁的弹性模量．结果表明，本文模型较好地体现了松质骨的各向异性力学性质．对骨小梁的计算结果表明，骨小梁的弹性模量在６７６ＧＰａ～１０９ＧＰａ范围内，平均意义上的结果是９２１ＧＰａ．     

碳纳米管在简单荷载作用下的变形屈曲行为研究

采用分子动力学方法研究了不同螺旋性的单壁碳纳米管在轴压、扭转和外压分别作用下的变形屈曲行为。针对每种加载形式，给出了宏观形式的内力一变形曲线来反映碳原子间相互作用力的变化与碳管变形量之间的关系。对碳管在这三种不同荷载作用下的失稳特性、微观结构和应变能的变化进行了详细的分析。结合经典弹性理论和本文的模拟结果，得出了单壁碳纳米管的等效壁厚、弹性模量和Poisson比。另外，通过对比不同螺旋性碳管的模拟结果，讨论了螺旋性对碳纳米管力学行为的影响。文中采用的势能函数为Brenner多体势。     

循环对称结构的多尺度拓扑优化方法

研究了循环对称结构多尺度拓扑优化问题,阐明了均匀化等效性能的基本性质,提出了循环对称单胞的特征参数概念,建立了均匀化映射计算方法,构造了等效性能的三元参数插值模型,有效简化了不同特征参数、微结构构型与体分比下的单胞均匀化等效过程.通过微结构的特征驱动建模与B样条参数化,克服了微结构拓扑优化变量多和变量离散难以保证其光滑连接的问题.给出了典型数值算例,比较了等比例排列单胞与等间距排列单胞对结构优化结果的影响,验证了多尺度优化方法的有效性.     

Tissue level microstructure and mechanical properties of the femoral head in the proximal femur of fracture patients

This study aims to investigate the regional variations of trabecular morphological parameters and mechanical parameters of the femoral head,as well as to determine the relationship between trabecular morphological and mechanical parameters.Seven femoral heads from patients with fractured proximal femur were scanned using a micro-CT system.Each femoral head was divided into 12 sub-regions according to the trabecular orientation.One 125 mm~3 trabecular cubic model was reconstructed from each sub-region.A total of 81 trabecular models were reconstructed,except three destroyed sub-regions from two femoral heads during the surgery.Trabecular morphological parameters,i.e.trabecular separation(Tb.Sp),trabecular thickness(Tb.Th),specific bone surface(BS/B V),bone volume fraction(BV/TV),structural model index(SMI),and degree of anisotropy(DA) were measured.Micro-finite element analyses were performed for each cube to obtain the apparent Young's modulus and tissue level von Mises stress distribution under 1%compressive strain along three orthogonal directions,respectively.Results revealed significant regional variations in the morphological parameters(P0.05).Young's moduli along the trabecular orientation were significantly higher than those along the other two directions.In general,trabecular mechanical properties in the medial region were lower than those in the lateral region.Trabecular mechanical parameters along the trabecular orientation were significantly correlated with BS/BV,BV/TV,Tb.Th,and DA.In this study,regional variations of microstructural features and mechanical properties in the femoral head of patients with proximal femur fracture were thoroughly investigated at the tissue level.The results of this study will help to elucidate the mechanism of femoral head fracture for reducing fracture risk and developing treatment strategies for the elderly.     

An Experimental–Numerical Methodology for a Rapid Prototyped Application Combined with Finite Element Models in Vertebral Trabecular Bone

In the present study, a novel evaluation method involving rapid prototyped (RP) technology and finite element (FE) analysis was used to study the elastic mechanical characteristics of human vertebral trabecular bone. Three-dimensional (3D) geometries of the RP and FE models were obtained from the central area of vertebral bones of female cadavers, age 70 and 85. RP and FE models were generated from the same high-resolution micro-computed tomography (μCT) scan data. We utilized RP technology along with FE analysis based on μCT for high-resolution vertebral trabecular bone specimens. RP models were used to fabricate complex 3D objects of vertebral trabecular bone that were created in a fused deposition modeling machine. RP models of vertebral trabecular bone are advantageous, particularly considering the repetition, risks, and ethical issues involved in using real bone from cadaveric specimens. A cubic specimen with a side length of 6.5 mm or a cylindrical specimen with a 7 mm diameter and 5 mm length proved better than a universal cubic specimen with a side length of 4 mm for the evaluation of elastic mechanical characteristics of vertebral trabecular bones through experimental and simulated compression tests. The results from the experimental compression tests of RP models closely matched those predicted by the FE models, and thus provided substantive corroboration of all three approaches (experimental tests using RP models and simulated tests using FE models with ABS and trabecular bone material properties). The RP technique combined with FE analysis has potential for widespread biomechanical use, such as the fabrication of dummy human skeleton systems for the investigation of elastic mechanical characteristics of various bones.     

基于近场动力学的单晶冰弹性各向异性数值模拟方法

天然冰材料在变形与破坏行为上的各向异性特征是冰与结构相互作用中产生复杂载荷过程的关键诱因, 而天然冰各向异性的根源则在于单晶冰的各向异性. 目前, 学术界针对单晶冰各向异性的数值模拟方法研究仍较为缺乏. 为了准确再现天然冰材料的特殊力学性质, 本文基于近场动力学理论, 提出了一种单晶冰弹性各向异性的数值模拟方法. 该方法的核心思想是将单晶冰杨氏模量沿不同加载方向的变化规律引入到近场动力学力密度向量的影响函数中. 以前人实验测试得到的杨氏模量值为参考, 通过开展与C轴呈0°, 45°和90°三个加载方向的单晶冰单轴压缩数值模拟实验, 提出了针对该影响函数的修正和辅助参数标定方法, 最终在15°, 30°, 60°和75°等其他四个加载方向进行了验证. 结果表明: 本文提出的针对影响函数的修正与参数标定方法, 能够较为便捷地找到数值模型杨氏模量与参考杨氏模量相一致的影响函数最优解, 即本文提出的基于影响函数的近场动力学数值模拟方法, 能够合理、准确地模拟单晶冰的弹性各向异性行为. 本文研究成果可为后续多晶冰各向异性数值模拟方法的建立提供基础性参考.      

3D FINITE ELEMENTANALYSIS OF THE DAMAGE EFFECTS ON THE DENTAL COMPOSITE SUBJECT TO WATER SORPTION

The damage effects of water sorption on the mechanical properties of the hydroxya-patite particle reinforced Bis-GMA/TEGDMA copolymer (HA/Bis-GMA/TEGDMA) have been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young's modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young's modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. The former results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.     

基于虚内键理论的材料多尺度力学模型

宏观上线弹性材料的力学属性只需杨氏模量和泊松比两个相互独立的参量来控制;相应地,微观上也需要两个相互独立的参量来控制.基于这个思想,在原VIB模型中引入了切向键,并提出了VMIB模型.该模型在材料的宏观力学属性与微观虚拟键力学属性之间建立起了一座桥梁.考虑到模型中能量密度函数含有坐标轴方向向量一项,该文对能量密度函数的张量性进行了严格的数学证明,并将VMIB模型初步应用到脆性材料的单轴受压破坏.     

Age- and direction-related adaptations of lumbar vertebral trabecular bone with respect to apparent stiffness and tissue level stress distribution

The objective of this study was to study the age-related adaptation of lumbar vertebral trabecular bone at the apparent level, as well as the tissue level in three orthogonal directions. Ninety trabecular specimens were obtained from six normal L4 vertebral bodies of six male cadavers in two age groups, three aged 62 years and three aged 69 years, and were scanned using a high-resolution micro-computed tomography （micro-CT） system, then converted to micro- finite element models to do micro-finite element analyses. The relationship between apparent stiffness and bone volume fraction, and the tissue level yon Mises stress distribution for each trabecular specimen when compressed separately in the longitudinal direction, medial-lateral and anterior-posterior directions （transverse directions） were derived and compared between two age groups. The results showed that at the apparent level, trabecular bones from 69-year group had stiffer bone structure relative to their volume fractions in all three directions, and in both age groups, changes in bone volume fraction could explain more variations in apparent stiffness in the longitudinal direction than the transverse directions; at the tissue level, aging had little effect on the tissue von Mises stress distributions for the compressions in all the three directions. The novelty of the present study was that it provided quantitative assessments on the age and direction- related adaptation of Chinese male lumbar vertebral trabecular bone from two different levels： stiffness at the apparent level and stress distribution at the tissue level. It may help to understand the failure mechanisms and fracture risks of vertebral body associated with aging and direction for the prevention of fracture risks in elder individuals.     

Theoretical investigation of the behavior of an acoustic metamaterial with extreme Young's modulus

A mechanical model with local resonators is proposed as an acoustic metamaterial that exhibits an unusual frequency-dependent effective stiffness. If treated as an equivalent elastic solid, its effective Young's modulus can become unbounded or vanishingly small at two respective frequencies. Moreover, in a certain frequency range, the effective Young's modulus would become negative, resulting in a band gap that coincides with this frequency range. The wave attenuation behavior and mechanism are studied through numerical simulations on the acoustic metamaterial model. The capability of the metamaterial to selectively block or filter unwanted waves is demonstrated by a numerical simulation example.     

纳米铜丝尺寸效应的分子动力学模拟

采用EAM原子势函数对不同截面尺寸纳米铜丝拉伸性能进行了零温分子动力学模拟。研究表明截面变化对纳米丝拉伸性能有明显影响。由于表明原子弛豫，纳米丝存在表面张应力以及与之相平衡的内部应力。这种本征应力的存在是纳米丝尺寸效应的根源。纳米丝截面减小，则拉伸强度提高、屈服推迟、初始拉伸模量的软化程度增加。     

不同模量弹性问题理论及有限元法研究进展

随着科学技术的日益发展，对材料力学性质的研究提出了更高的要求，研制新型的材料以及挖掘材料自身特性的潜力，已成为新的研究动向．简述了不同模量弹性问题理论及其有限元法的研究与发展．利用等效的概念，对Ambartsumyan有限元计算模型、Jones有限元计算模型、张允真等有限元计算模型、叶志明等计算模型进行改进和探讨．通过不同模量弹性理论及其有限元方法在实际工程构件问题分析中的应用表明，若沿用相同模量弹性理论或有限元对有关问题进行计算，其结果将与采用不同模量模型材料所得的刚度和强度有较大的偏差。     

Dependency of observation parameters on soil dynamic parameters in unconfined impact compression tests

Unconfined impact compression tests were used to measure dynamic parameters of soil samples in the domains of stress wave propagation and dynamic deformation. The high-frequency equivalent Young's modulus, which was determined by the reflection theory of one-dimensional elastic stress waves, increased exponentially with water content. The low-frequency equivalent Young's modulus, which was defined as the initial gradient of a stress-strain curve, decreased exponentially with increasing water content. The high-frequency equivalent Young's modulus was one to two orders of magnitude larger than the low-frequency one. Water content, height of soil samples and drop height affected the dynamic stress-strain relationships of the soil samples. The Maxwell's viscoelastic fluid model described the stress-strain relationships well, and the values of the model parameters varied with water content and frequency.