蠕变中文解释
蠕变
目前比较公认的是以位错理论对蠕变做出 的解释，但目前仍然停留在定性阶段。 位错理论可以用下图来简单表示：
施加应力 各晶粒内出现位错增殖 晶内加工硬化（低温时） 温度升高 热振动、原子扩散加剧
Balance
位错相消
回复（位错易移动）
3.2、对稳态蠕变的理论解释
当这种加工硬化与回复成平衡状态时就是 稳态蠕变。 所以实际上蠕变的位错理论可以总结为是 加工硬化产生的位错增殖与回复的竞争过 程。
蠕变断裂机理
晶界滑动机制 中等蠕变温度和较高应力水平。 空位聚集机制 较高温度和较低应力水平。
晶界滑动机制（V型裂纹形成）
空位聚集机制（O型裂纹形成）
4、蠕变强度及金属组织
在高温环境长期服役的构件通常会出现 蠕变现象，为此提出相应的性能指标以满足 设计的需求。 材料的蠕变强度目前尚未有一致的定义。 根据使用中的尺寸变化来规定设计条件时， 第一阶段和第二阶段蠕变应变或应变速率是 研究对象；根据到达断裂的耐用寿命来规定 设计条件时，断裂时间是研究对象。 目前常用的蠕变性能指标有：蠕变极限、 持久强度
5、蠕变试验方法
测定蠕变极限、持久强度的基本试验装置 多为一种杠杆式的静加载系统。 加载方法为：在杠杆上设有分载荷，随着试 样的伸长逐渐移动分载荷。 另外还有安德雷德的浮力法以及逐渐改 变杠杆有效长度的方法等。
拉伸蠕变试验机
6、实际中的蠕变断裂
6.1 焊接区的蠕变 随着焊接技术的 发展，在以发电用 锅炉为主的高温用 机器上，大量地采 用了焊接结构，焊 接区的蠕变强度， 实际上是一个极其 重要的问题。
4.1 蠕变极限
蠕变极限:高温长时载荷下材料对变形的抗 力指标。 表示方法（主要有以下两种）： ⑴在给定温度T(℃)下，使试样产生规定的恒 定蠕变速率的应力值， 。 ⑵在给定温度T(℃)和规定时间t(h)内，使试 样产生一定蠕变应变量的应力 值， T/ t , 1500 100N / mm2 。 /10
工业秤蠕变
工业秤蠕变标题：工业秤的蠕变之谜第一部分：工业秤的起源和功能工业秤，作为一种重要的计量工具，被广泛运用于各个行业。
它的主要功能是测量物体的质量，以便进行生产、贸易和研究等工作。
工业秤的设计和制造经历了长期的发展和改进，现如今已经成为一种高度精确的计量工具。
第二部分：工业秤的蠕变现象然而，近年来，一种令人困惑和让人难以理解的现象开始在工业秤中出现，那就是蠕变。
蠕变是指物体受到持续应力作用下，随时间发生的形变。
在工业秤中，蠕变现象意味着秤盘或秤座等部件会发生不可忽视的形变，从而导致测量结果的不准确。
第三部分：蠕变的原因及影响蠕变现象的出现往往是由于工业秤长时间承受着巨大的负荷和应力，这些应力会导致材料的分子结构发生变化，从而引起形变。
同时，环境的温度和湿度等因素也会加剧蠕变的发生。
蠕变对工业秤的影响是多方面的。
首先，蠕变会导致秤盘的变形，使得测量结果产生误差。
其次，蠕变还会对秤盘的稳定性和耐久性造成影响，从而缩短了工业秤的使用寿命。
此外，蠕变还可能对秤盘的结构强度产生不利影响，进一步降低了工业秤的可靠性和安全性。
第四部分：蠕变的解决方案为了解决工业秤中的蠕变问题，科学家和工程师们进行了大量的研究和探索。
他们发展出了一系列新型材料和设计理念，以提高工业秤的抗蠕变能力。
例如，采用高强度和耐热性能较好的材料制造秤盘和秤座，可以有效减少蠕变的发生。
此外，合理的结构设计和精确的制造工艺也可以降低蠕变的风险。
第五部分：工业秤的未来展望尽管蠕变问题对工业秤的准确度和可靠性带来了一定的挑战，但科学家和工程师们相信，通过不断的研究和创新，可以找到更好的解决方案。
他们希望能够开发出更稳定、更精确的工业秤，以满足不断发展的产业需求。
总结：工业秤的蠕变问题是一个复杂而又具有挑战性的课题。
通过深入研究和创新，我们可以找到解决方案，提高工业秤的准确度和可靠性，从而更好地服务于生产和贸易等领域。
希望未来能够有更多的科学家和工程师致力于解决这一问题，为人类创造更好的计量工具。
第1章第5节蠕变
28
48.2 45.1 41.8 38.2 34.3 30.3 26 21.5 16.7 11.7 6.6
表 1-2 叶片各截面σm 值
根据叶片的工作温度和应力选用温度为 550℃和应力为 98.07MN/m2 的蠕变曲线，其蠕 变速度 Vcr=1.2×10-7h-1,如图 1-20。 由已知各截面的应力值 σ ，计算各截面的σm 值，列入表 1-10 中，应用数值积分法得 到
550 1 × 10
0 − 5
= 90 MN/m2，是指当温度为 550℃时，蠕变速度为 1×10-5%/h 所对应
的应力为 90MN/m2。 表明材料高温强度特性的另一个重要指标时持久强度极限。 在一定温度下， 经过一定的 时间间隔后引起试件断裂（相当于图 1－ 19 上的 D 点）的应力叫做持久强度极限。例如 Cr11MoV 的持久强度极限为 σ
Ω(t ) 2 =
Vcr 1.2 × 10 −7 t = × 105 = 7.97 ×10 −8 ( MN / m 2 ) − m σm (98.07) 2.6
1 0
∆l2 = Ω(t ) 2 ∫ σ m dx = 7.97 × 10 −8 × 10.3 × 10 2 = 0.82 × 10 − 4 m
550 10 5
0
= 160 MN / m
2
，它表示温度为 550℃时，经过 105
小时造成断裂的应力为 160MN/m2。 为了进行蠕变计算，需要将蠕变引起的变形（ ε c⋅r = Vc⋅r t ）与有关因素的关系列为数学 表达式。从大量试验结果归纳出来的重要经验公式是
蠕变和疲劳的关系.docx
蠕变和疲劳的关系蠕变和疲劳的关系蠕变和疲劳是材料科学中两个重要的概念。
蠕变是材料在长时间高温条件下的塑性变形，而疲劳则是材料在循环加载下逐渐失效的现象。
这两个概念虽然在表面上看起来没有直接的联系，但实际上它们之间存在着一定的关系。
首先，我们来探讨蠕变和疲劳的定义和机制。
蠕变是指在高温下，材料受到持续的应力作用，逐渐发生塑性变形和材料结构的改变。
而疲劳则是指在循环应力作用下，材料在应力范围内反复加载和卸载，导致材料内部产生细小的裂纹，最终引发材料失效。
虽然蠕变和疲劳在机制上有所不同，但它们之间存在着共同点。
一方面，高温环境下的蠕变过程本身就是一个不断的应力加载和卸载的过程，因此可以说蠕变是一种特殊形式的疲劳。
另一方面，疲劳过程中也可能存在材料的温度升高，进而引发蠕变效应。
所以，可以说蠕变和疲劳在某种程度上是相互关联的。
其次，我们来讨论蠕变和疲劳对材料性能的影响。
蠕变会导致材料的塑性变形和结构的演变，破坏了材料的原始晶格结构，减弱材料的强度和刚度。
而疲劳则会在材料内部引发裂纹并逐渐扩展，最终导致材料的断裂。
因此，蠕变和疲劳都会对材料的可靠性和使用寿命产生重大影响。
此外，蠕变对疲劳行为的影响也值得关注。
蠕变过程中产生的微观结构演变可能会影响材料的疲劳寿命和疲劳裂纹扩展速率。
例如，蠕变过程中的晶粒生长和细小的裂纹形成，可能导致材料在疲劳加载下更容易发生裂纹扩展。
因此，蠕变可能会加速材料的疲劳失效。
为了克服蠕变和疲劳对材料性能的负面影响，科学家和工程师们采取了一系列的方法和措施。
例如，通过合理的材料选择和设计，可以选择具有良好耐蠕变和抗疲劳性能的材料。
此外，适当的热处理和表面处理也可以提高材料的蠕变和疲劳性能。
总之，蠕变和疲劳虽然在表面上看起来没有直接的联系，但实际上它们之间存在着一定的关系。
它们都会对材料的性能产生重大影响，并可能相互作用。
对于材料科学和工程领域的研究者和工程师们来说，深入理解蠕变和疲劳的关系，以及采取有效的措据研究表明，蠕变和疲劳之间的关系在材料的设计和应用中具有重要意义。
4.8 岩石的蠕变
弹性后效概念：
加载或卸载时，弹性应变滞后于应力的现象。
2. 影响蠕变的因素 岩性 应力 温度、湿度
长期强度概念：时间→∞的强度(σ∞)
/瞬时强度（σc ） σ∞/ σc ： 一般岩石为 0.4-0.8 中、软岩为 0.4-0.6 硬岩为 0.7-0.8
应力
时间
3. 蠕变模型及本构方程
（1）理想物体基本模型
4.8
岩石的蠕变
γ
2
蠕变模型
τ
（1）马克斯威尔(Maxwell)模型： 这种模型是用弹性单元和
γ
η
γ
粘性单元串联而成，见图4-11(a)。
当剪应力骤然施加并保持为常量时， τ 变形以常速率不断发展。这个模型用 两个常数G和来描述。
τ
b
γ
a b
a b
η
τ
η 1 图4-11 线性粘弹性模型及其蠕变曲线 c
弹性元件（弹簧）
粘性元件（阻尼器）
E
，如 G
塑性元件（摩擦片）
0 t
s时，＝0
s时，
4.8
岩石的蠕变
2 蠕变模型
2)粘性模型
或称粘性单元
这种模型完全服从牛顿粘性定律，它表示应力与 应变速率成比例，例如剪应力与剪应变速率的关系 为：
6
图3 试件P3剪切位移与剪切时间的关系曲
试件分别在0.3Mp、0.35Mp、0.45Mp剪切应力作用 下，剪切位移与剪切时间关系曲线发生转折点，由等 速蠕变转变为加速蠕变。因此，试件P1、P2、P3的长 期剪切蠕变强度分别为0.3Mp、0.35Mp、0.45Mp。
根据法向应力和对应的长期剪切蠕变强度的关系曲 线求出长期抗剪强度参数为71.2Kpa，为20.8。
蠕变变形机理
蠕变变形机理蠕变变形是一种材料在高温和恶劣环境下发生的塑性变形现象。
它是由材料内部微观结构的改变所引起的，具有一定的可逆性和不可逆性。
本文将从蠕变变形的定义、机理、影响因素以及应用等方面进行阐述。
蠕变变形是固体材料在高温下由于内部晶格结构的变化而引起的塑性变形。
蠕变变形机理主要有晶体滑移、晶粒边界滑移和扩散等。
晶体滑移是指晶体中的位错沿特定晶面和晶轴方向滑动，从而引起材料的塑性变形。
晶粒边界滑移是指晶体之间的位错滑动，晶粒与晶粒之间发生相对位移，导致材料的变形。
扩散是指材料中原子的相互迁移，在短时间内发生的微观结构的变化。
蠕变变形的机理主要与材料的温度、应力、时间和材料的微观结构等因素有关。
高温是引起蠕变变形的主要因素，因为高温会使材料的晶格结构发生变化，增加了位错的运动和晶粒边界的滑移。
应力是指施加在材料上的力，它会使位错发生运动，从而引起材料的塑性变形。
时间是影响蠕变变形的另一个重要因素，长时间的作用会使材料发生较大的塑性变形。
材料的微观结构也会影响蠕变变形的发生，晶粒的大小、晶界的特性以及杂质的存在都会影响蠕变变形的程度和速率。
蠕变变形在工程中具有重要的应用价值。
例如，在航空航天领域，高温合金材料可以在高温和高应力环境下保持较好的力学性能，从而保证航空发动机等关键部件的正常运行。
在能源领域，蠕变变形的研究可以帮助改进材料的耐热性能，提高能源装置的效率和寿命。
此外，蠕变变形还在材料加工和制备领域有着广泛的应用，例如高温变形、热处理和材料改性等。
蠕变变形是一种重要的材料塑性变形现象，其机理涉及晶体滑移、晶粒边界滑移和扩散等多个方面。
蠕变变形的发生受到温度、应力、时间和材料微观结构等因素的影响。
蠕变变形在航空航天、能源和材料加工等领域具有重要的应用价值。
随着科学技术的不断发展，对蠕变变形的研究将进一步深化，为材料设计和工程应用提供更好的支持和指导。
纤维的蠕变名词解释
纤维的蠕变名词解释纤维作为一种常见的材料，在我们日常生活中无处不在。
无论是衣物、被褥还是家具等物品，都离不开这种柔软而坚韧的材料。
纤维的蠕变是指其在使用过程中经历的形态或性质的变化。
本文将对纤维的蠕变进行解释和探讨。
一、纤维的蠕变现象纤维的蠕变是一种渐进且不可逆的变化。
这种变化可以是在纤维使用后，例如经过一段时间的拉伸或压力作用后引起的形状畸变，也可以是由于气候、湿度、温度等环境因素引起的纤维性能的变化。
在使用过程中，纤维常常会受到拉伸或压力的影响。
例如，一条弹性纤维拉伸后会出现延伸的现象，保持一段时间后，它可能不会完全恢复到最初的形状。
这是由于纤维分子在拉伸后发生了一定的内部重排，导致了物理结构的变化。
重复的拉伸过程会进一步加剧纤维的蠕变现象。
二、纤维的蠕变类型纤维的蠕变可以分为两种类型：弹性蠕变和塑性蠕变。
1. 弹性蠕变：这种蠕变是指纤维在受到外部拉伸或压力后能够部分恢复原来形态的现象。
弹性蠕变是可逆的，也就是说，一旦外部力消失，纤维就可以恢复到原来的形状。
这种蠕变通常发生在高弹性纤维中，如弹性织物中的弹性纤维。
2. 塑性蠕变：这种蠕变是指纤维在受到外部拉伸或压力后，无法完全恢复原来的形态，产生持久的形状变化。
塑性蠕变是不可逆的，一旦发生了蠕变，就无法恢复到原来的状态。
这种蠕变通常发生在非弹性纤维中，如棉纤维或人造纤维。
三、纤维蠕变的原因纤维蠕变的原因主要是由于纤维分子内部结构的变化，这些变化可以发生在微观和宏观层面。
1. 纤维分子的排列：纤维分子在受到拉伸或压力时会发生重排，导致纤维的形状和性质的变化。
拉伸时，纤维分子会被拉伸并平行排列，这会导致纤维的延伸。
而压力作用下，纤维分子会受到挤压，从而引起纤维的短缩。
2. 纤维结构的变化：纤维在使用过程中会受到多种环境因素的影响，如湿度、温度等。
这些因素会改变纤维的物理和化学性质，导致其结构的变化。
例如，高温可能使纤维分子与空气中的分子发生反应，从而导致纤维的硬化或变脆。
蠕变及其影响因素
蠕变及其影响因素1.蠕变钢在长时间的恒温、恒载作用下，发生缓慢的塑性变形现象称为蠕变。
蠕变是一种由热作用引起的，金属的晶格发生位错，点阵中的原子自身扩散（位置交换或位置迁移）形成了蠕变，蠕变可以在单一应力（拉力、压力或扭拒），也可在复合应力下产生。
蠕变曲线能够描述钢在一定温度、应力作用下蠕变整个变形过程，如图23-1所示。
Oa 为开始加载后所引起的瞬时变形；ab为蠕变的第一阶段，这个阶段蠕变速度随时间的增加而逐渐减少（称减速阶段或不稳定阶段）；cd为蠕变第三阶段，在这个阶段中蠕变加速进行（称加速阶段或最后阶段），直至d点断裂。
2.对蠕变行为的影响因素2.1通过固溶强化对晶格造成约束在热强钢中加入Mo、Mn、W、Cr、等元素实现固溶强化，增强了固溶体原子间结合力和晶格畸变，提高蠕变抗力和持久强度。
如低碳钢的工作温度一般在450~480℃，当加入0.5%Mo（0.5Mo钢），最高工作温度可达500℃左右。
为防止高温、长期运行会生产石墨化，加入Cr元素，同时也提高钢的抗氧化性。
2.2通过减少钢中的有害伴生元素净化晶界在高温、长时间承受应力时，晶界也参与变形，当变形速度越慢，晶界变形的比例越大。
这是由于晶界处原子排列不规则，位错和空位多，S、P等有害杂质易在晶界偏析聚集，造成晶界热强性很低，因此晶界是高温条件下的薄弱环节。
热强钢应严格限制杂质元素，选用优质钢和特殊优质钢实现净化晶界，同时还应加入B、Zr等微量元素减少晶界缺陷，提高晶界强度。
2.3通过正火处理使晶粒分布均匀对于常温力学性能来说，一般是晶粒越细小则强度和硬度越高，同时塑性和韧性也越好。
在高温条件下，原子沿晶界的扩散速度比晶粒内大得多，晶界成为最薄弱的部位，希望得到适中的晶粒度心减少晶界面积。
2.4通过弥散强化阻碍位错动力加入V、Ti、Nb等元素，形成（V4C3、TiC、NbC）高温时和稳定且不易聚集长大的碳化物相，析出的碳化物呈细小弥散状、均匀地分布在晶料的滑移面上，阻碍位错的运动，达到高温强化目的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
ANSYS提供了两个用户徐变方程：USERCR.F和USERCREEP.F。
其中：显式徐变用USERCR.F；前提是C6 = 100隐式徐变用USERCREEP.F，前提是TBOPT＝100（1）用户徐变子程序usercr，用于显式徐变subroutine usercr (elem,intpt,mat,ncomp,kfirst,kfsteq,e,posn,d,x proptb,timval,timinc,tem,dtem,toffst,fluen,dfluen,epel,epcrp,x statev,usvr,delcr)cc *** 基本功能: 允许用户写自己的徐变规律。
该逻辑仅在C6＝100时可用。
c *** 次要功能: 演示用户徐变方程的编写c *** 注意－本文件包含ANSYS 机要信息***c *** ansys(r) copyright(c) 2000c *** ansys, inc.cc 输入变量:c | （译者注）c |c | 类型：int－整型c | dp－双精度型c | 长度：sc－标量c | ar( , )－数组c | 目的：in－输入c | out－输出c | inout－输入输出c 变量（类型，长度，目的）－描述c elem (int,sc,in) －单元号(标识)c intpt (int,sc,in) －单元积分点数c mat (int,sc,in) －材料引用号c ncomp (int,sc,in) －应力/应变分量数(1,4 or 6)c 1 －xc 4 －x,y,z,xyc 6 －x,y,z,xy,yz,xzc kfirst (int,sc,in) －若是首次则值为1，否则为0c (对把状态变量初始化为非零值有用)c kfsteq (int,sc,in) －若是子步中首次平衡迭代则为1，否则为0 c e (dp,sc,in) －杨氏弹性模量c posn (dp,sc,in) －泊松比c d (dp,ar(ncomp,ncomp),in) －弹性应力－应变矩阵c proptb (dp,ar(72),in) －使用TB命令输入的材料属性值c (勿使用PROPTB(13), 因为它用在别处) c timval (dp,sc,in) －当前时间c timinc (dp,sc,in) －这一子步中的时间增量c tem (dp,sc,in) －这一子步结束时的温度c dtem (dp,sc,in) －这一子步结束时的温度增量c toffst (dp,sc,in) －相对于绝对零度的温度偏置c fluen (dp,sc,in) －这一子步结束时的流量c dfluen (dp,sc,in) －这一子步的流量增量cc epel (dp,ar(ncomp),inout) －弹性应变c epcrp (dp,ar(ncomp),inout) －前面子步的徐变应变c statev (dp,ar(ncomp*5+2),inout) －前面（已收敛）子步中的状态变量c 该变量仅用于显式徐变，并且是指一c 个内部变量，但不是用TB,stat命令定c 义的用于隐式徐变（用户徐变）以及用c 户材料的内部变量c usvr (dp,ar(nuval,nintp),inout) －附加状态变量，来自前面的迭代（使用c NSVR命令保存）cc 输出变量：c 变量（类型，长度，目的）－描述c epel (dp,ar(ncomp),inout) －考虑徐变增量后的弹性应变调整（elastic strain adjusted for creep increment）c epcrp (dp,ar(ncomp),inout) －更新后的徐变应变c statev (dp,ar(ncomp*5+2),inout) －更新后的各状态变量c usvr (dp,ar(nuval,nintp),inout) －更新后的附加状态变量c delcr (dp,sc,out) －等效徐变应变增量（用于徐变率计算）cc FORTRAN参数(由用户定义):c 变量（类型）－描述c nuval (int) －每个积分点处附加状态变量数目c nintp (int) －本子程序使用的单元最大积分点数目c （最大为14）c 注意: nuval×nintp = nstv （用nsvr 命令），不能超过840！cc 内部变量：c 变量(类型，长度) －描述c con (dp,sc) －临时变量c del (dp,ar(6)) －徐变应变增量c epet (dp,sc) －等效弹性应变(徐变前)c ept (dp,ar(6)) －总应变c eptot (dp,sc) －等效总应变，弹性+ 徐变c sigen (dp,sc) －等效应力(徐变前)c temabs (dp,sc) －绝对温度#Include "impcom.inc"external erhandlercc 用户定义的FORTRAN参数c --- usvr数据的容量integer nuval,nintpparameter (nuval=1,nintp=1)cc 外部子程序和函数external egendouble precision egenexternal vapb,vapb1,vamb1,vmultcc 整型变量integer elem,intpt,mat,ncomp,kfirst,kfsteqcc 双精度变量(&为续行标识符，原文为x，为清晰起见这里该为&。
译者注)double precision& e,posn,d(ncomp,ncomp),proptb(72),timval,timinc,tem,dtem,toffst,& fluen,dfluen,epel(ncomp),epcrp(ncomp),statev(ncomp*5+2),& usvr(nuval,nintp),delcr,& temabs,con,del(6),epet,sigen,eptot,ept(6)cc **** 当前版本的ANSYS程序不累计徐变效应（creep work）的影响，如需要，建议按如下步骤进行：c 1. 在本子程序中计算需要的徐变效应。
C 这可以在单元中累加实现。
可以这样做：c 在子程序的开头if (intpt .eq. 1) usvr(1,1) = 0.0，c 然后在所有积分点累加usvr(1,1) = usvr(1,1) + xxxxxc 注意这仅仅是徐变效应密度，如果需要总徐变效应，需要乘以体积。
这c 将在下面的步骤中阐述。
C 事实上，比这要稍微复杂一些，我们希望这个累加仅对收敛的平衡迭c 代进行，具体实现方法为：c if (kfstps .eq. 1) work = 0.0c if (kfsteq .eq. 1) workold = workc work = workold + delta(work)c 这样，在收敛的平衡迭代中，效应将是正确的c 2. 在userou中，将数据从usvr(1,1) 移到udbdat(1).c userou 仅在所有积分点都出来完毕后才被调用c 3. 在post1, 在nmisc record(各种各样的不可和记录).的末尾，c 获取udbdat(1)c 4. 乘以单元体积并累加.c 这个过程假设徐变效应密度在不规则单元区域内不迅速变化cc 这个子程序演示用基本徐变方程计算徐变应变。
更详细的资料参见ANSYS理论手册3.3。
c 在大应变分析中，所有应变都是hencky 应变. (即对数应变)c 所有量都是旋转中立（rotation-neutralized）并且是基于旋转单元坐标系统（不是全局）cc ***** 初始化*****cc ---初始检查c --- 当作者进行修改时，应删去下面的警告if (intpt.eq.1 .and. kfirst.eq.1)＆call erhandler('usercr',5000,2,＆'ANSYS,INC.- USERCR 的替代版已使用-'＆,0.0d0,' ')cc ---- 作者应删除下面一行，它仅仅是为分析者提供的con = d(1,1)+dfluen+dtem+fluen+statev(1)+usvr(1,1)+& dble(elem+kfsteq+mat)c ---检查作者定义的usvr 容量是否超界cif (nuval*nintp.gt.840)&call erhandler('usercr',5010,4,& ' USVR 超出允许的最大存储量.'& ,0.0d0,' ')cc --- 初始化徐变应变增量以防没有徐变delcr = 0.0d0c --- 没有徐变倘若时间没有变化if (timinc .le. 0.0d0 .and. timval .le. 0.0d0) go to 999c --- 没有徐变倘若温度未定义temabs = tem + toffstif (temabs .le. 0.0d0) thenccall erhandler('usercr',5020,3,& '徐变Temperature= %G 应大于绝对零度.'& ,temabs,' ')cgo to 999endifcc ***** 定义等效应变和应力*****cc --- 利用函数egen定义等效应变epet = egen (ncomp,epel(1),posn)c --- 若应变为0，则无徐变if (epet .eq. 0.0d0) go to 999c --- 定义应力sigen = e*epetccc ***** 定义徐变应变率*****cc ***** 用户修改正式开始c ***** 开始徐变方程的例子- same as c6 = 0cc --- 若c1为0则跳过if (proptb(1) .eq. 0.0d0) go to 999cc ---定义等效总应变call vapb (epel(1),epcrp(1),ept(1),ncomp)eptot = egen (ncomp,ept(1),posn)cc --- 定义徐变应变率if (eptot .gt. 0.0d0) thendelcr = delcr + exp( log(proptb(1)) + proptb(2)*log(sigen) + & proptb(3)*log(eptot) - proptb(4)/temabs)endifcc ***** 徐变方程的例子到此结束c ***** 用户修改正常结束cc ***** 为每一个分量计算徐变应变增量*****cc --- 根据徐变应变率计算徐变应变增量delcr = delcr*timinccc --- 应用prandtl-reuss 关系为每一分量计算增量if (ncomp .eq. 1) thendel(1) = delcr*epel(1)/epetelsecon = delcr/epetcon = con/(2.0d0*(1.0d0+posn))call vmult (epel(4),del(4),ncomp-3,3.0d0*con)del(1) = con*(2.0d0*epel(1) - epel(2) - epel(3))del(2) = con*(2.0d0*epel(2) - epel(3) - epel(1))del(3) = con*(2.0d0*epel(3) - epel(1) - epel(2))endifcc --- 更新应变call vapb1 (epcrp(1),del(1),ncomp)call vamb1 (epel(1),del(1),ncomp)c999 returnend＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝（2）用户徐变子程序usercreep，用于隐式徐变SUBROUTINE usercreep (impflg, ldstep, isubst, matId , elemId,& kDInPt, kLayer, kSecPt, nstatv, nprop,& prop , time , dtime , temp , dtemp ,& toffst, statev, creqv , pres , seqv ,& delcr , dcrda)c************************************************************************* c *** 基本功能***c 定义徐变规律，当TB,CREEP选项TBOPT=100时使用.c 演示如何补充用户徐变方程cc 徐变方程为c dotcreq := k0 * seqv ^ n * creqv ^ m * exp (-b/T)c 其中c seqv 等效有效应力(V on-Mises应力)c creqv 等效有效徐变c T 温度c k0, m, n, b 材料常数cc 该模型对应于基本徐变函数TBOPT = 1cc gal 10.01.1998cc************************************************************************* cc | （译者注）c |c | 类型：in－整型c | dp－双精度型c | 长度：sc－标量c | ar( , )－数组c | 目的：i－输入c | o－输出c | i/o－输入输出c | l－局部变量c impflg (in ,sc ,i) 显式/隐式积分标识(目前未使用)c ldstep (in ,sc ,i) 当前荷载步c isubst (in ,sc ,i) 当前子步c matId (in ,sc ,i) 材料号c elemId (in ,sc ,i) 单元号c kDInPt (in ,sc ,i) 材料积分点c kLayer (in ,sc ,i) 层号c kSecPt (in ,sc ,i) 断面点（Section point）c nstatv (in ,sc ,i) 状态变量数c nprop (in ,sc ,i) 材料属性数组长度c prop (dp ,ar(*),i) 材料属性数组c 该数组传递所有与TB,CREEP有关的TBDATA命令定义的在温度temp时的徐变常数（不要使用prop(13),因为它被用在别处）c time 当前时间c dtime 当前时间增量c temp 当前温度c dtemp 当前温度增量c toffst (dp, sc, i) 温度相对于绝对零度的偏置值c seqv (dp ,sc , i) 等效有效应力c creqv (dp ,sc , i) 等效有效徐变应变c pres (dp ,sc , i) 静水压力应力, (Sxx+Syy+Szz)/3ccc 输入输出参数c statev (dp,ar(*), i/o) 用户定义的在时间't' / 't+dt'时的内部状态变量.c 该数组将在时间增量的开始，传递这些变量的值。