大学生物医学工程专业三年级：肾单位滤过动力学数学建模与病理仿真教案

大学生物医学工程专业三年级：肾单位滤过动力学数学建模与病理仿真教案

一、课程基本信息

（一）课程名称：肾单位滤过动力学数学建模与病理仿真

（二）授课对象：大学本科三年级生物医学工程专业（医用仪器与系统方向、生物力学与建模方向）

（三）学时安排：总学时4学时（每学时50分钟），其中课堂讲授与互动研讨2学时，计算机仿真实验2学时（连贯进行）

（四）教学环境与资源：具备交互式白板的智慧教室；每位学生配备安装有MATLAB/SimulinkR2022b及以上版本或Python3.10+（Anaconda发行版）的笔记本电脑；课程专用云端仿真代码仓库；人体生理学虚拟仿真平台（肯德尔泌尿系统模块）；数字化切片库（肾小球超微结构）。

二、教学目标

（一）知识与技能目标

1.精准阐述肾单位滤过功能的核心解剖生理学基础，包括滤过膜的三层结构、滤过屏障的电荷与尺寸选择性、肾小球滤过率（GFR）的定义及生理波动范围。【基础】

2.运用Starling力的概念，完整推导肾小球滤过率经典公式GFR=Kf×(Pgc–Pbc–πgc+πbc)，并对其中各静水压、胶体渗透压、滤过系数Kf进行量纲分析与生理意义阐释。【非常重要】【高频考点】

3.识别常见肾小球疾病（如抗基底膜肾炎、糖尿病肾病、局灶节段性肾小球硬化）导致GFR改变的病理生理环节，并能将病理损伤映射至数学模型的参数扰动上。【难点】【热点】

4.构建并调试基于房室概念的肾单位滤过动力学集中参数模型，通过编写MATLAB脚本或Simulink框图，对正常生理状态及至少两种病理状态（如滤过膜面积减少、滤过膜通透性下降）进行数值仿真，输出GFR、滤过分数、肾小管流量等关键指标的动态变化曲线。【非常重要】

5.对比仿真结果与临床文献报道的典型数据，开展初步的模型验证与误差分析，并能从生物医学工程角度提出模型优化的可能路径。

（二）过程与方法目标

1.通过“生理机制→物理规律→数学方程→计算仿真”的转化链条，经历工程视角下复杂生命系统建模的全过程，强化跨学科类比思维（电路与流体网络类比）。【重要】

2.在小组协作完成病理参数扰动实验中，体验控制变量法、敏感性分析等系统辨识基本方法，能够基于仿真异常结果反推生理环节的故障定位。

3.利用云端代码仓库进行模型的版本管理与互评，初步养成开放科学（OpenScience）与可重复性研究的职业习惯。

（三）情感态度与价值观目标

1.感悟定量生理学对于精准医学的基石作用，树立将工程技术服务于生命健康的价值取向。

2.通过对肾脏代偿功能与失代偿临界点的仿真模拟，深化对机体稳态机制复杂精妙的敬畏，增强对慢性肾病患者早期干预重要性的社会责任感。

三、教学重点与难点

（一）教学重点

1.肾小球滤过动力学Starling公式的物理意义及其在GFR建模中的核心地位。【非常重要】【高频考点】

2.从集总参数房室模型到常微分方程组的状态空间描述，并利用数值积分求解动态过程。

3.病理生理变化与模型参数扰动之间的映射关系（例如：足细胞损伤→Kf降低；出球小动脉收缩→Pgc升高）。

（二）教学难点

1.理解沿肾小球毛细血管血流途径，胶体渗透压πgc并非定值，而是从入球端到出球端逐渐升高，从而导致有效滤过压净力沿程下降。经典集总参数模型如何简化处理这一空间分布效应，以及这种简化带来的误差与适用边界。【难点】

2.将抽象的病理学描述（如免疫复合物沉积、基底膜增厚、电荷屏障丢失）量化为具体的模型参数变化范围，并理解参数变化的高度非线性耦合。

3.仿真结果的临床转译：如何从模型预测的GFR下降趋势区分急性肾损伤与慢性肾脏病的不同代偿模式。

四、教学方法与策略

（一）教法：采用基于问题链的启发式讲授与翻转课堂相结合。核心推进路径为“临床悖论引出问题—生理机制还原—物理抽象—数学表达—计算实验验证—临床回归”。全程融合数字化切片与实时手绘板书，动态呈现多尺度关联。

（二）学法：采用拼图式协作学习（JigsawMethod）。将全班分为六个基础组，每个组内负责四个专家模块之一：模块A（超微结构与流体力学）、模块B（电类比与RC网络）、模块C（病理参数映射）、模块D（MATLAB/SimPy仿真实现）。专家交叉研讨后返回基础组完成完整模型建构与汇报。

