一种预测周期性材料性能的不确定性分析方法技术

本发明专利技术涉及一种预测周期性材料性能的不确定性分析方法，该方法先利用理论推导设计了一种细观微结构胞元加载边界条件，并证明该条件下通过有限元仿真计算所求得的等效性能可以代表整个材料的等效性能。本发明专利技术在对材料典型微细单胞尺度构型分析获取材料等效弹性模量过程中，充分考虑胞元尺寸与基体材料性能的不确定性，分析不确定性在宏细观等效过程中的传播，得到材料性能的实际波动范围，确保结果的安全性。

An uncertainty analysis method for predicting the properties of periodic materials

The invention relates to a method for predicting periodic material performance uncertainty analysis method. The method uses theory to design a mesoscopic microstructure cell loading boundary conditions, and prove the equivalent performance equivalent performance obtained by the finite element simulation under the condition of calculation can represent the entire material. The invention of materials in typical micro cell scale configuration analysis to obtain material equivalent elastic modulus in the process of fully considering the cell size and the properties of the matrix material uncertainty, uncertainty analysis in macro and micro propagation equivalent process, get the actual fluctuation of material properties of the range, to ensure the safety of the.

【技术实现步骤摘要】

本专利技术涉及周期性材料不确定等效分析方法领域，特别涉及一种预测周期性材料性能的不确定性分析方法；
技术介绍
工程实际常用含微细观结构的材料一般为有序或无序的轻质多孔金属或非金属材料，如类桁架材料、蜂窝材料与泡沫材料等。这类材料因其具有的轻质、高比强度、高比刚度与多功能性，逐渐在实际工程中不断推广应用。有序轻质多孔材料的特点是可以被表示为典型微结构胞元在不同方向的周期性排布，通过对微结构胞元的分析可获取整个材料的宏观力学性能如弹性模量。通常这种分析可利用代表体元法进行，即对满足周期性的典型胞元加上一定的周期性条件与位移边界条件，利用有限元计算的方法获取典型胞元的响应值，通过支反力与位移计算得到周期性材料的等效弹性模量，从而实现对周期性材料性能的预测。工程物理系统多同时受到不同种类不确定性的影响，如材料性质、几何特性、边界条件和载荷分布等，可使用随机变量或者区间变量对不确定性进行量化，而系统的响应可表示为多元随机变量或区间变量的函数，在周期性材料宏微观性能等效分析过程中，不确定性波动将对等效后材料弹性性能造成影响，真实弹性性能将在一定区间内随机分布。
技术实现思路
本专利技术要解决的技术问题为：建立一种预测周期性材料性能的不确定性分析方法，周期性材料的弹性模量由典型微结构胞元的有限元计算得到，在分析中施加合适的周期性条件与边界条件以模拟周期性材料内部位移与变形关系。通过对微结构胞元基体材料属性与几何尺寸加以一定不确定性波动，利用宏微观等效过程中的不确定传播可得到材料弹性性能的不确定性结果。本专利技术解决上述技术问题采用的技术方案为：一种预测周期性材料性能的不确定性分析方法，包括以下步骤：步骤(1)针对实际目标材料内部细观尺度的周期性排布规律，选取有效的立方形微结构胞元，该胞元在各个方向进行阵列能够重组复现材料内部构型分布规律，且满足对称性要求，记为典型胞元。步骤(2)对于上一步得到的典型胞元利用工程软件的CAD工具建立数字化几何模型，以关键点确定胞元的关键位置坐标，以线连结点确定胞元的关键结构，从而得到简化的胞元几何模型，上述即为几何建模过程。步骤(3)将典型胞元几何模型进一步转化为可用于数值计算的有限元模型，基于CAE工具，读入基体材料的基本性能数据与胞元尺寸数据，对所述典型胞元几何模型进行参数设置、材料属性设置、网格划分的有限元建模过程，以梁单元模拟胞元关键结构，各梁即为前述几何模型中关键点连线，各梁端点即为关键点，最终建立梁单元模拟下的典型胞元有限元模型，上述即为有限元建模过程。步骤(4)得到有限元模型后，施加特定周期性条件以模拟典型胞元在材料中的相互作用，在CAE工具内体现为关键点之间的位移关系方程，以实现不同关键点之间位移耦合的效果。步骤(5)进一步施加边界条件以模拟受到外载荷后通过材料分布在单个典型胞元上的位移与力的作用，在CAE工具内体现为典型胞元立方体顶点的强制位移与固定。进一步的，边界条件基于有效介质中总应变能U与非均匀代表胞元中的应变能U′相等，通过体积均匀化代表胞元的应力和应变张量得到宏观应力和宏观应变：其中，为体积均匀化代表胞元的宏观应力，为体积均匀化代表胞元的宏观应变，V为体积，σij为胞元内各点应力，εij为胞元内各点应变。有效介质中总应变能U与非均匀代表胞元中的应变能U′分别为：其中，U为有效介质中总应变能，U′为非均匀代表胞元中，为体积均匀化代表胞元的宏观应力，为体积均匀化代表胞元的宏观应变，V为体积，σij为胞元内各点应力，εij为胞元内各点应变，ui为位移，ui为位移。二者相减得到：即在该定义下的平均应力与平均应变可以满足宏微观等效前后材料内应变能相等，即等效介质可以用来替代非均匀的真实介质描述材料的弹性性能。步骤(6)施加周期性条件与边界条件后，在CAE工具内进行仿真计算，并利用后处理功能提取得到各关键点支反力、位移等结果数据，根据支反力与位移计算相应的胞元等效模量，将模量以结果文件的形式输出，上述操作即为计算与输出过程。步骤(7)将整个流程中所包含的几何建模过程、有限元建模过程、施加周期性条件过程、施加边界条件过程和计算与输出过程编写为APDL命令语句的形式，将语句整合为自动处理命令流文件。利用CAE软件调用命令流文件可以自动读入基体材料属性与胞元尺寸参数，生成相应几何模型与有限元模型，进行计算并输出相应结果文件，为不确定性分析做准备。步骤(8)根据一定间隔与范围，对基体材料属性与几何尺寸的平均值进行扰动，分别得到一定间隔的基体材料属性与几何尺寸的可能取值，对两组取值进行两两匹配，得到基体材料属性与胞元尺寸分别独立变化下的组合取值方案合集。步骤(9)利用上一步得到的自动处理命令流文件对合集中每一种方案进行分析，所有方案结果组合即为在典型微结构胞元性能等效中的不确定性传播结果变化情况，将所有结果中最大值与最小值分别取出即可得到典型微结构胞元性能等效结果的不确定性变化范围。本专利技术与现有技术相比的优点在于：本专利技术以简便方法建立有限元模型，通过周期性条件与边界条件模拟出典型微结构胞元在材料内部的变形与位移特点，并可方便计算出在特定基体材料属性与几何尺寸下周期性材料的等效弹性性能。本专利技术中的几何建模过程、有限元建模过程、施加周期性条件过程、施加边界条件过程和计算与输出过程都是通过自编语句实现，可针对特定输入参数自动计算并输出结果。同时，利用初始设计参数的不确定性波动，实现独立参数的不确定性方案组合，利用自动处理语句计算大量设计点，产生一系列弹性模量计算响应结果，以获取不确定性传播下周期性材料性能预测结果的波动范围，拓展了确定性下弹性性能唯一值，更加符合真实加工误差下所得材料的性能分布情况，比起常规方法，本专利技术所涉及方法具有更大的安全性与便捷性优势。附图说明图1为本专利技术中周期性材料阵列取典型胞元示意图；图2为本专利技术中典型胞元代表立方体关键位置编号示意图；图3为本专利技术典型胞元有限元模型示意图；图4为本专利技术中不确定性扰动后计算结果变化范围示意图；图5为本专利技术的方法实现流程图。具体实施方式下面结合附图以及具体实施方式进一步说明本专利技术。如图5所示，本专利技术预测周期性材料性能的不确定性传播分析方法，包括以下步骤：(1)针对实际目标材料，根据其内部细观尺度下微结构排布规律，选取合适的立方体块作为有效胞元，以如图1所示构型为例，所选取典型胞元在三轴方向进行阵列可重组复现材料内部周期性构型分布规律，同时，该胞元应满足三个方向对称性要求，以便施加后续条件。(2)对于上一步得到的典型胞元，在ANSY软件前处理器中建立数字化几何模型，以立方体顶点与体心作为关键点，关键点之间用线连结，以得到简化的胞元几何模型。(3)将胞元几何模型进一步转化为可用于数值计算的有限元模型，在ANSYS软件中，利用命令流读入基体材料基本性能数据(包括弹性模量与泊松比)和胞元尺寸数据(包括胞元内部连结梁的直径)，对所述典型胞元几何模型进行参数设置、材料属性设置、网格划分的有限元建模过程，以梁单元模拟胞元内部连接梁结构，最终建立梁单元模拟下的典型胞元有限元模型，如图3所示。(4)在典型胞元有限元模型中，对立方体关键位置节点施加周期性条件，各顶点与面心点编号如图2所示，对立方体顶点处节点施加周期性条件如下：其中，为各编号节点的位本文档来自技高网...

