基于神經(jīng)網(wǎng)絡(luò)算法的中藥制劑處方流動性預(yù)測模型研究＊

杜若飛，趙立杰，馮 怡

（上海中醫(yī)藥大學(xué) 中藥現(xiàn)代制劑技術(shù)教育部工程研究中心 上海 201203）

基于神經(jīng)網(wǎng)絡(luò)算法的中藥制劑處方流動性預(yù)測模型研究＊

杜若飛＊＊，趙立杰，馮 怡

（上海中醫(yī)藥大學(xué) 中藥現(xiàn)代制劑技術(shù)教育部工程研究中心 上海 201203）

目的：通過建立混合過程中物料流動性變化規(guī)律的神經(jīng)網(wǎng)絡(luò)模型，預(yù)測任意以微晶纖維素為稀釋劑的中藥處方流動性。方法：以五倍子、金錢草浸膏粉為模型藥物，將其與4種不同型號的微晶纖維素按不同比例均勻混合，分別測定混合前后各粉體的松密度、振實密度，并計算豪斯納比率和壓縮度以表征流動性，在此基礎(chǔ)上運用神經(jīng)網(wǎng)絡(luò)算法建立混合過程中流動性變化規(guī)律的數(shù)學(xué)模型，進而預(yù)測任意以微晶纖維素為稀釋劑的中藥處方流動性。結(jié)果：通過預(yù)測結(jié)果與實驗結(jié)果的比較發(fā)現(xiàn)，神經(jīng)網(wǎng)絡(luò)預(yù)測模型具有較高的擬合度和良好的預(yù)測效果。

中藥處方 微晶纖維素 松密度 振實密度 流動性 神經(jīng)網(wǎng)絡(luò)

在中藥固體制劑的生產(chǎn)中，制劑原料的流動性對產(chǎn)品質(zhì)量有重要影響[1]，如流動性差的物料，在膠囊劑填充時造成裝量差異；在干法制粒時，原料的流動性差易造成胚片密度不均，從而降低藥物的穩(wěn)定性和有效性[2，3]。因此，表征制劑原料的流動性對篩選處方，保證產(chǎn)品質(zhì)量和穩(wěn)定性具有指導(dǎo)意義。然而，表征流動性需要經(jīng)過篩分、混合、測量松密度和振實密度等操作過程，篩選制劑處方時每選擇一種輔料都需要重復(fù)一次上述過程，不僅費時費力，而且會造成原料的浪費。如果建立一個能夠預(yù)測處方流動性的數(shù)學(xué)模型，對添加不同種類、不同比例輔料制劑處方的流動性進行準確預(yù)測，則可以快速篩選出最優(yōu)處方，不僅可以節(jié)省原輔料，同時也可大大縮短處方藥的研發(fā)時間。這種通過建立數(shù)學(xué)模型預(yù)測和設(shè)計藥物處方和工藝的思路與國際上流行的QbD理念完全吻合。在該方向上，國內(nèi)上海中醫(yī)藥大學(xué)的馮怡教授的研究團隊在中藥固體制劑的各個領(lǐng)域均進行了較為系統(tǒng)的研究[4，5]，如探討中藥吸濕性與制劑工藝相關(guān)性的研究[6-8]，中藥物理屬性與干法制粒[9-11]、壓片工藝相關(guān)性[12-17]等，并初步建立了可運行的專家系統(tǒng)，為中藥制劑實現(xiàn)QbD提供了新思路。本文以五倍子和金錢草浸膏粉為制劑原料，以4種常用的微晶纖維素為稀釋劑，測量其按照不同比例混合后的流動性數(shù)據(jù)，運用神經(jīng)網(wǎng)絡(luò)算法建立松密度、振實密度變化規(guī)律的數(shù)學(xué)模型，以此為基礎(chǔ)計算豪斯納比率和壓縮度，從而預(yù)測混合粉體的流動性。

1 實驗材料

1.1 藥品與試劑

五倍子提取物（寧波立華制藥有限公司，批號：20120501）；金錢草提取物（寧波立華制藥有限公司，批號：20120201）；AvicelDG（美國FMC公司惠贈）；CeolusKG802（日本旭化成公司惠贈）；CeolusPH301（日本旭化成公司惠贈）；MCCWJ101（安徽山河藥用輔料有限公司，批號：111201）。

