简介本资源是一套面向生物力学与连续介质力学研究者的ABAQUS用户子程序UMAT通用开发框架专为实现大变形本构模型的三维Fortran编码而设计适用于软组织、细胞等复杂生物材料的本构建模需求。压缩包共82个文件涵盖49个核心Fortran源码.for/.f90、7个ABAQUS输入文件.inp、4个参数/宏定义头文件.inc、3个Mathematica符号推导脚本.mw以及测试脚本.sh、Python后处理工具.py、文档说明.md/.txt/.rpt和可视化图示.png整体仅800KB轻量但结构完整。已有239人学习下载。用户可直接复用模块化子程序如invariants.for、spectral.for、proj_lagrangian.for、umat_general.for等快速集成各向异性、体积响应、纤维方向演化、应力更新算法等关键功能并通过配套inp算例与sta/pes输出文件验证模型行为显著降低UMAT开发门槛与调试成本。 如果你已经在ABAQUS里跑过几个分析迟早会碰上这样的事内置本构模型要么调不出想要的硬化曲线要么面对实验数据对不上。这时候周围人就会意味深长地丢给你一句去写UMAT啊。可等你真打开Fortran文件面对那一长串变量和让人头大的环境配置往往一卡就是好几天。这篇东西我不讲高深的本构理论推导只聊清楚一件事怎么把一套能用的UMAT开发环境搭起来再照着通用框架把材料模型写进去最后让它真的在ABAQUS里跑出结果。顺带把那些高频出现的环境错误比如-97、FlexNet版本问题、ABAQUS和UG许可证打架一次性掰开揉碎。1. UMAT在求解器里到底扮演什么角色1.1 材料子程序的核心任务更新应力和给刚度先明确一个基本认知UMAT在ABAQUS/Standard隐式分析中做的事情本质上就两件——根据当前应变增量更新应力以及给出材料雅可比矩阵DDSDDE。说得直白一点求解器在每个增量步把应变增量DSTRAN扔给你UMAT算完把新的应力STRESS和切线刚度DDSDDE还回去。这个DDSDDE决定了Newton-Raphson迭代的收敛速度如果给的刚度方向不对整个增量步就可能要反复迭代甚至直接发散。很多第一次写UMAT的人会把精力全放在推导公式上结果算出来总是不收敛。原因不是本构模型错了而是DDSDDE和应力更新算法不自洽。比如做弹塑性的时候应力返回映射写得没问题但一致切线模量图省事用连续切线模量代替大增量步下收敛性就会明显变差。记住一个原则UMAT的返回值不只是材料行为更是在告诉求解器这个材料在当前状态下往哪个方向走、走多快后者直接决定了数值稳定性。1.2 内置材料模型不能满足的那部分需求ABAQUS内置的材料库其实相当庞大从线弹性、金属塑性到混凝土损伤、Hyperelastic覆盖了大多数常规工程场景。但实际工程项目里总会出现一些内置模型覆盖不到的情况。举几个最常见的例子做金属热加工仿真需要位错密度作为内变量、考虑动态回复和再结晶的率相关本构做橡胶和生物软组织多场耦合分析需要把纤维分散度写进应变能密度做岩土本构研发要用非关联流动法则配合双屈服面、考虑剪胀角随剪切带演化。这些模型在内置库里找不到现成方案而UMAT就是官方留出来的自定义接口。ABAQUS在每个积分点调用UMAT用户完全可以决定状态变量的演化、应力的更新方式以及刚度矩阵的组装形式。1.3 UMAT与VUMAT怎么选正式动手之前先分清UMAT和VUMAT的应用边界。UMAT用于ABAQUS/Standard隐式分析最大的特点是每一步都要组装整体刚度并迭代对DDSDDE的准确性要求极高也因此很适合那些增量步大、对精度要求高的静力、蠕变、疲劳问题。VUMAT则跑在ABAQUS/Explicit显式分析里面采用的是中心差分格式不需要迭代收敛但为了保证稳定时间步会被临界稳定步长卡住。我个人的选型习惯是如果问题本身是准静态或者需要用隐式方法处理强接触优先走UMAT如果问题是冲击、碰撞、爆炸这类高度非线性动态过程或者模型单元数量太大显式配合VUMAT更合适。千万不要把UMAT代码直接改改关键字硬挪到Explicit里用UMAT和VUMAT的变量约定、切向刚度计算方法完全不同。2. 开发环境搭建Fortran编译器、Visual Studio与ABAQUS的版本匹配2.1 版本匹配是第一条生死线UMAT开发绕不开Fortran而ABAQUS在Windows平台上的子程序编译依赖的是一整套Intel Fortran编译器 Microsoft Visual Studio组合。版本匹配这个事说多了都是泪。不少人的UMAT代码本身写得没问题结果因为编译器版本不对卡在abaqus verify阶段过不去。以ABAQUS 2026为代表的较新版本官方对Intel oneAPI工具链支持得比较好。历史版本里常见的组合是ABAQUS 2020配合Intel Parallel Studio XE 2020 Visual Studio 2019ABAQUS 2023之后推荐迁移到Intel oneAPI Fortran编译器。安装的时候建议先装Visual Studio再装Intel编译器最后装或者配置ABAQUS。这个顺序不是随便写的因为Intel编译器安装时会自动探测Visual Studio并写入相关环境配置如果顺序反了后面就得手动去设置IFORT_COMPILER和VS安装路径非常折腾。安装完成后别急着写代码。先打开ABAQUS Command执行abaqus verify -user_std这一步会实际编译一个自带测试子程序确认编译器链路是通的。假如输出里出现类似Problem during compilation - ifort.exe not found或者Unable to locate the compiler九成是环境变量或者版本匹配出了问题。这时检查两个地方一是IFORT_COMPILER环境变量是否指向Intel编译器的安装根目录二是ABAQUS安装目录下的custom_v6.env里是否能看到compile_fortran [ifort]这样的配置。2.2 环境变量与PATH配置的细节很多人以为装完软件就能直接用实际上ABAQUS调用Fortran编译器的时候依赖的是当前终端环境下能访问的ifort或ifx命令。ABAQUS图形界面启动的子程序任务有时读不到你手动设置的临时环境变量所以最稳妥的做法是把编译器目录和Visual Studio的vcvars64批处理路径固定写进系统环境变量。以Intel oneAPI为例常见路径是C:\Program Files (x86)\Intel\oneAPI\compiler\latest\env\vars.bat和C:\Program Files (x86)\Intel\oneAPI\setvars.bat。在运行ABAQUS之前确保这台机器的系统环境变量里已经有Intel编译器相关的注册信息。如果你用的是Visual Studio Build Tools而不是完整版VS要尤其小心因为ABAQUS在检测VC编译器时经常因为找不到cl.exe而报错这时需要在custom_v6.env里明确设置VS的安装路径或者手动把cl.exe所在目录加入PATH。提示abaqus verify -user_std验证失败时别急着重装。先在命令行里手动执行一次ifort如果系统提示不是内部或外部命令说明编译器环境根本没生效这是环境变量问题不是ABAQUS问题。2.3 VSCode做Fortran子程序编辑的配置思路新版本ABAQUS装好之后很多人喜欢用VSCode写Fortran子程序但会遇到VSCode提示找不到ifort这类问题。原因很简单VSCode内置终端不会自动加载Intel编译器环境。解决办法是在VSCode的settings.json里配置终端启动参数让每开一个终端就执行一次Intel的vars.batterminal.integrated.profiles.windows: { Intel Fortran: { path: C:\\Windows\\System32\\cmd.exe, args: [ /k, \C:\\Program Files (x86)\\Intel\\oneAPI\\setvars.bat\ intel64 vs2019 ] } }这样配置之后VSCode里打开的终端就天然具备ifort命令。同时建议安装Modern Fortran扩展它能识别Fortran源文件的语法、跳转和错误提示。不过注意VSCode里的lint提示只是辅助最终编译和链接仍然是通过ABAQUS调用编译器完成的VSCode本身不会也不能替你编译UMAT。如果你想从头排查、看详细编译输出可以在提交job的时候在命令行里加-verbose参数。abaqus jobjobname userumat.for -verbose会输出完整的编译链接过程哪里找不到头文件、哪里链接失败一目了然。这个做法比在CAE界面里反复看一遍又一遍模糊的报错弹窗高效得多。3. 从零写一个UMAT接口结构和通用代码骨架3.1 先认识调用接口不是所有参数都要用到UMAT子程序的固定开头是长长的一行SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL, 1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, 2 TIME, DTIME, TEMP, DTEMP, PREDEF, DPRED, 3 CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, 4 NPROPS, COORDS, DROT, PNEWDT, CELENT, 5 DFGRD0, DFGRD1, NOEL, NPT, LAYER, KSPT, 6 KSTEP, KINC)新手很容易被这一长串参数吓住其实每个版本都大同小异关键是先把高频使用的那几个弄明白。STRESS进入子程序时是增量步开始时的应力退出时必须更新为增量步结束时的应力。STATEV状态变量数组适合存放塑性应变、内变量、损伤变量等。DDSDDE当前切线刚度矩阵NTENS×NTENS大小直接影响收敛。STRAN和DSTRAN总应变和当前增量步的应变增量。DTIME当前增量步的时间步长做率相关材料的时候必用。PROPS和NPROPS材料参数数组及数量对应CAE里输入的材料常数。NDI和NSHR正常应力分量个数和剪应力分量个数三维问题NDI3、NSHR3平面应变NDI2、NSHR1。COORDS积分点坐标处理非均匀材料或梯度材料时非常有用。NOEL和NPT单元编号和积分点编号用于调试。3.2 一个最简线性弹性UMAT示例先别急着上超弹性、晶体塑性我建议所有人写第一个UMAT都从线性弹性开始。一方面能验证环境是否完全打通另一方面也方便和内置弹性模型的结果做对照。SUBROUTINE UMAT(STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL, 1 DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, 2 TIME, DTIME, TEMP, DTEMP, PREDEF, DPRED, 3 CMNAME, NDI, NSHR, NTENS, NSTATV, PROPS, 4 NPROPS, COORDS, DROT, PNEWDT, CELENT, 5 DFGRD0, DFGRD1, NOEL, NPT, LAYER, KSPT, 6 KSTEP, KINC) C IMPLICIT NONE C CHARACTER*80 CMNAME DIMENSION STRESS(NTENS), STATEV(NSTATV), DDSDDE(NTENS, NTENS) DIMENSION STRAN(NTENS), DSTRAN(NTENS), TIME(2), PREDEF(1), DPRED(1) DIMENSION PROPS(NPROPS), COORDS(3), DROT(3, 3) DIMENSION DFGRD0(3, 3), DFGRD1(3, 3) C INTEGER NDI, NSHR, NTENS, NSTATV, NPROPS INTEGER NOEL, NPT, LAYER, KSPT, KSTEP, KINC DOUBLE PRECISION STRESS, STATEV, DDSDDE, SSE, SPD, SCD, RPL DOUBLE PRECISION DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN DOUBLE PRECISION TIME, DTIME, TEMP, DTEMP, PREDEF, DPRED DOUBLE PRECISION PROPS, COORDS, DROT, PNEWDT, CELENT DOUBLE PRECISION DFGRD0, DFGRD1 C C PROPS(1) E, PROPS(2) NU C DOUBLE PRECISION E, NU, LAMBDA, MU DOUBLE PRECISION C11, C12, C44 INTEGER I, J C E PROPS(1) NU PROPS(2) LAMBDA E * NU / ((1.0D0 NU) * (1.0D0 - 2.0D0 * NU)) MU E / (2.0D0 * (1.0D0 NU)) C C 初始化切线刚度矩阵 DO I 1, NTENS DO J 1, NTENS DDSDDE(I, J) 0.0D0 END DO END DO C C11 LAMBDA 2.0D0 * MU C12 LAMBDA C44 MU C DO I 1, NDI DDSDDE(I, I) C11 DO J 1, NDI IF (J .NE. I) THEN DDSDDE(I, J) C12 END IF END DO END DO DO I NDI 1, NTENS DDSDDE(I, I) C44 END DO C C 更新应力 DO I 1, NTENS DO J 1, NTENS STRESS(I) STRESS(I) DDSDDE(I, J) * DSTRAN(J) END DO END DO C C 记录等效弹性应变能 SSE 0.0D0 DO I 1, NTENS SSE SSE STRESS(I) * STRAN(I) * 0.5D0 END DO C RETURN END这个骨架里最有价值的地方在于它展示了DDSDDE的组装方式和应力更新的基本流程。平面应力和轴对称情况的DDSDDE需要特殊处理因为存在约束关系不过那是后续再说的事。把这版跑通、和内置弹性模型对比完全重合就可以确认编译器链路、UMAT调用逻辑、数组维度全部正确。3.3 UMAT通用代码结构初始化、计算、输出三段式随着材料模型复杂度上升代码逻辑需要有一个稳定的分层结构。我的习惯是把UMAT写成三个逻辑段。第一段是初始化。把PROPS里的材料参数读出来命名要直观不要直接PROPS(1)、PROPS(2)满天飞。工程上常见参数有上百个变量名就是维护成本。第二段是核心计算。根据STATEV里存放的历史状态依据当前应变增量更新应力和状态变量。这一步是本构模型的主战场。写代码前先在草稿纸上把增量步的更新逻辑写清楚哪个量是隐式更新、哪个量需要迭代求解、状态变量在何时更新都要有明确的分工。第三段是收尾输出。把需要输出的物理量写入STATEV这些变量可以在后续的CAE后处理中用SDV1、SDV2等输出非常方便可视化。比如等效塑性应变、塑性耗散、损伤量、硬化参数尽量单独分配到独立的STATEV槽位里不要一段代码无数个魔法数字穿插。另外强烈建议在代码里加一个底层的收敛判断变量PNEWDT。当本构迭代不收敛或预测当前增量步无法准确计算时给PNEWDT赋一个小于1的值比如0.5ABAQUS就会自动折减时间步重新尝试。这是UMAT中非常实用的稳定性控制手段很多本构代码分支复杂没有这层保护一个错误的参数组合就可能让整个job陷入无穷的增量步缩小循环。4. 高频环境与许可证错误排查-97、FlexNet、UG冲突4.1 错误-97的本质和排查链路ABAQUS报错-97在子程序开发和安装调试时非常常见又极其劝退新人。所谓-97很多时候是FlexNet许可证通信层返回的错误码它本身并不告诉你具体的故障点是证书过期、服务器不可达还是版本不兼容需要按链路一层层排出。我的排查顺序是固定的。先打开ABAQUS License环境工具检查当前激活的许可证服务器地址是否可达然后确认系统环境变量LM_LICENSE_FILE没有被其他软件的变量覆盖掉接着检查许可证日志看服务器是否把feature发放给了当前客户端最后再考虑重装或更新License Manager。注意多个CAD/CAE软件共存时LM_LICENSE_FILE被后安装软件覆盖是特别常见的状态这也是后面要说的ABAQUS和UG冲突的根源之一。4.2 Your ABAQUS license server is running with an unsupported version of FlexNet的成因与处理这个错误几乎出现在某个新版本ABAQUS比如2026配旧许可证服务的场景里。ABAQUS新版安装包对FlexNet版本的依赖是上行的旧版本FlexNet做不到让新版本license feature正常激活。报错信息说得很直白你的许可证服务器正在用一个不支持的FlexNet版本跑。处理思路就两条升级许可证管理器或者让新ABAQUS指向一台已经升级过FlexNet的服务器。关键是别只在客户端这边折腾因为错误描述里的license server是宿主服务不是CAE安装目录所以要找到实际跑着许可证服务的那个机器更新其FlexNet版本。然后在客户端重新运行许可证连接测试确认不会再把旧版本识别为首选服务。4.3 ABAQUS和UG许可证冲突根源在LM_LICENSE_FILE很多工程师的电脑上同时装着ABAQUS和NXUG安装顺序一乱或者两个软件共用一套FlexNet/FlexLM服务就可能出现ABAQUS无法启动或者UG打开正常ABAQUS启动就闪退的怪象。这里真正的根源绝大多数是环境变量冲突。ABAQUS默认用LM_LICENSE_FILE或ABAQUSLM_LICENSE_FILE定位许可证UG的许可证管理也使用FlexLM体系安装时往往会修改LM_LICENSE_FILE指向自己的许可证文件。一旦ABAQUS被写到这个系统变量后面就可能找错服务造成启动即失败。处理办法是在系统环境变量里为ABAQUS单独指定许可证变量并确保其优先级高于通用变量。打开高级系统设置 - 环境变量检查是否有LM_LICENSE_FILE变量确认它指向的端口主机地址到底是谁的许可证服务器。如果是UG的就新建或修改ABAQUSLM_LICENSE_FILE显式指向ABAQUS的那台license server。这样ABAQUS启动时会优先使用专用的变量再也不会被UG干扰。4.4 ABAQUS无法启动的常见隐藏原因除了许可证冲突ABAQUS无法启动还有一个非常隐蔽的原因当前Windows用户名或者工作目录包含中文字符和空格。UMAT编译过程涉及临时文件生成和路径传递有些编译器和批处理脚本对中文路径兼容性很差。典型症状是GUI能打开但提交standard求解时提示找不到临时目录或权限不足。这种问题排查起来很简单新建一个纯英文的Windows用户或者把工作目录改到C:\temp这种纯ASCII路径下再试一次。不少用户折腾了半天许可证最后发现是当前用户名里的中文导致link阶段失败白白浪费一晚上。另外安装完新版本ABAQUS后建议重新启动一次电脑确保Intel编译器环境变量、FlexNet服务路径设置全部生效。有的同事装完不重启就急着验证结果遇到ABAQUS无法启动其实是服务还没注册完全。5. 调试UMAT的实用技巧从print埋点到单单元验证5.1 在UMAT里埋print输出这是最快的内因定位方式写UMAT最难受的处境是job提交了分析报错Too many attempts made for this increment你完全不知道是本构方程求解失败还是切线刚度给错。这时候print输出是救命稻草。在Fortran子程序里WRITE(6,*)会输出到job的标准输出也就是正常情况下可以在ABAQUS交互界面或者命令行看到。WRITE(7,*)则会输出到.msg文件适合输出大量增量步过程信息。你在关键变量更新后加一行IF (KINC .EQ. 1 .AND. NPT .EQ. 1) THEN WRITE(7,*) KSTEP, KSTEP, KINC, KINC, NPT, NPT WRITE(7,*) DSTRAN, DSTRAN(1), DSTRAN(2), DSTRAN(3) WRITE(7,*) STRESS_NEW, STRESS(1), STRESS(2), STRESS(3) END IF用KSTEP、KINC、NOEL、NPT几个变量做过滤条件就能只输出你关心的那个单元、那个积分点、那个增量步的信息而不是被海量数据淹没。注意调试完成后正式跑大规模计算前务必把这些调试输出注释干净。密集的WRITE语句会严重拖慢计算速度几百个增量步的job可能因为输出I/O多跑好几倍时间。5.2 单单元单增量步验证法调试UMAT科学与否决定了你后面整条开发线会不会崩。我的习惯是把一切复杂条件都降维建一个单单元模型一个增量步加载然后手动对比结果。这一步在umat验证里的价值怎么强调都不过分。比如你要验证弹塑性模型就建一个C3D8R单元所有面施加约束只留一个方向加载时间步设成1然后通过*USER MATERIAL输入参数。跑完以后把UMAT算出来的应力-应变数据导出和同一问题用内置塑性模型跑出来的结果对比。如果曲线重合UMAT的算法和内置模型等价如果偏了就用5.1节的print埋点检查每个增量步结束时STATEV的演化是否合理。单单元验证之后的第二步是做一个非均匀应力状态测试比如含孔板单拉用来检验UMAT在大变形、剪切和复杂应力状态下的稳定性。这样两步走完UMAT才算真正可以交到工程项目上使用。5.3 利用ABAQUS内置模型做背靠背对照很多本构参数是从内置模型调试出来的所以背靠背对照非常有用。比如你想写一个更复杂一点的Ramberg-Osgood硬化本构可以用ABAQUS内置的塑性模型配上同样参数跑一遍再把两边的位移场、应力场、支反力做全局对比。如果UMAT写错了单单元可能碰巧能过但整体结构响应会在某几个单元先发散背靠背对照能把这种分布性错误揪出来。再推荐一个冷门但有奇效的手段在后处理里输出单元积分点的SDV再叠加输出S、LE这些内置场变量。一旦发现SDV值异常到离谱第一时间在UMAT里用5.1的过滤条件锁定该单元定位是材料参数读取错误还是状态变量更新逻辑出问题。5.4 编译期与运行期报错的区分UMAT报错分两类处理思路完全不同。编译期报错比如Problem during compilation说明Fortran代码本身有语法错误或者编译器环境没打通要看VSCode或命令行输出的错误行号直接对着代码改。运行期报错比如Standard aborted with error code或计算中途不收敛多半是本构逻辑错误比如状态变量未初始化、DDSDDE零矩阵、除零、数组越界。有一个细节我见过很多人踩坑UMAT里的隐式DOUBLE PRECISION问题。用IMPLICIT REAL*8 (A-H, O-Z)可以省去一长串类型声明但也容易掩盖问题。比如循环变量I、J会被隐式声明成INTEGER而PROPS这种变量名以P开头会被隐式声明成REAL*8一旦你在PROPS传递时不小心用错类型就可能出现数组越界或者NaN。最稳的方法还是像我前面代码里那样用IMPLICIT NONE强制声明所有变量初期麻烦一点后期调试能省十倍时间。6. 进阶从UMAT框架到可维护的材料库6.1 UMAT封装与模块化设计UMAT本质上是ABAQUS在每一个积分点上调用的一段Fortran过程。如果材料本构逻辑复杂把所有代码都堆在SUBROUTINE UMAT里最后会变成几万行的大泥球完全没法维护。我推荐把核心本构放到独立的Fortran模块或子函数里UMAT只做接口转换。比如建立material_lib.f90里面包含弹性子程序、塑性本构子程序、损伤演化子程序每个子程序的入参明确返回新的应力和一致的切线刚度。UMAT里只负责把ABAQUS传来的STATEV数组拆包、解码成实际问题中的物理量调用子程序得到结果后再装包回去。这样每个材料子程序可以单独单元测试后续增加新材料也只需在UMAT里加一个CASE分支。6.2 从单元测试到完整验证的流程建议框架搭好后我把整个UMAT开发流程分成四层验证。第一层是单单元单轴加载看基本应力应变响应第二层是单单元多轴加载验证不同应力路径下状态变量演化是否合理第三层是含缺口的代表性结构模型看整体响应和收敛性第四层是真实零件模型加入接触、边界条件、多载荷步看UMAT在复杂工况里是否稳定。这四层下来本构算法问题基本都能暴露出来。6.3 常见UMAT错误速查表这里整理一份高频报错的对照表基本覆盖我踩过和帮朋友排过的大多数坑。错误表现常见原因处理建议Problem during compilation - ifort.exe not found编译器环境没生效或版本不配套执行setvars.bat后重新验证Unable to open file: umat.obj源代码有语法错误或模块引用路径不对检查Fortran代码行号License error -97FlexNet服务异常或版本不匹配按4.1节逐层排查Unsupported version of FlexNet许可证服务器FlexNet版本过旧升级服务器端FlexNetABAQUS startup fails without clear reason环境变量LM_LICENSE_FILE被其他软件覆盖设置ABAQUSLM_LICENSE_FILEToo many attempts for increment本构迭代不收敛或DDSDDE给错检查PNEWDT和切线刚度NaN in stress/strain变量未初始化或除零用PRINT定位首个异常增量步SDV output is zeroSTATEV赋值位置错误检查是否在最后一步被覆盖这六类问题的共性根子基本都能往前回溯到两条要么环境没配对要么本构状态更新和切线刚度不一致。环境问题按前面章节的方法处理本构问题就得回到算法本身把每一个状态变量的更新式和DDSDDE的推导过关。做UMAT这件事外表看着是一大堆Fortran代码和ABAQUS界面里的报错弹窗底子里其实是数值算法、编程规范、软件工程、求解器机制四样东西拧在一起。大多数人卡住缺的往往不是最前沿的本构理论而是从理论公式到稳定数值实现之间那一层调试能力。能坚持把环境搭建、普通跑通、对比验证这三个阶段走完后面再换什么材料本构主动权基本就在你手里了。我自己的习惯是每套UMAT都保留一张单单元验证模型和一份背靠背对照数据换机器、换版本、换编译器之后先跑一遍回归省下的排查时间比写代码的时间多得多。本文还有配套的精品资源点击获取