五、教学实施过程（核心环节，完整呈现四个课时的无缝衔接）

（一）课前预置阶段（学习通平台发布，课前48小时完成）

1.微课视频1（12分钟）：回顾肾脏大体解剖与肾单位微观结构，重点强调入球小动脉、出球小动脉、肾小球毛细血管网、鲍曼氏囊的毗邻关系。视频末尾设置选择题：肾小球滤过膜损伤时，下列哪项最先出现异常？A.静水压B.胶体渗透压C.滤过系数KfD.血浆流量。【基础】【高频考点】学生需提交答案并生成词云，教师根据错误率调整课中讲解权重。

2.阅读材料：GuytonHall医学生理学第14版第26章“肾小球滤过率、肾血流量及其调控”节选（英文原版+中文核心段落翻译对照）。要求学生用思维导图工具绘制GFR调控的三层网络：局部反馈、神经调节、激素调节。此思维导图将在课中进行二次深化。

3.代码预置：云端仓库已提供空白的Simulink模型框架（.slx文件）和一份带有函数占位符的Python脚本model_placeholder.py。学生需提前完成开发环境配置并运行测试文件test_env.py，输出“KidenymodelingenvironmentOK”。此项作为形成性评价的签到依据。

（二）课中第一学时：生理机制向数学模型的第一次抽象——从结构到变量

1.临床案例引爆（7分钟）：展示一位32岁男性患者的临床资料——尿常规示大量蛋白尿（+++），血肌酐132μmol/L，eGFR56mL/min/1.73m²。同时展示同一患者肾活检电镜图片：弥漫性足细胞足突融合。教师设问：足细胞损伤，肉眼看不到摸不着，工程师如何用数学语言描述它对肾脏功能的影响？学生迅速进入问题情境，动机被激发。

2.板书重构滤过屏障（13分钟）：教师手绘（或利用数字笔在白板上绘制）肾小球滤过膜三层超微结构图，并同步投射电镜原图。依次引出内皮细胞窗孔（直径70-100nm）、基底膜（富含负电荷蛋白聚糖）、足细胞裂孔（直径约40nm）。指出决定滤过性能的三个独立子因素：有效孔隙面积、等效厚度、电荷密度。此处进行重要等级标注：【基础】【非常重要】——此三因素为后续Kf参数分解的核心。

3.物理量提取与命名（15分钟）：引导学生将上述结构特征转化为可测或可估的物理量。师生对话形成如下映射：有效滤过压差ΔP→实际驱动压力；滤过膜总有效通透面积×单位面积水力传导系数→集总参数Kf。教师强调：Kf并不是一个基本物理常数，而是强烈依赖于生理状态的集总参数，其量纲为mL/(min·mmHg)，数值约为12.5mL/(min·mmHg)正常成人双肾。此处引出Kf的两种组分：Kf=k×S，其中k为单位面积水力传导系数（反映膜孔隙率与厚度），S为总滤过面积（约0.8㎡）。【重要】【高频考点】

4.构建基本稳态方程（15分钟）：基于Starling假设，列出肾小球毛细血管段某一点的有效滤过压Puf=(Pgc–Pbc)–(πgc–πbc)。由于πbc在鲍曼囊中极低且接近恒定（约0mmHg），πbc常忽略或设为固定小值。关键在于πgc的沿程变化——教师指出，这是经典模型与工程学生直觉最容易冲突之处。教师并不直接给出简化，而是提出思考题：若πgc从入球端28mmHg升至出球端36mmHg，那么整个毛细血管网的平均有效滤过压应当如何计算？学生短暂讨论后，教师引出“算术平均胶渗压”近似法，并给出临床常用简化公式：GFR=Kf×(ΔP–πgc平均)。此环节标注【难点】，并要求学生在笔记中用五角星标出。

（三）课中第二学时：从静态方程到动态房室模型——状态空间的建立

1.类比迁移（8分钟）：回顾电路理论中的一阶RC电路充放电过程。电荷q(t)对应于肾小囊中滤液容积V(t)；电阻R对应于滤过阻力（1/Kf）；电容C对应于肾小囊顺应性（ΔV/ΔPbc）。教师瞬时画出类比图，指出这种类比并非完美但极具启发性，使学生快速接受用微分方程描述肾小囊容积变化的思维。【重要】

2.房室划分与流率平衡（17分钟）：严格定义两个房室：房室1——肾小球毛细血管网（视为压力源，压力Pgc受入球/出球小动脉张力调节，本课暂时简化为常数值50mmHg）；房室2——鲍曼氏囊（容积Vbc，压力Pbc是Vbc的函数，与近端小管流出阻力有关）。根据质量守恒：dVbc/dt=GFR–Qout。其中Qout为流入近端小管的流量，可建模为Pbc与近端小管阻力Rpt的函数：Qout=(Pbc–Ppt0)/Rpt，Ppt0为近端小管远端参考压（约10mmHg）。此为本课的核心状态方程，教师逐项板书，强调单位统一（mmHg转换为cmH₂O需补充说明，但课内直接使用mmHg，以与临床指标对接）。【非常重要】【高频考点】