1.2 儀器設(shè)備

FA1004N電子天平（上海精密科學(xué)儀器有限公司）；RPM mini2三維混合機（德國J. Engelsmann AG公司）；BT-1000型粉體綜合特性測試儀（丹東市百特儀器有限公司）；R2013a版MATLAB/SIMULINK軟件（美國MathWorks公司）。

2 實驗方法

2.1 原料混合

采用三維混合機將五倍子提取物、金錢草提取物與4種微晶纖維素分別按照質(zhì)量比5、10、20、40、60、80%的比例混合，混合轉(zhuǎn)速為40 rpm，混合時間為10 min，共得到48個處方。

2.2 松密度（ρb）的測定

預(yù)先精密稱定密度金屬容器（容積用100 cm3）的質(zhì)量，計算質(zhì)量G0，再通過振動篩加入樣品，當(dāng)粉末充滿金屬密度容器后停止加料，用刮板將多余的粉體刮出，并用毛刷清除容器外部粉末，精密稱定，計算總質(zhì)量G1，按如下公式計算：

2.3 振實密度（ρt）的測定

通過振動篩加入樣品，開啟震動5 min至容器內(nèi)樣品高度不發(fā)生變化，精密稱定總質(zhì)量G2。按如下公式計算：

2.4 豪斯納比率（HR）

豪斯納比率定義為粉體振實密度和松密度之比，反映了粉體壓縮性和流動性，同樣通過大小來評價粉體流動性的方法，即粉體越大，可壓縮性越強，流動性越差。按如下公式計算：

2.5 壓縮度（Cp）的計算

Cp使用BT-1 000型粉體綜合特性測試儀測試，計算松密度ρb與振實密度ρt。計算公式如下：

3 實驗結(jié)果

表1 中藥浸膏粉和微晶纖維素的ρb、ρt、HR及Cp（±s，n=3）

表1 中藥浸膏粉和微晶纖維素的ρb、ρt、HR及Cp（±s，n=3）

制劑原料ρb/g·cm-3ρt/g·cm-3HR Cp/%五倍子浸膏粉0.442 6±0.010 3 0.750 3±0.000 3 1.70 41.01金錢草浸膏粉0.454 5±0.002 5 0.709 2±0.001 4 1.56 35.91 AvicelDG 0.354 6±0.001 80.523 7±0.001 7 1.48 32.29 CeolusKG802 0.254 2±0.008 0 0.394 0±0.002 6 1.55 35.48 CeolusPH301 0.409 7±0.001 4 0.586 5±0.002 0 1.43 30.15 MCCWJ101 0.447 1±0.003 8 0.678 4±0.002 5 1.52 34.09

表2 不同比例AvicelDG與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n= 3）

表2 不同比例AvicelDG與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n= 3）

五倍子+AvicelDG ρb/g·cm-3ρt/g·cm-3HR Cp/%五倍子+5%AvicelDG 0.481 5±0.002 8 0.797 6±0.003 6 1.66 39.63五倍子+10%AvicelDG 0.484 5±0.003 8 0.801 0±0.003 8 1.65 39.52五倍子+20%AvicelDG 0.463 7±0.013 40.769 5±0.007 5 1.66 39.74五倍子+40%AvicelDG 0.456 1±.007 3 0.734 8±0.001 9 1.61 37.92五倍子+60%AvicelDG 0.429 0±0.010 4 0.664 3±0.002 1 1.55 35.42五倍子+80%AvicelDG 0.393 3±0.004 9 0.592 2±0.001 3 1.51 33.58

表3 不同比例CeolusKG802與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

表3 不同比例CeolusKG802與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

五倍子+CeolusKG802 ρb/g·cm-3ρt/g·cm-3HR Cp/%五倍子+5%CeolusKG802 0.479 2±0.008 2 0.804 2±0.001 2 1.68 40.41五倍子+10%CeolusKG802 0.463 6±0.011 4 0.781 3±0.000 9 1.69 40.66五倍子+20%CeolusKG8020.4433±0.0009 0.743 5±0.000 9 1.68 40.38五倍子+40%CeolusKG802 0.400 1±0.005 3 0.653 5±0.003 4 1.63 38.77五倍子+60%CeolusKG802 0.354 7±0.001 6 0.555 2±0.002 9 1.57 36.12五倍子+80%CeolusKG802 0.300 9±0.002 0 0.463 7±0.002 2 1.54 35.10