【技术保护点】
一种预测周期性材料性能的不确定性分析方法，其特征在于实现步骤如下：步骤(1)针对实际目标材料内部细观尺度的周期性排布规律，选取有效的立方形微结构胞元，该胞元在各个方向进行阵列能够重组复现材料内部构性分布规律，且满足对称性要求，记为典型胞元；步骤(2)对于上一步得到的典型胞元利用工程软件的CAD(计算机辅助设计)工具建立数字化几何模型，以关键点确定胞元的关键位置坐标，以线连结点确定胞元的关键结构，从而得到简化的胞元几何模型，上述即为几何建模过程；步骤(3)将典型胞元几何模型进一步转化为可用于数值计算的有限元模型，基于CAE(计算机辅助工程)工具，读入基体材料的基本性能数据与胞元尺寸数据，对所述典型胞元几何模型进行参数设置、材料属性设置、网格划分的有限元建模过程，以梁单元模拟胞元关键结构，各梁即为前述几何模型中关键点连线，各梁端点即为关键点，最终建立梁单元模拟下的典型胞元有限元模型，上述即为有限元建模过程；步骤(4)得到有限元模型后，施加特定周期性条件以模拟典型胞元在材料中的相互作用，在CAE工具内体现为关键点之间的位移关系方程，以实现不同关键点之间位移耦合的效果；步骤(5)进一步施加边界条件以模拟受到外载荷后通过材料分布在单个典型胞元上的位移与力的作用，在CAE工具内体现为典型胞元立方体顶点的强制位移与固定；步骤(6)施加周期性条件与边界条件后，在CAE工具内进行仿真计算，并利用后处理功能提取得到各关键点支反力、位移等结果数据，根据支反力与位移计算相应的胞元等效模量，将模量以结果文件的形式输出，上述操作即为计算与输出过程；步骤(7)将整个流程中所包含的几何建模过程、有限元建模过程、施加周期性条件过程、施加边界条件过程和计算与输出过程编写为APDL命令语句的形式，将语句整合为自动处理命令流文件，利用CAE软件调用命令流文件可以自动读入基体材料属性与胞元尺寸参数，生成相应几何模型与有限元模型，进行计算并输出相应结果文件，为不确定性分析做准备；步骤(8)根据一定间隔与范围，对基体材料属性与几何尺寸的平均值进行扰动，分别得到一定间隔的基体材料属性与几何尺寸的可能取值，对两组取值进行两两匹配，得到基体材料属性与胞元尺寸分别独立变化下的组合取值方案合集；步骤(9)利用上一步得到的自动处理命令流文件对合集中每一种方案进行分析，所有方案结果组合即为在典型微结构胞元性能等效中的不确定性传播结果变化情况，将所有结果中最大值与最小值分别取出即可得到典型微结构胞元性能等效结果的不确定性变化范围。...