3.非线性环节处理（15分钟）：指出GFR本身也依赖于Pbc，因为有效滤过压中减去了囊内压。因此状态方程变为：dVbc/dt=Kf×(Pgc–Pbc(t)–πgc平均+πbc)–(Pbc(t)–Ppt0)/Rpt。这是一个关于Vbc（或Pbc）的一阶非线性常微分方程。教师指出，由于πgc平均实际上依赖于血浆流量和滤过分数，严格模型应引入第四个变量——毛细血管内蛋白浓度，但本课聚焦核心思维，先假设πgc平均为定值（32mmHg），后续课程扩展。此处标注【难点】，但要求学生务必掌握该方程各变量的生理对应物。

4.参数生理边界确定（10分钟）：以表格形式（但此处用段落描述）逐一给出典型生理参数范围：Kf4.8~12.5mL/(min·mmHg)（年龄、性别差异）；Pgc45~60mmHg；Pbc10~15mmHg；πgc25~35mmHg；Rpt0.5~2.0mmHg·min/mL。强调此为群体统计值，个体差异显著，且测量手段不同（微穿刺、间接估算）结果不同。此处进行【基础】标注，并指出在后续仿真作业中，学生必须从文献中自行查找并引用一组自洽的参数。

（四）第三、四学时（连堂）：计算机仿真实验——从参数扰动到病理表型

本环节在计算机机房进行，2学时连续，保证思维不中断。教学形态为“半自主仿真工作坊”，教师巡回指导，每15分钟进行一次阶段性发现共享。

1.模型原型实现（25分钟）：学生基于课前代码框架，在教师分步提示下完成核心微分方程的数值积分。MATLAB学生可使用ode45，Python学生使用egrate.solve_ivp。教师投影演示如何将状态方程写为函数文件，特别注意Pbc与Vbc的关系——此处简化假设肾小囊为弹性囊，顺应性Cbc=dV/dP≈常数0.5mL/mmHg。这一简化虽然粗糙，但足以反映基本动态。教师强调这是模型假设的一部分，鼓励学生在拓展作业中替换为更真实的非线性顺应曲线。

2.正常生理稳态仿真（18分钟）：运行代码，设定典型参数（Kf=10.0,Pgc=50,πgc平均=30,Pbc初始=10,Rpt=1.2）。观察Pbc(t)与GFR(t)从初始值过渡到稳态的过程。学生记录稳态GFR数值（应接近临床正常值125mL/min）。教师引导发现：由于负反馈（Pbc升高抑制GFR），系统是自稳定的。此处可略微调高Rpt，模拟近端小管梗阻，可见Pbc上升、GFR下降——第一次病理扰动。【重要】

3.病理参数映射专项研讨（22分钟）：此为整个课时的巅峰环节。教师展示四种典型肾小球疾病的病理生理改变卡片（电子随机分发至各组），要求各专家组完成从病变到参数变化的映射，并在组内用自然语言描述模型应修改哪个参数。例如：

1.4.糖尿病肾病早期：基底膜增厚，但滤过面积代偿性增加→Kf可能先升高后降低；同时因高滤过，Pgc升高。学生需决定模型中是改Kf还是改Pgc，并陈述理由。

2.5.新月体肾炎：细胞增生压迫毛细血管袢，导致有效滤过面积S急剧下降→Kf显著降低。且常伴有入球小动脉收缩→Pgc降低。此为双重打击。

3.6.膜性肾病：免疫复合物沉积于基底膜上皮下，膜本身通透性下降（k降低），但面积S变化不大→Kf降低，且由于电荷屏障丢失，大分子蛋白滤过增加——但本模型不涉及溶质，故只改Kf。

各专家完成映射后回到基础组，选择一种病理场景，修改仿真代码中的参数（例如Kf乘以0.4，或Pgc降低15%），重新运行仿真，记录GFR稳态值下降百分比及达到新稳态的时间常数。此环节教师需反复强调：参数变化必须依据病理生理逻辑，不可随意捏造数值。并提供参考文献截图作为参考依据（课上口头说明来源，不显示网址）。【非常重要】【热点】【难点】