表4 不同比例CeolusPH301與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

表4 不同比例CeolusPH301與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

五倍子+CeolusPH301 ρb/g·cm-3ρt/g·cm-3HR Cp/%五倍子+5%CeolusPH301 0.491 3±0.006 8 0.823 1±0.001 5 1.68 40.31五倍子+10%CeolusPH301 0.497 5±0.005 4 0.821 4±0.001 6 1.65 39.44五倍子+20%CeolusPH3010.494 7±0.004 6 0.814 2±0.001 7 1.65 39.24五倍子+40%CeolusPH301 0.498 3±0.004 9 0.795 0±0.000 2 1.60 37.32五倍子+60%CeolusPH301 0.486 4±0.000 9 0.753 6±0.001 2 1.55 35.46五倍子+80%CeolusPH301 0.473 8±0.005 6 0.680 8±0.000 4 1.44 30.40

表5 不同比例的MCCWJ101與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

表5 不同比例的MCCWJ101與五倍子浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

五倍子+MCCWJ101 ρb/g·cm-3ρt/g·cm-3HR Cp/%五倍子+5%MCCWJ101 0.485 8±0.000 5 0.820 8±0.000 6 1.69 40.81五倍子+10%MCCWJ101 0.476 2±0.007 6 0.802 2±0.004 3 1.68 40.64五倍子+20%MCCWJ101 0.467 6±0.007 8 0.785 4±0.005 4 1.68 40.47五倍子+40%MCCWJ101 0.449 3±0.005 1 0.741 0±0.002 9 1.65 39.37五倍子+60%MCCWJ101 0.410 5±0.003 9 0.660 7±0.000 2 1.61 37.87五倍子+80%MCCWJ101 0.368 3±0.001 0 0.583 5±0.001 1 1.58 36.87

表6 不同比例的AvicelDG與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

表6 不同比例的AvicelDG與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

金錢草+AVICELDG ρb/g·cm-3ρt/g·cm-3HR Cp/%金錢草+5%AvicelDG 0.451 9±0.008 3 0.704 7±0.009 2 1.56 35.88金錢草+10%AvicelDG 0.444 4±0.000 7 0.692 2±0.000 9 1.56 35.79金錢草+20%AvicelDG0.443 0±0.005 3 0.674 8±0.002 3 1.52 34.35金錢草+40%AvicelDG 0.418 3±0.000 3 0.634 6±0.001 5 1.52 34.09金錢草+60%AvicelDG 0.402 2±0.003 2 0.596 8±0.002 0 1.48 32.60金錢草+80%AvicelDG 0.374 3±0.000 6 0.553 2±0.002 4 1.48 32.34

表7 不同比例的CeolusKG802與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

表7 不同比例的CeolusKG802與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

金錢草+CeolusKG802 ρb/g·cm-3ρt/g·cm-3HR Cp/%金錢草+5%CeolusKG802 0.443 8±0.000 2 0.692 9±0.000 7 1.56 35.96金錢草+10%CeolusKG802 0.434 4±0.004 1 0.664 6±0.000 9 1.53 34.64金錢草+20%CeolusKG8020.409 3±0.003 8 0.625 9±0.001 6 1.53 34.60金錢草+40%CeolusKG802 0.363 6±0.000 2 0.554 1±0.000 1 1.52 34.38金錢草+60%CeolusKG802 0.320 7±0.000 4 0.486 9±0.003 0 1.52 34.13金錢草+80%CeolusKG802 0.279 0±0.000 2 0.423 3±0.001 6 1.52 34.10

表8 不同比例的CeolusPH301與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

表8 不同比例的CeolusPH301與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp（±s，n=3）

金錢草+CeolusPH301 ρb/g·cm-3ρt/g·cm-3HR Cp/%金錢草+5%CeolusPH301 0.457 8±0.001 1 0.697 0±0.002 7 1.52 34.32金錢草+10%CeolusPH301 0.461 1±0.001 1 0.695 7±0.003 6 1.51 33.72金錢草+20%CeolusPH3010.458 3±0.005 6 0.684 4±0.002 0 1.49 33.04金錢草+40%CeolusPH301 0.450 7±0.006 6 0.661 8±0.001 4 1.47 31.90金錢草+60%CeolusPH301 0.438 1±0.002 4 0.625 2±0.001 9 1.43 29.93金錢草+80%CeolusPH301 0.421 1±0.002 6 0.596 4±0.001 9 1.42 29.40