【技术特征摘要】
1.一种预测周期性材料性能的不确定性分析方法，其特征在于实现步骤如下：步骤(1)针对实际目标材料内部细观尺度的周期性排布规律，选取有效的立方形微结构胞元，该胞元在各个方向进行阵列能够重组复现材料内部构性分布规律，且满足对称性要求，记为典型胞元；步骤(2)对于上一步得到的典型胞元利用工程软件的CAD(计算机辅助设计)工具建立数字化几何模型，以关键点确定胞元的关键位置坐标，以线连结点确定胞元的关键结构，从而得到简化的胞元几何模型，上述即为几何建模过程；步骤(3)将典型胞元几何模型进一步转化为可用于数值计算的有限元模型，基于CAE(计算机辅助工程)工具，读入基体材料的基本性能数据与胞元尺寸数据，对所述典型胞元几何模型进行参数设置、材料属性设置、网格划分的有限元建模过程，以梁单元模拟胞元关键结构，各梁即为前述几何模型中关键点连线，各梁端点即为关键点，最终建立梁单元模拟下的典型胞元有限元模型，上述即为有限元建模过程；步骤(4)得到有限元模型后，施加特定周期性条件以模拟典型胞元在材料中的相互作用，在CAE工具内体现为关键点之间的位移关系方程，以实现不同关键点之间位移耦合的效果；步骤(5)进一步施加边界条件以模拟受到外载荷后通过材料分布在单个典型胞元上的位移与力的作用，在CAE工具内体现为典型胞元立方体顶点的强制位移与固定；步骤(6)施加周期性条件与边界条件后，在CAE工具内进行仿真计算，并利用后处理功能提取得到各关键点支反力、位移等结果数据，根据支反力与位移计算相应的胞元等效模量，将模量以结果文件的形式输出，上述操作即为计算与输出过程；步骤(7)将整个流程中所包含的几何建模过程、有限元建模过程、施加周期性...

【专利技术属性】
技术研发人员：王磊，蔡逸如，王晓军，刘易斯，刘东亮，夏海军，耿新宇，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11