7.仿真结果临床关联分析（20分钟）：每组派代表展示病理仿真曲线，并尝试回答“如果该患者肾小球滤过率下降至60mL/min，根据模型推测，其主要损伤因素是Kf还是血管压力？”这一推理题。通过对比Pgc下降与Kf下降两种情景下Pbc的代偿性变化差异，引导学生识别不同病因的血流动力学指纹。例如，单纯Kf下降时，Pbc往往略低于正常（因为滤液生成少，囊内压难以上升）；而单纯Pgc下降时，Pbc变化不明显。此分析将抽象数学模型与临床鉴别诊断思维紧密联结。【高阶能力目标】

8.模型敏感性初步分析（10分钟）：教师演示如何通过循环脚本，在Kf正常值上下30%范围内扫描，观察GFR变化率，并计算敏感度系数。学生直观看到GFR对Kf高度敏感，提示临床监测Kf相关生物标志物（如足细胞脱落数目）对于早期诊断的价值。同时，扫描Rpt发现，只有当Rpt增大至3倍以上时GFR才显著下降，这与慢性肾病后期肾小管间质纤维化的临床进程相符。此部分渗透系统生理学的工程思维，并为后续课程的多尺度建模埋下伏笔。

（五）课堂总结与高阶思维提升（10分钟，穿插在第四学时后半段）

1.教师以概念图形式（语言叙述）总结从“足细胞”到“微分方程”再到“仿真曲线”的完整转化路径，强调生物医学工程师的核心竞争力正是搭建这座桥梁。

2.提出三个开放式延伸问题，作为课后的思辨任务，不计入分数但纳入课堂表现记录：[1]本模型假设πgc平均恒定，请查阅文献设计一个迭代算法，使πgc随滤过分数动态更新，并评估新模型与现有模型的差异大小。【难点】【重要】[2]如果考虑肾小管-肾小球反馈（TGF），模型需要增加什么环节？请绘制信号流图并用语言描述。[3]肿瘤患者在化疗期间使用顺铂，易发生急性肾损伤。根据顺铂主要损伤近端小管上皮细胞的病理知识，应修改本模型中哪个参数？你的修改依据是什么？

六、教学评价与反馈

（一）形成性评价

1.课前预习题准确率：关注对Kf含义的理解，错误率高于30%则在课中增加Kf组分拆解互动。

2.课中代码调试速度：教师通过局域网监控工具观察学生模型运行成功的时间节点，对20分钟内未完成ode45基础框架的学生进行邻近同伴互助。

3.病理映射研讨记录单：各专家组需提交一份电子便签（50-100字），描述“病变-参数-改动幅度”三元组，教师根据逻辑一致性评定为A/B/C等级，计入平时成绩30%。

4.随堂测验：第四学时最后5分钟发布一道选择题——“在肾小球滤过动力学模型中，若将鲍曼囊顺应性改为非线性（高压下顺应性降低），预测在急性肾梗阻时，GFR的下降速度会比线性模型更快还是更慢？”意在检测学生对负反馈环路动态特性的理解深度。

（二）终结性评价

课后大作业（三人小组）：自选一种非典型肾小球疾病（如Alport综合征、淀粉样变性、Fabry病），完成从文献综述、模型参数扰动设计、仿真实现到结果临床意义解读的完整报告，篇幅限2000字以内，附源代码及注释。评价量规涵盖：病理机制与参数的映射合理性（40%）、仿真代码正确性与可重复性（30%）、结果讨论与临床数据对比（20%）、参考文献规范（10%）。此作业提交至课程GitHubClassroom仓库，并开启同伴互评（PeerReview）环节。【高频考点】【热点】

七、教学资源与支持系统

（一）数字化资源库：PhysioNet开源肾脏生理数据集；美国生理学会教学资源中心“RenalExcretion”交互模块；3D解剖软件CompleteAnatomy的肾脏滤过单元。

（二）核心参考文献（课上口述或PPT展示，不印发链接）：[1]GuytonAC,HallJE.bookofMedicalPhysiology.14thed.Chapter26.[2]WeinsteinAM.Amathematicalmodeloftheratkidney:Kimmel–Kaufman–Knoxequations.AmJPhysiolRenalPhysiol.2020.[3]EdwardsA,LaytonAT.Mathematicalmodelingofrenalphysiology.ComprPhysiol.2014.（此处仅提供文献描述，无链接无网页序号）

（三）硬件保障：教室配备高分辨率数位屏，教师板书与学生代码界面可同时无线投射至大屏，实现讲、练、评一体化。

八、教学预设与应变策略

（一）若多数学生对微分方程数值解法感到困难，则立即暂停编程推进，用5分钟类比匀加速运动位移计算（欧拉法思想），强调“虽不精确但可逼近”的工程实用主义，降低认知负荷。

（二）若仿真实验中出现参数设置合理但结果发散（例如Pbc无限增长），引导学生检查符号错误：GFR表达式中Pbc前应为负号。这种错误是建模学习中极有价值的“认知冲突”，教师借此