表9 不同比例的MCCWJ101與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

表9 不同比例的MCCWJ101與金錢草浸膏粉混合粉體的ρb、ρt、HR及Cp （±s，n=3）

金錢草+MCCWJ101 ρb/g·cm-3ρt/g·cm-3HR Cp/%金錢草+5%MCCWJ101 0.455 8±0.003 0 0.703 0±0.001 2 1.54 35.17金錢草+10%MCCWJ101 0.447 7±0.002 3 0.690 5±0.000 5 1.54 35.16金錢草+20%MCCWJ1010.441 2±0.001 8 0.676 5±0.000 3 1.53 34.78金錢草+40%MCCWJ101 0.417 1±0.001 1 0.638 5±0.0005 1.53 34.68金錢草+60%MCCWJ101 0.391 6±0.001 1 0.599 1±0.000 1 1.53 34.64金錢草+80%MCCWJ101 0.358 1±0.001 7 0.547 5±0.000 5 1.53 34.59

4 數(shù)據(jù)分析

編程前將實驗數(shù)據(jù)保存為SMDZSMDYSD. xlsx格式的Excel文件。文件中按照AvicelDG、CeolusKG802、CeolusPH301、MCCWJ101的順序?qū)⑺泻?、松密度、振實密度的?shù)據(jù)依次排列。運用的建模軟件為Matlab 2013，建模方法為BP神經(jīng)網(wǎng)絡(luò)算法。經(jīng)系統(tǒng)編程后，由計算機完成神經(jīng)網(wǎng)絡(luò)建模。如圖1所示，神經(jīng)網(wǎng)絡(luò)在預(yù)測松密度的過程中設(shè)定輸入層為3組參數(shù)，分別為中藥浸膏粉的松密度、輔料所占百分比、混合粉體的松密度，經(jīng)過多次嘗試后隱藏層節(jié)點數(shù)設(shè)為15，輸出層為混合粉體的松密度。如圖5所示，在預(yù)測振實密度的過程中設(shè)定輸入層也為3組參數(shù)，分別為中藥浸膏粉的振實密度、輔料所占百分比、混合粉體的振實密度，經(jīng)過多次嘗試后隱藏層節(jié)點數(shù)設(shè)為20，輸出層為混合粉體的振實密度。

建模結(jié)果如下：

圖2、6為神經(jīng)網(wǎng)絡(luò)模型的回歸分析，包括訓(xùn)練、測試、驗證以及所有數(shù)據(jù)的回歸分析，結(jié)果顯示：松密度預(yù)測模型的相關(guān)性均大于0.97，振實密度預(yù)測模型的相關(guān)性均大于0.96，相關(guān)性非常好。

圖3、7是誤差直方圖，是實驗值和預(yù)測值的差值，從圖中可以看出：大部分數(shù)據(jù)分布在0誤差周圍，這說明預(yù)測值和實驗值之間的差值較小。

圖4中紅色線條為松密度神經(jīng)網(wǎng)絡(luò)模型的預(yù)測值，藍色點代表松密度的實測值，圖8中紅色線條為振實密度神經(jīng)網(wǎng)絡(luò)模型的預(yù)測值，藍色三角代表振實密度的實測值，結(jié)果顯示：松密度、振實密度的預(yù)測值與真實值具有較高重合度，說明兩個神經(jīng)網(wǎng)絡(luò)模型擬合效果良好。

4.1 神經(jīng)網(wǎng)絡(luò)模型預(yù)測松密度的結(jié)果

4.2 神經(jīng)網(wǎng)絡(luò)模型預(yù)測振實密度的結(jié)果

圖1 松密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型的網(wǎng)絡(luò)結(jié)構(gòu)

圖2 松密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型的回歸分析

圖3 松密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型誤差直方圖

圖4 神經(jīng)網(wǎng)絡(luò)預(yù)測的松密度與實驗松密度比較圖

