CN115859715A - 一种周期性多孔结构材料等效性能的计算方法 - Google Patents
一种周期性多孔结构材料等效性能的计算方法 Download PDFInfo
- Publication number
- CN115859715A CN115859715A CN202211475546.3A CN202211475546A CN115859715A CN 115859715 A CN115859715 A CN 115859715A CN 202211475546 A CN202211475546 A CN 202211475546A CN 115859715 A CN115859715 A CN 115859715A
- Authority
- CN
- China
- Prior art keywords
- unit cell
- cell model
- porous structure
- structure material
- periodic
- 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.)
- Pending
Links
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于多孔材料的等效力学性能分析与计算技术领域,公开了一种周期性多孔结构材料等效性能的计算方法,通过代表性胞元模型的建立、周期性边界条件、位移载荷的施加、计算宏观应力四个步骤,可有效实现不同孔隙率且孔分布形式为正六边形的材料等效弹塑性力学性能的求解,不仅计算精度高,而且能给出材料在细观尺度下的具体细节,有效分析孔隙的形状、体积分数对材料等效性能的影响规律,揭示材料变形和破坏的物理机制。
Description
技术领域
本发明属于多孔材料的等效力学性能分析与计算技术领域,具体来说涉及一种周期性多孔结构材料等效性能的计算方法。
背景技术
多孔材料具有质轻、多功能化应用等特点而被广泛应用于各种领域。如离子推力器中的栅极,其在结构上呈现周期性多孔薄壁型设计以满足减重、工作稳定性等方面的综合要求。目前,多孔材料等效性能的研究多限于弹性性能,而多孔材料在服役过程中通常会呈现出明显的弹塑性行为,同时经典的解析法虽能给出联系宏、细观量的具体列式,但难以揭示材料变形和破坏的物理机制。因此,利用有限元法研究多孔材料的宏观弹塑性力学行为具有重要意义。
但是由于多孔材料细观结构的特殊性,其周期性单元结构与宏观尺度相比存在数量级的差异,导致难以直接利用常规的有限元方法对其宏观弹塑性力学性能进行计算。
发明内容
本发明的目的是提供一种周期性多孔结构材料等效性能的计算方法,可有效实现不同孔隙率且孔分布形式为正六边形的材料等效弹塑性力学性能的求解,不仅计算精度高,而且能给出材料在细观尺度下的具体细节,有效分析孔隙的形状、体积分数对材料等效性能的影响规律。
本发明采取的技术方案是周期性多孔结构材料等效性能的计算方法,包括以下步骤:
步骤一、矩形单胞模型的建立:将孔洞进行填充处理,划分网格,赋予材料属性;
步骤二、周期性边界条件的施加:在矩形单胞模型平行相对面上相应网格节点处建立线性约束方程;
进一步地,所述步骤一中,周期性多孔结构材料由基体和孔洞两部分组成,填充材料的弹性模量比基体材料小10-3数量级。
进一步地,所述材料属性包括弹性模量、屈服强度、真实应力应变。
进一步地,所述步骤二中,矩形单胞平行相对面上的节点位移满足如下关系:
式中,ui、vi分别表示单胞模型上任意一点的横坐标和纵坐标,ui‘、vi‘是指单胞平行相对面上对应的节点横坐标和纵坐标;ua,va、ub,vb、ud,vd分别表示单胞模型上顶点a、b、d处的横坐标和纵坐标。
式中H为单胞模型宽度,uH和uL分别为沿单胞模型宽度和长度方向的真实位移。
进一步地,所述步骤四中根据均匀化方法得到的公式为:
本发明的有益效果在于:
1、本方法提出了将孔洞看作特殊的夹杂相,并进行填充处理。建立代表性胞元,通过施加简单的周期性和位移边界条件,根据均匀化方法,就可计算得出宏观应力;
2、该方法能给出材料在细观尺度下的具体细节,有效分析孔隙的形状、体积分数对材料等效性能的影响规律,揭示材料变形和破坏的物理机制;
3、该方法经过实例验证,计算精度高,为其他形状复杂的周期性多孔结构等效性能的求解提供了可靠的参考和借鉴。
附图说明
图1为周期性多孔结构材料等效性能的计算方法步骤流程图;
图2为实验得到的不锈钢平板试样的真实应力应变曲线;
图3中(a)为矩形单胞模型,(b)为施加周期性边界条件后的单胞模型;(c)为宏观应变为0.07时的有限元单胞应力分布云图;
图4为矩形单胞模型网格节点示意图;
图5为试验得到的周期性多孔不锈钢的应力应变曲线与模拟值的对比。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的周期性多孔结构材料等效性能的计算方法,参见图1至图4,包括单胞模型的建立、周期性边界条件的施加、边界约束条件的施加、宏观应力的计算四个步骤。该方法从具有周期性多孔结构中选取矩形单胞模型,模型包括基体和孔洞(夹杂)。周期性边界条件的施加是通过在单胞平行相对面上施加线性约束方程实现。再在单胞模型上施加位移载荷,求解模型,根据均匀化方法,求得宏观应力。
步骤一、矩形单胞模型的建立:周期性多孔结构材料由基体和孔洞两部分组成,将孔洞当作等效细观力学模型中的夹杂相,并进行填充处理,即设置填充材料的弹性模量比基体材料小10-3数量级。划分网格,赋予材料属性;
其中,材料属性包括弹性模量、屈服强度、真实应力应变。
步骤二、周期性边界条件的施加:参看图4,在矩形单胞模型平行相对面上相应网格节点处建立线性约束方程,即满足如下关系:
式中,ui、vi分别表示单胞模型上任意一点的横坐标和纵坐标,ui‘、Vi‘是指单胞平行相对面上对应的节点横坐标和纵坐标;ua,va分别表示单胞模型上顶点a处的横坐标和纵坐标;ub,vb分别表示单胞模型上顶点b处的横坐标和纵坐标;ud,vd分别表示单胞模型上顶点d处的横坐标和纵坐标。
步骤三中边界约束条件的施加的具体过程为:
在利用细观力学有限元法研究材料宏观力学行为时,应使得单胞应变能与相应的匀质材料单胞应变能相等。
均匀位移边界载荷条件表示为:
均匀力边界载荷条件表示为:
在匀质等效体内,总应变能为:
在单胞模型体内,总应变能为:
为便于描述,建立一原点固定在模型的笛卡尔坐标系。在所建立的笛卡尔坐标系中,沿单胞模型宽度和长度方向的真实位移可表示为:
式中,H、L分别为单胞模型的宽度和长度。
步骤四、宏观应力计算:求解矩形单胞模型,得到单胞模型的细观应力场,根据均匀化方法,公式表达如下:
采用如下算例,对本方法进行验证说明:
制备不锈钢平板拉伸试样并进行拉伸实验,得到弹性模量、屈服强度及真实应力应变曲线等参数,具体数值见图2和表1:
表1不锈钢基体材料力学性能
制备具有周期性多孔结构的不锈钢拉伸试样并进行DIC拉伸实验(DIC,DigitalImage Correlation,数字图像相关法,是一种测量物体表面应变和变形的方法),得到真实应力应变曲线见图5所示。将图2和表1数据输入有限元模拟软件中,建立矩形单胞有限元模型,施加周期性边界条件和位移载荷进行计算,得到多孔不锈钢真实应力应变曲线的模拟值见图5所示。
由上面的算例可知,应力应变曲线的模拟值与实验值在弹性阶段和加工硬化阶段都具有很好的相符性,总的来说吻合良好。所以该计算方法不仅计算精度高,而且能给出材料在细观尺度下的具体细节,有效分析孔隙的形状、体积分数对材料等效性能的影响规律,揭示材料变形和破坏的物理机制。
Claims (6)
2.根据权利要求1所述的周期性多孔结构材料等效性能的计算方法,其特征在于,所述步骤一中,周期性多孔结构材料由基体和孔洞组成,填充材料的弹性模量比基体材料小10-3数量级。
3.根据权利要求1所述的周期性多孔结构材料等效性能的计算方法,其特征在于,所述材料属性包括弹性模量、屈服强度、真实应力应变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211475546.3A CN115859715A (zh) | 2022-11-23 | 2022-11-23 | 一种周期性多孔结构材料等效性能的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211475546.3A CN115859715A (zh) | 2022-11-23 | 2022-11-23 | 一种周期性多孔结构材料等效性能的计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115859715A true CN115859715A (zh) | 2023-03-28 |
Family
ID=85665422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211475546.3A Pending CN115859715A (zh) | 2022-11-23 | 2022-11-23 | 一种周期性多孔结构材料等效性能的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115859715A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116738707A (zh) * | 2023-06-08 | 2023-09-12 | 西安交通大学 | 部分周期性换热器通道的等效力学性能预测方法及*** |
-
2022
- 2022-11-23 CN CN202211475546.3A patent/CN115859715A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116738707A (zh) * | 2023-06-08 | 2023-09-12 | 西安交通大学 | 部分周期性换热器通道的等效力学性能预测方法及*** |
CN116738707B (zh) * | 2023-06-08 | 2024-05-17 | 西安交通大学 | 部分周期性换热器通道的等效力学性能预测方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Karamooz Ravari et al. | A computationally efficient modeling approach for predicting mechanical behavior of cellular lattice structures | |
Camas et al. | Numerical study of the thickness transition in bi-dimensional specimen cracks | |
Camas et al. | Crack front curvature: Influence and effects on the crack tip fields in bi-dimensional specimens | |
Vizotto | Computational generation of free-form shells in architectural design and civil engineering | |
Francesconi et al. | An investigation of the enhanced fatigue performance of low-porosity auxetic metamaterials | |
Zghal et al. | Buckling responses of porous structural components with gradient power-based and sigmoid material variations under different types of compression loads | |
Zhuang et al. | Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow | |
CN110688789B (zh) | 一种碳纤维复合单向板偏载冲击的多尺度响应分析方法 | |
CN114970260B (zh) | 一种用于模拟复合材料破坏的格构相场方法 | |
Njim et al. | Analytical and numerical free vibration analysis of porous functionally graded materials (FGPMs) sandwich plate using Rayleigh-Ritz method | |
Lignola et al. | Exact stiffness–matrix of two nodes timoshenko beam on elastic medium. an analogy with eringen model of nonlocal euler–bernoulli nanobeams | |
CN115859715A (zh) | 一种周期性多孔结构材料等效性能的计算方法 | |
Deng et al. | An improved interpolating complex variable meshless method for bending problem of Kirchhoff plates | |
Patni et al. | Efficient modelling of beam-like structures with general non-prismatic, curved geometry | |
Cai et al. | Nonlinear stability analysis of hybrid grid shells | |
CN110321571B (zh) | 一种蜂窝板壳结构的力学参数数值提取方法 | |
Zhang et al. | Efficient isogeometric formulation for vibration analysis of complex spatial beam structures | |
CN106503385B (zh) | 一种点阵夹层材料等效弹性模量的计算方法 | |
Zengin et al. | Investigation of force and moment effect of St 37 and St 70 roof lattice steels in Ansys program | |
CN109252441B (zh) | 一种变截面箱梁剪力滞效应的分析方法 | |
CN115906333A (zh) | 一种桁架结构的几何非线性等效板动力学建模与响应分析方法 | |
Wang et al. | The study on bending performance of aluminum alloy honeycomb panel-beam composite grid structure | |
Park et al. | Equivalent static wind loads vs. database-assisted design of tall buildings: An assessment | |
CN111353229B (zh) | 一种实体结构光滑粒子动力学建模方法 | |
CN103399974A (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 |