Preparation and properties of microporous scaffolds by 3D printing of bio-based polylactic acid composites
-
摘要:目的
探讨生物基聚乳酸复合材料用于3D打印直接构建含微孔支架的可行性。
方法采用热重分析仪和差示扫描量热仪探究聚乳酸复合材料的热性能,扫描电镜表征支架的微观形貌,活/死细胞染色检测支架的细胞黏附情况。
结果所制备的0.6%ADC-PHAP和40%NaCl-PHAP复合材料具有良好的热稳定性和加工性,适用于熔融沉积3D打印。当压缩应变为80%时,0.6%ADC-PHAP和40%NaCl-PHAP支架相应的压缩应力分别为45.27和52.11 MPa;0.6%ADC-PHAP复合材料的初始分解温度比40%NaCl-PHAP复合材料低19.5 ℃;0.6%ADC-PHAP支架的孔隙率达到63.33%,有利于细胞黏附,且细胞相容性比40%NaCl-PHAP支架更好。
结论生物基聚乳酸复合材料可通过熔融沉积型3D打印直接构建含微孔支架,所制备的0.6%ADC-PHAP支架具有一定的应用潜力。
Abstract:ObjectiveTo explore the feasibility of using bio-based polylactic acid composites to directly construct microporous scaffolds by 3D printing.
MethodThermogravimetric analyzer and differential scanning calorimeter were used to explore the thermal properties of bio-based polylactic acid composites, scanning electron microscopy was used to characterize the microscopic morphology of the bio-scaffold, and the live/dead cell staining was used for detecting cell adhesion of the scaffold.
ResultThe prepared 0.6%ADC-PHAP and 40%NaCl-PHAP composites had good thermal stability and processability, and were suitable for the fused deposition modeling 3D printing process. When the compressive strain was 80%, the corresponding compressive stresses of the 0.6%ADC-PHAP and 40%NaCl-PHAP scaffolds were 45.27 and 52.11 MPa, respectively. The initial decomposition temperature of the 0.6%ADC-PHAP composite was 19.5 ℃ lower than that of the 40%NaCl-PHAP composite. The porosity of the 0.6%ADC-PHAP scaffold reached 63.33% which was conducive to cell adhesion, and the cell compatibility was better than that of the 40%NaCl-PHAP scaffold.
ConclusionThe bio-based polylactic acid composites can be used to directly construct microporous bio-scaffolds through fused deposition modeling 3D printing, and the prepared 0.6%ADC-PHAP bio-scaffold has certain application potential.
-
Keywords:
- Melt blending /
- 3D printing /
- Scaffold /
- Polylactic acid /
- Bio-based material
-
近年来,随着我国农业劳动力的下降,收获高峰期找不到劳动力[1-3],农业采摘机器人成为农业工程领域的研究热点[4-9]。夹持机构作为与农业生物对象直接接触的末端执行器是农业采摘机械装备的关键机构。目前国内外文献报道的夹持机构大多数针对体积较小、质量较轻的串形或单粒水果,如苹果、荔枝、葡萄和草莓等。尹建军等[10]设计的葡萄采摘机构中,用电机通过连杆机构驱动四指,从果实外包络抓取葡萄。苗玉彬等[11]设计的苹果夹持机构中,利用压力传感器测得所需要的夹持力,并将数值输给夹持末端,使得夹持力的数值比较精确。张凯良等[12]设计的草莓夹持器中,先利用摄像头定位,再通过2个光纤传感器探测机械爪中是否存在果柄,并进行切割,在切割的同时机械手对草莓进行夹持,一次性完成草莓的采摘。陈燕等[13]设计了具有4个自由度的结构,根据D-H坐标系建立了荔枝采摘机械手运动学模型并用于采摘荔枝。比利时的Baeten等[14]研制出一种苹果采摘机器人,该机器人的手臂具有6个自由度,由液压升降装置驱动机械手臂作业。澳大利亚团队研制的香蕉采收机,以液压伸缩臂和切割机进行采收,节省了人力物力[15]。日本研制了一种由多自由度的机械手臂和视觉控制系统构建而成的葡萄采收器[16]。
中国是香蕉生产大国,与其他果蔬相比,香蕉属于串形水果,质量和体积大、果柄粗,难以用常规夹持机构夹持。人工采摘时,通常一人肩托着香蕉串、一人用刀砍断果柄[17]。本研究模拟人抱紧香蕉串的动作行为,设计了一种夹持蕉串果柄的末端执行机构,即内掌带钉抱掌夹持装置;通过搭建夹持试验平台,选择影响夹持效果的气缸压力、钉子数目以及钉子排列方式3个主要因素,采用单因素试验和响应曲面法进行了抱掌夹持力不同组合因素的效果测试,取得了较为理想的设计参数和试验数据,为夹持机构的进一步优化设计提供了理论和实践参考。
1. 带钉抱掌夹持机构设计
模拟人工收获时手掌抱握对象的夹持方式[18],本夹持机构设计圆弧形抱掌,通过抱掌内表面设置钉子增大抱掌与蕉柄间的附着力,设计的香蕉果柄夹持机构如图1a所示。电机通过联轴器驱动丝杆转动,螺旋传动带动螺套往复运动,与螺套固连的抱掌在四杆机构驱动下实现张开与闭合。夹持机构推动抱掌闭合时,抱掌在抱紧蕉柄时的受力状态如图1c所示。
机构静力平衡关系如下:
$$f_1 = f_2\sin \beta ,$$ (1) $$f_2\cos \alpha = {f / 2},$$ (2) 式中,f1为钉子对蕉柄的挤压力,f2为连接杆对夹持杆的拉力,
$\alpha $ 为连接杆和丝杆之间的夹角,$\beta $ 为连接杆和夹持杆之间的夹角,f为夹持机构对蕉柄的挤压力。由式(1)和式(2)可知,夹持机构对蕉柄的挤压力(f)是影响抱掌与蕉柄之间静摩擦力的主要因素,f是选择电机功率和螺旋机构设计的重要参数。为了轻量化夹持机构并获得可靠的夹持力,需要研究夹持因素对夹持力的影响规律,以获得优化的机构组合参数。
2. 试验方法
2.1 夹持试验
2.1.1 夹持力分析
夹持机构由螺旋机构和四杆机构组成。夹持过程中,抱掌通过螺旋机构带动四杆机构运动,获取作用在蕉柄上的压力,当蕉串从蕉柄上切割分离后,蕉串的质量(m)由抱掌与蕉柄之间产生的摩擦力(
$f' $ )承担(图1b),抱掌夹持状态下与蕉柄截面受力示意图如图1d。由图1d可知,夹持机构作用在抱掌上的挤压力(f1)、抱掌内表面钉子的排列方式和钉子数量是影响抱掌和蕉柄之间静摩擦力的主要因素,通过试验测试夹持各因素的影响规律及夹持效果,获取最佳组合参数。2.1.2 试验平台及试验材料
试验平台如图2a所示,主要包括:万能试验机、抱掌、夹持元件、动力系统。抱掌由半圆形夹套和钉子(直径5 mm)构成,夹持元件由支撑平台、挡板和支柱组成,动力系统由空气压缩机、MGPM40-50Z气缸、手动换向阀、气源两联件和气管等组成。试验材料选用落梳后的果柄,切割成长度约为150 mm的蕉柄段,方便放入抱掌中进行夹持试验。
图 2 试验平台和原理图1:万能试验机,2:STC-100 kg拉力传感器,3:手动换向阀,4:钉子,5:半圆形夹套,6:调速器,7:支撑平台,8:MGPM40-50Z气缸,9:挡板Figure 2. Schematic diagrams of experimental platform and principle1: Universal testing machine, 2: STC-100 kg tension sensor, 3: Manual reversing valve, 4: Nail, 5: Semicircular jacket, 6: Governor, 7: Support platform, 8: MGPM40-50Z cylinder, 9: Baffle2.1.3 试验过程
试验场景如图2b。测试过程中,由空气压缩机提供动力,经过气源两联件的减压阀调节输入气缸的气压大小,手动换向阀可调节气缸前进、后退和暂停3个档位。各元件安装好后,打开手动换向阀,抱掌在两侧气缸的推动下夹紧蕉柄,然后启动万能试验机。根据拉拔试验原理向上拉动蕉柄,模拟蕉柄切断后香蕉串的重力作用,测试抱掌对蕉柄的夹持效果,得出拉动过程中拉力和位移的数据信息,分析得到各因素对夹持效果的影响。
2.2 单因素试验
选取刚采收下来的香蕉果柄,裁截为便于夹持的蕉段,香蕉柄外形近似为圆柱形,蕉柄直径为60~70 mm,成熟度为7~8成熟。为获得气缸压强(即挤压力)、钉子数量和钉子排列方式对香蕉柄的夹持影响规律,分别对这3个因素进行单因素试验,每组试验分别进行3次拉拔试验,取平均值作为试验结果进行分析。试验取香蕉串质量30 kg为载荷参照,安全系数取1.5,摩擦夹持情况下,夹持机构和蕉柄的最大静摩擦力为500 N左右,可保持夹持机构能有效夹持蕉串。
单因素试验设计见表1。为测试气缸压强对夹持效果的影响规律,选取钉子的排列方式为X型、钉子数量为10个、气缸直径为40 mm,拉拔速度5 mm/min。在气缸压强对夹持的影响研究中,确定了气缸压强至少0.40 MPa可以稳定夹持,为此,控制钉子数量为10个、气缸压强0.40 MPa、拉拔速度5 mm/min;参考加工的方便性,分别选择钉子的排列方式为X型、矩形和三角形,测试钉子的最优排列方式;选择气缸压强0.40 MPa和矩形排列方式不变,钉子数量分别为4、6、8、10、12和14个,测试钉子数量的影响效果。
表 1 气缸压强、钉子数量和排列方式单因素试验设计Table 1. Single factor test designs of cylinder pressure, number and arrangement of nails因素
Factor试验编号
Test number气缸压强/MPa
Cylinder
pressure钉子数量/个
Number of
nails排列方式
Nail
arrangement气缸压强
Cylinder pressure1 0.30 10 X型 2 0.35 10 X型 3 0.40 10 X型 4 0.45 10 X型 5 0.50 10 X型 排列方式
Nail arrangement1 0.40 10 X型 2 0.40 10 三角形 3 0.40 10 矩形 钉子数量
Number of nails1 0.40 4 矩形 2 0.40 6 矩形 3 0.40 8 矩形 4 0.40 10 矩形 5 0.40 12 矩形 6 0.40 14 矩形 2.3 响应曲面法试验
响应曲面法的设计模型主要是中心复合设计(Central composite design,CCD)[19-20]。为验证3个因素的组合夹持效果,本文采用CCD中的三因素三水平试验进行试验设计,测试气缸压强(x1)、钉子排列方式(x2)和钉子数量(x3)对夹持效果的影响,测试指标为拉力(即夹持机构承受的重力,y),表2给出了因素编码及试验设计水平。利用Design-Expert 8.0.6软件对试验数据进行分析,采用线性和二次多项式拟合得出线性回归方程,确定最优的数学模型。
表 2 因素编码及试验水平1)Table 2. Factor coding and test level编码 Coding x1 x2 x3 −1 0.3 X型 4 0 0.4 矩形 10 1 0.5 三角形 16 1) x1:气缸压强,MPa;x2:钉子排列方式;x3:钉子数量
1) x1: Cylinder pressure, MPa; x2: Arrangement of nails; x3: Number of nails3. 结果与分析
3.1 单因素试验结果
3.1.1 气缸压强对夹持的影响
拉拔过程中蕉柄的滑移受力曲线如图3a。由图3a可以看出,拉力从零点至最高点的曲线段表示试验过程中蕉柄受力压缩阶段,此时拉力并未达到蕉柄和抱掌之间的最大静摩擦力;最高点为蕉柄和夹抱掌之间的最大静摩擦力点;最高点之后的曲线段表示蕉柄在拉力作用下平缓滑动,由于钉子的滑动使蕉柄表面破裂,蕉柄和抱掌之间的滑动摩擦力减小,蕉柄出现滑移,夹持失效。试验曲线表明,当气缸压强为0.35 MPa、拉力还未达到500 N时,蕉柄在抱掌中产生滑移的现象,夹持失效;当气缸压强为0.40 MPa、拉力为500 N左右时,最接近蕉柄和抱掌之间的最大静摩擦力,可以稳定夹持;当气缸压强为0.45或0.50 MPa、拉力达到600或800 N时,蕉柄在抱掌中出现滑移现象。
上述结果表明,气缸压强约为0.4 MPa时,夹持状态最佳。查阅MGPM40−50气缸手册,气缸压强为0.4 MPa时,气缸输出压力为500 N,即挤压力(f1)=500 N[21]。
3.1.2 钉子排列方式对夹持的影响
图3b为不同钉子排列方式下,拉拔过程中香蕉柄的滑移受力曲线。从图3b可以看出,零点至最高点的曲线表示拉拔过程中蕉柄的受力压缩阶段,最高点之后是蕉柄在拉力作用下的滑动阶段。在气缸压强0.40 MPa、钉子数量为10个时,矩形排列方式拉力最大,最高点拉力为800 N,说明在滑动位移相同时,矩形可承受的摩擦力最大且最大值(fmax)为800 N,大于蕉柄和抱掌之间的fmax,因此,选择矩形的排列方式可以获得稳定的夹持效果。
3.1.3 钉子数量对夹持的影响
图3c是蕉柄在不同钉子数量情况下可承受的最大拉力。由图3c可知,钉子数量分别为6、8和10个时,可承受的拉力均为800 N左右,大于蕉柄和抱掌之间的fmax,说明这3个水平均可获得稳定的夹持。钉子数量为8和10个时获得的拉力比较接近,并未因钉子数量增加而拉力增加,其原因可能是香蕉柄不是规则的圆柱形,导致圆形抱掌上的钉子并未全部扎进蕉柄内,钉子数量为8或10个时,只有6个钉子扎进蕉柄进行有效夹持。钉子数量超过10个后,拉力不升反降,出现该现象的原因可能是钉子越多,出现互相干扰的现象越严重,成为有效钉子的数量不升反降。试验结束后,观察蕉柄表面的滑痕,进一步印证了该结论。
3.2 响应曲面法的结果与方差分析
用Design-Expert8.0.6软件对试验因素进行分析,结果见表3。对表3中的数据进行回归拟合,保留显著性因素,对拉力进行二次回归方程方差分析,结果见表4。表4的结果表明,模型的P=0.0002<0.01,达极显著水平,拟合系数(R2)为0.9702;拉力失拟项的P=0.109 0>0.05,不显著,这主要是因为蕉柄的表面不规则导致钉子有部分没有刺入,造成了试验误差。由于R2大于0.80,根据响应曲面法拟合系数的有效性参数原则,评价该系数有效,说明本试验采用的数学模型具有较好的拟合度。在设计的因素水平范围内,由各因素的F值大小得出各因素对拉力的影响,由表4可见,排列方式的二次项有较好的拟合度,各因素对拉力的影响大小依次为:排列方式的二次项效应(x22)>气缸压强(x1)>气缸压强与钉子数量的相互作用(x1x3)>钉子排列方式(x2)>钉子数量与排列方式的相互作用(x2x3)>钉子数量(x3)>气缸压强与排列方式的相互作用(x1x2)>钉子数量的二次效应(x32)>气缸压强的二次效应(x12)。
表 3 响应曲面法的试验结果1)Table 3. Experimental result by the response surface method运行顺序
Running order试验序号
Test numberx1 x2 x3 f 5 1 0.3 0 4 564 6 2 0.5 0 4 582 7 3 0.3 0 16 387 4 4 0.5 1 10 401 3 5 0.3 1 10 374 17 6 0.4 0 10 637 13 7 0.4 0 10 590 2 8 0.5 −1 10 411 1 9 0.3 −1 10 190 15 10 0.4 0 10 605 8 11 0.5 0 16 888 16 12 0.4 0 10 617 11 13 0.4 −1 16 307 10 14 0.4 1 4 336 9 15 0.4 −1 4 349 14 16 0.4 0 10 560 12 17 0.4 1 16 556 1) x1:气缸压强,MPa;x2:钉子排列方式;x3:钉子数量;f:拉力,N
1) x1: Cylinder pressure, MPa; x2: Nail arrangement; x3: Number of nails; f: Tension, N表 4 拉力方差分析Table 4. Variance analysis of tension方差来源1)
Source of variance平方和
Square sumF P 模型 Model 435 000 25.36 0.0002 ${x_1}$ 73 536 38.62 0.0004 ${x_2}$ 21 012 11.04 0.0127 ${x_3}$ 11 781 6.19 0.0417 ${x_1}{x_2}$ 9 409 4.94 0.0616 ${x_1}{x_3}$ 58 322 30.63 0.0009 ${x_2}{x_3}$ 17 161 9.01 0.0199 $x_1^2$ 1 646 0.86 0.3833 $x_2^2$ 239 000 125.29 <0.0001 $x_3^2$ 2 271 1.19 0.3109 残差 Residual 13 327 失拟项 Lack of fit 9 961 3.94 0.1090 纯误差 Pure error 3367 总离差 Total deviation 448 000 1) x1:气缸压强,MPa;x2:钉子排列方式;x3:钉子数量
1) x1: Cylinder pressure, MPa; x2: Nail arrangement; x3: Number of nails3.3 两因素相互作用对拉力的影响
两因素相互作用对于拉力的影响和单因素的影响大致相同(图4)。钉子数量为10个时,气缸压强与钉子排列方式相互作用对拉力影响的3D响应曲面和等高线图见图4a1和4a2。由图4a可以看出,气缸压强对拉力的影响更大,排列方式次之。在排列方式为矩形时,随着气缸压强的增大,拉力增大;在气缸压强一定时,排列方式对拉力的影响是先增大后减小。
气缸压强为0.4 MPa时,钉子数量与排列方式相互作用对拉力影响的3D响应曲面和等高线图见图4b1和4b2。由图4b可以看出,钉子数量对拉力的影响更大,排列方式次之。在排列方式一定时,钉子数量增多,拉力增大;在钉子数量一定时,排列方式对拉力的影响是先增大后减小。
钉子排列方式为矩形时,气缸压强与钉子数量相互作用对拉力影响的3D响应曲面和等高线图见图4c1和4c2。由图4c可以看出,钉子数量和气缸压强对拉力影响均较大。在气缸压强一定时,随着钉子数量增多,拉力增大;在钉子数量一定时,随着气缸压强增大,压力增大。
3.4 参数优化与验证结果
根据以上的分析,得到最优的参数组合,利用Design-Expert软件的目标优化法,根据以下优化方程式:
$$\left\{ \begin{array}{l} 0.3 \leqslant {x_1} \leqslant 0.5 \\ - 1 \leqslant {x_2} \leqslant 1 \\ 4 \leqslant {x_3} \leqslant 16 \\ y = f({x_1},{x_2},{x_3},y) \\ {\rm{Max}} \; y \end{array} \right.,$$ (3) 得出预测和实测的优化值,结果见表5。由表5可知,优化条件下的拉力预测值和实测值误差为4.6%,该误差在考虑范围内,即优化结果具有可参考性。拉力预测值和实测值的散点值差别如图5所示,预测拟合度高。表明气缸压强0.4 MPa、钉子数量10个、排列方式为矩形是抱掌夹持机构设计最佳组合。
表 5 优化条件下预测和实测的优化值Table 5. Predicted and actual optimization values under optimized conditions项目
Item挤压力/MPa
Extrusion钉子排列方式
Nail arrangement钉子数量/个
Number of nails拉力/N
Tension最佳条件(预测) Optimum condition (forecast) 0.38 0.22 10.27 500 修改条件(实测) Modify condition (actual) 0.40 0 10 524 4. 结论
对于香蕉这类大质量串形水果的机械化采摘,可靠稳定的夹持方式是夹持机构设计的关键。本文通过搭建试验平台,进行夹持影响因素的单因素和正交试验测试,分析得如下结论:1)针对本次夹持装置中的研究对象,影响夹持效果的主次因素顺序为:气缸压强>钉子排列方式>钉子数量,最优组合为气缸压强0.4 MPa、钉子数量10个、钉子排列方式为矩形,此时,实测拉力524 N,大于最大静摩擦力,可以稳定夹持;2)在进行钉子数量的最优化试验过程中发现,当钉子数量超过一定数目时,实测拉力直线下降,甚至低于最大静摩擦力,原因是钉子数量过多,使蕉柄横截面破坏过大,蕉柄部分断裂;3)对于香蕉这种体积大的果蔬收获,传统摩擦式夹持机构难以保证稳定的夹持,为了保证香蕉从采摘到运输至指定地点的平稳,采用钉刺夹持方式具有可行性。
-
表 1 PLA复合材料热降解数值统计
Table 1 Numerical statistics of thermal degradation of PLA composites
样品 Sample θ0/℃ θmax1(θmax2)/℃ θend/℃ 残余量/% Residue PHAP 344.4 363.4 377.5 6.03 10%NaCl-PHAP 336.1 355.4(399.9) 374.9 17.26 40%NaCl-PHAP 328.5 348.1(397.8) 368.5 37.96 0.2%ADC-PHAP 310.5 332.9(390.0) 367.4 10.72 0.6%ADC-PHAP 309.0 336.7(390.3) 373.4 5.82 表 2 PLA复合材料的DSC测试结果
Table 2 DSC test results of PLA composites
样品
Sample玻璃化温度(θg)/℃
Glass transition temperature熔融温度(θm)/℃
Melting temperature熔融焓(ΔΗm)/(J·g−1)
Melting enthalpy结晶度(Xc)/%
CrystallinityPHAP 56.23 167.53 27.06 37.75 10%NaCl-PHAP 56.04 167.66 26.44 40.98 40%NaCl-PHAP 56.79 167.50 22.28 51.80 0.2%ADC-PHAP 60.89 167.77 22.43 31.35 0.6%ADC-PHAP 56.16 167.89 21.69 30.44 -
[1] OSTAFINSKA A, FORTELNÝ I, HODAN J, et al. Strong synergistic effects in PLA/PCL blends: Impact of PLA matrix viscosity[J]. Journal of the Mechanical Behavior of Biomedical Materials, 2017, 69(5): 229-241.
[2] 宋中波, 甄卫军. PCL增韧PLA材料的流变行为及其热降解动力学[J]. 塑料, 2017, 46(4): 32-38. [3] BAI Z F, DOU Q. Rheology, morphology, crystallization behaviors, mechanical and thermal Properties of poly(lactic acid)/polypropylene/maleic anhydride-grafted polypropylene blends[J]. Journal of Polymers and the Environment, 2018, 26(3): 959-969. doi: 10.1007/s10924-017-1006-5
[4] 程思敏, 陈丽杰, 洪阳阳, 等. 羟基磷灰石的表面改性及其对聚乳酸基多孔支架性能的影响[J]. 复合材料学报, 2018, 35(5): 1087-1094. [5] ZHOU C, YANG K, WANG K, et al. Combination of fused deposition modeling and gas foaming technique to fabricated hierarchical macro/microporous polymer scaffolds[J]. Materials & Design, 2016, 109(21): 415-424.
[6] SCOTT G D, KILGOUR D M. The density of random close packing of spheres[J]. Journal of Physics D: Applied Physics, 1969, 2(6): 863-866. doi: 10.1088/0022-3727/2/6/311
[7] ZHANG J, XIAO D, HE X, et al. A novel porous bioceramic scaffold by accumulating hydroxyapatite spheres for large bone tissue engineering: III: Characterization of porous structure[J]. Materials Science & Engineering : C, 2018, 89(8): 223-229.
[8] KAKUMANU V, SRINIVAS SUNDARRAM S. Dual pore network polymer foams for biomedical applications via combined solid state foaming and additive manufacturing[J]. Materials Letters, 2018, 213(4): 366-369.
[9] SONG P, ZHOU C, FAN H, et al. Novel 3D porous biocomposite scaffolds fabricated by fused deposition modeling and gas foaming combined technology[J]. Composites Part B:Engineering, 2018, 152(21): 151-159.
[10] LU T, LI Y,CHEN T. Techniques for fabrication and construction of three-dimensional scaffolds for tissue engineering[J]. International Journal of Nanomedicine, 2013, 8(1): 337-350.
[11] HAIDER A, HAIDER S, KUMMARA M R, et al. Advances in the scaffolds fabrication techniques using biocompatible polymers and their biomedical application: A technical and statistical review[J]. Journal of Saudi Chemical Society, 2020, 24(2): 186-215. doi: 10.1016/j.jscs.2020.01.002
[12] LI X, ZHANG S J, ZHANG X, et al. Biocompatibility and physicochemical characteristics of poly(Ɛ-caprolactone)/poly(lactide-co-glycolide)/nano-hydroxyapatite composite scaffolds for bone tissue engineering[J]. Materials & Design, 2017, 114(2): 149-160.
[13] PATI F, SONG T H, RIJAL G, et al. Ornamenting 3D printed scaffolds with cell-laid extracellular matrix for bone tissue regeneration[J]. Biomaterials, 2015, 37(2): 230-241.
[14] FERRI J M, JORDÁ J, MONTANES N, et al. Manufacturing and characterization of poly(lactic acid) composites with hydroxyapatite[J]. Journal of Thermoplastic Composite Materials, 2017, 31(7): 865-881.
[15] ZHAO H, ZHAO G. Mechanical and thermal properties of conventional and microcellular injection molded poly (lactic acid)/poly (ε-caprolactone) blends[J]. Journal of the Mechanical Behavior of Biomedical Materials, 2016, 53(1): 59-67.
[16] HO M P, LAU K T, WANG H, et al. Improvement on the properties of polylactic acid (PLA) using bamboo charcoal particles[J]. Composites Part B:Engineering, 2015, 81(14): 14-25.
[17] ABBASI H, ANTUNES M, VELASCO J I. Graphene nanoplatelets-reinforced polyetherimide foams prepared by water vapor-induced phase separation[J]. Express Polymer Letters, 2015, 9(5): 412-423. doi: 10.3144/expresspolymlett.2015.40
[18] ABU HASSAN N A, AHMAD S, CHEN R S, et al. Cells analyses, mechanical and thermal stability of extruded polylactic acid/kenaf bio-composite foams[J]. Construction and Building Materials, 2020, 240(11): 117884.
[19] FRACKOWIAK S, LUDWICZAK J, LELUK K, et al. Foamed poly(lactic acid) composites with carbonaceous fillers for electromagnetic shielding[J]. Materials and Design, 2015, 65(1): 749-756.
[20] 张学盈, 崔永岩. AC发泡剂与增塑剂对PVC发泡材料性能的影响[J]. 塑料, 2016, 45(1): 32-34.