圖5 振實密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型的網(wǎng)絡(luò)結(jié)構(gòu)

5 討論

根據(jù)模型擬合所得到的松密度與振實密度，即可計算處方的壓縮度與豪斯納比率，從而評價處方的流動性。應(yīng)用該模型預(yù)測新的中藥制劑處方的流動性，只需輸入中藥制劑原料（中藥浸膏粉）的松密度與振實密度，設(shè)定由低到高一系列目標(biāo)載藥量，即可得到以“AvicelDG、CeolusKG802、CeolusPH301、MCCWJ101”4種微晶纖維素為充填劑的處方流動性，進而篩選最佳處方。但是由于本研究所選用的輔料僅限于微晶纖維素類充填劑，若要實現(xiàn)預(yù)測中藥浸膏粉和任意輔料混合物的松密度、振實密度、豪斯納比率及壓縮度，則還需要選擇多種其他類型輔料進一步增加數(shù)據(jù)量，構(gòu)建新的模型。

運用人工神經(jīng)網(wǎng)絡(luò)建立復(fù)雜數(shù)學(xué)模型具有方法簡單、擬合度高的優(yōu)勢，但該方法也同樣存在一定缺陷，即神經(jīng)網(wǎng)絡(luò)模型無法提供具體解析式。如果能夠在今后的研究中以該模型為基礎(chǔ)，建立微分方程，給出具體解析式，并賦予解析式中各系數(shù)的物理意義，將會使擬合結(jié)果更有說服力。

圖6 振實密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型的回歸分析

圖7 振實密度神經(jīng)網(wǎng)絡(luò)預(yù)測模型誤差直方圖

圖8 神經(jīng)網(wǎng)絡(luò)預(yù)測的振實密度與實驗振實密度比較圖

本研究將神經(jīng)網(wǎng)絡(luò)算法和計算機編程模擬引入中藥處方篩選過程，通過神經(jīng)網(wǎng)絡(luò)算法搭建了實驗數(shù)據(jù)之間的相關(guān)性，建立了基于流動性的中藥充填劑篩選模型，初步實現(xiàn)了合理化、智能化篩選、設(shè)計中藥處方。本研究不僅實現(xiàn)了學(xué)科交叉，而且避免了傳統(tǒng)工藝中輔料篩選的盲目性，為中藥制劑現(xiàn)代化開辟了一個新思路，同時也為中藥制劑實現(xiàn)QbD起到一定的示范作用。

1 崔福德.藥劑學(xué)（5版）. 北京:人民衛(wèi)生出版社, 2003: 298.

2 曹韓韓,杜若飛,馮怡,等.干法制粒技術(shù)在中藥研究者的應(yīng)用進展.中草藥, 2013, 44(19): 2772-2776.

3 Miguélez-Morán A M, Wu C Y, Seville J P. The effect of lubrication on density distributions of roller compacted ribbons. Int J Pharm, 2008, 362(1-2): 52-59.

4 馮怡,杜若飛,洪燕龍,等. 關(guān)于構(gòu)建中藥制劑工藝設(shè)計專家系統(tǒng)的思考. 世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化, 2013, 15(1): 25-28.

5 馮怡,洪燕龍,鮮潔晨,等. 基于QbD理念的中藥新藥成型工藝研發(fā)模式的探討.中國中藥雜志, 2014, 39(17): 3404-3408.

6 趙立杰,馮怡,徐德生,等.基于多元數(shù)據(jù)分析的中藥制劑原料吸濕性與其他物理特性的相關(guān)性研究.藥學(xué)學(xué)報, 2012, 47(4): 517-521.

7 趙立杰,馮怡,沈嵐,等.相分散法降低中藥制劑原料吸濕性機理初步探討.世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化, 2011, 13(6): 1057-1060.

8 趙立杰,馮怡,徐德生,等.工藝與處方因素對中藥制劑原料防潮性能的影響.中國中藥雜志, 2009, 34(1): 35-38.

9 曹韓韓,楊嘉寧,杜若飛,等. 基于多元回歸分析干法制粒工藝參數(shù)與顆粒質(zhì)量的相關(guān)性.藥學(xué)學(xué)報, 2014, 49 (3): 406-410.

