第一章 引言
药物释放技术属于现代药剂学的关键领域，它的主要目的就是借助科学的载体设计，准确地控制药物进入体内的速度和部位，进而保持稳定的血药浓度，并且减小毒副作用。传统制剂在实际使用中常常会遇到释放突释或者释放不完全的问题，不能满足复杂的临床治疗需要。基于拓扑优化的药物释放模型，就是针对以上问题提出的一种新的设计方法。该技术把数学计算模型同材料工程学原理融合起来，在给定的设计区域内，按照预先设定的药物释放目标，智能算法会自动找到并改善材料的微观分布结构，剔除多余的材料，重新塑造出拥有最佳流体通道和扩散路径的药物载体结构。该过程可以明显提高药物释放的可控性、重现性，在保证药效的基础上减少载体材料的用量，有重要的科研价值和临床应用前景[1]。
第二章 拓扑优化与药物释放模型的基础理论及研究进展
2.1 药物控释系统的核心作用机制与现有构型局限
药物控释系统的关键之处在于借助物理化学手段来控制药物分子由载体内部往外部体液环境迁移的速度。主流机制主要有三种，分别是扩散控释依靠浓度梯度推动药物通过多孔骨架或者膜层释放出来，适合于水溶性药物，溶蚀控释是聚合物材料在体内逐渐降解而释放出药物，常用来制备长效制剂，渗透压控释则是利用半透膜内外的渗透压差把药物以恒定的速率泵出，适合于对释药精度要求很高的治疗场景。但是目前常规构型如实心载药块、规则通孔骨架和微球阵列等存在着明显的不足。实心结构不能维持长周期的零级释药，容易因为表面积减小而使速率衰减，规则通孔骨架在追求高载药量的时候常常会牺牲结构稳定性，释放曲线的调整空间也比较小，微球阵列虽然可以实现局部释放，但是很难准确地满足复杂的个体化释药需求。这些固有的缺陷说明传统的设计方法已经达到了性能的天花板，必须用新的设计思想来突破现有的限制。
表1 药物控释系统典型作用机制与现有构型局限对比
2.2 面向可控释药的拓扑优化核心算法逻辑与研究现状
面向可控释药的拓扑优化技术，就是用数学算法在一定的设计空间里寻找材料的最佳分布，从而达到对药物释放动力学精确控制的目的。SIMP（固体各向同性材料惩罚模型）是目前应用最广的一个框架，它的基本思想就是用惩罚因子来惩罚材料的中间密度，使设计变量趋近于0或者1，从而得到清晰的多孔释药结构。对于药物释放场景，该算法要兼顾药物传质效率和结构力学支撑两方面的要求，构造出包含扩散方程和结构刚度方程的耦合目标函数。算法根据灵敏度信息不断地对材料分布进行更新，直到满足设定的释放速率和力学强度收敛条件为止。回顾近十年国内外研究，虽然已经取得了很大的进展，但是仍然存在着明显的不足，即现有的研究目标函数对于特定的释药曲线拟合精度不高，多物理场耦合计算在处理复杂的生物环境时存在误差，并且缺少系统的体外实验和仿真结果的完整性检验。因此本文主要从修正目标函数的针对性、提高耦合场计算精度和完善仿真校验体系三个方面来达到技术上的创新。
第三章 基于拓扑优化的药物释放构型模型构建与仿真验证
3.1 药物释放过程的多物理场耦合控制方程搭建
图 1 药物释放过程多物理场耦合控制方程逻辑
药物释放过程是复杂的多物理场相互作用，本节主要建立扩散、溶蚀和流体流动的控制方程组。根据Fick定律描述药物分子在多孔骨架中的扩散行为，得到有效扩散系数与骨架孔隙率、曲折度的关系；根据一级降解动力学推导载体材料溶蚀方程，得到溶蚀速率常数和环境依赖性；根据Darcy定律建立渗透压驱动下流体流动方程，得到流体黏度和渗透率参数。核心难点就是确定物理场之间的耦合关系，骨架局部密度的实时变化会同步改变药物扩散系数和流体渗透系数，从而达到质量传递和骨架结构演化的动态耦合。另外根据释放环境来设定浓度梯度和流体压力边界条件，保证数值求解的稳定性以及准确性，为之后拓扑优化运算和仿真验证打下坚实的数理基础。
3.2 面向零级释药目标的拓扑优化目标函数与约束条件设置
为了达到零级释药的目的，首先要把药物释放速率的恒定性当作主要的性能指标，确定单位时间内释放量的波动幅度控制阈值以及预设周期内总释放率的要求。构建目标函数时，以最小化不同时间节点的实际释药速率与目标恒定速率的偏差为核心，数学表达为
，其中 为 时刻的释放速率， 为目标速率。然后设置多个约束条件，即最低力学刚度约束来保证植入时结构的完整性，最大孔隙率约束来保证骨架的可制造性，局部材料密度最小阈值来防止出现非物理的完全镂空区域。同时用SIMP插值模型对密度相关的传质系数和弹性模量做正则化处理，保证优化过程收敛到具有良好释药性能和工程可行性的多孔构型上。
3.3 优化后多孔释药骨架结构的药物释放行为仿真校验
3.3.1 优化骨架的释药速率随时间演化规律分析
用仿真全周期输出的数据，提取出各个时间节点上优化多孔释药骨架的累计释药率和瞬时释药速率，绘制出释药速率随时间变化的动态演化曲线，分阶段分析释药初期、中期和末期的速率波动情况。对瞬时释药速率和目标零级速率的偏差百分比进行量化统计，就可以准确地评价模型的释药精度。同时结合不同时刻骨架内部药物浓度分布云图，对释药孔隙随着传质过程不断扩大的动态作用机理进行详细的拆解释药孔隙的动态作用机理，明确拓扑优化构型通过改变孔隙分布来调节扩散阻力的途径，阐明拓扑优化构型能长时间保持稳定释药速率的原因，给以后的制剂开发提供理论基础。
3.3.2 与传统规则孔道骨架的释药均一性对比验证
选择临床常用的等直径平行规则孔道释药骨架作为对照样本，在保证两组样本总载药量、整体外形尺寸和骨架基体材料属性完全一致的情况下，同时进行相同的边界条件下药物释放仿真实验。校验过程中主要检测释药初期的突释量、全周期释药速率的标准差和各个批次样本释药曲线的重合度，以此来作为量化评价释药均一性的指标。用直观的对照数据和误差分析可知，拓扑优化得到的多孔释药骨架可以有效地减小药物释放的波动，提高零级释药的维持能力，说明该模型对于改善药物释放均一性和控释精度有实际应用的价值。
第四章 结论
本文利用拓扑优化的方法建立药物释放模型，对药物递送系统的设计进行了系统的验证。核心原理就是采用变密度法，在给定的设计域内用数学算法自动寻找最优的材料分布路径，从而形成具有特定释放动力学的载体微观结构。实现过程包含建立有限元模型、定义目标函数和灵敏度分析等主要环节，实现了由理论设计到结构优化的闭环。实际使用中发现，该模型可以很好地控制药物的扩散速度，提高药物载量的同时又能很好地抑制突释现象，使释放曲线更加平稳。以拓扑优化为基础的设计方法可以克服传统经验设计的不足，给缓控释制剂的结构创新提供标准化、可量化的方法，有重要的临床应用价值和学术推广意义。
参考文献
\[1\]高扬, 孟繁舒, 闵浩楠. PCL/ASA药物控释系统结构拓扑优化及瞬态释药性能[J]. 工程塑料应用, 2024, 52(2): 87-93.
\[2\]机械. 水凝胶/静电纺丝膜复合敷料的结构设计与性能研究[D]. 2023.