CN116882239A - 一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 - Google Patents
一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 Download PDFInfo
- Publication number
- CN116882239A CN116882239A CN202310837525.XA CN202310837525A CN116882239A CN 116882239 A CN116882239 A CN 116882239A CN 202310837525 A CN202310837525 A CN 202310837525A CN 116882239 A CN116882239 A CN 116882239A
- Authority
- CN
- China
- Prior art keywords
- composite material
- component
- representative volume
- mesoscopic
- volume unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 113
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 49
- 239000011159 matrix material Substances 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 239000000463 material Substances 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 3
- 229910052782 aluminium Inorganic materials 0.000 claims description 3
- 239000004568 cement Substances 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000452 restraining effect Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 4
- 239000004567 concrete Substances 0.000 description 7
- 239000011083 cement mortar Substances 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000000919 ceramic Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质,属于复合材料仿真技术领域。为解决有效节省运算资源、缩短复合材料仿真计算时间。本发明根据有限元分析方法构建复合材料的代表性体积单元RVE,设置复合材料的代表性体积单元刚度矩阵,将复合材料的代表性体积单元的内部结构划分为三种细观组分,划分RVE的背景网格为细观网格,根据有限元分析方法设置第i种细观组分所有单元的材料属性矩阵为[mi],设置第i种细观组分中第e个单元出现在空间中位置x的概率矩阵为[pi,(e)(x)],则得到第e个单元的单元刚度矩阵[k(e)],基于随机分形理论测算组分空间分布的概率,进而推算复合材料的均质化等效材料属性,可有效节省运算资源、缩短计算时间。
Description
技术领域
本发明属于复合材料仿真技术领域,具体涉及一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质。
背景技术
受组分在空间中的随机分布影响,复合材料的宏观模量和强度呈现明显的非均质特性。为了方便应用,通常采用均质化方法,将非均匀的复合材料分析问题转化为一个对等效均质材料分析的问题。
按基本原理,均质化方法可分为两类:解析方法适用分析夹杂相对规则、空间分布均匀的复合材料,其理论相对完备,但是应用范围较小;数值方法适用范围较广,但是需要针对不同的细观结构分别建模,运算量较大。
发明内容
本发明要解决的问题是为了有效节省运算资源、缩短复合材料仿真计算时间,提出一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质。
为实现上述目的,本发明通过以下技术方案实现:
一种基于组分空间分布的复合材料属性分析方法,其特征在于,包括如下步骤:
S1、根据有限元分析方法构建复合材料的代表性体积单元RVE,设置复合材料的代表性体积单元刚度矩阵为[K];
S2、将复合材料的代表性体积单元的内部结构划分为三种细观组分,包括夹杂、基质、夹杂和基质间界面;
S3、划分RVE的背景网格为细观网格;
S4、根据有限元分析方法设置第i种细观组分所有单元的材料属性矩阵为[mi],设置第i种细观组分中第e个单元出现在空间中位置x的概率矩阵为[pi,(e)(x)],则得到第e个单元的单元刚度矩阵[k(e)],计算表达式为:
[k(e)]=[mi][p(e)(x)];
S5、利用随机分形方法,将[pi,(e)(x)]表达为复合材料的细观组分在三维空间里的分形维数D3和复合材料的细观组分的多重分形测度μ3的函数f(D3,μ3),则得到第e个单元的单元刚度矩阵[k(e)]的计算表达式为:
S6、基于步骤S5得到的第e个单元的单元刚度矩阵计算复合材料的代表性体积单元刚度矩阵[K],计算表达式为:
S7、对步骤S6得到的复合材料的代表性体积单元刚度矩阵施加力边界条件,分析基于组分空间分布的复合材料属性。
进一步的,步骤S1中的复合材料包括水泥基复合材料、陶瓷增强铝基复合材料中的一种。
进一步的,步骤S3的具体实现方法为以dw为横向间距、以dh为竖向间距,在复合材料的代表性体积单元上布置纵横向直线,划分为背景网格。
进一步的,步骤S5中利用计盒方法计算复合材料的细观组分在三维空间里的分形维数,具体实现方法包括如下步骤:
S5.1、在在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的边长分别为δ1,δ2,…,δn;
S5.2、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);
S5.3、以δ1,δ2,…,δn为横坐标,以Ni(δ1),Ni(δ2),…,Ni(δn)为纵坐标,在双对数坐标轴上绘制散点,计算函数lnN(δ)和-lnδ的斜率为复合材料第i种细观组分在三维空间里的分形维数。
进一步的,步骤S5中利用粗线条方法计算复合材料的细观组分的多重分形测度,具体实现方法包括如下步骤:
S5.4、在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的边长分别为δ1,δ2,…,δn;
S5.5、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);根据相交盒子数目与总盒子数目之比,统计所有盒子与第i种细观组分相交的概率为Pi(δ1),Pi(δ2),…,Pi(δn);
S5.6、取q为一定范围内的实数,设置与组分空间分布概率相关的计算过程函数χq(δ)的计算表达式为:
然后计算复合材料的细观组分的多重分形测度,计算公式为:
进一步的,步骤S7的具体实现方法包括如下步骤:
S7.1、对步骤S6得到的复合材料的代表性体积单元刚度矩阵上边界的节点施加力边界条件[F],同时约束下边界节点;
S7.2、由矩阵方程[K][U]=[F]计算位移值[U],进而由力的值与相应位移值之比,计算复合材料的均质化等效材料属性。
电子设备,包括存储器和处理器,存储器存储有计算机程序,所述的处理器执行所述计算机程序时实现所述的一种基于组分空间分布的复合材料属性分析方法的步骤。
计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述的一种基于组分空间分布的复合材料属性分析方法。
本发明的有益效果:
本发明所述的一种基于组分空间分布的复合材料属性分析方法,基于随机分形理论测算组分空间分布的概率,进而推算复合材料的均质化等效材料属性,可有效节省运算资源、缩短计算时间。
附图说明
图1为本发明所述的一种基于组分空间分布的复合材料属性分析方法的流程图;
图2为本发明所述的一种基于组分空间分布的复合材料属性分析方法构建的混凝土细观模型;
图3为本发明所述的一种基于组分空间分布的复合材料属性分析方法构建混凝土细观模型网格示意图;
图4为本发明所述的一种基于组分空间分布的复合材料属性分析方法构建混凝土细观模型单元刚度矩阵示意图;
图5为本发明所述的一种基于组分空间分布的复合材料属性分析方法混凝土弹性模量统计分布直方图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施方式,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅用以解释本发明,并不用于限定本发明,即所描述的具体实施方式仅仅是本发明一部分实施方式,而不是全部的具体实施方式。通常在此处附图中描述和展示的本发明具体实施方式的组件可以以各种不同的配置来布置和设计,本发明还可以具有其他实施方式。
因此,以下对在附图中提供的本发明的具体实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定具体实施方式。基于本发明的具体实施方式,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他具体实施方式,都属于本发明保护的范围。
为能进一步了解本发明的发明内容、特点及功效,兹例举以下具体实施方式,并配合附图1-附图5详细说明如下:
具体实施方式一:
一种基于组分空间分布的复合材料属性分析方法,包括如下步骤:
S1、根据有限元分析方法构建复合材料的代表性体积单元RVE,设置复合材料的代表性体积单元刚度矩阵为[K];
进一步的,步骤S1中的复合材料包括水泥基复合材料、陶瓷增强铝基复合材料中的一种;
S2、将复合材料的代表性体积单元的内部结构划分为三种细观组分,包括夹杂、基质、夹杂和基质间界面;
S3、划分RVE的背景网格为细观网格;
进一步的,步骤S3的具体实现方法为以dw为横向间距、以dh为竖向间距,在复合材料的代表性体积单元上布置纵横向直线,划分为背景网格;
S4、根据有限元分析方法设置第i种细观组分所有单元的材料属性矩阵为[mi],设置第i种细观组分中第e个单元出现在空间中位置x的概率矩阵为[pi,(e)(x)],则得到第e个单元的单元刚度矩阵[k(e)],计算表达式为:
[k(e)]=[mi][p(e)(x)];
S5、利用随机分形方法,将[pi,(e)(x)]表达为复合材料的细观组分在三维空间里的分形维数D3和复合材料的细观组分的多重分形测度μ3的函数f(D3,μ3),则得到第e个单元的单元刚度矩阵[k(e)]的计算表达式为:
进一步的,步骤S5中利用计盒方法计算复合材料的细观组分在三维空间里的分形维数,具体实现方法包括如下步骤:
S5.1、在在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的形状为正方形,设置盒子的边长分别为δ1,δ2,…,δn;
S5.2、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);
S5.3、以δ1,δ2,…,δn为横坐标,以Ni(δ1),Ni(δ2),…,Ni(δn)为纵坐标,在双对数坐标轴上绘制散点,计算函数lnN(δ)和-lnδ的斜率为复合材料第i种细观组分在三维空间里的分形维数;
步骤S5中利用粗线条方法计算复合材料的细观组分的多重分形测度,具体实现方法包括如下步骤:
S5.4、在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的形状为正方形,设置盒子的边长分别为δ1,δ2,…,δn;
S5.5、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);根据相交盒子数目与总盒子数目之比,统计所有盒子与第i种细观组分相交的概率为Pi(δ1),Pi(δ2),…,Pi(δn);
S5.6、取q为一定范围内的实数,设置与组分空间分布概率相关的计算过程函数χq(δ)的计算表达式为:
然后计算复合材料的细观组分的多重分形测度,计算公式为:
S6、基于步骤S5得到的第e个单元的单元刚度矩阵计算复合材料的代表性体积单元刚度矩阵[K],计算表达式为:
S7、对步骤S6得到的复合材料的代表性体积单元刚度矩阵施加力边界条件,分析基于组分空间分布的复合材料属性。
进一步的,步骤S7的具体实现方法包括如下步骤:
S7.1、对步骤S6得到的复合材料的代表性体积单元刚度矩阵上边界的节点施加力边界条件[F],同时约束下边界节点;
S7.2、由矩阵方程[K][U]=[F]计算位移值[U],进而由力的值与相应位移值之比,计算复合材料的均质化等效材料属性。
以复合材料为混凝土为例:将复合材料的代表性体积单元的内部结构划分为三种细观组分,包括夹杂、基质、夹杂和基质间界面如图2所示;
将横向间距dw和竖向间距dh均取为2mm,可划分出如图3红线所示的正方形背景网格;
先按不考虑单元出现在空间中相应位置的概率考虑,按有限单元法,按每一个网格的节点位置,并代入如下由力学试验结果统计所得数据:混凝土骨料弹性模量平均值4.7×104MPa、标准差为0.72×104MPa,泊松比平均值0.18、标准差为0.029;PO 42.5水泥砂浆弹性模量平均值2.01×104MPa、标准差为0.31×104MPa,泊松比平均值0.23、标准差为0.035,均服从正态分布。ITZ弹性模量按水泥砂浆的80%计。可得骨料、水泥砂浆和ITZ的单元刚度矩阵,分别记为[m(1)]、[m(2)]和[m(3)],与普通正方形网格的单元刚度矩阵相同。
先测得每个组分的分形维数,骨料、水泥砂浆和ITZ的分形维数分别为3.0001、2.1450、2.9998,进而计算其多重分形测度,可以测得每个组分在单元格出现在的概率。以图4所示单元为例,在该单元处出现骨料、水泥砂浆和ITZ的概率分别为p(1)=0.48、p(2)=0.51和p(3)=0.01。则此单元刚度矩阵可表示为:
[k(e)]=0.48[m(1)]+0.51[m(2)]+0.01[m(3)]
用此方式表示所有单元的单元刚度矩阵,进而通过直接叠加法计算RVE的刚度矩阵[K]。对RVE边界的节点上施加力边界条件[F],同时约束下边界节点,可列出矩阵方程[K][U]=[F]计算出其位移值[U],进而由力的值与相应位移值之比,推算此混凝土材料的均质化等效材料属性,均值为5.869×104MPa,标准差为1.179×104MPa,偏度系数为-0.0095,其统计分布如图5所示,近似服从正态分布。所得结果可用于后续可靠度分析。
具体实施方式二:
电子设备,包括存储器和处理器,存储器存储有计算机程序,所述的处理器执行所述计算机程序时实现具体实施方式一所述的一种基于组分空间分布的复合材料属性分析方法的步骤。
本发明的计算机装置可以是包括有处理器以及存储器等装置,例如包含中央处理器的单片机等。并且,处理器用于执行存储器中存储的计算机程序时实现上述的一种基于组分空间分布的复合材料属性分析方法的步骤。
所称处理器可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等;存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外,存储器可以包括高速随机存取存储器,还可以包括非易失性存储器,例如硬盘、内存、插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
具体实施方式三:
计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述的一种基于组分空间分布的复合材料属性分析方法。
本发明的计算机可读存储介质可以是被计算机装置的处理器所读取的任何形式的存储介质,包括但不限于非易失性存储器、易失性存储器、铁电存储器等,计算机可读存储介质上存储有计算机程序,当计算机装置的处理器读取并执行存储器中所存储的计算机程序时,可以实现上述的一种基于组分空间分布的复合材料属性分析方法的步骤。
所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括电载波信号和电信信号。
需要说明的是,术语“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
虽然在上文中已经参考具体实施方式对本申请进行了描述,然而在不脱离本申请的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,本申请所披露的具体实施方式中的各项特征均可通过任意方式相互结合起来使用,在本说明书中未对这些组合的情况进行穷举性的描述仅仅是出于省略篇幅和节约资源的考虑。因此,本申请并不局限于文中公开的特定具体实施方式,而是包括落入权利要求的范围内的所有技术方案。
Claims (8)
1.一种基于组分空间分布的复合材料属性分析方法,其特征在于,包括如下步骤:
S1、根据有限元分析方法构建复合材料的代表性体积单元RVE,设置复合材料的代表性体积单元刚度矩阵为[K];
S2、将复合材料的代表性体积单元的内部结构划分为三种细观组分,包括夹杂、基质、夹杂和基质间界面;
S3、划分RVE的背景网格为细观网格;
S4、根据有限元分析方法设置第i种细观组分所有单元的材料属性矩阵为[mi],设置第i种细观组分中第e个单元出现在空间中位置x的概率矩阵为[pi,(e)(x)],则得到第e个单元的单元刚度矩阵[k(e)],计算表达式为:
[k(e)]=[mi][p(e)(x)];
S5、利用随机分形方法,将[pi,(e)(x)]表达为复合材料的细观组分在三维空间里的分形维数D3和复合材料的细观组分的多重分形测度μ3的函数f(D3,μ3),则得到第e个单元的单元刚度矩阵[k(e)]的计算表达式为:
S6、基于步骤S5得到的第e个单元的单元刚度矩阵计算复合材料的代表性体积单元刚度矩阵[K],计算表达式为:
S7、对步骤S6得到的复合材料的代表性体积单元刚度矩阵施加力边界条件,分析基于组分空间分布的复合材料属性。
2.根据权利要求1所述的一种基于组分空间分布的复合材料属性分析方法,其特征在于,步骤S1中的复合材料包括水泥基复合材料、陶瓷增强铝基复合材料中的一种。
3.根据权利要求1或2所述的一种基于组分空间分布的复合材料属性分析方法,其特征在于,步骤S3的具体实现方法为以dw为横向间距、以dh为竖向间距,在复合材料的代表性体积单元上布置纵横向直线,划分为背景网格。
4.根据权利要求3所述的一种基于组分空间分布的复合材料属性分析方法,其特征在于,步骤S5中利用计盒方法计算复合材料的细观组分在三维空间里的分形维数,具体实现方法包括如下步骤:
S5.1、在在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的边长分别为δ1,δ2,…,δn;
S5.2、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);
S5.3、以δ1,δ2,…,δn为横坐标,以Ni(δ1),Ni(δ2),…,Ni(δn)为纵坐标,在双对数坐标轴上绘制散点,计算函数lnN(δ)和-lnδ的斜率为复合材料第i种细观组分在三维空间里的分形维数。
5.根据权利要求4所述的一种基于组分空间分布的复合材料属性分析方法,其特征在于,步骤S5中利用粗线条方法计算复合材料的细观组分的多重分形测度,具体实现方法包括如下步骤:
S5.4、在复合材料的代表性体积单元上构造一系列边长不同的盒子,设置盒子的边长分别为δ1,δ2,…,δn;
S5.5、判断盒子所含的区域内是否含有指定的复合材料的代表性体积单元的内部结构的细观组分,判断为是则盒子与指定的复合材料的代表性体积单元的内部结构的细观组分相交,统计所有盒子与第i种细观组分相交的个数为Ni(δ1),Ni(δ2),…,Ni(δn);根据相交盒子数目与总盒子数目之比,统计所有盒子与第i种细观组分相交的概率为Pi(δ1),Pi(δ2),…,Pi(δn);
S5.6、取q为一定范围内的实数,设置与组分空间分布概率相关的计算过程函数χq(δ)的计算表达式为:
然后计算复合材料的细观组分的多重分形测度,计算公式为:
6.根据权利要求5所述的一种基于组分空间分布的复合材料属性分析方法,其特征在于,步骤S7的具体实现方法包括如下步骤:
S7.1、对步骤S6得到的复合材料的代表性体积单元刚度矩阵上边界的节点施加力边界条件[F],同时约束下边界节点;
S7.2、由矩阵方程[K][U]=[F]计算位移值[U],进而由力的值与相应位移值之比,计算复合材料的均质化等效材料属性。
7.电子设备,其特征在于,包括存储器和处理器,存储器存储有计算机程序,所述的处理器执行所述计算机程序时实现权利要求1-6任一项所述的一种基于组分空间分布的复合材料属性分析方法的步骤。
8.计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-6任一项所述的一种基于组分空间分布的复合材料属性分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310837525.XA CN116882239B (zh) | 2023-07-10 | 2023-07-10 | 一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310837525.XA CN116882239B (zh) | 2023-07-10 | 2023-07-10 | 一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116882239A true CN116882239A (zh) | 2023-10-13 |
CN116882239B CN116882239B (zh) | 2024-01-23 |
Family
ID=88263747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310837525.XA Active CN116882239B (zh) | 2023-07-10 | 2023-07-10 | 一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116882239B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016177357A1 (de) * | 2015-05-04 | 2016-11-10 | Skl Schwergewebekonfektion Lichtenstein Gmbh | Hybrider, textilverstärkter multi-material-verbundwerkstoff und verfahren zu seiner herstellung |
US20190124314A1 (en) * | 2016-04-14 | 2019-04-25 | National Institute Of Advanced Industrial Science And Technology | Three-dimensional state estimation device, three-dimensional state estimation program, and three-dimensional state estimation method |
CN113239583A (zh) * | 2021-04-19 | 2021-08-10 | 北京理工大学 | 一种复合材料本构模型的构建方法及系统 |
CN113722879A (zh) * | 2021-07-16 | 2021-11-30 | 山东大学 | 一种基于均质化理论的微结构力学性能确定方法与系统 |
CN113761763A (zh) * | 2021-08-06 | 2021-12-07 | 上海索辰信息科技股份有限公司 | 微观和细观尺度结构rve多尺度宏观材料性质分析方法 |
CN116227276A (zh) * | 2023-01-16 | 2023-06-06 | 西安交通大学 | 多重材料分布情况下电场电位多尺度有限元计算方法、设备及存储介质 |
-
2023
- 2023-07-10 CN CN202310837525.XA patent/CN116882239B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016177357A1 (de) * | 2015-05-04 | 2016-11-10 | Skl Schwergewebekonfektion Lichtenstein Gmbh | Hybrider, textilverstärkter multi-material-verbundwerkstoff und verfahren zu seiner herstellung |
US20190124314A1 (en) * | 2016-04-14 | 2019-04-25 | National Institute Of Advanced Industrial Science And Technology | Three-dimensional state estimation device, three-dimensional state estimation program, and three-dimensional state estimation method |
CN113239583A (zh) * | 2021-04-19 | 2021-08-10 | 北京理工大学 | 一种复合材料本构模型的构建方法及系统 |
CN113722879A (zh) * | 2021-07-16 | 2021-11-30 | 山东大学 | 一种基于均质化理论的微结构力学性能确定方法与系统 |
CN113761763A (zh) * | 2021-08-06 | 2021-12-07 | 上海索辰信息科技股份有限公司 | 微观和细观尺度结构rve多尺度宏观材料性质分析方法 |
CN116227276A (zh) * | 2023-01-16 | 2023-06-06 | 西安交通大学 | 多重材料分布情况下电场电位多尺度有限元计算方法、设备及存储介质 |
Non-Patent Citations (4)
Title |
---|
XU YANG 等: "Quantity and shape modification for random-fractal-based 3D concrete meso-simulation", POWDER TECHNOLOGY, vol. 320 * |
李伟;宋卫东;宁建国;: "考虑粒子分布特征的复合材料细观力学方法", 固体力学学报, no. 04 * |
杨旭: ""基于分形和混沌理论混凝土细观分析方法的研究"", 中国博士学位论文全文数据库 工程科技Ⅱ辑 * |
覃立宁;戴自航;周瑞忠;: "基于平均值定理和点积分方案的自然单元法及其程序实现", 计算力学学报, no. 05 * |
Also Published As
Publication number | Publication date |
---|---|
CN116882239B (zh) | 2024-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7379780B2 (en) | Equivalent material constant calculation system, storage medium storing an equivalent material constant calculation program, equivalent material constant calculation method, design system, and structure manufacturing method | |
CN116757146B (zh) | 分布式随机行走寄生电容提取方法、装置、设备及介质 | |
CN116882239B (zh) | 一种基于组分空间分布的复合材料属性分析方法、电子设备及存储介质 | |
DE102018119919A1 (de) | Nicht-flüchtiges computerlesbares Medium, das eine Standard-Zell-Library speichert, integrierter Schaltkreis mit synchronem Schaltkreis, und Computersystem zum Designen des integrierten Schaltkreises | |
CN108122267A (zh) | 一种冗余金属的填充方法及装置 | |
US20240126950A1 (en) | Thermal analysis method, thermal analysis apparatus, and computer program | |
Wheeler | Problems with skewness and kurtosis, Part One | |
CN111158732B (zh) | 访问数据处理方法、装置、计算机设备及存储介质 | |
JP2003347301A (ja) | 均質化法を用いた構造解析方法及び構造設計方法 | |
CN113506600B (zh) | 颗粒均匀性评价方法、装置及设备 | |
CN116070324A (zh) | 一种含预制裂缝的混凝土梁细观模型生成方法及系统 | |
CN114492272A (zh) | Fpga芯片及其顶层网表的搭建方法、装置、设备、存储介质 | |
CN113157806B (zh) | 网格数据分布式存储服务系统、方法、装置、设备及介质 | |
CN115422882A (zh) | 集成电路产品布线区域等效导热系数的批量化计算方法 | |
CN114463477A (zh) | 模型贴图方法、装置和电子设备 | |
US11408730B2 (en) | Stress measuring device and stress measuring method | |
CN111597096A (zh) | 一种基准测试方法、系统及终端设备 | |
DE112019007256T5 (de) | Adaptiver Leistungsmessungsakkumulator mit sich verändernder Abtastfrequenz | |
DE102014116418A1 (de) | Bereitstellung von, zumindest teilweise, mindestens einer Indikation, dass zumindest ein Teil der Daten zur Verarbeitung verfügbar sind | |
CN115062439A (zh) | 简化模型构建方法、装置、设备及可读存储介质 | |
CN116385571B (zh) | 一种基于多维化动态可变分辨率的点云压缩方法及系统 | |
Dutta et al. | Some studies on distribution parameters for facilities design problems | |
CN103823789B (zh) | 一种低复杂度的通用混合基fft设计方法 | |
CN113777547B (zh) | 在片s参数测量系统校准判断方法、装置及终端 | |
CN108573012B (zh) | 一种数据处理方法、装置、设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |