1. 引言
多孔材料是一种新型复合材料,在现代生活中广泛存在。此材料在结构、减震、隔热、渗透等多方面发挥着重要的作用。多孔材料存在孔隙的变化,孔隙的存在会使多孔材料具备良好的力学性能,如相对密度低、高比表面积、高冲击韧性、良好的传导性等,多孔材料的应用十分广泛,对其弯曲性能的研究非常重要。
为了解多孔材料梁的力学性能,研究者们运用了多种方法理论 [1] [2] 。在经典力学理论的背景下,研究人员的观点是材料的弹性系数在变形过程中不改变,也就是说,应力–应变本构关系在弹性范围内是线性关系。但事实上许多材料的本构关系并非是线性关系, 多孔材料的本构关系就具有较强的非线性特性。1982年,Ambartsumyan [3] 在研究这种非线性本构关系时,发现多数材料的应力应变关系在受压和受拉时分别接近线性,提出了双模量理论。很多学者利用该理论研究了构件的弯曲变形、屈曲、过屈曲问题以及振动问题。张鹏 [4] 等人利用双模量理论对四边固支的板进行了研究,将所计算的数据与有限元软件对比得出结论,双模量理论计算的数值更加精确。王康健和赵永刚等人 [5] 采用双模量理论对杆的过屈曲进行了分析。Fraternali [6] 等人对曲梁进行了研究,曲梁的后屈曲响应会受到拉压不同模量的强烈影响。穆琳 [7] 等人对简支梁的振动问题进行了研究,采用打靶法求解了弯曲挠度和固定频率。还有许多研究结果表明采用双模量理论的计算精度得到了很大提高 [8] [9] 。此外,还有许多的研究者对新型材料梁板的力学行为进行了研究 [10] [11] [12] 。然而,在材料的本构关系非线性特征较强时,即受拉和受压时的应力应变关系分别线性化显著不合适时,双模量理论的研究结果就还不够精确。
本文思考非线性本构关系是连续变化状态,认为材料的弹性系数是应变的线性函数 [13] ,对多孔梁两端固定边界条件的弯曲问题进行了研究。经过与文献的比较验证了本文方程和数据的准确性,并基于所推导的基本方程研究了孔隙分布类型、孔隙率、非线性本构参数和载荷对多孔梁的弯曲变形的影响。
2. 基本方程
文章对如图1所示的具有非线性应力–应变本构关系的矩形截面多孔梁的弯曲变形问题进行了研究。多孔梁的长、宽、高分别为l、b和h,构件上加载了一个铅垂方向分布的作用载荷q。梁的厚度方向的孔隙分布有三种典型的分布形式,如图2所示。其中,类型1 (T1)分布表示孔隙在梁的厚度方向均匀分布,类型2 (T2)分布表示梁的杨氏模量
在上下表面处最大,向几何中面递减,分布类型3 (T3),其杨氏模量的变化趋势与类型2 (T2)的孔隙分布相反。三种情况下杨氏模量的表达式如下所示 [14] [15] :
类型1 (T1):
(1)
类型2 (T2):
(2)
类型3 (T3):
(3)
上式中,
表示含孔隙的梁的有效杨氏模量,
为无孔(
)时材料的杨氏模量,
是孔隙度系数。
![](//html.hanspub.org/file/1-2650287x15_hanspub.png?20240315081413193)
Figure 1. The geometric structure of multi hole beams
图1. 多孔梁的几何结构
![](//html.hanspub.org/file/1-2650287x18_hanspub.png?20240315081413193)
Figure 2. Three types of porosity distributions for beams
图2. 梁的三种孔隙度分布类型
2.1. 几何方程
基于几何结构,对多孔梁建立笛卡尔坐标系O-xyz (见图1),x、y、z轴分别为梁的轴线,横截面的水平对称轴和铅垂对称轴。多孔梁在经历大挠度弯曲变形时,其几何中面上的轴向应变
为:
(4)
其中:u为多孔梁几何中面上的轴向位移,w多孔梁几何中面上的铅垂位移。关于多孔梁的非几何中面层上的其他处的弯曲变形情况,其应变
为:
(5)
2.2. 物理方程
考虑到材料的本构关系的非线性特征,本文假设材料的弹性模量E与应变是线性函数的关系(线性弹性模量理论),即:
(6)
其中,
和
分别是与多孔梁孔隙分布类型相对应的弹性常数,它们与孔隙率的函数关系具有相同形式,即可写为
。
是反映该材料本构关系非线性特点的常数。则梁弯曲时,截面上的正应力表达式为:
(7)
多孔梁横截面上的轴力FN和弯矩M的表达式为:
(8)
其中,对于T1类型的梁:
;
对于T2类型的梁:
;
对于T3类型的梁:
。
2.3. 静力学方程
当多孔梁上受到了如图1所示的载荷分布时,在大挠度弯曲变形时,对构件的微段平衡进行了研究,得到了内力的平衡方程为:
(9)
两端固定边界条件为:
(10)
3. 对基本方程的无量纲处理
将相关方程以及边界条件进行了无量纲化以达到计算结果的准确性和求解的便利性要求,相关无量纲参数如下所示:
(11)
将相关参数代入(9)~(10)式,得到无量纲内力表达式为:
(12)
无量纲控制方程为:
(13)
无量纲边界条件(固定端)为:在X = 0,1处,
。
4. 弯曲问题的数值分析
利用打靶法对两端固定边界条件下的无量纲方程进行求解。令:
(14)
无量纲化后的基本式子可转化为如下的一阶非线性常微分方程:
(15)
其中:
相应的边界条件(固定端):
。
当三种孔隙分布类型(T1, T2, T3)梁的孔隙率参数
时,多孔梁此时是一种无孔的均匀材料梁。针对大挠度弯曲问题,此时对长细比
,非线性本构参数
不同时的两端固定梁的中点挠度W (0.5)与载荷Q变化关系进行数据计算,如图3所示。将计算数据与文献 [13] 中的相关数据图进行了比较,其变形结果高度吻合,以此验证所推导方程的正确性和计算得到的数据的可靠性,然后研究了两端固定梁在不同孔隙分布类型、物理非线性本构参数、载荷和孔隙率下的大挠度弯曲问题。
![](//html.hanspub.org/file/1-2650287x47_hanspub.png?20240315081413193)
Figure 3. Equilibrium path of large deflection problem of T1, T2, T3 beams (
)
图3. T1,T2,T3梁大挠度问题的平衡路径(
)
图4所示为
,孔隙率
,非线性本构参数Er不同数值下,三种孔隙分布情况的两端固定多孔梁中点挠度W (0.5)与载荷Q变化的平衡路径。观察图可知,Er取值的变化对中点挠度与载荷的非线性特征无影响,因为该模型的几何关系具有非线性特性。并且当Er > 0和Er < 0时,其相关数据是分居在Er = 0这条线的两侧。Er取正值时,曲线在下侧,W (0.5)变小,抗弯曲能力增强;当Er取负值时,曲线在上侧,W (0.5)变大,抗弯曲能力减弱。在处于不同孔隙分布类型下时,三种孔隙分布的梁抗弯能力按T1、T3、T2的顺序依次增强,且T2梁的抗弯曲能力明显强于T1和T3梁。所以在现实工况中,对于多孔材料分布类型的选择,应该尽量选择孔隙度为T2分布形式的梁。
![](//html.hanspub.org/file/1-2650287x52_hanspub.png?20240315081413193)
Figure 4. Equilibrium path of large deflection problem fixed at both ends of T1, T2 and T3 beams
图4. T1,T2,T3梁两端固定大挠度问题的平衡路径
图5所示为
,载荷Q = 130,非线性本构参数Er变化下,大挠度问题下的不同孔隙分布类型的梁在中点挠度W (0.5)随孔隙率
的趋势变化图。观察图可以发现,随着孔隙率
的增加,T1,T2,T3梁中点挠度也是逐渐增大的,但是从变化的缓急程度来看,T2梁变化的最缓和,其次是T3梁,最后是T1梁,即T1梁受孔隙率
的影响最大。从曲线的变化节点来看,大约
> 0.6后,T1梁的中点挠度数值随着孔隙率
的增大而急剧上升,T2、T3梁也有类似的情形,这说明在选择多孔材料时,面对孔隙率的选择,在保持材料刚度能力最优时,孔隙率
的大小最好不要超过0.6。并且,在相同的外部条件下,T2梁的变化趋势远弱于T1,T3梁,这也证实了T2梁的抗弯能力优于另外两种梁。
![](//html.hanspub.org/file/1-2650287x60_hanspub.png?20240315081413193)
Figure 5. Diagram of the large deflection problem W at both ends of T1, T2 and T3 beams with porosity
图5. T1,T2,T3梁两端固定大挠度问题W随孔隙率
的变化图
![](//html.hanspub.org/file/1-2650287x63_hanspub.png?20240315081413193)
Figure 6. Equilibrium configuration of fixed large deflection problem at both ends of T1, T2, T3 beams
图6. T1,T2,T3梁两端固定大挠度问题的平衡构形
图6所示为
,Q = 130,
,非线性本构参数Er变化时,T1,T2,T3梁的大挠度弯曲情况下的弯曲构形。看图可知,梁的弯曲变形程度与Er的变化是呈相反趋势的,对于相同的非线性参数Er时,三种孔隙分布类型梁的抗弯性能均随着Er的增大而减小。而面对不同孔隙分布类型,三种孔隙分布类型的抗弯能力也是按T1,T3,T2逐次增强的,且T2的抗弯能力明显优于另外两种。
图7给出了
,Q = 130,非线性本构参数Er = 0时,两端固定T1,T2,T3梁的大挠度的弯曲变形与孔隙率
的变化关系。看图可知,随着孔隙率的增加,T1,T2,T3梁的弯曲变形逐渐增大,这是因为孔隙率的增加会让梁的孔隙部分越来越多,梁的抗弯能力减弱,因此在受到相同的外载荷的情况下,多孔梁发生的弯曲变形就会随着孔隙率的增加而增加。且对于相同的孔隙率
,T1梁的挠度最大。当三种梁的孔隙率
> 0.5后,孔隙率
的变化对弯曲变形的影响愈加明显。
![](//html.hanspub.org/file/1-2650287x71_hanspub.png?20240315081413193)
Figure 7. Equilibrium configuration of fixed large deflection problem at both ends of T1, T2, T3 beams with respect to the change of porosity
图7. T1,T2,T3梁两端固定大挠度问题的平衡构形关于孔隙率
的变化
5. 结论
本文以非线性本构关系为基础,设多孔材料梁的弹性模量与应变的关系是线性函数,对多孔梁不同孔隙分布类型下的弯曲变形问题的基本方程进行了推导,采用打靶法求解了两端固定多孔梁在大挠度弯曲问题下受横向均布载荷的数值解。基于文中所研究的结果可得到以下结论:
(1) 多孔梁大挠度弯曲变形的数值会随着非线性本构参数Er数值的增大而减小。
(2) 孔隙率
存在的情况下,T1,T2,T3梁的弯曲变形数值随着孔隙率
的增加而增加,其中三种梁的抗弯能力按T2,T3,T1的顺序依次减弱。
(3) 不同的孔隙分布类型下,面对相同的外界条件,T2梁受到外载荷Q、孔隙率
和非线性本构参数Er的变化影响小于T1和T3梁,即T2梁的综合抗弯能力最强。
基金项目
本文得到国家自然科学基金11862012,12062010的资助。