10 曹韓韓,杜若飛,楊嘉寧,等. 微晶纖維素對五倍子提取物可壓縮性和成型性的影響.中草藥, 2014, 45(8): 1072-1077.

11 曹韓韓,楊嘉寧,杜若飛,等. 中藥粉體的物理性質(zhì)與屈服壓力的相關(guān)性分析.中國實驗方劑學(xué)雜志, 2014, 20(5):14-17.

12 李曉海,趙立杰,阮克萍,等. 應(yīng)用感應(yīng)壓片機研究直接壓片輔料的形變行為和成型性.中國藥學(xué)雜志, 2014, 10: 849-857.

13 李曉海,李美黛,杜焰,等. 應(yīng)用因子分析法評價常用片劑輔料的形變行為.藥學(xué)學(xué)報, 2013, 48(9): 1475-1483.

14 李曉海,李美黛,杜閻,等. 感應(yīng)壓片機Korsch XP1數(shù)據(jù)采集方法研究.中國新藥雜志, 2013, 22(19): 2301-2308.

15 Li Xiao Hai, Zhao LiJie, Ruan Ke Ping. The application of factor analysis to evaluate deforming behaviors of directly compressed powders. Powder Technology, 2013, 247: 47-54.

16 李曉海,趙立杰,馮怡,等. 阮克物理性質(zhì)對微晶纖維素可壓縮性和成型性的影響.中國藥學(xué)雜志, 2013, 48(2): 116-122.

17 李曉海,趙立杰,李美黛,等. 藥物粉體壓縮研究數(shù)學(xué)模型的意義與應(yīng)用特點.中國新藥雜志, 2012, 21(12): 1362-1366.

Abstract:This study aimed to forecast the flowability of formulations in traditional Chinese medicine (TCM) regarding microcrystalline cellulose (MCC) as diluent through establishing the artificial neural networks (ANN) model on disciplines of flowability changes in the process of mixing materials. In this study, the Galla chinensis extract and the Herba lysimachiae extract were taken as model drugs which mixed with four kinds of MCC in different proportions, respectively. Then the bulk density and tap density of the extracts were measured before and after mixture. The hausner ratio and compressibility index were calculated to represent their flowability. On the basis of the preparation above, the ANN model involving the discipline of flowability changes was established during the mixing process to forecast the flowability of any TCM formulation with MCC diluent. As a result, the ANN model involving the discipline of the flowability changes had its share of a higher degree of fitting and better forecasting effectiveness when comparing the forecast results with the actual experimental outcomes .

Determining the Flowability of Formulations in Traditional Chinese Medicine Using the Forecasting Model Established by Artificial Neural Networks

Du Ruofei, Zhao Lijie, Feng Yi
(Engineering Research Center of Modern Preparation Technology of Traditional Chinese Medicine, Ministry of Education, Shanghai University of Traditional Chinese Medicine, Shanghai 201203, China)

Traditional Chinese medicine formulations, microcrystalline cellulose, bulk density, tap density, flow ability, artificial neural networks

10.11842/wst.2016.04.010

R283

A

（責(zé)任編輯：馬雅靜，責(zé)任譯審：朱黎婷）

2015-11-10

修回日期：2015-12-23

＊ 上海市衛(wèi)計委科研項目（ZY3-CCCX-3-5001）：中藥制劑關(guān)鍵技術(shù)研究與應(yīng)用，負責(zé)人：馮怡；上海市科委平臺資助項目（15DZ2292000）：上海市中藥固體制劑專業(yè)技術(shù)服務(wù)平臺，負責(zé)人：馮怡；國家自然科學(xué)基金委XX項目（81403110）：基于人工神經(jīng)網(wǎng)絡(luò)算法的中藥硬膠囊劑設(shè)計原理研究，負責(zé)人：杜若飛；上海市教育委員會預(yù)算內(nèi)項目（2013JW27）：羥丙基甲基纖維素改善中藥噴霧干燥熱熔型黏壁及其機理研究，負責(zé)人：王優(yōu)杰；上海市教育委員會預(yù)算內(nèi)項目（2014YSN23）：固體中藥制劑原料吸濕行為模擬與分析，負責(zé)人：趙立杰。

＊＊ 通訊作者：杜若飛，助理研究員，博士，主要研究方向：中藥制劑關(guān)鍵技術(shù)研究。