{"id":261,"date":"2026-06-29T21:14:44","date_gmt":"2026-06-29T13:14:44","guid":{"rendered":"https:\/\/numsimlab.com\/?p=261"},"modified":"2026-06-29T21:18:53","modified_gmt":"2026-06-29T13:18:53","slug":"%e6%99%b6%e4%bd%93%e5%a1%91%e6%80%a7%e5%8a%9b%e5%ad%a6huang-umat%e5%ad%90%e7%a8%8b%e5%ba%8f%e8%af%a6%e8%a7%a3","status":"publish","type":"post","link":"https:\/\/numsimlab.com\/?p=261","title":{"rendered":"\u6676\u4f53\u5851\u6027\u529b\u5b66Huang UMAT\u5b50\u7a0b\u5e8f\u8be6\u89e3"},"content":{"rendered":"\n<!DOCTYPE html>\n<html lang=\"zh\">\n<head>\n<meta charset=\"UTF-8\">\n<meta name=\"viewport\" content=\"width=device-width, initial-scale=1.0\">\n<meta name=\"description\" content=\"HuangUMAT\u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u8be6\u89e3 - Detailed Explanation of Huang's UMAT Crystal Plasticity Subroutine\">\n<title data-lang=\"zh\">HuangUMAT\u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u8be6\u89e3<\/title>\n<title data-lang=\"en\" style=\"display:none\">Huang&#8217;s UMAT Crystal Plasticity Subroutine<\/title>\n<style>\n:root {\n  --bg: #f5f7fa;\n  --bg2: #eef1f6;\n  --ink: #1a2332;\n  --muted: #5a6a7e;\n  --rule: #d0d7e2;\n  --accent: #2563eb;\n  --accent2: #0ea5e9;\n  --code-bg: #1a2332;\n  --code-fg: #e2e8f0;\n}\n* { box-sizing: border-box; margin: 0; padding: 0; }\nhtml { scroll-behavior: smooth; }\nbody {\n  font-family: \"Microsoft YaHei\", \"PingFang SC\", system-ui, -apple-system, sans-serif;\n  background: var(--bg);\n  color: var(--ink);\n  line-height: 1.75;\n  font-size: 16px;\n}\nbody.en-mode {\n  font-family: system-ui, -apple-system, \"Segoe UI\", sans-serif;\n}\n.lang-bar {\n  position: sticky;\n  top: 0;\n  z-index: 1000;\n  background: linear-gradient(90deg, #1a2332, #2c3e5a);\n  padding: 10px 20px;\n  display: flex;\n  justify-content: center;\n  gap: 12px;\n  box-shadow: 0 2px 8px rgba(0,0,0,0.15);\n}\n.lang-bar button {\n  background: rgba(255,255,255,0.12);\n  border: 1px solid rgba(255,255,255,0.25);\n  color: #fff;\n  padding: 6px 18px;\n  border-radius: 20px;\n  cursor: pointer;\n  font-size: 14px;\n  transition: all 0.25s ease;\n}\n.lang-bar button:hover {\n  background: rgba(255,255,255,0.25);\n  transform: translateY(-1px);\n}\n.lang-bar button.active {\n  background: var(--accent);\n  border-color: var(--accent);\n  box-shadow: 0 2px 8px rgba(37,99,235,0.35);\n}\n.container {\n  max-width: 1100px;\n  margin: 0 auto;\n  padding: 40px 24px;\n}\n.cover {\n  text-align: center;\n  padding: 60px 20px 50px;\n  background: linear-gradient(135deg, var(--bg2), var(--bg));\n  border-bottom: 1px solid var(--rule);\n  margin-bottom: 40px;\n}\n.cover h1 {\n  font-size: 2.4rem;\n  font-weight: 700;\n  color: var(--ink);\n  margin-bottom: 12px;\n  letter-spacing: -0.5px;\n}\n.cover h2 {\n  font-size: 1.5rem;\n  font-weight: 400;\n  color: var(--muted);\n  margin-bottom: 8px;\n}\n.cover .subtitle {\n  font-size: 1.05rem;\n  color: var(--muted);\n  margin-top: 16px;\n  font-style: italic;\n}\n.cover .meta {\n  margin-top: 24px;\n  font-size: 0.95rem;\n  color: var(--muted);\n}\n.toc {\n  background: #fff;\n  border: 1px solid var(--rule);\n  border-radius: 12px;\n  padding: 32px;\n  margin-bottom: 40px;\n  box-shadow: 0 2px 12px rgba(0,0,0,0.04);\n}\n.toc h2 {\n  font-size: 1.4rem;\n  margin-bottom: 20px;\n  color: var(--ink);\n  border-bottom: 2px solid var(--accent);\n  padding-bottom: 8px;\n  display: inline-block;\n}\n.toc ul {\n  list-style: none;\n  columns: 2;\n  column-gap: 40px;\n}\n.toc li {\n  margin-bottom: 8px;\n  break-inside: avoid;\n}\n.toc a {\n  color: var(--accent);\n  text-decoration: none;\n  font-size: 0.95rem;\n}\n.toc a:hover { text-decoration: underline; }\n.chapter {\n  background: #fff;\n  border: 1px solid var(--rule);\n  border-radius: 12px;\n  padding: 36px;\n  margin-bottom: 32px;\n  box-shadow: 0 2px 12px rgba(0,0,0,0.04);\n}\n.chapter h2 {\n  font-size: 1.6rem;\n  color: var(--ink);\n  margin-bottom: 6px;\n}\n.chapter .chap-num {\n  display: inline-block;\n  background: var(--accent);\n  color: #fff;\n  padding: 2px 10px;\n  border-radius: 4px;\n  font-size: 0.85rem;\n  margin-right: 8px;\n  vertical-align: middle;\n}\n.chapter h3 {\n  font-size: 1.25rem;\n  color: var(--ink);\n  margin: 28px 0 14px;\n  padding-left: 12px;\n  border-left: 4px solid var(--accent2);\n}\n.chapter h4 {\n  font-size: 1.1rem;\n  color: var(--ink);\n  margin: 20px 0 10px;\n}\n.chapter p {\n  margin-bottom: 12px;\n  text-align: justify;\n}\n.chapter ul, .chapter ol {\n  margin: 10px 0 10px 24px;\n}\n.chapter li {\n  margin-bottom: 6px;\n}\ntable {\n  width: 100%;\n  border-collapse: collapse;\n  margin: 16px 0;\n  font-size: 0.92rem;\n  background: #fff;\n}\nth, td {\n  border: 1px solid var(--rule);\n  padding: 10px 12px;\n  text-align: left;\n}\nth {\n  background: var(--bg2);\n  font-weight: 600;\n  color: var(--ink);\n}\ntr:nth-child(even) { background: var(--bg); }\ntr:hover { background: #e8edf5; }\npre {\n  background: var(--code-bg);\n  color: var(--code-fg);\n  padding: 18px;\n  border-radius: 8px;\n  overflow-x: auto;\n  font-family: Consolas, \"Courier New\", monospace;\n  font-size: 0.88rem;\n  line-height: 1.6;\n  margin: 16px 0;\n  position: relative;\n}\ncode {\n  font-family: Consolas, \"Courier New\", monospace;\n  background: #e8edf5;\n  padding: 2px 6px;\n  border-radius: 4px;\n  font-size: 0.9em;\n  color: var(--ink);\n}\npre code {\n  background: none;\n  padding: 0;\n  color: var(--code-fg);\n  font-size: inherit;\n}\n.formula {\n  background: var(--bg2);\n  border-left: 4px solid var(--accent);\n  padding: 14px 18px;\n  margin: 16px 0;\n  border-radius: 0 8px 8px 0;\n  font-family: \"Times New Roman\", serif;\n  font-size: 1.05rem;\n  overflow-x: auto;\n}\n.mermaid {\n  background: #fff;\n  border: 1px solid var(--rule);\n  border-radius: 8px;\n  padding: 20px;\n  margin: 16px 0;\n  text-align: center;\n}\n.note {\n  background: #fff9e6;\n  border-left: 4px solid #f59e0b;\n  padding: 14px 18px;\n  margin: 16px 0;\n  border-radius: 0 8px 8px 0;\n  color: #713f12;\n}\n.warning {\n  background: #fef2f2;\n  border-left: 4px solid #ef4444;\n  padding: 14px 18px;\n  margin: 16px 0;\n  border-radius: 0 8px 8px 0;\n  color: #7f1d1d;\n}\n.ref-item {\n  margin-bottom: 10px;\n  padding-left: 24px;\n  text-indent: -24px;\n}\n@media (max-width: 768px) {\n  .container { padding: 20px 14px; }\n  .cover h1 { font-size: 1.6rem; }\n  .cover h2 { font-size: 1.15rem; }\n  .chapter { padding: 22px; }\n  .chapter h2 { font-size: 1.3rem; }\n  .toc ul { columns: 1; }\n  table { font-size: 0.82rem; }\n  th, td { padding: 6px 8px; }\n}\n.hljs-keyword { color: #c678dd; }\n.hljs-string { color: #98c379; }\n.hljs-number { color: #d19a66; }\n.hljs-comment { color: #5c6370; font-style: italic; }\n.hljs-function { color: #61afef; }\n.hljs-params { color: #e06c75; }\n.hljs-built_in { color: #e5c07b; }\n<\/style>\n<script>\nfunction switchLanguage(lang) {\n  document.documentElement.lang = lang;\n  document.querySelectorAll('[data-lang]').forEach(function(el) {\n    if (el.tagName.toLowerCase() === 'title') {\n      if (el.getAttribute('data-lang') === lang) el.style.display = '';\n      else el.style.display = 'none';\n      return;\n    }\n    if (el.getAttribute('data-lang') === lang) {\n      el.style.display = '';\n    } else {\n      el.style.display = 'none';\n    }\n  });\n  document.querySelectorAll('.lang-bar button').forEach(function(btn) {\n    btn.classList.toggle('active', btn.getAttribute('data-set-lang') === lang);\n  });\n  if (lang === 'en') document.body.classList.add('en-mode');\n  else document.body.classList.remove('en-mode');\n  if (window.mermaid) {\n    mermaid.initialize({ startOnLoad: true, theme: 'neutral', securityLevel: 'loose' });\n  }\n}\ndocument.addEventListener('DOMContentLoaded', function() {\n  switchLanguage('zh');\n});\n<\/script>\n<\/head>\n<body>\n\n<div class=\"lang-bar\">\n  <button data-set-lang=\"zh\" onclick=\"switchLanguage('zh')\" class=\"active\">\n    <span data-lang=\"zh\">\u4e2d\u6587<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Chinese<\/span>\n  <\/button>\n  <button data-set-lang=\"en\" onclick=\"switchLanguage('en')\">\n    <span data-lang=\"zh\">English<\/span>\n    <span data-lang=\"en\" style=\"display:none\">English<\/span>\n  <\/button>\n<\/div>\n\n<div class=\"container\">\n\n<!-- ==================== COVER ==================== -->\n<div class=\"cover\">\n  <h1>\n    <span data-lang=\"zh\">Huang UMAT \u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u8be6\u89e3<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Detailed Explanation of Huang&#8217;s UMAT Crystal Plasticity Subroutine<\/span>\n  <\/h1>\n  <h2>\n    <span data-lang=\"zh\">Crystal Plasticity Finite Element Method (CPFEM) \/ \u6676\u4f53\u5851\u6027\u6709\u9650\u5143\u65b9\u6cd5<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Crystal Plasticity Finite Element Method (CPFEM)<\/span>\n  <\/h2>\n  <div class=\"subtitle\">\n    <span data-lang=\"zh\">\u57fa\u4e8e Abaqus\/Standard \u9690\u5f0f\u6c42\u89e3\u5668<\/span>\n    <span data-lang=\"en\" style=\"display:none\">For Abaqus\/Standard Implicit Solver<\/span>\n  <\/div>\n  <div class=\"meta\">\n    <span data-lang=\"zh\">\u5b8c\u6574\u672c\u6784\u7406\u8bba\u3001UMAT\u63a5\u53e3\u3001\u4ee3\u7801\u7ed3\u6784\u4e0e\u5b9e\u73b0\u7ec6\u8282<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Complete Constitutive Theory, UMAT Interface, Code Structure and Implementation Details<\/span>\n  <\/div>\n<\/div>\n\n\n\n<!-- ==================== TOC ==================== -->\n<div class=\"toc\" id=\"toc\">\n  <h2>\n    <span data-lang=\"zh\">\u76ee\u5f55 \/ Table of Contents<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Table of Contents<\/span>\n  <\/h2>\n  <ul>\n    <li><a href=\"#ch1\">\n      <span data-lang=\"zh\">Chapter 1 \u6982\u8ff0 \/ Overview<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 1 Overview<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch2\">\n      <span data-lang=\"zh\">Chapter 2 \u672c\u6784\u7406\u8bba \/ Constitutive Theory<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 2 Constitutive Theory<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch3\">\n      <span data-lang=\"zh\">Chapter 3 Abaqus UMAT\u63a5\u53e3 \/ UMAT Interface<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 3 Abaqus UMAT Interface<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch4\">\n      <span data-lang=\"zh\">Chapter 4 \u6750\u6599\u53c2\u6570 \/ Material Parameters<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 4 Material Parameters<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch5\">\n      <span data-lang=\"zh\">Chapter 5 \u4ee3\u7801\u7ed3\u6784 \/ Code Structure<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 5 Code Structure<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch6\">\n      <span data-lang=\"zh\">Chapter 6 Jacobian\u77e9\u9635 \/ Jacobian Matrix<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 6 Jacobian Matrix<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch7\">\n      <span data-lang=\"zh\">Chapter 7 \u5178\u578b\u53c2\u6570 \/ Typical Parameters<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 7 Typical Parameters<\/span>\n    <\/a><\/li>\n    <li><a href=\"#ch8\">\n      <span data-lang=\"zh\">Chapter 8 \u4f7f\u7528\u6ce8\u610f\u4e8b\u9879 \/ Usage Notes<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Chapter 8 Usage Notes<\/span>\n    <\/a><\/li>\n    <li><a href=\"#refs\">\n      <span data-lang=\"zh\">\u53c2\u8003\u6587\u732e \/ References<\/span>\n      <span data-lang=\"en\" style=\"display:none\">References<\/span>\n    <\/a><\/li>\n  <\/ul>\n<\/div>\n\n<!-- ==================== CHAPTER 1 ==================== -->\n<div class=\"chapter\" id=\"ch1\">\n  <h2>\n    <span class=\"chap-num\">01<\/span>\n    <span data-lang=\"zh\">\u6982\u8ff0 \/ Overview<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Overview<\/span>\n  <\/h2>\n\n  <h3 id=\"ch1-1\">\n    <span data-lang=\"zh\">1.1 \u672c\u6784\u6a21\u578b\u7b80\u4ecb \/ Constitutive Model Introduction<\/span>\n    <span data-lang=\"en\" style=\"display:none\">1.1 Constitutive Model Introduction<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    Huang\uff08Yonggang Huang\uff09\u6559\u6388\u5f00\u53d1\u7684UMAT\u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u662f\u6676\u4f53\u5851\u6027\u6709\u9650\u5143\u65b9\u6cd5\uff08CPFEM\uff09\u9886\u57df\u6700\u5e7f\u6cdb\u4f7f\u7528\u7684\u5f00\u6e90\u5b9e\u73b0\u4e4b\u4e00\u3002\n    \u8be5\u5b50\u7a0b\u5e8f\u57fa\u4e8eAbaqus\/Standard\u7528\u6237\u6750\u6599\u5b50\u7a0b\u5e8f\uff08UMAT\uff09\u63a5\u53e3\uff0c\u91c7\u7528\u9690\u5f0f\u65f6\u95f4\u79ef\u5206\u65b9\u6848\uff0c\n    \u80fd\u591f\u51c6\u786e\u6a21\u62df\u5355\u6676\u53ca\u591a\u6676\u6750\u6599\u5728\u590d\u6742\u52a0\u8f7d\u6761\u4ef6\u4e0b\u7684\u5851\u6027\u53d8\u5f62\u884c\u4e3a\u3002\n    \u8be5\u6a21\u578b\u8003\u8651\u4e86\u6676\u4f53\u7684\u5404\u5411\u5f02\u6027\u5f39\u6027\u3001\u6676\u4f53\u5b66\u6ed1\u79fb\u7cfb\u4e0a\u7684\u526a\u5207\u53d8\u5f62\uff0c\u4ee5\u53ca\u4f4d\u9519\u4ea4\u4e92\u4f5c\u7528\u5bfc\u81f4\u7684\u786c\u5316\u6548\u5e94\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The UMAT single-crystal plasticity subroutine developed by Professor Yonggang Huang\n    is one of the most widely used open-source implementations in the field of Crystal Plasticity Finite Element Method (CPFEM).\n    Based on the Abaqus\/Standard User Material (UMAT) interface with an implicit time integration scheme,\n    it accurately simulates the plastic deformation behavior of single-crystal and polycrystalline materials under complex loading conditions.\n    The model accounts for anisotropic elasticity, crystallographic slip system shear deformation,\n    and hardening effects due to dislocation interactions.\n  <\/p>\n  <p data-lang=\"zh\">\n    \u672c\u7a0b\u5e8f\u7684\u6838\u5fc3\u7279\u5f81\u5305\u62ec\uff1a\u57fa\u4e8e\u4e2d\u95f4\u6784\u578b\u7684\u4e58\u6cd5\u5206\u89e3\u8fd0\u52a8\u5b66\u6846\u67b6\u3001Schmid\u5b9a\u5f8b\u9a71\u52a8\u7684\u6ed1\u79fb\u7cfb\u6fc0\u6d3b\u5224\u636e\u3001\n    \u7387\u76f8\u5173\u5e42\u5f8b\u6d41\u52a8\u6cd5\u5219\u3001\u53ef\u9009\u62e9\u7684\u786c\u5316\u6a21\u578b\uff08\u53cc\u66f2\u6b63\u5272\u786c\u5316\u6216Bassani-Wu\u786c\u5316\uff09\uff0c\n    \u4ee5\u53ca\u5b8c\u6574\u7684\u9690\u5f0f\u5e94\u529b\u66f4\u65b0\u7b97\u6cd5\u548c\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635\uff08Jacobian\uff09\u8ba1\u7b97\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Key features of this subroutine include: multiplicative decomposition kinematics based on the intermediate configuration,\n    Schmid&#8217;s law driven slip system activation criterion, rate-dependent power-law flow rule,\n    selectable hardening models (hyperbolic secant hardening or Bassani-Wu hardening),\n    and a complete implicit stress update algorithm with consistent tangent stiffness matrix (Jacobian) computation.\n  <\/p>\n\n  <h3 id=\"ch1-2\">\n    <span data-lang=\"zh\">1.2 \u7406\u8bba\u57fa\u7840 \/ Theoretical Foundation<\/span>\n    <span data-lang=\"en\" style=\"display:none\">1.2 Theoretical Foundation<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u9ec4\u6c38\u521aUMAT\u7684\u7406\u8bba\u57fa\u7840\u4e3b\u8981\u6765\u6e90\u4e8ePeirce\u3001Asaro\u548cNeedleman\uff081983\uff09\u4ee5\u53caPeirce\u3001Shih\u548cNeedleman\uff081984\uff09\u7684\u7ecf\u5178\u5de5\u4f5c\u3002\n    \u8fd9\u4e9b\u6587\u732e\u5efa\u7acb\u4e86\u6676\u4f53\u5851\u6027\u7406\u8bba\u5728\u73b0\u4ee3\u6709\u9650\u5143\u6846\u67b6\u4e0b\u7684\u6570\u5b66\u8868\u8ff0\uff0c\u5305\u62ec\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The theoretical foundation of Huang&#8217;s UMAT is primarily derived from the seminal works of\n    Peirce, Asaro and Needleman (1983) and Peirce, Shih and Needleman (1984).\n    These papers established the mathematical formulation of crystal plasticity theory within the modern finite element framework, including:\n  <\/p>\n  <ul>\n    <li>\n      <span data-lang=\"zh\">\u603b\u53d8\u5f62\u68af\u5ea6F\u7684\u4e58\u6cd5\u5206\u89e3\u4e3a\u5f39\u6027\u90e8\u5206Fe\u548c\u5851\u6027\u90e8\u5206Fp<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Multiplicative decomposition of the total deformation gradient F into elastic Fe and plastic Fp parts<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u57fa\u4e8eSchmid\u5206\u89e3\u7684\u5206\u89e3\u5e94\u529b\uff08resolved shear stress\uff09\u6982\u5ff5<\/span>\n      <span data-lang=\"en\" style=\"display:none\">The concept of resolved shear stress based on Schmid decomposition<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u7387\u76f8\u5173\u6d41\u52a8\u6cd5\u5219\uff0c\u907f\u514d\u7387\u65e0\u5173\u6a21\u578b\u7684\u975e\u552f\u4e00\u6027\u548c\u6570\u503c\u4e0d\u7a33\u5b9a\u6027<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Rate-dependent flow rule to avoid non-uniqueness and numerical instability of rate-independent models<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u81ea\u786c\u5316\uff08self-hardening\uff09\u548c\u6f5c\u786c\u5316\uff08latent hardening\uff09\u7684\u7edf\u4e00\u63cf\u8ff0<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Unified description of self-hardening and latent hardening<\/span>\n    <\/li>\n  <\/ul>\n  <p data-lang=\"zh\">\n    \u540e\u7eedAsaro\u548cNeedleman\uff081985\uff09\u4ee5\u53caHarren\u3001Deve\u548cAsaro\uff081988\uff09\u7684\u5de5\u4f5c\u8fdb\u4e00\u6b65\u5b8c\u5584\u4e86\u6676\u4f53\u5851\u6027\u7406\u8bba\u5728\u6709\u9650\u53d8\u5f62\u6761\u4ef6\u4e0b\u7684\u6570\u503c\u5b9e\u73b0\u3002\n    \u9ec4\u6c38\u521a\u5728\u5176\u5b50\u7a0b\u5e8f\u4e2d\u6574\u5408\u4e86\u8fd9\u4e9b\u7406\u8bba\u6210\u679c\uff0c\u5e76\u9488\u5bf9Abaqus\/Standard\u7684\u9690\u5f0f\u6c42\u89e3\u5668\u8fdb\u884c\u4e86\u7b97\u6cd5\u4f18\u5316\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Subsequent works by Asaro and Needleman (1985) and Harren, Deve and Asaro (1988) further refined\n    the numerical implementation of crystal plasticity theory under finite deformation conditions.\n    Huang integrated these theoretical achievements into his subroutine and optimized the algorithm\n    specifically for the Abaqus\/Standard implicit solver.\n  <\/p>\n\n  <h3 id=\"ch1-3\">\n    <span data-lang=\"zh\">1.3 \u9002\u7528\u8303\u56f4 \/ Scope<\/span>\n    <span data-lang=\"en\" style=\"display:none\">1.3 Scope<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u672c\u5b50\u7a0b\u5e8f\u9002\u7528\u4e8e\u4ee5\u4e0b\u7814\u7a76\u573a\u666f\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    This subroutine is applicable to the following research scenarios:\n  <\/p>\n  <ul>\n    <li>\n      <span data-lang=\"zh\">FCC\u3001BCC\u548cHCP\u6676\u4f53\u7ed3\u6784\u7684\u5355\u6676\u548c\u591a\u6676\u91d1\u5c5e\u6750\u6599\u7684\u529b\u5b66\u884c\u4e3a\u6a21\u62df<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Mechanical behavior simulation of single-crystal and polycrystalline metals with FCC, BCC, and HCP crystal structures<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u5355\u8c03\u52a0\u8f7d\u3001\u5faa\u73af\u52a0\u8f7d\u53ca\u5e94\u53d8\u8def\u5f84\u53d8\u5316\u6761\u4ef6\u4e0b\u7684\u5851\u6027\u53d8\u5f62\u5206\u6790<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Plastic deformation analysis under monotonic, cyclic, and strain-path-change loading conditions<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u7ec7\u6784\u6f14\u5316\uff08texture evolution\uff09\u548c\u5851\u6027\u5404\u5411\u5f02\u6027\u53d1\u5c55\u7684\u5b9a\u91cf\u9884\u6d4b<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Quantitative prediction of texture evolution and plastic anisotropy development<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u5c40\u90e8\u5316\u53d8\u5f62\u3001\u526a\u5207\u5e26\u5f62\u6210\u53ca\u7ec6\u89c2\u5c3a\u5ea6\u635f\u4f24\u8d77\u59cb\u7814\u7a76<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Studies of localized deformation, shear band formation, and mesoscale damage initiation<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u5fae\u67f1\u538b\u7f29\u3001\u7eb3\u7c73\u538b\u75d5\u7b49\u5c0f\u5c3a\u5ea6\u5355\u6676\u5b9e\u9a8c\u7684\u6709\u9650\u5143\u5efa\u6a21<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Finite element modeling of micro-pillar compression, nanoindentation, and other small-scale single-crystal experiments<\/span>\n    <\/li>\n  <\/ul>\n  <div class=\"note\">\n    <span data-lang=\"zh\"><strong>\u6ce8\u610f\uff1a<\/strong>\u672c\u5b50\u7a0b\u5e8f\u57fa\u4e8e\u5c0f\u5f39-\u5927\u5851\u5047\u8bbe\uff0c\u5373\u5f39\u6027\u53d8\u5f62\u68af\u5ea6Fe\u63a5\u8fd1\u5355\u4f4d\u5f20\u91cf\uff08||Fe-I|| << 1\uff09\uff0c\u4f46\u5851\u6027\u53d8\u5f62Fp\u53ef\u4ee5\u4efb\u610f\u5927\u3002\u8be5\u5047\u8bbe\u5bf9\u4e8e\u5927\u591a\u6570\u91d1\u5c5e\u6676\u4f53\u5728\u4e2d\u7b49\u5e94\u53d8\u8303\u56f4\u5185\u662f\u5408\u7406\u7684\u3002<\/span>\n    <span data-lang=\"en\" style=\"display:none\"><strong>Note:<\/strong> This subroutine is based on the small-elastic-large-plastic assumption, where the elastic deformation gradient Fe is close to the identity tensor (||Fe-I|| << 1), but the plastic deformation Fp can be arbitrarily large. This assumption is reasonable for most metallic crystals at moderate strain ranges.<\/span>\n  <\/div>\n<\/div>\n\n<!-- ==================== CHAPTER 2 ==================== -->\n<div class=\"chapter\" id=\"ch2\">\n  <h2>\n    <span class=\"chap-num\">02<\/span>\n    <span data-lang=\"zh\">\u672c\u6784\u7406\u8bba \/ Constitutive Theory<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Constitutive Theory<\/span>\n  <\/h2>\n\n  <h3 id=\"ch2-1\">\n    <span data-lang=\"zh\">2.1 \u5206\u89e3\u8fd0\u52a8\u5b66 \/ Kinematics of Decomposition<\/span>\n    <span data-lang=\"en\" style=\"display:none\">2.1 Kinematics of Decomposition<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u6676\u4f53\u5851\u6027\u7406\u8bba\u91c7\u7528\u4e58\u6cd5\u5206\u89e3\uff08multiplicative decomposition\uff09\u5c06\u603b\u53d8\u5f62\u68af\u5ea6F\u5206\u89e3\u4e3a\u5f39\u6027\u90e8\u5206Fe\u548c\u5851\u6027\u90e8\u5206Fp\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Crystal plasticity theory employs multiplicative decomposition to separate the total deformation gradient F into elastic Fe and plastic Fp parts:\n  <\/p>\n  <div class=\"formula\">\n    <b>F = F<sup>e<\/sup> \u00b7 F<sup>p<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5176\u4e2d\uff0c\u5851\u6027\u53d8\u5f62\u68af\u5ea6Fp\u7531\u6240\u6709\u6fc0\u6d3b\u6ed1\u79fb\u7cfb\u4e0a\u7684\u526a\u5207\u5e94\u53d8\u7d2f\u79ef\u4ea7\u751f\uff0c\u4fdd\u6301\u6676\u683c\u7ed3\u6784\u4e0d\u53d8\uff08\u7b49\u5bb9\u5851\u6027\u53d8\u5f62\uff09\uff1b\n    \u5f39\u6027\u53d8\u5f62\u68af\u5ea6Fe\u63cf\u8ff0\u6676\u683c\u7684\u5f39\u6027\u7578\u53d8\u548c\u521a\u6027\u8f6c\u52a8\u3002\n    \u8be5\u5206\u89e3\u5c06\u53c2\u8003\u6784\u578b\uff08reference configuration\uff09\u6620\u5c04\u5230\u4e2d\u95f4\u6784\u578b\uff08intermediate configuration\uff09\uff0c\n    \u518d\u6620\u5c04\u5230\u5f53\u524d\u6784\u578b\uff08current configuration\uff09\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Here, the plastic deformation gradient Fp is produced by the accumulation of shear strains on all active slip systems,\n    preserving the crystal lattice structure (isochoric plastic deformation);\n    the elastic deformation gradient Fe describes the elastic distortion and rigid rotation of the lattice.\n    This decomposition maps from the reference configuration to the intermediate configuration,\n    and then to the current configuration.\n  <\/p>\n  <p data-lang=\"zh\">\n    \u5851\u6027\u53d8\u5f62\u68af\u5ea6\u7684\u6f14\u5316\u65b9\u7a0b\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The evolution equation for the plastic deformation gradient is:\n  <\/p>\n  <div class=\"formula\">\n    <b>F&#775;<sup>p<\/sup> = L<sup>p<\/sup> \u00b7 F<sup>p<\/sup><\/b>&nbsp;&nbsp;&nbsp;&nbsp;\n    <span data-lang=\"zh\">\u5176\u4e2d\u5851\u6027\u901f\u5ea6\u68af\u5ea6<\/span>\n    <span data-lang=\"en\" style=\"display:none\">where the plastic velocity gradient<\/span>\n    <b>L<sup>p<\/sup> = &Sigma;<sub>&alpha;<\/sub> &gamma;&#775;<sup>(&alpha;)<\/sup> s<sup>(&alpha;)<\/sup> &otimes; m<sup>(&alpha;)<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u4e0a\u5f0f\u4e2d\uff0cs<sup>(&alpha;)<\/sup>\u548cm<sup>(&alpha;)<\/sup>\u5206\u522b\u662f\u7b2c&alpha;\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u65b9\u5411\uff08slip direction\uff09\u548c\u6ed1\u79fb\u9762\u6cd5\u5411\uff08slip plane normal\uff09\uff0c\n    \u4e8c\u8005\u5b9a\u4e49\u5728\u4e2d\u95f4\u6784\u578b\u4e2d\u4e14\u4fdd\u6301\u6b63\u4ea4\u5f52\u4e00\u3002&gamma;&#775;<sup>(&alpha;)<\/sup>\u4e3a\u5bf9\u5e94\u7684\u526a\u5207\u5e94\u53d8\u7387\u3002\n    \u6ed1\u79fb\u65b9\u5411s\u548c\u6cd5\u5411m\u5728\u5f53\u524d\u6784\u578b\u4e2d\u7684\u6620\u50cf\u5206\u522b\u4e3aF<sup>e<\/sup>s\u548c(F<sup>e<\/sup>)<sup>-T<\/sup>m\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In the above equation, s<sup>(&alpha;)<\/sup> and m<sup>(&alpha;)<\/sup> are the slip direction and slip plane normal of the &alpha;-th slip system, respectively,\n    defined in the intermediate configuration and maintained orthonormal.\n    &gamma;&#775;<sup>(&alpha;)<\/sup> is the corresponding shear strain rate.\n    The images of slip direction s and normal m in the current configuration are F<sup>e<\/sup>s and (F<sup>e<\/sup>)<sup>-T<\/sup>m, respectively.\n  <\/p>\n  <p data-lang=\"zh\">\n    \u5f39\u6027\u53d8\u5f62\u68af\u5ea6F<sup>e<\/sup>\u7684\u6781\u5206\u89e3\u4e3aF<sup>e<\/sup> = R<sup>e<\/sup>U<sup>e<\/sup> = V<sup>e<\/sup>R<sup>e<\/sup>\uff0c\n    \u5176\u4e2dU<sup>e<\/sup>\u548cV<sup>e<\/sup>\u5206\u522b\u4e3a\u53f3\u3001\u5de6\u5f39\u6027\u62c9\u4f38\u5f20\u91cf\uff0cR<sup>e<\/sup>\u4e3a\u5f39\u6027\u8f6c\u52a8\u5f20\u91cf\u3002\n    Green\u5f39\u6027\u5e94\u53d8\u5f20\u91cf\u5b9a\u4e49\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The polar decomposition of the elastic deformation gradient is F<sup>e<\/sup> = R<sup>e<\/sup>U<sup>e<\/sup> = V<sup>e<\/sup>R<sup>e<\/sup>,\n    where U<sup>e<\/sup> and V<sup>e<\/sup> are the right and left elastic stretch tensors, and R<sup>e<\/sup> is the elastic rotation tensor.\n    The Green elastic strain tensor is defined as:\n  <\/p>\n  <div class=\"formula\">\n    <b>E<sup>e<\/sup> = (1\/2)(U<sup>e<\/sup><sup>2<\/sup> &#8211; I) = (1\/2)(F<sup>eT<\/sup>F<sup>e<\/sup> &#8211; I)<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5f39\u6027\u672c\u6784\u5173\u7cfb\u91c7\u7528\u5e7f\u4e49\u80e1\u514b\u5b9a\u5f8b\uff0c\u4ee5\u7b2c\u4e8c\u7c7bPiola-Kirchhoff\u5e94\u529bS<sup>e<\/sup>\uff08\u5b9a\u4e49\u5728\u4e2d\u95f4\u6784\u578b\uff09\u548cGreen\u5f39\u6027\u5e94\u53d8E<sup>e<\/sup>\u8868\u793a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The elastic constitutive relation adopts the generalized Hooke&#8217;s law, expressed in terms of the second Piola-Kirchhoff stress S<sup>e<\/sup>\n    (defined on the intermediate configuration) and Green elastic strain E<sup>e<\/sup>:\n  <\/p>\n  <div class=\"formula\">\n    <b>S<sup>e<\/sup> = C : E<sup>e<\/sup><\/b>&nbsp;&nbsp;&nbsp;&nbsp;\n    <span data-lang=\"zh\">\u5176\u4e2dC\u4e3a\u56db\u9636\u5f39\u6027\u521a\u5ea6\u5f20\u91cf<\/span>\n    <span data-lang=\"en\" style=\"display:none\">where C is the fourth-order elastic stiffness tensor<\/span>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5f53\u524d\u6784\u578b\u4e2d\u7684Kirchhoff\u5e94\u529b&tau;\u548cCauchy\u5e94\u529b&sigma;\u5206\u522b\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The Kirchhoff stress &tau; and Cauchy stress &sigma; in the current configuration are:\n  <\/p>\n  <div class=\"formula\">\n    <b>&tau; = F<sup>e<\/sup> S<sup>e<\/sup> F<sup>eT<\/sup> = J &sigma;<\/b>&nbsp;&nbsp;&nbsp;&nbsp;\n    <span data-lang=\"zh\">\u5176\u4e2dJ = det(F)\u4e3a\u4f53\u79ef\u6bd4<\/span>\n    <span data-lang=\"en\" style=\"display:none\">where J = det(F) is the volume ratio<\/span>\n  <\/div>\n\n  <h3 id=\"ch2-2\">\n    <span data-lang=\"zh\">2.2 Schmid\u5b9a\u5f8b \/ Schmid&#8217;s Law<\/span>\n    <span data-lang=\"en\" style=\"display:none\">2.2 Schmid&#8217;s Law<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    Schmid\u5b9a\u5f8b\u662f\u6676\u4f53\u5851\u6027\u7406\u8bba\u7684\u6838\u5fc3\uff0c\u5b83\u6307\u51fa\uff1a\u6ed1\u79fb\u7cfb&alpha;\u4e0a\u7684\u5206\u89e3\u526a\u5e94\u529b\uff08resolved shear stress\uff09&tau;<sup>(&alpha;)<\/sup>\u7b49\u4e8eCauchy\u5e94\u529b&sigma;\u5728\u6ed1\u79fb\u65b9\u5411-\u6ed1\u79fb\u9762\u6cd5\u5411\u6295\u5f71\u4e0a\u7684\u5206\u91cf\u3002\n    \u5728\u7a0b\u5e8f\u5b9e\u73b0\u4e2d\uff0c\u66f4\u5e38\u7528\u7684\u662f\u5bf9\u79f0\u5316\u7684Schmid\u5f20\u91cfP<sup>(&alpha;)<\/sup>\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Schmid&#8217;s law is the core of crystal plasticity theory, stating that the resolved shear stress &tau;<sup>(&alpha;)<\/sup> on slip system &alpha;\n    equals the component of Cauchy stress &sigma; projected onto the slip direction-slip plane normal pair.\n    In program implementation, the symmetrized Schmid tensor P<sup>(&alpha;)<\/sup> is more commonly used:\n  <\/p>\n  <div class=\"formula\">\n    <b>P<sup>(&alpha;)<\/sup> = (1\/2)( s<sup>(&alpha;)<\/sup> &otimes; m<sup>(&alpha;)<\/sup> + m<sup>(&alpha;)<\/sup> &otimes; s<sup>(&alpha;)<\/sup> )<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5206\u89e3\u526a\u5e94\u529b\u53ef\u8868\u793a\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The resolved shear stress can be expressed as:\n  <\/p>\n  <div class=\"formula\">\n    <b>&tau;<sup>(&alpha;)<\/sup> = &sigma; : P<sup>(&alpha;)<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u8fd9\u91cc&#8221;:&#8221;\u8868\u793a\u53cc\u70b9\u79ef\uff08\u5f20\u91cf\u7f29\u5e76\uff09\uff0c\u5373&sigma;<sub>ij<\/sub>P<sub>ij<\/sub>\u3002Schmid\u5f20\u91cfP<sup>(&alpha;)<\/sup>\u662f\u5bf9\u79f0\u5f20\u91cf\uff0c\n    \u5b83\u5c06\u5e94\u529b\u7a7a\u95f4\u4e2d\u9a71\u52a8\u7b2c&alpha;\u4e2a\u6ed1\u79fb\u7cfb\u526a\u5207\u53d8\u5f62\u7684\u65b9\u5411\u63d0\u53d6\u51fa\u6765\u3002\n    \u5f53&tau;<sup>(&alpha;)<\/sup>\u8fbe\u5230\u8be5\u6ed1\u79fb\u7cfb\u7684\u4e34\u754c\u5206\u5207\u5e94\u529b\uff08critical resolved shear stress, CRSS\uff09g<sup>(&alpha;)<\/sup>\u65f6\uff0c\u6ed1\u79fb\u5f00\u59cb\u6fc0\u6d3b\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Here &#8220;:&#8221; denotes the double dot product (tensor contraction), i.e., &sigma;<sub>ij<\/sub>P<sub>ij<\/sub>.\n    The Schmid tensor P<sup>(&alpha;)<\/sup> is symmetric and extracts the direction in stress space that drives shear deformation on the &alpha;-th slip system.\n    Slip activates when &tau;<sup>(&alpha;)<\/sup> reaches the critical resolved shear stress (CRSS) g<sup>(&alpha;)<\/sup> of that slip system.\n  <\/p>\n  <div class=\"note\">\n    <span data-lang=\"zh\"><strong>\u7269\u7406\u610f\u4e49\uff1a<\/strong>Schmid\u5b9a\u5f8b\u8868\u660e\uff0c\u65e0\u8bba\u5e94\u529b\u72b6\u6001\u591a\u4e48\u590d\u6742\uff0c\u53ea\u6709\u4f5c\u7528\u5728\u6ed1\u79fb\u9762\u548c\u6ed1\u79fb\u65b9\u5411\u4e0a\u7684\u526a\u5e94\u529b\u5206\u91cf\u624d\u80fd\u9a71\u52a8\u4f4d\u9519\u6ed1\u79fb\u3002\u8fd9\u662f\u6676\u4f53\u5851\u6027\u5404\u5411\u5f02\u6027\u7684\u6839\u672c\u6765\u6e90\u2014\u2014\u4e0d\u540c\u6ed1\u79fb\u7cfb\u5728\u4e0d\u540c\u5e94\u529b\u72b6\u6001\u4e0b\u7684\u6fc0\u6d3b\u7a0b\u5ea6\u4e0d\u540c\u3002<\/span>\n    <span data-lang=\"en\" style=\"display:none\"><strong>Physical meaning:<\/strong> Schmid&#8217;s law indicates that regardless of how complex the stress state is, only the shear stress component acting on the slip plane and in the slip direction can drive dislocation glide. This is the fundamental source of plastic anisotropy in crystals\u2014the activation degree of different slip systems varies under different stress states.<\/span>\n  <\/div>\n\n  <h3 id=\"ch2-3\">\n    <span data-lang=\"zh\">2.3 \u7387\u76f8\u5173\u6d41\u52a8\u5f8b \/ Rate-Dependent Flow Law<\/span>\n    <span data-lang=\"en\" style=\"display:none\">2.3 Rate-Dependent Flow Law<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u9ec4\u6c38\u521aUMAT\u91c7\u7528\u7387\u76f8\u5173\uff08viscoplastic\uff09\u6d41\u52a8\u6cd5\u5219\uff0c\u907f\u514d\u4e86\u7387\u65e0\u5173\u6a21\u578b\u4e2d\u6ed1\u79fb\u7cfb\u6fc0\u6d3b-\u975e\u6fc0\u6d3b\u5207\u6362\u5e26\u6765\u7684\u6570\u503c\u56f0\u96be\u3002\n    \u6bcf\u4e2a\u6ed1\u79fb\u7cfb\u4e0a\u7684\u526a\u5207\u5e94\u53d8\u7387\u7531\u4ee5\u4e0b\u5e42\u5f8b\uff08power law\uff09\u63cf\u8ff0\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Huang&#8217;s UMAT employs a rate-dependent (viscoplastic) flow rule, avoiding numerical difficulties\n    associated with the activation-deactivation switching of slip systems in rate-independent models.\n    The shear strain rate on each slip system is described by the following power law:\n  <\/p>\n  <div class=\"formula\">\n    <b>&gamma;&#775;<sup>(&alpha;)<\/sup> = &gamma;&#775;<sub>0<\/sub> &middot; | &tau;<sup>(&alpha;)<\/sup> \/ g<sup>(&alpha;)<\/sup> |<sup>n<\/sup> &middot; sign( &tau;<sup>(&alpha;)<\/sup> \/ g<sup>(&alpha;)<\/sup> )<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5f0f\u4e2d\u5404\u53c2\u6570\u542b\u4e49\u5982\u4e0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The parameters in the equation are defined as follows:\n  <\/p>\n  <ul>\n    <li>\n      <b>&gamma;&#775;<sub>0<\/sub><\/b>:\n      <span data-lang=\"zh\">\u53c2\u8003\u526a\u5207\u7387\uff08reference shear rate\uff09\uff0c\u5355\u4f4d\u4e3as<sup>-1<\/sup>\uff0c\u901a\u5e38\u53d6\u4e0e\u5b9e\u9a8c\u5e94\u53d8\u901f\u7387\u76f8\u540c\u91cf\u7ea7<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Reference shear rate, in s<sup>-1<\/sup>, typically of the same order as the experimental strain rate<\/span>\n    <\/li>\n    <li>\n      <b>n<\/b>:\n      <span data-lang=\"zh\">\u7387\u654f\u611f\u6307\u6570\uff08rate sensitivity exponent\uff09\uff0c\u65e0\u91cf\u7eb2\u3002n\u503c\u8d8a\u5927\uff0c\u6750\u6599\u884c\u4e3a\u8d8a\u63a5\u8fd1\u7387\u65e0\u5173\uff1bn&rarr;&infin;\u65f6\u9000\u5316\u4e3a\u7387\u65e0\u5173\u6a21\u578b<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Rate sensitivity exponent, dimensionless. Larger n means behavior closer to rate-independent; n&rarr;&infin; recovers the rate-independent model<\/span>\n    <\/li>\n    <li>\n      <b>&tau;<sup>(&alpha;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u7b2c&alpha;\u4e2a\u6ed1\u79fb\u7cfb\u4e0a\u7684\u5206\u89e3\u526a\u5e94\u529b<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Resolved shear stress on the &alpha;-th slip system<\/span>\n    <\/li>\n    <li>\n      <b>g<sup>(&alpha;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u7b2c&alpha;\u4e2a\u6ed1\u79fb\u7cfb\u7684\u5f53\u524d\u4e34\u754c\u5206\u5207\u5e94\u529b\uff08\u5373\u5f53\u524d\u5f3a\u5ea6\/\u786c\u5316\u72b6\u6001\uff09<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Current critical resolved shear stress (i.e., current strength\/hardening state) of the &alpha;-th slip system<\/span>\n    <\/li>\n    <li>\n      <b>sign(x)<\/b>:\n      <span data-lang=\"zh\">\u7b26\u53f7\u51fd\u6570\uff0c\u53d6\u503c\u4e3a-1\u30010\u62161<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Sign function, taking values -1, 0, or 1<\/span>\n    <\/li>\n  <\/ul>\n  <p data-lang=\"zh\">\n    \u7387\u76f8\u5173\u6a21\u578b\u7684\u4f18\u52bf\u5728\u4e8e\uff1a\u6240\u6709\u6ed1\u79fb\u7cfb\u59cb\u7ec8\u5904\u4e8e&#8221;\u90e8\u5206\u6fc0\u6d3b&#8221;\u72b6\u6001\uff0c\u4e0d\u9700\u8981\u663e\u5f0f\u7684\u6fc0\u6d3b\/\u975e\u6fc0\u6d3b\u5207\u6362\u6761\u4ef6\uff1b\n    \u8fd9\u5927\u5927\u63d0\u9ad8\u4e86\u6570\u503c\u8ba1\u7b97\u7684\u7a33\u5b9a\u6027\u548c\u6536\u655b\u6027\u3002\n    \u5728\u5178\u578b\u91d1\u5c5e\u4e2d\uff0cn\u503c\u8303\u56f4\u7ea6\u4e3a20-200\uff1aFCC\u91d1\u5c5en&asymp;20-100\uff0cBCC\u91d1\u5c5en&asymp;50-200\uff0cHCP\u91d1\u5c5en&asymp;20-80\u3002\n    n=1\u5bf9\u5e94\u7ebf\u6027\u7c98\u6027\uff08\u725b\u987f\u6d41\u4f53\uff09\uff0cn\u5f88\u5927\u65f6\u8fd1\u4f3c\u4e8e\u7406\u60f3\u5851\u6027\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The advantage of the rate-dependent model is that all slip systems are always in a &#8220;partially active&#8221; state,\n    without requiring explicit active\/inactive switching conditions;\n    this greatly improves numerical stability and convergence.\n    In typical metals, n ranges from about 20-200: FCC metals n&asymp;20-100, BCC metals n&asymp;50-200, HCP metals n&asymp;20-80.\n    n=1 corresponds to linear viscosity (Newtonian fluid), and large n approximates ideal plasticity.\n  <\/p>\n\n  <h3 id=\"ch2-4\">\n    <span data-lang=\"zh\">2.4 \u786c\u5316\u5f8b \/ Hardening Laws<\/span>\n    <span data-lang=\"en\" style=\"display:none\">2.4 Hardening Laws<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u4e34\u754c\u5206\u5207\u5e94\u529bg<sup>(&alpha;)<\/sup>\u968f\u7d2f\u79ef\u5851\u6027\u526a\u5207\u5e94\u53d8&gamma; = &Sigma;<sub>&alpha;<\/sub> &int;|&gamma;&#775;<sup>(&alpha;)<\/sup>|dt\u800c\u6f14\u5316\uff0c\n    \u6f14\u5316\u65b9\u7a0b\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The critical resolved shear stress g<sup>(&alpha;)<\/sup> evolves with the accumulated plastic shear strain\n    &gamma; = &Sigma;<sub>&alpha;<\/sub> &int;|&gamma;&#775;<sup>(&alpha;)<\/sup>|dt, according to:\n  <\/p>\n  <div class=\"formula\">\n    <b>&#287;<sup>(&alpha;)<\/sup> = &Sigma;<sub>&beta;<\/sub> h<sub>&alpha;&beta;<\/sub> |&gamma;&#775;<sup>(&beta;)<\/sup>|<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5176\u4e2dh<sub>&alpha;&beta;<\/sub>\u4e3a\u786c\u5316\u77e9\u9635\uff0c\u5176\u5bf9\u89d2\u5143\u7d20h<sub>&alpha;&alpha;<\/sub>\u8868\u793a\u81ea\u786c\u5316\uff08self-hardening\uff09\uff0c\n    \u975e\u5bf9\u89d2\u5143\u7d20h<sub>&alpha;&beta;<\/sub>(&alpha;&ne;&beta;)\u8868\u793a\u6f5c\u786c\u5316\uff08latent hardening\uff09\uff0c\n    \u63cf\u8ff0\u4e0d\u540c\u6ed1\u79fb\u7cfb\u4e4b\u95f4\u7684\u4ea4\u4e92\u5f3a\u5316\u6548\u5e94\u3002\n    \u9ec4\u6c38\u521aUMAT\u63d0\u4f9b\u4e24\u79cd\u786c\u5316\u6a21\u578b\u9009\u9879\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    where h<sub>&alpha;&beta;<\/sub> is the hardening matrix. The diagonal elements h<sub>&alpha;&alpha;<\/sub> represent self-hardening,\n    and the off-diagonal elements h<sub>&alpha;&beta;<\/sub> (&alpha;&ne;&beta;) represent latent hardening,\n    describing the interaction strengthening effect between different slip systems.\n    Huang&#8217;s UMAT provides two hardening model options:\n  <\/p>\n\n  <h4>(a) \u53cc\u66f2\u6b63\u5272\u786c\u5316\u6a21\u578b \/ Hyperbolic Secant Hardening Model<\/h4>\n  <p data-lang=\"zh\">\n    \u81ea\u786c\u5316\u6a21\u91cf\u91c7\u7528\u53cc\u66f2\u6b63\u5272\uff08sech\u00b2\uff09\u5f62\u5f0f\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The self-hardening modulus adopts a hyperbolic secant (sech\u00b2) form:\n  <\/p>\n  <div class=\"formula\">\n    <b>h<sub>&alpha;&alpha;<\/sub> = h(&gamma;) = H<sub>0<\/sub> &middot; sech\u00b2( H<sub>0<\/sub> &gamma; \/ (&tau;<sub>s<\/sub> &#8211; &tau;<sub>0<\/sub>) )<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u7d2f\u79ef\u526a\u5207\u5e94\u53d8&gamma;\u4ece0\u589e\u52a0\u5230&infin;\u65f6\uff0c\u786c\u5316\u6a21\u91cfh\u4eceH<sub>0<\/sub>\u5355\u8c03\u9012\u51cf\u52300\uff0c\n    \u4e34\u754c\u5206\u5207\u5e94\u529b\u4ece\u521d\u59cb\u503c&tau;<sub>0<\/sub>\u6e10\u8fd1\u8d8b\u8fd1\u4e8e\u9971\u548c\u503c&tau;<sub>s<\/sub>\u3002\n    \u79ef\u5206\u540e\u53ef\u5f97g(&gamma;)\u7684\u663e\u5f0f\u8868\u8fbe\u5f0f\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    As the accumulated shear strain &gamma; increases from 0 to &infin;, the hardening modulus h decreases monotonically from H<sub>0<\/sub> to 0,\n    and the critical resolved shear stress asymptotically approaches the saturation value &tau;<sub>s<\/sub> from the initial value &tau;<sub>0<\/sub>.\n    Integration yields the explicit expression for g(&gamma;):\n  <\/p>\n  <div class=\"formula\">\n    <b>g(&gamma;) = &tau;<sub>s<\/sub> + (&tau;<sub>0<\/sub> &#8211; &tau;<sub>s<\/sub>) &middot; tanh( H<sub>0<\/sub> &gamma; \/ (&tau;<sub>s<\/sub> &#8211; &tau;<sub>0<\/sub>) )<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u53c2\u6570\u8bf4\u660e\uff1aH<sub>0<\/sub>\u4e3a\u521d\u59cb\u786c\u5316\u6a21\u91cf\uff08MPa\uff09\uff0c&tau;<sub>0<\/sub>\u4e3a\u521d\u59cbCRSS\uff08MPa\uff09\uff0c&tau;<sub>s<\/sub>\u4e3a\u9971\u548cCRSS\uff08MPa\uff09\u3002\n    \u8be5\u6a21\u578b\u7b80\u6d01\u4e14\u7269\u7406\u610f\u4e49\u660e\u786e\uff0c\u9002\u7528\u4e8e\u5927\u591a\u6570FCC\u91d1\u5c5e\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Parameter definitions: H<sub>0<\/sub> is the initial hardening modulus (MPa), &tau;<sub>0<\/sub> is the initial CRSS (MPa), and &tau;<sub>s<\/sub> is the saturation CRSS (MPa).\n    This model is concise with clear physical meaning and is suitable for most FCC metals.\n  <\/p>\n\n  <h4>(b) Bassani-Wu\u786c\u5316\u6a21\u578b \/ Bassani-Wu Hardening Model<\/h4>\n  <p data-lang=\"zh\">\n    Bassani\u548cWu\uff081991\uff09\u63d0\u51fa\u7684\u786c\u5316\u6a21\u578b\u8003\u8651\u4e86\u6ed1\u79fb\u7cfb\u95f4\u4ea4\u4e92\u4f5c\u7528\u7684\u5404\u5411\u5f02\u6027\u7279\u5f81\u3002\n    \u81ea\u786c\u5316\u6a21\u91cfF(&gamma;)\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The hardening model proposed by Bassani and Wu (1991) accounts for the anisotropic characteristics of interaction between slip systems.\n    The self-hardening modulus F(&gamma;) is:\n  <\/p>\n  <div class=\"formula\">\n    <b>F(&gamma;) = (H<sub>0<\/sub> &#8211; H<sub>s<\/sub>) &middot; sech\u00b2( (&tau;<sub>s<\/sub> &#8211; &tau;<sub>0<\/sub>) &gamma; \/ (&gamma;<sub>0<\/sub> + &gamma;) ) + H<sub>s<\/sub><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u6f5c\u786c\u5316\u7cfb\u6570q<sub>&alpha;&beta;<\/sub>\u53d6\u51b3\u4e8e\u6ed1\u79fb\u7cfb\u5bf9(&alpha;,&beta;)\u7684\u4ea4\u4e92\u7c7b\u578b\uff08\u5171\u9762\/\u5171\u5411\/\u6b63\u4ea4\/Lomer-Cottrell\u9501\u7b49\uff09\uff0c\n    \u786c\u5316\u77e9\u9635\u4e00\u822c\u5199\u4f5ch<sub>&alpha;&beta;<\/sub> = q<sub>&alpha;&beta;<\/sub> &middot; h<sub>&beta;&beta;<\/sub>\u3002\n    \u5b8c\u6574\u7684\u6f5c\u786c\u5316\u77e9\u9635\u662f12&times;12\u7684\u5bf9\u79f0\u77e9\u9635\uff0c\u5305\u542b7\u79cd\u4e0d\u540c\u7684\u4ea4\u4e92\u7c7b\u578b\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The latent hardening coefficient q<sub>&alpha;&beta;<\/sub> depends on the interaction type of the slip system pair (&alpha;,&beta;)\n    (coplanar\/codirectional\/orthogonal\/Lomer-Cottrell lock, etc.).\n    The hardening matrix is generally written as h<sub>&alpha;&beta;<\/sub> = q<sub>&alpha;&beta;<\/sub> &middot; h<sub>&beta;&beta;<\/sub>.\n    The full latent hardening matrix is a 12&times;12 symmetric matrix containing 7 different interaction types.\n  <\/p>\n  <p data-lang=\"zh\">\n    Bassani\u6a21\u578b\u7684\u4f18\u52bf\u5728\u4e8e\u80fd\u591f\u66f4\u51c6\u786e\u5730\u63cf\u8ff0FCC\u6676\u4f53\u4e2d\u4e0d\u540c\u6ed1\u79fb\u7cfb\u4e4b\u95f4\u7684\u4ea4\u4e92\u786c\u5316\u6548\u5e94\uff0c\n    \u7279\u522b\u662f\u975e\u5171\u9762\u6ed1\u79fb\u7cfb\u4e4b\u95f4\u7684\u5f3a\u76f8\u4e92\u4f5c\u7528\uff08\u5982Lomer-Cottrell\u9501\uff09\u3002\n    \u4f46\u8be5\u6a21\u578b\u9700\u8981\u989d\u5916\u7684\u6f5c\u786c\u5316\u53c2\u6570\uff0c\u4e14\u4ea4\u4e92\u7c7b\u578b\u5206\u7c7b\u8f83\u4e3a\u590d\u6742\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The advantage of the Bassani model is its ability to more accurately describe the interaction hardening effects between different slip systems in FCC crystals,\n    especially the strong interactions between non-coplanar slip systems (such as Lomer-Cottrell locks).\n    However, this model requires additional latent hardening parameters, and the interaction type classification is relatively complex.\n  <\/p>\n\n  <h3 id=\"ch2-5\">\n    <span data-lang=\"zh\">2.5 \u5e94\u529b\u66f4\u65b0 \/ Stress Update (Jaumann Rate)<\/span>\n    <span data-lang=\"en\" style=\"display:none\">2.5 Stress Update (Jaumann Rate)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u5728\u6709\u9650\u53d8\u5f62\u6846\u67b6\u4e0b\uff0c\u672c\u6784\u5173\u7cfb\u5fc5\u987b\u4ee5\u5ba2\u89c2\uff08objective\uff09\u7684\u5e94\u529b\u7387\u5f62\u5f0f\u8868\u8ff0\uff0c\u4ee5\u6d88\u9664\u521a\u6027\u8f6c\u52a8\u5e26\u6765\u7684\u865a\u5047\u5e94\u529b\u53d8\u5316\u3002\n    \u9ec4\u6c38\u521aUMAT\u91c7\u7528Jaumann\u5e94\u529b\u7387\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In the finite deformation framework, the constitutive relation must be expressed in an objective stress rate form\n    to eliminate spurious stress changes due to rigid body rotation.\n    Huang&#8217;s UMAT adopts the Jaumann stress rate:\n  <\/p>\n  <div class=\"formula\">\n    <b>&nabla;&sigma; = C<sup>ep<\/sup> : D<sup>e<\/sup> = C<sup>ep<\/sup> : (D &#8211; D<sup>p<\/sup>)<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5176\u4e2d&nabla;&sigma; = &sigma;&#775; &#8211; W&middot;&sigma; + &sigma;&middot;W\u4e3aJaumann\u5e94\u529b\u7387\uff0c\n    D\u4e3a\u53d8\u5f62\u7387\uff08symmetric part of velocity gradient L\uff09\uff0c\n    W\u4e3a\u65cb\u7387\uff08skew-symmetric part of L\uff09\uff0c\n    C<sup>ep<\/sup>\u4e3a\u5f39\u5851\u6027\u521a\u5ea6\u5f20\u91cf\uff0cD<sup>p<\/sup>\u4e3a\u5851\u6027\u53d8\u5f62\u7387\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    where &nabla;&sigma; = &sigma;&#775; &#8211; W&middot;&sigma; + &sigma;&middot;W is the Jaumann stress rate,\n    D is the rate of deformation (symmetric part of velocity gradient L),\n    W is the spin (skew-symmetric part of L),\n    C<sup>ep<\/sup> is the elastoplastic stiffness tensor, and D<sup>p<\/sup> is the plastic rate of deformation:\n  <\/p>\n  <div class=\"formula\">\n    <b>D<sup>p<\/sup> = &Sigma;<sub>&alpha;<\/sub> &gamma;&#775;<sup>(&alpha;)<\/sup> P<sup>(&alpha;)<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u9690\u5f0f\u65f6\u95f4\u79ef\u5206\u91c7\u7528Newton-Raphson\u8fed\u4ee3\u6c42\u89e3\u6bcf\u4e2a\u65f6\u95f4\u6b65\u672b\u7684\u5e94\u529b\u548c\u72b6\u6001\u53d8\u91cf\u3002\n    \u7ebf\u6027\u5316\u540e\u5f97\u5230\u4ee5&Delta;&gamma;<sup>(&alpha;)<\/sup>\u4e3a\u672a\u77e5\u91cf\u7684\u7ebf\u6027\u65b9\u7a0b\u7ec4\uff0c\u901a\u8fc7LU\u5206\u89e3\u6c42\u89e3\u3002\n    \u6536\u655b\u540e\u66f4\u65b0STATEV\u6570\u7ec4\u4e2d\u7684\u5e94\u529b\u3001\u5e94\u53d8\u548c\u5404\u6ed1\u79fb\u7cfb\u53d8\u91cf\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Implicit time integration employs Newton-Raphson iteration to solve for the stress and state variables at the end of each time increment.\n    Linearization yields a linear system of equations with &Delta;&gamma;<sup>(&alpha;)<\/sup> as unknowns, solved by LU decomposition.\n    After convergence, the STATEV array is updated with stress, strain, and slip system variables.\n  <\/p>\n<\/div>\n\n<!-- ==================== CHAPTER 3 ==================== -->\n<div class=\"chapter\" id=\"ch3\">\n  <h2>\n    <span class=\"chap-num\">03<\/span>\n    <span data-lang=\"zh\">Abaqus UMAT\u63a5\u53e3 \/ UMAT Interface<\/span>\n    <span data-lang=\"en\" style=\"display:none\">UMAT Interface<\/span>\n  <\/h2>\n\n  <h3 id=\"ch3-1\">\n    <span data-lang=\"zh\">3.1 \u5b8c\u6574Fortran\u7b7e\u540d \/ Full Fortran Signature<\/span>\n    <span data-lang=\"en\" style=\"display:none\">3.1 Full Fortran Signature<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    Abaqus\/Standard\u7684\u7528\u6237\u6750\u6599\u5b50\u7a0b\u5e8fUMAT\u5177\u6709\u56fa\u5b9a\u7684\u63a5\u53e3\u7b7e\u540d\u3002\u9ec4\u6c38\u521a\u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u5728\u8be5\u6846\u67b6\u4e0b\u5b9e\u73b0\uff0c\u5b8c\u6574\u58f0\u660e\u5982\u4e0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The Abaqus\/Standard User Material subroutine (UMAT) has a fixed interface signature. Huang&#8217;s single-crystal plasticity subroutine is implemented within this framework. The complete declaration is:\n  <\/p>\n  <pre><code>      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,\n     1 RPL,DDSDDT,DRPLDE,DRPLDT,\n     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,\n     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,\n     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)\nC-----------------------------------------------------------------------\nC     UMAT for single crystal plasticity\nC     Developed by Yonggang Huang\nC-----------------------------------------------------------------------\n      INCLUDE 'ABA_PARAM.INC'\n      CHARACTER*80 CMNAME\n      DIMENSION STRESS(NTENS),STATEV(NSTATV),\n     1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),\n     2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),\n     3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)\nC     ... user coding ...\n      RETURN\n      END<\/code><\/pre>\n\n  <h3 id=\"ch3-2\">\n    <span data-lang=\"zh\">3.2 \u53d8\u91cf\u8bf4\u660e\u8868 \/ Variable Specification Table<\/span>\n    <span data-lang=\"en\" style=\"display:none\">3.2 Variable Specification Table<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u4e0b\u8868\u5217\u51fa\u4e86UMAT\u63a5\u53e3\u4e2d\u4e0e\u5355\u6676\u5851\u6027\u5b9e\u73b0\u5bc6\u5207\u76f8\u5173\u7684\u8f93\u5165\/\u8f93\u51fa\u53d8\u91cf\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The following table lists the input\/output variables in the UMAT interface most relevant to the single-crystal plasticity implementation:\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>\n          <span data-lang=\"zh\">\u53d8\u91cf\u540d \/ Variable<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Variable<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u7c7b\u578b \/ Type<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Type<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr>\n        <td><code>STRESS(NTENS)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6570\u7ec4 \/ Array<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Array<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5e94\u529b\u5f20\u91cf\uff08\u8f93\u5165\u4e3a\u589e\u91cf\u5f00\u59cb\u65f6\u7684\u5e94\u529b\uff0c\u8f93\u51fa\u4e3a\u589e\u91cf\u7ed3\u675f\u65f6\u7684\u5e94\u529b\uff09\uff0c\u6309Abaqus\u7ea6\u5b9a\u5b58\u50a8\u4e3aVoigt\u5411\u91cf\uff1a&sigma;11, &sigma;22, &sigma;33, &sigma;12, &sigma;13, &sigma;23<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Stress tensor (input: stress at increment start; output: stress at increment end), stored as Voigt vector per Abaqus convention: &sigma;11, &sigma;22, &sigma;33, &sigma;12, &sigma;13, &sigma;23<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>STATEV(NSTATV)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6570\u7ec4 \/ Array<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Array<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u72b6\u6001\u53d8\u91cf\u6570\u7ec4\uff08solution-dependent state variables\uff09\uff0c\u7528\u4e8e\u5b58\u50a8\u5851\u6027\u5e94\u53d8\u3001\u7d2f\u79ef\u526a\u5207\u3001\u6ed1\u79fb\u7cfb\u5f3a\u5ea6\u7b49\u5386\u53f2\u4fe1\u606f\u3002\u5fc5\u987b\u5728\u672c\u5b50\u7a0b\u5e8f\u4e2d\u5b9a\u4e49\u66f4\u65b0\u89c4\u5219\u3002<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Solution-dependent state variable array for storing plastic strain, accumulated shear, slip system strength, and other history information. Update rules must be defined in this subroutine.<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DDSDDE(NTENS,NTENS)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u4e8c\u7ef4\u6570\u7ec4 \/ Matrix<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Matrix<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u6750\u6599Jacobian\u77e9\u9635\uff08\u4e00\u81f4\u5207\u7ebf\u6a21\u91cf\uff09\uff0c&part;&Delta;&sigma;\/&part;&Delta;&epsilon;\u3002\u8fd9\u662f\u9690\u5f0f\u6c42\u89e3\u5668\u4fdd\u8bc1\u4e8c\u9636\u6536\u655b\u7684\u5173\u952e\uff0c\u5fc5\u987b\u63d0\u4f9b\u51c6\u786e\u7684\u5f39\u5851\u6027\u5207\u7ebf\u521a\u5ea6\u3002<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Material Jacobian matrix (consistent tangent modulus), &part;&Delta;&sigma;\/&part;&Delta;&epsilon;. This is essential for the implicit solver to achieve quadratic convergence; accurate elastoplastic tangent stiffness must be provided.<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>STRAN(NTENS)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6570\u7ec4 \/ Array<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Array<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u589e\u91cf\u5f00\u59cb\u65f6\u7684\u603b\u5e94\u53d8\uff08engineering strain\uff0c\u526a\u5e94\u53d8\u4e0d\u9664\u4ee52\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Total strain at increment start (engineering strain; shear strains are not divided by 2)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DSTRAN(NTENS)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6570\u7ec4 \/ Array<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Array<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5e94\u53d8\u589e\u91cf\uff08engineering strain increment\uff09\uff0c\u7531Abaqus\u6839\u636e\u6574\u4f53\u5e73\u8861\u8fed\u4ee3\u63d0\u4f9b<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Strain increment (engineering strain increment), provided by Abaqus based on global equilibrium iteration<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DFGRD0(3,3)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u77e9\u9635 \/ Matrix<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Matrix<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u589e\u91cf\u5f00\u59cb\u65f6\u7684\u53d8\u5f62\u68af\u5ea6F\uff083&times;3\u77e9\u9635\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Deformation gradient F at increment start (3&times;3 matrix)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DFGRD1(3,3)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u77e9\u9635 \/ Matrix<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Matrix<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u589e\u91cf\u7ed3\u675f\u65f6\u7684\u53d8\u5f62\u68af\u5ea6F\uff08\u7531Abaqus\u57fa\u4e8e\u6574\u4f53\u4f4d\u79fb\u573a\u8ba1\u7b97\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Deformation gradient F at increment end (computed by Abaqus based on global displacement field)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>PROPS(NPROPS)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6570\u7ec4 \/ Array<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Array<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u6750\u6599\u53c2\u6570\u6570\u7ec4\uff0c\u5728Abaqus\u8f93\u5165\u6587\u4ef6\u4e2d\u901a\u8fc7*Material, *User Material\u5173\u952e\u5b57\u5b9a\u4e49<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Material property array, defined in the Abaqus input file via *Material, *User Material keywords<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>NPROPS<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6574\u6570 \/ Integer<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Integer<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">PROPS\u6570\u7ec4\u957f\u5ea6\uff0c\u9ec4\u6c38\u521a\u7a0b\u5e8f\u901a\u5e38\u8981\u6c42NPROPS &ge; 160<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Length of PROPS array; Huang&#8217;s program typically requires NPROPS &ge; 160<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>NSTATV<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6574\u6570 \/ Integer<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Integer<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">STATEV\u6570\u7ec4\u957f\u5ea6\uff0c\u53d6\u51b3\u4e8e\u6ed1\u79fb\u7cfb\u6570\u91cf\u548c\u5b58\u50a8\u9700\u6c42\uff0c\u901a\u5e38NSTATV = 6 + NSLIP + NSLIP + &#8230;<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Length of STATEV array, depending on number of slip systems and storage requirements; typically NSTATV = 6 + NSLIP + NSLIP + &#8230;<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DROT(3,3)<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u77e9\u9635 \/ Matrix<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Matrix<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u589e\u91cf\u65cb\u8f6c\u5f20\u91cf\uff08\u7531\u6781\u5206\u89e3\u5f97\u5230\uff09\uff0cAbaqus\u81ea\u52a8\u8ba1\u7b97\u5e76\u63d0\u4f9b<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Incremental rotation tensor (from polar decomposition), automatically computed and provided by Abaqus<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>DTIME<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u53cc\u7cbe\u5ea6 \/ Double<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Double<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u65f6\u95f4\u589e\u91cf\u5927\u5c0f\uff08\u5728\u7387\u76f8\u5173\u6750\u6599\u4e2d\u7528\u4e8e\u8ba1\u7b97&gamma;&#775;\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Time increment size (used to compute &gamma;&#775; in rate-dependent materials)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>PNEWDT<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u53cc\u7cbe\u5ea6 \/ Double<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Double<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5efa\u8bae\u7684\u4e0b\u4e00\u589e\u91cf\u65f6\u95f4\u6bd4\u4f8b\u56e0\u5b50\uff080 &lt; PNEWDT &lt; 1\uff09\u3002\u82e5\u5c40\u90e8\u8fed\u4ee3\u4e0d\u6536\u655b\uff0c\u53ef\u8bbe\u7f6e\u6b64\u503c\u4f7fAbaqus\u81ea\u52a8\u51cf\u5c0f\u65f6\u95f4\u6b65\u3002<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Suggested ratio for next increment time step (0 &lt; PNEWDT &lt; 1). If local iteration does not converge, set this to let Abaqus automatically reduce the time step.<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td><code>NOEL, NPT<\/code><\/td>\n        <td>\n          <span data-lang=\"zh\">\u6574\u6570 \/ Integer<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Integer<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5355\u5143\u7f16\u53f7\u548c\u79ef\u5206\u70b9\u7f16\u53f7\uff0c\u7528\u4e8e\u8c03\u8bd5\u548c\u5b9a\u4f4d<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Element number and integration point number, for debugging and localization<\/span>\n        <\/td>\n      <\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch3-3\">\n    <span data-lang=\"zh\">3.3 STATEV\u5e03\u5c40 \/ State Variable Layout Table<\/span>\n    <span data-lang=\"en\" style=\"display:none\">3.3 State Variable Layout Table<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    STATEV\u6570\u7ec4\u662fUMAT\u5b50\u7a0b\u5e8f\u7684\u6838\u5fc3\u6570\u636e\u7ed3\u6784\uff0c\u7528\u4e8e\u5728\u589e\u91cf\u6b65\u4e4b\u95f4\u4f20\u9012\u548c\u4fdd\u5b58\u6750\u6599\u7684\u5386\u53f2\u72b6\u6001\u3002\n    \u9ec4\u6c38\u521aUMAT\u7684STATEV\u5e03\u5c40\u5982\u4e0b\uff08\u4ee5FCC 12\u6ed1\u79fb\u7cfb\u4e3a\u4f8b\uff09\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The STATEV array is the core data structure of the UMAT subroutine, used to pass and preserve material history states between increments.\n    The STATEV layout in Huang&#8217;s UMAT is as follows (for FCC 12 slip systems as an example):\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>\n          <span data-lang=\"zh\">\u7d22\u5f15\u8303\u56f4 \/ Index Range<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Index Range<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u5185\u5bb9 \/ Content<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Content<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr>\n        <td>1 &#8211; 6<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5851\u6027\u5e94\u53d8 \/ Plastic strain<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Plastic strain<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u7d2f\u79ef\u5851\u6027\u5e94\u53d8\u5206\u91cf &epsilon;<sup>p<\/sup><sub>11<\/sub>, &epsilon;<sup>p<\/sup><sub>22<\/sub>, &epsilon;<sup>p<\/sup><sub>33<\/sub>, &gamma;<sup>p<\/sup><sub>12<\/sub>, &gamma;<sup>p<\/sup><sub>13<\/sub>, &gamma;<sup>p<\/sup><sub>23<\/sub>\uff08\u5de5\u7a0b\u526a\u5e94\u53d8\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Accumulated plastic strain components &epsilon;<sup>p<\/sup><sub>11<\/sub>, &epsilon;<sup>p<\/sup><sub>22<\/sub>, &epsilon;<sup>p<\/sup><sub>33<\/sub>, &gamma;<sup>p<\/sup><sub>12<\/sub>, &gamma;<sup>p<\/sup><sub>13<\/sub>, &gamma;<sup>p<\/sup><sub>23<\/sub> (engineering shear strains)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>7 &#8211; 18<\/td>\n        <td>\n          <span data-lang=\"zh\">\u7d2f\u79ef\u526a\u5207\u5e94\u53d8 \/ Accumulated shear<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Accumulated shear<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5404\u6ed1\u79fb\u7cfb\u7d2f\u79ef\u526a\u5207\u5e94\u53d8 &gamma;<sup>(&alpha;)<\/sup> = &int;|&gamma;&#775;<sup>(&alpha;)<\/sup>|dt\uff0c&alpha; = 1,&#8230;,12<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Accumulated shear strain on each slip system &gamma;<sup>(&alpha;)<\/sup> = &int;|&gamma;&#775;<sup>(&alpha;)<\/sup>|dt, &alpha; = 1,&#8230;,12<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>19 &#8211; 30<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5f53\u524d\u6ed1\u79fb\u7cfb\u5f3a\u5ea6 \/ Current slip strength<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Current slip strength<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5404\u6ed1\u79fb\u7cfb\u5f53\u524d\u4e34\u754c\u5206\u5207\u5e94\u529b g<sup>(&alpha;)<\/sup>\uff0c&alpha; = 1,&#8230;,12<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Current critical resolved shear stress g<sup>(&alpha;)<\/sup> for each slip system, &alpha; = 1,&#8230;,12<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>31 &#8211; 39<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5f39\u6027\u53d8\u5f62\u68af\u5ea6 \/ Elastic def. grad.<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Elastic def. grad.<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">F<sup>e<\/sup>\u76849\u4e2a\u5206\u91cf\uff083&times;3\u77e9\u9635\u6309\u5217\u4f18\u5148\u5b58\u50a8\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">9 components of F<sup>e<\/sup> (3&times;3 matrix stored column-major)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>40 &#8211; 48<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5851\u6027\u53d8\u5f62\u68af\u5ea6 \/ Plastic def. grad.<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Plastic def. grad.<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">F<sup>p<\/sup>\u76849\u4e2a\u5206\u91cf\uff083&times;3\u77e9\u9635\u6309\u5217\u4f18\u5148\u5b58\u50a8\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">9 components of F<sup>p<\/sup> (3&times;3 matrix stored column-major)<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>49 &#8211; 54<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5e94\u529b \/ Stress<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Stress<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u5b58\u50a8Cauchy\u5e94\u529b\u5206\u91cf\uff0c\u4fbf\u4e8e\u540e\u5904\u7406\u8f93\u51fa<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Cauchy stress components stored for post-processing output<\/span>\n        <\/td>\n      <\/tr>\n      <tr>\n        <td>55+<\/td>\n        <td>\n          <span data-lang=\"zh\">\u5176\u4ed6\u8f85\u52a9\u53d8\u91cf \/ Aux. variables<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Aux. variables<\/span>\n        <\/td>\n        <td>\n          <span data-lang=\"zh\">\u6676\u4f53\u53d6\u5411\u3001Euler\u89d2\u3001\u65cb\u8f6c\u77e9\u9635\u5206\u91cf\u7b49\uff08\u5177\u4f53\u53d6\u51b3\u4e8ePROPS\u914d\u7f6e\uff09<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Crystal orientation, Euler angles, rotation matrix components, etc. (depending on PROPS configuration)<\/span>\n        <\/td>\n      <\/tr>\n    <\/tbody>\n  <\/table>\n\n  <div class=\"note\">\n    <span data-lang=\"zh\"><strong>\u91cd\u8981\u63d0\u793a\uff1a<\/strong>STATEV\u6570\u7ec4\u5728Abaqus\u4e2d\u901a\u8fc7*Depvar\u5173\u952e\u5b57\u5b9a\u4e49\u5927\u5c0f\u3002\u5fc5\u987b\u786e\u4fddNSTATV\u8db3\u591f\u5bb9\u7eb3\u6240\u6709\u9700\u8981\u5b58\u50a8\u7684\u72b6\u6001\u53d8\u91cf\uff0c\u5426\u5219\u4f1a\u53d1\u751f\u6570\u7ec4\u8d8a\u754c\u9519\u8bef\u6216\u72b6\u6001\u4fe1\u606f\u4e22\u5931\u3002<\/span>\n    <span data-lang=\"en\" style=\"display:none\"><strong>Important:<\/strong> The STATEV array size is defined in Abaqus via the *Depvar keyword. Ensure NSTATV is large enough to hold all state variables; otherwise array overflow or state loss will occur.<\/span>\n  <\/div>\n<\/div>\n\n<!-- ==================== CHAPTER 4 ==================== -->\n<div class=\"chapter\" id=\"ch4\">\n  <h2>\n    <span class=\"chap-num\">04<\/span>\n    <span data-lang=\"zh\">\u6750\u6599\u53c2\u6570 \/ Material Parameters<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Material Parameters<\/span>\n  <\/h2>\n\n  <p data-lang=\"zh\">\n    \u9ec4\u6c38\u521aUMAT\u901a\u8fc7PROPS\u6570\u7ec4\u4f20\u5165\u5168\u90e8\u6750\u6599\u53c2\u6570\u3002PROPS\u7684\u7d22\u5f15\u5e03\u5c40\u7ecf\u8fc7\u7cbe\u5fc3\u8bbe\u8ba1\uff0c\n    \u5c06\u5f39\u6027\u5e38\u6570\u3001\u6ed1\u79fb\u7cfb\u5b9a\u4e49\u3001\u6676\u4f53\u53d6\u5411\u3001\u7387\u76f8\u5173\u53c2\u6570\u3001\u786c\u5316\u53c2\u6570\u548c\u79ef\u5206\u63a7\u5236\u53c2\u6570\u5206\u533a\u5b58\u50a8\u3002\n    \u4ee5\u4e0b\u5404\u8282\u8be6\u7ec6\u5217\u51fa\u6bcf\u4e2aPROPS\u7d22\u5f15\u5bf9\u5e94\u7684\u53c2\u6570\u542b\u4e49\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Huang&#8217;s UMAT receives all material parameters through the PROPS array. The PROPS index layout is carefully designed,\n    partitioning elastic constants, slip system definitions, crystal orientation, rate-dependent parameters, hardening parameters, and integration control parameters.\n    The following sections detail the meaning of each PROPS index.\n  <\/p>\n\n  <h3 id=\"ch4-1\">\n    <span data-lang=\"zh\">4.1 \u5f39\u6027\u5e38\u6570 PROPS(1-21) \/ Elastic Constants<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.1 Elastic Constants PROPS(1-21)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u5f39\u6027\u5e38\u6570\u533a\u5b9a\u4e49\u6676\u4f53\u6750\u6599\u7684\u5404\u5411\u5f02\u6027\u5f39\u6027\u884c\u4e3a\u3002\u5bf9\u4e8e\u7acb\u65b9\u6676\u7cfb\uff0c\u53ea\u9700\u89813\u4e2a\u72ec\u7acb\u5e38\u6570\uff08C11, C12, C44\uff09\uff0c\n    \u4f46\u7a0b\u5e8f\u63a5\u53e3\u9884\u7559\u4e86\u5b8c\u657421\u4e2a\u53c2\u6570\u4ee5\u652f\u6301\u4e00\u822c\u5404\u5411\u5f02\u6027\u6750\u6599\uff08\u4e09\u659c\u6676\u7cfb\uff09\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The elastic constant section defines the anisotropic elastic behavior of the crystal material.\n    For cubic crystals, only 3 independent constants (C11, C12, C44) are needed,\n    but the program interface reserves a full 21 parameters to support general anisotropic materials (triclinic).\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>1<\/td><td>C11<\/td><td>C11<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC11\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C11 (MPa)<\/span><\/td><\/tr>\n      <tr><td>2<\/td><td>C12<\/td><td>C12<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC12\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C12 (MPa)<\/span><\/td><\/tr>\n      <tr><td>3<\/td><td>C13<\/td><td>C13<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC13\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC13=C12<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C13 (MPa); C13=C12 for cubic<\/span><\/td><\/tr>\n      <tr><td>4<\/td><td>C14<\/td><td>C14<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC14\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC14=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C14 (MPa); C14=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>5<\/td><td>C15<\/td><td>C15<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC15\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC15=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C15 (MPa); C15=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>6<\/td><td>C16<\/td><td>C16<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC16\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC16=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C16 (MPa); C16=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>7<\/td><td>C22<\/td><td>C22<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC22\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC22=C11<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C22 (MPa); C22=C11 for cubic<\/span><\/td><\/tr>\n      <tr><td>8<\/td><td>C23<\/td><td>C23<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC23\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC23=C12<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C23 (MPa); C23=C12 for cubic<\/span><\/td><\/tr>\n      <tr><td>9<\/td><td>C24<\/td><td>C24<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC24\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC24=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C24 (MPa); C24=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>10<\/td><td>C25<\/td><td>C25<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC25\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC25=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C25 (MPa); C25=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>11<\/td><td>C26<\/td><td>C26<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC26\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC26=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C26 (MPa); C26=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>12<\/td><td>C33<\/td><td>C33<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC33\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC33=C11<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C33 (MPa); C33=C11 for cubic<\/span><\/td><\/tr>\n      <tr><td>13<\/td><td>C34<\/td><td>C34<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC34\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC34=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C34 (MPa); C34=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>14<\/td><td>C35<\/td><td>C35<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC35\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC35=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C35 (MPa); C35=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>15<\/td><td>C36<\/td><td>C36<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC36\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC36=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C36 (MPa); C36=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>16<\/td><td>C44<\/td><td>C44<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC44\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C44 (MPa)<\/span><\/td><\/tr>\n      <tr><td>17<\/td><td>C45<\/td><td>C45<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC45\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC45=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C45 (MPa); C45=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>18<\/td><td>C46<\/td><td>C46<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC46\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC46=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C46 (MPa); C46=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>19<\/td><td>C55<\/td><td>C55<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC55\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC55=C44<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C55 (MPa); C55=C44 for cubic<\/span><\/td><\/tr>\n      <tr><td>20<\/td><td>C56<\/td><td>C56<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC56\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC56=0<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C56 (MPa); C56=0 for cubic<\/span><\/td><\/tr>\n      <tr><td>21<\/td><td>C66<\/td><td>C66<\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u5206\u91cfC66\uff08MPa\uff09\uff0c\u7acb\u65b9\u6676\u7cfb\u4e0bC66=C44<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness matrix component C66 (MPa); C66=C44 for cubic<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch4-2\">\n    <span data-lang=\"zh\">4.2 \u6ed1\u79fb\u7cfb\u5b9a\u4e49 PROPS(25-56) \/ Slip System Definition<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.2 Slip System Definition PROPS(25-56)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u6ed1\u79fb\u7cfb\u5b9a\u4e49\u533a\u7528\u4e8e\u6307\u5b9a\u6bcf\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u65b9\u5411s\u548c\u6ed1\u79fb\u9762\u6cd5\u5411m\uff08\u5728\u6676\u4f53\u5750\u6807\u7cfb\u4e2d\uff09\u3002\n    \u6bcf\u7ec4\u6ed1\u79fb\u7cfb\u9700\u89813\u4e2a\u65b9\u5411\u5206\u91cf\u548c3\u4e2a\u6cd5\u5411\u5206\u91cf\u3002\u5bf9\u4e8eFCC\u7ed3\u6784\uff0c\u901a\u5e38\u5b9a\u4e4912\u4e2a\u6ed1\u79fb\u7cfb\uff08{111}&lt;110&gt;\uff09\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The slip system definition section specifies the slip direction s and slip plane normal m (in crystal coordinates) for each slip system.\n    Each slip system requires 3 direction components and 3 normal components. For FCC, 12 slip systems ({111}&lt;110&gt;) are typically defined.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>22<\/td><td>NSLIP<\/td><td>NSLIP<\/td><td><span data-lang=\"zh\">\u6ed1\u79fb\u7cfb\u603b\u6570\uff08\u6574\u6570\uff0c\u5982FCC=12\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Total number of slip systems (integer, e.g., FCC=12)<\/span><\/td><\/tr>\n      <tr><td>23<\/td><td>NSLIP1<\/td><td>NSLIP1<\/td><td><span data-lang=\"zh\">\u7b2c\u4e00\u65cf\u6ed1\u79fb\u7cfb\u6570\u91cf\uff08\u5982FCC\u4e2d\u4e3a12\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Number of first-family slip systems (e.g., 12 for FCC)<\/span><\/td><\/tr>\n      <tr><td>24<\/td><td>NSLIP2<\/td><td>NSLIP2<\/td><td><span data-lang=\"zh\">\u7b2c\u4e8c\u65cf\u6ed1\u79fb\u7cfb\u6570\u91cf\uff08\u5982\u6709\u9700\u8981\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Number of second-family slip systems (if needed)<\/span><\/td><\/tr>\n      <tr><td>25-27<\/td><td>s1(1-3)<\/td><td>Slip dir. 1<\/td><td><span data-lang=\"zh\">\u7b2c1\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u65b9\u5411\u5206\u91cf(sx,sy,sz)<\/span><span data-lang=\"en\" style=\"display:none\">Slip direction components (sx,sy,sz) for slip system 1<\/span><\/td><\/tr>\n      <tr><td>28-30<\/td><td>m1(1-3)<\/td><td>Slip normal 1<\/td><td><span data-lang=\"zh\">\u7b2c1\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u9762\u6cd5\u5411\u5206\u91cf(mx,my,mz)<\/span><span data-lang=\"en\" style=\"display:none\">Slip plane normal components (mx,my,mz) for slip system 1<\/span><\/td><\/tr>\n      <tr><td>31-33<\/td><td>s2(1-3)<\/td><td>Slip dir. 2<\/td><td><span data-lang=\"zh\">\u7b2c2\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u65b9\u5411\u5206\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Slip direction components for slip system 2<\/span><\/td><\/tr>\n      <tr><td>34-36<\/td><td>m2(1-3)<\/td><td>Slip normal 2<\/td><td><span data-lang=\"zh\">\u7b2c2\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u9762\u6cd5\u5411\u5206\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Slip plane normal components for slip system 2<\/span><\/td><\/tr>\n      <tr><td colspan=\"4\" style=\"text-align:center\">\n        <span data-lang=\"zh\">&#8230; \u4f9d\u6b64\u7c7b\u63a8\uff0c\u6bcf6\u4e2aPROPS\u5b9a\u4e49\u4e00\u4e2a\u6ed1\u79fb\u7cfb &#8230;<\/span>\n        <span data-lang=\"en\" style=\"display:none\">&#8230; and so on; each slip system occupies 6 PROPS entries &#8230;<\/span>\n      <\/td><\/tr>\n      <tr><td>55-57<\/td><td>s12(1-3)<\/td><td>Slip dir. 12<\/td><td><span data-lang=\"zh\">\u7b2c12\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u65b9\u5411\u5206\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Slip direction components for slip system 12<\/span><\/td><\/tr>\n      <tr><td>58-60<\/td><td>m12(1-3)<\/td><td>Slip normal 12<\/td><td><span data-lang=\"zh\">\u7b2c12\u4e2a\u6ed1\u79fb\u7cfb\u7684\u6ed1\u79fb\u9762\u6cd5\u5411\u5206\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Slip plane normal components for slip system 12<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n\n  <p data-lang=\"zh\">\n    FCC\u6676\u4f53\u5178\u578b\u768412\u4e2a\u6ed1\u79fb\u7cfb\u4e3a{111}&lt;110&gt;\u65cf\uff0c\u5206\u4e3a4\u4e2a\u6ed1\u79fb\u9762\uff0c\u6bcf\u4e2a\u9762\u67093\u4e2a\u6ed1\u79fb\u65b9\u5411\u3002\n    \u4ee5\u4e0b\u4e3a\u6676\u4f53\u5750\u6807\u7cfb\u4e2d\u7684\u6807\u51c6\u5b9a\u4e49\uff08\u65e0\u9700\u5f52\u4e00\u5316\uff0c\u7a0b\u5e8f\u5185\u90e8\u4f1a\u81ea\u52a8\u5904\u7406\uff09\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The typical 12 slip systems for FCC crystals belong to the {111}&lt;110&gt; family, divided into 4 slip planes with 3 slip directions each.\n    The following is the standard definition in crystal coordinates (normalization is handled internally by the program):\n  <\/p>\n\n  <pre><code>C     FCC {111}&lt;110&gt; slip systems (standard definition)\nC     Plane 1 (111): directions [1-10], [10-1], [01-1]\nC     Plane 2 (1-11): directions [110], [101], [0-11]\nC     Plane 3 (-111): directions [110], [10-1], [011]\nC     Plane 4 (11-1): directions [1-10], [101], [011]<\/code><\/pre>\n\n  <h3 id=\"ch4-3\">\n    <span data-lang=\"zh\">4.3 \u6676\u4f53\u53d6\u5411 PROPS(57-70) \/ Crystal Orientation<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.3 Crystal Orientation PROPS(57-70)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u6676\u4f53\u53d6\u5411\u533a\u5b9a\u4e49\u6676\u4f53\u5750\u6807\u7cfb\u4e0e\u5168\u5c40\uff08\u8bd5\u6837\uff09\u5750\u6807\u7cfb\u4e4b\u95f4\u7684\u65cb\u8f6c\u5173\u7cfb\u3002\n    \u901a\u5e38\u91c7\u7528Bunge Euler\u89d2(&phi;1, &Phi;, &phi;2)\u6216\u65b9\u5411\u4f59\u5f26\u77e9\u9635\u8868\u793a\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The crystal orientation section defines the rotation relationship between the crystal coordinate system and the global (sample) coordinate system.\n    Bunge Euler angles (&phi;1, &Phi;, &phi;2) or the direction cosine matrix are commonly used.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>57-59<\/td><td>Euler\u89d21-3<\/td><td>Euler angles<\/td><td><span data-lang=\"zh\">Bunge Euler\u89d2 &phi;1, &Phi;, &phi;2\uff08\u5355\u4f4d\uff1a\u5ea6\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Bunge Euler angles &phi;1, &Phi;, &phi;2 (in degrees)<\/span><\/td><\/tr>\n      <tr><td>60-62<\/td><td>\u6676\u4f53\u8f741<\/td><td>Crystal axis 1<\/td><td><span data-lang=\"zh\">\u5168\u5c40\u5750\u6807\u7cfb\u4e2d\u6676\u4f53[100]\u8f74\u7684\u65b9\u5411\u4f59\u5f26<\/span><span data-lang=\"en\" style=\"display:none\">Direction cosines of crystal [100] axis in global coordinates<\/span><\/td><\/tr>\n      <tr><td>63-65<\/td><td>\u6676\u4f53\u8f742<\/td><td>Crystal axis 2<\/td><td><span data-lang=\"zh\">\u5168\u5c40\u5750\u6807\u7cfb\u4e2d\u6676\u4f53[010]\u8f74\u7684\u65b9\u5411\u4f59\u5f26<\/span><span data-lang=\"en\" style=\"display:none\">Direction cosines of crystal [010] axis in global coordinates<\/span><\/td><\/tr>\n      <tr><td>66-68<\/td><td>\u6676\u4f53\u8f743<\/td><td>Crystal axis 3<\/td><td><span data-lang=\"zh\">\u5168\u5c40\u5750\u6807\u7cfb\u4e2d\u6676\u4f53[001]\u8f74\u7684\u65b9\u5411\u4f59\u5f26<\/span><span data-lang=\"en\" style=\"display:none\">Direction cosines of crystal [001] axis in global coordinates<\/span><\/td><\/tr>\n      <tr><td>69<\/td><td>\u53d6\u5411\u8f93\u5165\u6a21\u5f0f<\/td><td>Orientation mode<\/td><td><span data-lang=\"zh\">0=\u4f7f\u7528Euler\u89d2\uff0c1=\u4f7f\u7528\u65b9\u5411\u4f59\u5f26\u77e9\u9635<\/span><span data-lang=\"en\" style=\"display:none\">0=use Euler angles, 1=use direction cosine matrix<\/span><\/td><\/tr>\n      <tr><td>70<\/td><td>\u53d6\u5411\u66f4\u65b0\u6807\u5fd7<\/td><td>Orientation update<\/td><td><span data-lang=\"zh\">1=\u5728\u53d8\u5f62\u8fc7\u7a0b\u4e2d\u66f4\u65b0\u6676\u4f53\u53d6\u5411\uff08\u8003\u8651\u6676\u683c\u8f6c\u52a8\uff09\uff0c0=\u56fa\u5b9a\u53d6\u5411<\/span><span data-lang=\"en\" style=\"display:none\">1=update crystal orientation during deformation (lattice rotation), 0=fixed orientation<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch4-4\">\n    <span data-lang=\"zh\">4.4 \u7387\u76f8\u5173\u53c2\u6570 PROPS(73-96) \/ Rate-Dependent Parameters<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.4 Rate-Dependent Parameters PROPS(73-96)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u7387\u76f8\u5173\u53c2\u6570\u533a\u63a7\u5236\u5e42\u5f8b\u6d41\u52a8\u6cd5\u5219\u7684\u884c\u4e3a\u3002\u5bf9\u4e8e\u591a\u65cf\u6ed1\u79fb\u7cfb\uff0c\u6bcf\u65cf\u53ef\u4ee5\u72ec\u7acb\u6307\u5b9a\u7387\u76f8\u5173\u53c2\u6570\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The rate-dependent parameter section controls the behavior of the power-law flow rule. For multiple families of slip systems, each family can have independent rate-dependent parameters.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>71<\/td><td>GAMMA0_1<\/td><td>GAMMA0_1<\/td><td><span data-lang=\"zh\">\u7b2c1\u65cf\u6ed1\u79fb\u7cfb\u7684\u53c2\u8003\u526a\u5207\u7387 &gamma;&#775;<sub>0<\/sub>\uff08s<sup>-1<\/sup>\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Reference shear rate &gamma;&#775;<sub>0<\/sub> for 1st family (s<sup>-1<\/sup>)<\/span><\/td><\/tr>\n      <tr><td>72<\/td><td>AM_1<\/td><td>AM_1<\/td><td><span data-lang=\"zh\">\u7b2c1\u65cf\u6ed1\u79fb\u7cfb\u7684\u7387\u654f\u611f\u6307\u6570\u5012\u6570 1\/n<\/span><span data-lang=\"en\" style=\"display:none\">Inverse rate sensitivity exponent 1\/n for 1st family<\/span><\/td><\/tr>\n      <tr><td>73<\/td><td>N_1<\/td><td>N_1<\/td><td><span data-lang=\"zh\">\u7b2c1\u65cf\u6ed1\u79fb\u7cfb\u7684\u7387\u654f\u611f\u6307\u6570 n<\/span><span data-lang=\"en\" style=\"display:none\">Rate sensitivity exponent n for 1st family<\/span><\/td><\/tr>\n      <tr><td>74<\/td><td>GAMMA0_2<\/td><td>GAMMA0_2<\/td><td><span data-lang=\"zh\">\u7b2c2\u65cf\u6ed1\u79fb\u7cfb\u7684\u53c2\u8003\u526a\u5207\u7387\uff08\u5982\u6709\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Reference shear rate for 2nd family (if applicable)<\/span><\/td><\/tr>\n      <tr><td>75<\/td><td>AM_2<\/td><td>AM_2<\/td><td><span data-lang=\"zh\">\u7b2c2\u65cf\u6ed1\u79fb\u7cfb\u7684\u7387\u654f\u611f\u6307\u6570\u5012\u6570 1\/n<\/span><span data-lang=\"en\" style=\"display:none\">Inverse rate sensitivity exponent 1\/n for 2nd family<\/span><\/td><\/tr>\n      <tr><td>76<\/td><td>N_2<\/td><td>N_2<\/td><td><span data-lang=\"zh\">\u7b2c2\u65cf\u6ed1\u79fb\u7cfb\u7684\u7387\u654f\u611f\u6307\u6570 n<\/span><span data-lang=\"en\" style=\"display:none\">Rate sensitivity exponent n for 2nd family<\/span><\/td><\/tr>\n      <tr><td>77-84<\/td><td>\u6d41\u52a8\u5f8b\u53c2\u6570<\/td><td>Flow params<\/td><td><span data-lang=\"zh\">\u7b2c3-4\u65cf\u6ed1\u79fb\u7cfb\u7684\u7387\u76f8\u5173\u53c2\u6570\uff08\u9884\u7559\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Rate-dependent parameters for 3rd-4th families (reserved)<\/span><\/td><\/tr>\n      <tr><td>85<\/td><td>TEMP_REF<\/td><td>TEMP_REF<\/td><td><span data-lang=\"zh\">\u53c2\u8003\u6e29\u5ea6\uff08K\uff09\uff0c\u7528\u4e8e\u70ed\u6fc0\u6d3b\u6a21\u578b<\/span><span data-lang=\"en\" style=\"display:none\">Reference temperature (K) for thermal activation model<\/span><\/td><\/tr>\n      <tr><td>86<\/td><td>TEMP_CUR<\/td><td>TEMP_CUR<\/td><td><span data-lang=\"zh\">\u5f53\u524d\u6e29\u5ea6\uff08K\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Current temperature (K)<\/span><\/td><\/tr>\n      <tr><td>87-96<\/td><td>\u9884\u7559<\/td><td>Reserved<\/td><td><span data-lang=\"zh\">\u4e3a\u540e\u7eed\u6269\u5c55\u9884\u7559<\/span><span data-lang=\"en\" style=\"display:none\">Reserved for future extensions<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch4-5\">\n    <span data-lang=\"zh\">4.5 \u786c\u5316\u53c2\u6570 PROPS(97-144) \/ Hardening Parameters<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.5 Hardening Parameters PROPS(97-144)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u786c\u5316\u53c2\u6570\u533a\u5305\u542b\u81ea\u786c\u5316\u548c\u6f5c\u786c\u5316\u7684\u5168\u90e8\u53c2\u6570\u3002\u6839\u636e\u9009\u62e9\u7684\u786c\u5316\u6a21\u578b\u7c7b\u578b\uff08PROPS\u4e2d\u7684\u6807\u5fd7\uff09\uff0c\u8fd9\u4e9b\u53c2\u6570\u7684\u89e3\u91ca\u6709\u6240\u4e0d\u540c\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The hardening parameter section contains all parameters for self-hardening and latent hardening. Their interpretation depends on the selected hardening model type (flag in PROPS).\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>97<\/td><td>\u786c\u5316\u6a21\u578b\u7c7b\u578b<\/td><td>Hardening model<\/td><td><span data-lang=\"zh\">0=\u53cc\u66f2\u6b63\u5272\u786c\u5316\uff0c1=Bassani\u786c\u5316\uff0c2=\u5176\u4ed6<\/span><span data-lang=\"en\" style=\"display:none\">0=Hyperbolic secant, 1=Bassani, 2=Other<\/span><\/td><\/tr>\n      <tr><td>98<\/td><td>TAU0<\/td><td>TAU0<\/td><td><span data-lang=\"zh\">\u521d\u59cb\u4e34\u754c\u5206\u5207\u5e94\u529b &tau;<sub>0<\/sub>\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Initial critical resolved shear stress &tau;<sub>0<\/sub> (MPa)<\/span><\/td><\/tr>\n      <tr><td>99<\/td><td>TAUS<\/td><td>TAUS<\/td><td><span data-lang=\"zh\">\u9971\u548c\u4e34\u754c\u5206\u5207\u5e94\u529b &tau;<sub>s<\/sub>\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Saturation critical resolved shear stress &tau;<sub>s<\/sub> (MPa)<\/span><\/td><\/tr>\n      <tr><td>100<\/td><td>H0<\/td><td>H0<\/td><td><span data-lang=\"zh\">\u521d\u59cb\u786c\u5316\u6a21\u91cf H<sub>0<\/sub>\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Initial hardening modulus H<sub>0<\/sub> (MPa)<\/span><\/td><\/tr>\n      <tr><td>101<\/td><td>HS<\/td><td>HS<\/td><td><span data-lang=\"zh\">Bassani\u6a21\u578b\u4e2d\u7684\u9971\u548c\u786c\u5316\u6a21\u91cf H<sub>s<\/sub>\uff08MPa\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Saturation hardening modulus H<sub>s<\/sub> in Bassani model (MPa)<\/span><\/td><\/tr>\n      <tr><td>102<\/td><td>GAMMA0_HARD<\/td><td>GAMMA0_HARD<\/td><td><span data-lang=\"zh\">Bassani\u6a21\u578b\u4e2d\u7684\u7279\u5f81\u526a\u5207\u5e94\u53d8 &gamma;<sub>0<\/sub><\/span><span data-lang=\"en\" style=\"display:none\">Characteristic shear strain &gamma;<sub>0<\/sub> in Bassani model<\/span><\/td><\/tr>\n      <tr><td>103<\/td><td>Q1<\/td><td>Q1<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q1\uff08\u5171\u9762\u6ed1\u79fb\u7cfb\u95f4\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q1 (coplanar)<\/span><\/td><\/tr>\n      <tr><td>104<\/td><td>Q2<\/td><td>Q2<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q2\uff08\u5171\u5411\u6ed1\u79fb\u7cfb\u95f4\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q2 (codirectional)<\/span><\/td><\/tr>\n      <tr><td>105<\/td><td>Q3<\/td><td>Q3<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q3\uff08\u6b63\u4ea4\u6ed1\u79fb\u7cfb\u95f4\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q3 (orthogonal)<\/span><\/td><\/tr>\n      <tr><td>106<\/td><td>Q4<\/td><td>Q4<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q4\uff08Hirth\u9501\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q4 (Hirth lock)<\/span><\/td><\/tr>\n      <tr><td>107<\/td><td>Q5<\/td><td>Q5<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q5\uff08Lomer-Cottrell\u9501\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q5 (Lomer-Cottrell lock)<\/span><\/td><\/tr>\n      <tr><td>108<\/td><td>Q6<\/td><td>Q6<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q6\uff08\u4ea4\u53c9\u6ed1\u79fb\u76f8\u5173\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q6 (cross-slip related)<\/span><\/td><\/tr>\n      <tr><td>109<\/td><td>Q7<\/td><td>Q7<\/td><td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570q7\uff08glissile\u9501\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient q7 (glissile lock)<\/span><\/td><\/tr>\n      <tr><td>110-144<\/td><td>\u9884\u7559\/\u6269\u5c55<\/td><td>Reserved<\/td><td><span data-lang=\"zh\">\u4e3a\u591a\u79cd\u6ed1\u79fb\u7cfb\u3001\u591a\u9636\u6bb5\u786c\u5316\u3001\u6e29\u5ea6\u6548\u5e94\u7b49\u9884\u7559<\/span><span data-lang=\"en\" style=\"display:none\">Reserved for multiple slip systems, multi-stage hardening, thermal effects, etc.<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch4-6\">\n    <span data-lang=\"zh\">4.6 \u79ef\u5206\u63a7\u5236 PROPS(145-160) \/ Integration Control<\/span>\n    <span data-lang=\"en\" style=\"display:none\">4.6 Integration Control PROPS(145-160)<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u79ef\u5206\u63a7\u5236\u53c2\u6570\u5f71\u54cd\u9690\u5f0fNewton-Raphson\u8fed\u4ee3\u7684\u6570\u503c\u884c\u4e3a\uff0c\u5305\u62ec\u6536\u655b\u5bb9\u5dee\u3001\u6700\u5927\u8fed\u4ee3\u6b21\u6570\u3001\n    \u81ea\u52a8\u65f6\u95f4\u6b65\u957f\u63a7\u5236\u7b56\u7565\u7b49\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Integration control parameters affect the numerical behavior of the implicit Newton-Raphson iteration,\n    including convergence tolerance, maximum iteration count, and automatic time step control strategy.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>PROPS Index<\/th>\n        <th>\n          <span data-lang=\"zh\">\u4e2d\u6587\u540d\u79f0<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Name (Chinese)<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u82f1\u6587\u540d\u79f0 \/ English Name<\/span>\n          <span data-lang=\"en\" style=\"display:none\">English Name<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td>145<\/td><td>ITMAX<\/td><td>ITMAX<\/td><td><span data-lang=\"zh\">\u6700\u5927Newton-Raphson\u8fed\u4ee3\u6b21\u6570\uff08\u9ed8\u8ba450\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Maximum Newton-Raphson iterations (default 50)<\/span><\/td><\/tr>\n      <tr><td>146<\/td><td>TOL<\/td><td>TOL<\/td><td><span data-lang=\"zh\">\u76f8\u5bf9\u6536\u655b\u5bb9\u5dee\uff08\u9ed8\u8ba41e-6\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Relative convergence tolerance (default 1e-6)<\/span><\/td><\/tr>\n      <tr><td>147<\/td><td>TOL2<\/td><td>TOL2<\/td><td><span data-lang=\"zh\">\u7edd\u5bf9\u6536\u655b\u5bb9\u5dee\uff08\u9ed8\u8ba41e-10\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Absolute convergence tolerance (default 1e-10)<\/span><\/td><\/tr>\n      <tr><td>148<\/td><td>DGMAX<\/td><td>DGMAX<\/td><td><span data-lang=\"zh\">\u5355\u6b65\u6700\u5927\u5141\u8bb8\u526a\u5207\u5e94\u53d8\u589e\u91cf\u9650\u5236<\/span><span data-lang=\"en\" style=\"display:none\">Maximum allowable shear strain increment per step<\/span><\/td><\/tr>\n      <tr><td>149<\/td><td>ITERMOD<\/td><td>ITERMOD<\/td><td><span data-lang=\"zh\">\u8fed\u4ee3\u6a21\u5f0f\uff1a0=\u6807\u51c6NR\uff0c1=\u7ebf\u641c\u7d22\u589e\u5f3a<\/span><span data-lang=\"en\" style=\"display:none\">Iteration mode: 0=standard NR, 1=line-search enhanced<\/span><\/td><\/tr>\n      <tr><td>150<\/td><td>OUTPUT<\/td><td>OUTPUT<\/td><td><span data-lang=\"zh\">\u8f93\u51fa\u63a7\u5236\u6807\u5fd7\uff1a0=\u6700\u5c0f\u8f93\u51fa\uff0c1=\u6807\u51c6\uff0c2=\u8be6\u7ec6\u8c03\u8bd5<\/span><span data-lang=\"en\" style=\"display:none\">Output control flag: 0=minimal, 1=standard, 2=verbose debug<\/span><\/td><\/tr>\n      <tr><td>151-160<\/td><td>\u9884\u7559<\/td><td>Reserved<\/td><td><span data-lang=\"zh\">\u4e3a\u81ea\u9002\u5e94\u6b65\u957f\u63a7\u5236\u3001\u5b50\u6b65\u79ef\u5206\u7b49\u9884\u7559<\/span><span data-lang=\"en\" style=\"display:none\">Reserved for adaptive step control, sub-stepping, etc.<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n<\/div>\n\n<!-- ==================== CHAPTER 5 ==================== -->\n<div class=\"chapter\" id=\"ch5\">\n  <h2>\n    <span class=\"chap-num\">05<\/span>\n    <span data-lang=\"zh\">\u4ee3\u7801\u7ed3\u6784 \/ Code Structure<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Code Structure<\/span>\n  <\/h2>\n\n  <h3 id=\"ch5-1\">\n    <span data-lang=\"zh\">5.1 \u5b50\u7a0b\u5e8f\u8c03\u7528\u5173\u7cfb\u56fe \/ Subroutine Call Graph<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.1 Subroutine Call Graph<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u9ec4\u6c38\u521aUMAT\u7531\u4e3b\u7a0b\u5e8fUMAT\u548c\u591a\u4e2a\u8f85\u52a9\u5b50\u7a0b\u5e8f\u7ec4\u6210\uff0c\u5404\u5b50\u7a0b\u5e8f\u627f\u62c5\u660e\u786e\u7684\u804c\u8d23\u3002\n    \u4ee5\u4e0b\u4e3a\u8c03\u7528\u5173\u7cfb\u56fe\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Huang&#8217;s UMAT consists of the main UMAT subroutine and multiple auxiliary subroutines, each with a well-defined responsibility.\n    The call graph is as follows:\n  <\/p>\n\n  <div class=\"mermaid\">\ngraph TD\n    A[UMAT<br\/>\u4e3b\u7a0b\u5e8f \/ Main Routine] &#8211;> B[ROTATION<br\/>\u8ba1\u7b97\u65cb\u8f6c\u77e9\u9635 \/ Compute Rotation]\n    A &#8211;> C[SLIPSYS<br\/>\u6ed1\u79fb\u7cfb\u521d\u59cb\u5316 \/ Slip System Init]\n    C &#8211;> D[GSLPINIT<br\/>\u521d\u59cb\u5316\u6ed1\u79fb\u5f3a\u5ea6 \/ Init Slip Strength]\n    D &#8211;> E[GSLP0<br\/>\u521d\u59cbCRSS \/ Initial CRSS]\n    A &#8211;> F[STRAINRATE<br\/>\u8ba1\u7b97\u526a\u5207\u7387 \/ Compute Shear Rate]\n    F &#8211;> G[F<br\/>\u6d41\u52a8\u51fd\u6570 \/ Flow Function]\n    F &#8211;> H[DFDX<br\/>\u6d41\u52a8\u51fd\u6570\u5bfc\u6570 \/ Flow Derivative]\n    A &#8211;> I[LATENTHARDEN<br\/>\u6f5c\u786c\u5316\u8ba1\u7b97 \/ Latent Hardening]\n    I &#8211;> J[HSELF<br\/>\u81ea\u786c\u5316 \/ Self-Hardening]\n    I &#8211;> K[HLATNT<br\/>\u6f5c\u786c\u5316\u7cfb\u6570 \/ Latent Coeffs]\n    A &#8211;> L[ITERATION<br\/>Newton-Raphson\u8fed\u4ee3 \/ NR Iteration]\n    L &#8211;> M[DHSELF<br\/>\u81ea\u786c\u5316\u5bfc\u6570 \/ Self-Hardening Deriv]\n    L &#8211;> N[DHLATN<br\/>\u6f5c\u786c\u5316\u5bfc\u6570 \/ Latent Hardening Deriv]\n    L &#8211;> O[LUDCMP<br\/>LU\u5206\u89e3 \/ LU Decomposition]\n    L &#8211;> P[LUBKSB<br\/>LU\u56de\u4ee3 \/ LU Back-Subst]\n    A &#8211;> Q[CROSS<br\/>\u5411\u91cf\u53c9\u4e58 \/ Vector Cross Product]\n    style A fill:#2563eb,stroke:#1a2332,color:#fff\n    style L fill:#0ea5e9,stroke:#1a2332,color:#fff\n  <\/div>\n\n  <h3 id=\"ch5-2\">\n    <span data-lang=\"zh\">5.2 \u4e3b\u7b97\u6cd5\u6d41\u7a0b\u56fe \/ Main Algorithm Flowchart<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.2 Main Algorithm Flowchart<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u6bcf\u4e2a\u589e\u91cf\u6b65\u4e2d\uff0cUMAT\u4e3b\u7a0b\u5e8f\u6267\u884c\u4ee5\u4e0b\u7b97\u6cd5\u6d41\u7a0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In each increment, the main UMAT routine executes the following algorithm:\n  <\/p>\n\n  <div class=\"mermaid\">\nflowchart TD\n    START([\u5f00\u59cb \/ Start]) &#8211;> A[\u8bfb\u53d6\u6750\u6599\u53c2\u6570<br\/>Read PROPS]\n    A &#8211;> B[\u8bfb\u53d6\u72b6\u6001\u53d8\u91cf<br\/>Read STATEV]\n    B &#8211;> C[\u8ba1\u7b97\u5f39\u6027\u8bd5\u5e94\u529b<br\/>Elastic Trial Stress]\n    C &#8211;> D{Newton-Raphson<br\/>\u8fed\u4ee3\u6536\u655b?}\n    D &#8212; \u5426 \/ No &#8211;> E[\u8ba1\u7b97\u6d41\u52a8\u5f8b\u548c\u786c\u5316<br\/>Flow &#038; Hardening]\n    E &#8211;> F[\u7ec4\u88c5Jacobian<br\/>Assemble Jacobian]\n    F &#8211;> G[LU\u5206\u89e3\u6c42\u89e3<br\/>LU Decomposition]\n    G &#8211;> H[\u66f4\u65b0\u526a\u5207\u589e\u91cf<br\/>Update &Delta;&gamma;]\n    H &#8211;> D\n    D &#8212; \u662f \/ Yes &#8211;> I[\u66f4\u65b0\u5e94\u529b<br\/>Update Stress]\n    I &#8211;> J[\u66f4\u65b0STATEV<br\/>Update STATEV]\n    J &#8211;> K[\u8ba1\u7b97DDSDDE<br\/>Compute Consistent Tangent]\n    K &#8211;> END([\u7ed3\u675f \/ End])\n    style D fill:#f59e0b,stroke:#1a2332,color:#fff\n    style START fill:#2563eb,stroke:#1a2332,color:#fff\n    style END fill:#ef4444,stroke:#1a2332,color:#fff\n  <\/div>\n\n  <h3 id=\"ch5-3\">\n    <span data-lang=\"zh\">5.3 \u5173\u952e\u4ee3\u7801\u6bb5 \/ Key Code Segments<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3 Key Code Segments<\/span>\n  <\/h3>\n\n  <h4>\n    <span data-lang=\"zh\">5.3.1 \u5f39\u6027\u8bd5\u5e94\u529b\u8ba1\u7b97 \/ Elastic Trial Stress<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3.1 Elastic Trial Stress<\/span>\n  <\/h4>\n  <pre><code>C-----------------------------------------------------------------------\nC     \u5f39\u6027\u8bd5\u5e94\u529b\u8ba1\u7b97 \/ Elastic trial stress computation\nC     \u57fa\u4e8e\u589e\u91cf\u5f00\u59cb\u65f6\u7684Fe\u548c\u6574\u4e2a\u5e94\u53d8\u589e\u91cf\u8ba1\u7b97\u8bd5\u5e94\u529b\nC-----------------------------------------------------------------------\n      DO I=1,3\n         DO J=1,3\n            FE_INV(I,J)=DFGRD1(I,J)   ! \u521d\u59cb\u4f30\u8ba1: F^e = F (\u5ffd\u7565\u5851\u6027)\n         ENDDO\n      ENDDO\nC\nC     \u4eceSTATEV\u6062\u590d\u4e0a\u4e00\u589e\u91cf\u6b65\u7684F^p\u548cF^e\n      CALL KINTOTAL(STATEV,NSTATV,FE_FP,NSLIP)\nC\nC     \u8ba1\u7b97\u8bd5\u5f39\u6027\u53d8\u5f62\u68af\u5ea6: F^{e,trial} = F_{n+1} * (F^p_n)^{-1}\n      CALL LUDCMP(FP_TRY,3,3,INDX,D)\n      CALL LUBKSB(FP_TRY,3,3,INDX,FE_TRY_COL1)\n      CALL LUBKSB(FP_TRY,3,3,INDX,FE_TRY_COL2)\n      CALL LUBKSB(FP_TRY,3,3,INDX,FE_TRY_COL3)\nC\nC     \u8ba1\u7b97Green\u5f39\u6027\u5e94\u53d8 E^e = 0.5*(F^{eT}F^e - I)\n      DO I=1,3\n         DO J=1,3\n            FE_TF_FE(I,J)=0.0\n            DO K=1,3\n               FE_TF_FE(I,J)=FE_TF_FE(I,J)+FE_TRY(K,I)*FE_TRY(K,J)\n            ENDDO\n         ENDDO\n      ENDDO\n      DO I=1,3\n         DO J=1,3\n            EE(I,J)=0.5*(FE_TF_FE(I,J)-DELTA(I,J))\n         ENDDO\n      ENDDO<\/code><\/pre>\n\n  <h4>\n    <span data-lang=\"zh\">5.3.2 \u6d41\u52a8\u5f8b\u5b9e\u73b0 \/ Flow Law Implementation<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3.2 Flow Law Implementation<\/span>\n  <\/h4>\n  <pre><code>C-----------------------------------------------------------------------\nC     \u7387\u76f8\u5173\u6d41\u52a8\u5f8b\u8ba1\u7b97 \/ Rate-dependent flow rule computation\nC     gamma_dot = gamma0 * |tau\/g|^n * sign(tau\/g)\nC-----------------------------------------------------------------------\n      SUBROUTINE STRAINRATE(TAU, G, GAMMA0, AM, NSLIP, DGAMMA)\n      IMPLICIT NONE\n      INTEGER NSLIP, I\n      REAL*8 TAU(NSLIP), G(NSLIP), GAMMA0, AM\n      REAL*8 DGAMMA(NSLIP), RATIO, ABSRAT, POWER\nC\n      DO I=1,NSLIP\n         IF (G(I).LE.0.0) THEN\n            WRITE(*,*) 'ERROR: g(',I,') =',G(I)\n            STOP\n         ENDIF\n         RATIO = TAU(I)\/G(I)\n         ABSRAT = DABS(RATIO)\nC        \u907f\u514d\u6570\u503c\u6ea2\u51fa \/ Avoid numerical overflow\n         IF (ABSRAT.GT.1.0D+10) THEN\n            POWER = 1.0D+10**AM\n         ELSE\n            POWER = ABSRAT**AM\n         ENDIF\nC        \u526a\u5207\u5e94\u53d8\u7387 \/ Shear strain rate\n         DGAMMA(I) = GAMMA0 * POWER\n         IF (RATIO.LT.0.0) DGAMMA(I) = -DGAMMA(I)\n      ENDDO\n      RETURN\n      END<\/code><\/pre>\n\n  <h4>\n    <span data-lang=\"zh\">5.3.3 \u81ea\u786c\u5316\u5b9e\u73b0 \/ Self-Hardening Implementation<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3.3 Self-Hardening Implementation<\/span>\n  <\/h4>\n  <pre><code>C-----------------------------------------------------------------------\nC     \u53cc\u66f2\u6b63\u5272\u81ea\u786c\u5316\u6a21\u578b \/ Hyperbolic secant self-hardening model\nC     h(gamma) = H0 * sech^2( H0*gamma \/ (tau_s - tau_0) )\nC-----------------------------------------------------------------------\n      FUNCTION HSELF(GAMMA, H0, TAUS, TAU0)\n      IMPLICIT NONE\n      REAL*8 HSELF, GAMMA, H0, TAUS, TAU0\n      REAL*8 ARG, SECH2\nC\n      IF (TAUS.LE.TAU0) THEN\n         HSELF = 0.0\n         RETURN\n      ENDIF\nC\n      ARG = H0 * GAMMA \/ (TAUS - TAU0)\nC     sech^2(x) = 1\/cosh^2(x)\n      SECH2 = 1.0 \/ (DCOSH(ARG)**2)\n      HSELF = H0 * SECH2\n      RETURN\n      END\nC\nC-----------------------------------------------------------------------\nC     \u5f53\u524d\u6ed1\u79fb\u7cfb\u5f3a\u5ea6 \/ Current slip system strength\nC     g(gamma) = tau_s + (tau_0 - tau_s) * tanh( H0*gamma\/(tau_s-tau_0) )\nC-----------------------------------------------------------------------\n      FUNCTION GSLP0(GAMMA, H0, TAUS, TAU0)\n      IMPLICIT NONE\n      REAL*8 GSLP0, GAMMA, H0, TAUS, TAU0\n      REAL*8 ARG\nC\n      IF (TAUS.LE.TAU0) THEN\n         GSLP0 = TAU0\n         RETURN\n      ENDIF\nC\n      ARG = H0 * GAMMA \/ (TAUS - TAU0)\n      GSLP0 = TAUS + (TAU0 - TAUS) * DTANH(ARG)\n      RETURN\n      END<\/code><\/pre>\n\n  <h4>\n    <span data-lang=\"zh\">5.3.4 Newton-Raphson\u8fed\u4ee3\u6838\u5fc3 \/ NR Iteration Core<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3.4 NR Iteration Core<\/span>\n  <\/h4>\n  <pre><code>C-----------------------------------------------------------------------\nC     Newton-Raphson\u8fed\u4ee3\u6c42\u89e3\u526a\u5207\u5e94\u53d8\u589e\u91cf \/ NR iteration for shear increment\nC     \u672a\u77e5\u91cf: DGA(NSLIP) = &Delta;&gamma;^(&alpha;)\nC     \u6b8b\u5dee:   R(&alpha;) = &Delta;&gamma;^(&alpha;) - &Delta;t * gamma_dot^{n+1}_(&alpha;) = 0\nC-----------------------------------------------------------------------\n      SUBROUTINE ITERATION(...)\n      IMPLICIT NONE\n      INTEGER NSLIP, ITMAX, ITER, I, J, INFO\n      REAL*8 TOL, DTIME, DGA(NSLIP), GA(NSLIP), TAU(NSLIP)\n      REAL*8 G(NSLIP), DG(NSLIP), GAMMA0, AM\n      REAL*8 JAC(NSLIP,NSLIP), RES(NSLIP), DGA_INC(NSLIP)\n      REAL*8 DFDX(NSLIP), DHDG(NSLIP,NSLIP)\n      REAL*8 DET, RESNORM, RESNORM0\n      INTEGER INDX(NSLIP)\nC\n      RESNORM0 = 1.0\n      DO ITER=1,ITMAX\nC        \u8ba1\u7b97\u5f53\u524d\u5e94\u529b\u72b6\u6001 \/ Compute current stress state\n         CALL STRESSUPDATE(DGA, TAU, SIGMA, ...)\nC\nC        \u8ba1\u7b97\u6d41\u52a8\u5f8b \/ Compute flow rule\n         CALL STRAINRATE(TAU, G, GAMMA0, AM, NSLIP, DG)\nC\nC        \u8ba1\u7b97\u6b8b\u5dee \/ Compute residual\n         DO I=1,NSLIP\n            RES(I) = DGA(I) - DTIME * DG(I)\n         ENDDO\nC\nC        \u8ba1\u7b97\u6b8b\u5dee\u8303\u6570 \/ Compute residual norm\n         RESNORM = 0.0\n         DO I=1,NSLIP\n            RESNORM = RESNORM + RES(I)*RES(I)\n         ENDDO\n         RESNORM = DSQRT(RESNORM)\n         IF (ITER.EQ.1) RESNORM0 = RESNORM\nC\nC        \u68c0\u67e5\u6536\u655b \/ Check convergence\n         IF (RESNORM.LT.TOL*RESNORM0 .OR. RESNORM.LT.1.0D-12) THEN\n            RETURN   ! \u6536\u655b \/ Converged\n         ENDIF\nC\nC        \u8ba1\u7b97Jacobian\u77e9\u9635 \/ Compute Jacobian matrix\nC        J_{&alpha;&beta;} = &delta;_{&alpha;&beta;} - &Delta;t * (&part;gamma_dot_&alpha;\/&part;DGA_&beta;)\n         DO I=1,NSLIP\n            DO J=1,NSLIP\n               JAC(I,J) = -DTIME * DHDG(I,J)\n            ENDDO\n            JAC(I,I) = JAC(I,I) + 1.0\n         ENDDO\nC\nC        LU\u5206\u89e3\u5e76\u6c42\u89e3 \/ LU decomposition and solve\n         CALL LUDCMP(JAC, NSLIP, NSLIP, INDX, DET)\n         DO I=1,NSLIP\n            DGA_INC(I) = -RES(I)\n         ENDDO\n         CALL LUBKSB(JAC, NSLIP, NSLIP, INDX, DGA_INC)\nC\nC        \u66f4\u65b0\u672a\u77e5\u91cf \/ Update unknowns\n         DO I=1,NSLIP\n            DGA(I) = DGA(I) + DGA_INC(I)\n         ENDDO\n      ENDDO\nC\n      WRITE(*,*) 'WARNING: NR iteration not converged in UMAT'\n      RETURN\n      END<\/code><\/pre>\n\n  <h4>\n    <span data-lang=\"zh\">5.3.5 DDSDDE\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6 \/ Consistent Tangent DDSDDE<\/span>\n    <span data-lang=\"en\" style=\"display:none\">5.3.5 Consistent Tangent DDSDDE<\/span>\n  <\/h4>\n  <pre><code>C-----------------------------------------------------------------------\nC     \u8ba1\u7b97\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635 DDSDDE \/ Compute consistent tangent modulus\nC     DDSDDE = &part;&Delta;&sigma; \/ &part;&Delta;&epsilon;\nC-----------------------------------------------------------------------\n      SUBROUTINE CONSISTENTTANGENT(...)\n      IMPLICIT NONE\n      INTEGER NTENS, NSLIP, I, J, K, L\n      REAL*8 DDSDDE(NTENS,NTENS), CEL(6,6)\n      REAL*8 P_ALPHA(3,3,NSLIP), DGA(NSLIP)\n      REAL*8 DTAU_DSIG(NTENS,NSLIP), DHDG(NSLIP,NSLIP)\n      REAL*8 JAC(NSLIP,NSLIP), INV_JAC(NSLIP,NSLIP)\n      REAL*8 TERM(NTENS,NSLIP), SUM\n      INTEGER INDX(NSLIP)\n      REAL*8 DET\nC\nC     \u5f39\u6027\u521a\u5ea6\u77e9\u9635 \/ Elastic stiffness matrix\n      DO I=1,NTENS\n         DO J=1,NTENS\n            DDSDDE(I,J) = CEL(I,J)\n         ENDDO\n      ENDDO\nC\nC     \u8ba1\u7b97 &part;tau_&alpha;\/&part;&sigma; \u5e76\u7ec4\u88c5\u4fee\u6b63\u9879\nC     Compute &part;tau_&alpha;\/&part;&sigma; and assemble correction\n      DO K=1,NSLIP\n         DO I=1,NTENS\n            TERM(I,K) = 0.0\n            DO J=1,NTENS\n               TERM(I,K) = TERM(I,K) + CEL(I,J)*P_ALPHA_VOIGT(J,K)\n            ENDDO\n         ENDDO\n      ENDDO\nC\nC     \u8ba1\u7b97 (I - &Delta;t * &part;dg\/&part;g * &part;g\/&part;DGA)^{-1}\n      CALL LUDCMP(JAC, NSLIP, NSLIP, INDX, DET)\n      CALL MATINV(JAC, NSLIP, NSLIP, INDX, INV_JAC)\nC\nC     DDSDDE = C^e - C^e:P * INV_JAC * DTIME * dg\/dtau * d tau\/d eps\n      DO I=1,NTENS\n         DO J=1,NTENS\n            SUM = 0.0\n            DO K=1,NSLIP\n               DO L=1,NSLIP\n                  SUM = SUM + TERM(I,K)*INV_JAC(K,L)*DTAU_DSIG(L,J)\n               ENDDO\n            ENDDO\n            DDSDDE(I,J) = DDSDDE(I,J) - SUM * DTIME\n         ENDDO\n      ENDDO\n      RETURN\n      END<\/code><\/pre>\n<\/div>\n\n<!-- ==================== CHAPTER 6 ==================== -->\n<div class=\"chapter\" id=\"ch6\">\n  <h2>\n    <span class=\"chap-num\">06<\/span>\n    <span data-lang=\"zh\">Jacobian\u77e9\u9635 \/ Jacobian Matrix<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Jacobian Matrix<\/span>\n  <\/h2>\n\n  <h3 id=\"ch6-1\">\n    <span data-lang=\"zh\">6.1 \u5e94\u53d8\u5206\u89e3 \/ Strain Decomposition<\/span>\n    <span data-lang=\"en\" style=\"display:none\">6.1 Strain Decomposition<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u5728Abaqus\/Standard\u7684\u9690\u5f0f\u6c42\u89e3\u6846\u67b6\u4e2d\uff0cUMAT\u5fc5\u987b\u63d0\u4f9b\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635DDSDDE\uff0c\n    \u5373\u5e94\u529b\u589e\u91cf\u5bf9\u5e94\u53d8\u589e\u91cf\u7684\u504f\u5bfc\u6570DDSDDE = &part;&Delta;&sigma; \/ &part;&Delta;&epsilon;\u3002\n    \u5bf9\u4e8e\u5f39\u5851\u6027\u6750\u6599\uff0c\u603b\u5e94\u53d8\u589e\u91cf\u5206\u89e3\u4e3a\u5f39\u6027\u90e8\u5206\u548c\u5851\u6027\u90e8\u5206\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In the Abaqus\/Standard implicit solver framework, UMAT must provide the consistent tangent stiffness matrix DDSDDE,\n    i.e., the partial derivative of stress increment with respect to strain increment, DDSDDE = &part;&Delta;&sigma; \/ &part;&Delta;&epsilon;.\n    For elastoplastic materials, the total strain increment is decomposed into elastic and plastic parts:\n  <\/p>\n  <div class=\"formula\">\n    <b>&Delta;&epsilon; = &Delta;&epsilon;<sup>e<\/sup> + &Delta;&epsilon;<sup>p<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5728\u5c0f\u53d8\u5f62\u5047\u8bbe\u4e0b\uff0c\u5f39\u6027\u5e94\u53d8\u589e\u91cf\u4e0e\u5e94\u529b\u589e\u91cf\u7684\u5173\u7cfb\u7531\u5e7f\u4e49\u80e1\u514b\u5b9a\u5f8b\u7ed9\u51fa\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Under the small-deformation assumption, the relationship between elastic strain increment and stress increment is given by the generalized Hooke&#8217;s law:\n  <\/p>\n  <div class=\"formula\">\n    <b>&Delta;&sigma; = C<sup>e<\/sup> : &Delta;&epsilon;<sup>e<\/sup> = C<sup>e<\/sup> : (&Delta;&epsilon; &#8211; &Delta;&epsilon;<sup>p<\/sup>)<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5851\u6027\u5e94\u53d8\u589e\u91cf\u7531\u5404\u6ed1\u79fb\u7cfb\u7684\u526a\u5207\u5e94\u53d8\u589e\u91cf\u8d21\u732e\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The plastic strain increment is contributed by the shear strain increments on all slip systems:\n  <\/p>\n  <div class=\"formula\">\n    <b>&Delta;&epsilon;<sup>p<\/sup> = &Sigma;<sub>&alpha;<\/sub> &Delta;&gamma;<sup>(&alpha;)<\/sup> P<sup>(&alpha;)<\/sup><\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u56e0\u6b64\uff0c\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635\u7684\u5173\u952e\u5728\u4e8e\u8ba1\u7b97 &part;&Delta;&gamma;<sup>(&alpha;)<\/sup> \/ &part;&Delta;&epsilon;\uff0c\n    \u5373\u526a\u5207\u5e94\u53d8\u589e\u91cf\u5bf9\u5e94\u53d8\u589e\u91cf\u7684\u654f\u611f\u5ea6\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Therefore, the key to the consistent tangent stiffness matrix is computing &part;&Delta;&gamma;<sup>(&alpha;)<\/sup> \/ &part;&Delta;&epsilon;,\n    i.e., the sensitivity of shear strain increment to strain increment.\n  <\/p>\n\n  <h3 id=\"ch6-2\">\n    <span data-lang=\"zh\">6.2 \u7ebf\u6027\u65b9\u7a0b\u7ec4 \/ Linear System<\/span>\n    <span data-lang=\"en\" style=\"display:none\">6.2 Linear System<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u5728Newton-Raphson\u8fed\u4ee3\u4e2d\uff0c\u6bcf\u4e2a\u6ed1\u79fb\u7cfb\u7684\u526a\u5207\u5e94\u53d8\u589e\u91cf\u6ee1\u8db3\u4ee5\u4e0b\u975e\u7ebf\u6027\u65b9\u7a0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In the Newton-Raphson iteration, the shear strain increment for each slip system satisfies the following nonlinear equation:\n  <\/p>\n  <div class=\"formula\">\n    <b>R<sup>(&alpha;)<\/sup> = &Delta;&gamma;<sup>(&alpha;)<\/sup> &#8211; &Delta;t &middot; &gamma;&#775;<sub>0<\/sub> &middot; | &tau;<sup>(&alpha;)<\/sup> \/ g<sup>(&alpha;)<\/sup> |<sup>n<\/sup> &middot; sign(&tau;<sup>(&alpha;)<\/sup>\/g<sup>(&alpha;)<\/sup>) = 0<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5bf9\u4e0a\u8ff0\u65b9\u7a0b\u5173\u4e8e&Delta;&gamma;<sup>(&beta;)<\/sup>\u6c42\u504f\u5bfc\uff0c\u5f97\u5230Newton-Raphson\u8fed\u4ee3\u7684Jacobian\u77e9\u9635\uff08\u5c40\u90e8Newton\u8fed\u4ee3\u77e9\u9635\uff09\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Taking the partial derivative of the above equation with respect to &Delta;&gamma;<sup>(&beta;)<\/sup> yields the Jacobian matrix for Newton-Raphson iteration (local Newton iteration matrix):\n  <\/p>\n  <div class=\"formula\">\n    <b>J<sub>&alpha;&beta;<\/sub> = &part;R<sup>(&alpha;)<\/sup> \/ &part;&Delta;&gamma;<sup>(&beta;)<\/sup> = &delta;<sub>&alpha;&beta;<\/sub> &#8211; &Delta;t &middot; (&part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;&tau;<sup>(&alpha;)<\/sup>) &middot; (&part;&tau;<sup>(&alpha;)<\/sup> \/ &part;&Delta;&gamma;<sup>(&beta;)<\/sup>) &#8211; &Delta;t &middot; (&part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;g<sup>(&alpha;)<\/sup>) &middot; (&part;g<sup>(&alpha;)<\/sup> \/ &part;&Delta;&gamma;<sup>(&beta;)<\/sup>)<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u5404\u9879\u504f\u5bfc\u6570\u7684\u7269\u7406\u542b\u4e49\u548c\u8868\u8fbe\u5f0f\u5982\u4e0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The physical meaning and expressions of each partial derivative are:\n  <\/p>\n  <ul>\n    <li>\n      <b>&part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;&tau;<sup>(&alpha;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u6d41\u52a8\u5f8b\u5bf9\u5206\u89e3\u526a\u5e94\u529b\u7684\u654f\u611f\u5ea6\uff0c\u7b49\u4e8e n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ &tau;<sup>(&alpha;)<\/sup><\/span>\n      <span data-lang=\"en\" style=\"display:none\">Sensitivity of flow rule to resolved shear stress, equal to n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ &tau;<sup>(&alpha;)<\/sup><\/span>\n    <\/li>\n    <li>\n      <b>&part;&tau;<sup>(&alpha;)<\/sup> \/ &part;&Delta;&gamma;<sup>(&beta;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u5206\u89e3\u526a\u5e94\u529b\u5bf9\u526a\u5207\u589e\u91cf\u7684\u654f\u611f\u5ea6\uff0c\u7b49\u4e8e -C<sup>e<\/sup> : P<sup>(&beta;)<\/sup> : P<sup>(&alpha;)<\/sup>\uff08\u5e94\u529b\u677e\u5f1b\u6548\u5e94\uff09<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Sensitivity of resolved shear stress to shear increment, equal to -C<sup>e<\/sup> : P<sup>(&beta;)<\/sup> : P<sup>(&alpha;)<\/sup> (stress relaxation effect)<\/span>\n    <\/li>\n    <li>\n      <b>&part;g<sup>(&alpha;)<\/sup> \/ &part;&Delta;&gamma;<sup>(&beta;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u6ed1\u79fb\u7cfb\u5f3a\u5ea6\u5bf9\u526a\u5207\u589e\u91cf\u7684\u654f\u611f\u5ea6\uff0c\u7b49\u4e8e h<sub>&alpha;&beta;<\/sub> &middot; sign(&tau;<sup>(&beta;)<\/sup>)<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Sensitivity of slip system strength to shear increment, equal to h<sub>&alpha;&beta;<\/sub> &middot; sign(&tau;<sup>(&beta;)<\/sup>)<\/span>\n    <\/li>\n    <li>\n      <b>&part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;g<sup>(&alpha;)<\/sup><\/b>:\n      <span data-lang=\"zh\">\u6d41\u52a8\u5f8b\u5bf9\u5f3a\u5ea6\u7684\u654f\u611f\u5ea6\uff0c\u7b49\u4e8e -n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ g<sup>(&alpha;)<\/sup><\/span>\n      <span data-lang=\"en\" style=\"display:none\">Sensitivity of flow rule to strength, equal to -n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ g<sup>(&alpha;)<\/sup><\/span>\n    <\/li>\n  <\/ul>\n  <p data-lang=\"zh\">\n    \u5c06\u6240\u6709\u504f\u5bfc\u6570\u4ee3\u5165\u540e\uff0c\u5c40\u90e8Newton\u77e9\u9635\u53ef\u663e\u5f0f\u5199\u51fa\u3002\u8be5\u77e9\u9635\u7684\u5927\u5c0f\u4e3aNSLIP &times; NSLIP\uff08\u4f8b\u5982FCC\u4e3a12&times;12\uff09\u3002\n    \u5728\u6bcf\u6b21Newton\u8fed\u4ee3\u4e2d\uff0c\u5bf9\u8be5\u77e9\u9635\u8fdb\u884cLU\u5206\u89e3\u5e76\u56de\u4ee3\u6c42\u89e3\u526a\u5207\u5e94\u53d8\u589e\u91cf\u4fee\u6b63\u503c&Delta;DGA\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    After substituting all partial derivatives, the local Newton matrix can be written explicitly. Its size is NSLIP &times; NSLIP (e.g., 12&times;12 for FCC).\n    In each Newton iteration, this matrix undergoes LU decomposition and back-substitution to solve for the shear strain increment correction &Delta;DGA.\n  <\/p>\n\n  <h3 id=\"ch6-3\">\n    <span data-lang=\"zh\">6.3 DDSDDE\u6784\u9020 \/ DDSDDE Construction<\/span>\n    <span data-lang=\"en\" style=\"display:none\">6.3 DDSDDE Construction<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635\u7684\u6784\u9020\u516c\u5f0f\u4e3a\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The consistent tangent stiffness matrix is constructed as:\n  <\/p>\n  <div class=\"formula\">\n    <b>DDSDDE = C<sup>e<\/sup> &#8211; C<sup>e<\/sup> : P<sup>(&alpha;)<\/sup> &otimes; [ J<sup>-1<\/sup><sub>&alpha;&beta;<\/sub> &middot; &Delta;t &middot; (&part;&gamma;&#775;<sup>(&beta;)<\/sup> \/ &part;&tau;<sup>(&beta;)<\/sup>) &middot; P<sup>(&beta;)<\/sup> : C<sup>e<\/sup> ]<\/b>\n  <\/div>\n  <p data-lang=\"zh\">\n    \u4e0a\u8ff0\u516c\u5f0f\u4e2d\uff0c\u7b2c\u4e00\u9879C<sup>e<\/sup>\u662f\u7eaf\u5f39\u6027\u521a\u5ea6\uff0c\u7b2c\u4e8c\u9879\u662f\u5851\u6027\u4fee\u6b63\u3002\n    \u5f53\u6750\u6599\u5904\u4e8e\u7eaf\u5f39\u6027\u72b6\u6001\uff08\u65e0\u6ed1\u79fb\u6fc0\u6d3b\uff09\u65f6\uff0c\u7b2c\u4e8c\u9879\u4e3a\u96f6\uff0cDDSDDE\u9000\u5316\u4e3a\u5f39\u6027\u521a\u5ea6\u77e9\u9635\u3002\n    \u5f53\u5927\u91cf\u6ed1\u79fb\u7cfb\u6fc0\u6d3b\u65f6\uff0c\u7b2c\u4e8c\u9879\u663e\u8457\u964d\u4f4e\u6750\u6599\u7684\u8868\u89c2\u521a\u5ea6\uff0c\u53cd\u6620\u5851\u6027\u6d41\u52a8\u5f15\u8d77\u7684\u5e94\u529b\u677e\u5f1b\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    In the above formula, the first term C<sup>e<\/sup> is the pure elastic stiffness, and the second term is the plastic correction.\n    When the material is in a purely elastic state (no slip activated), the second term vanishes and DDSDDE reduces to the elastic stiffness matrix.\n    When many slip systems are active, the second term significantly reduces the apparent material stiffness, reflecting stress relaxation due to plastic flow.\n  <\/p>\n  <p data-lang=\"zh\">\n    \u5177\u4f53\u8ba1\u7b97\u6b65\u9aa4\u5982\u4e0b\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The specific calculation steps are:\n  <\/p>\n  <ol>\n    <li>\n      <span data-lang=\"zh\">\u8ba1\u7b97\u5f39\u6027\u521a\u5ea6\u77e9\u9635C<sup>e<\/sup>\u7684Voigt\u5f62\u5f0f\uff086&times;6\u77e9\u9635\uff09<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Compute the Voigt form of elastic stiffness matrix C<sup>e<\/sup> (6&times;6 matrix)<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u8ba1\u7b97\u6bcf\u4e2a\u6ed1\u79fb\u7cfb\u7684Schmid\u5f20\u91cfP<sup>(&alpha;)<\/sup>\u7684Voigt\u5411\u91cf\u5f62\u5f0f<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Compute the Voigt vector form of Schmid tensor P<sup>(&alpha;)<\/sup> for each slip system<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u8ba1\u7b97\u6d41\u52a8\u5f8b\u5bfc\u6570 &part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;&tau;<sup>(&alpha;)<\/sup> = n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ &tau;<sup>(&alpha;)<\/sup><\/span>\n      <span data-lang=\"en\" style=\"display:none\">Compute flow rule derivative &part;&gamma;&#775;<sup>(&alpha;)<\/sup> \/ &part;&tau;<sup>(&alpha;)<\/sup> = n &middot; &gamma;&#775;<sup>(&alpha;)<\/sup> \/ &tau;<sup>(&alpha;)<\/sup><\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u7ec4\u88c5\u5c40\u90e8Newton\u77e9\u9635J<sub>&alpha;&beta;<\/sub>\u5e76\u8fdb\u884cLU\u5206\u89e3\u6c42\u9006<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Assemble local Newton matrix J<sub>&alpha;&beta;<\/sub> and perform LU decomposition for inversion<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u8ba1\u7b97\u4e2d\u95f4\u5f20\u91cf\uff1a\u5bf9\u6bcf\u4e2a\u6ed1\u79fb\u7cfb&alpha;\uff0c\u8ba1\u7b97 Q<sup>(&alpha;)<\/sup><sub>ij<\/sub> = (J<sup>-1<\/sup><sub>&alpha;&beta;<\/sub> &middot; DFLOW<sub>&beta;<\/sub>) &middot; (P<sup>(&beta;)<\/sup><sub>kl<\/sub> C<sup>e<\/sup><sub>klmn<\/sub>)<\/span>\n      <span data-lang=\"en\" style=\"display:none\">Compute intermediate tensor: for each slip system &alpha;, compute Q<sup>(&alpha;)<\/sup><sub>ij<\/sub> = (J<sup>-1<\/sup><sub>&alpha;&beta;<\/sub> &middot; DFLOW<sub>&beta;<\/sub>) &middot; (P<sup>(&beta;)<\/sup><sub>kl<\/sub> C<sup>e<\/sup><sub>klmn<\/sub>)<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\">\u7ec4\u88c5DDSDDE\uff1aDDSDDE<sub>ij<\/sub> = C<sup>e<\/sup><sub>ij<\/sub> &#8211; &Sigma;<sub>&alpha;<\/sub> (C<sup>e<\/sup><sub>ijkl<\/sub> P<sup>(&alpha;)<\/sup><sub>kl<\/sub>) Q<sup>(&alpha;)<\/sup><sub>mn<\/sub><\/span>\n      <span data-lang=\"en\" style=\"display:none\">Assemble DDSDDE: DDSDDE<sub>ij<\/sub> = C<sup>e<\/sup><sub>ij<\/sub> &#8211; &Sigma;<sub>&alpha;<\/sub> (C<sup>e<\/sup><sub>ijkl<\/sub> P<sup>(&alpha;)<\/sup><sub>kl<\/sub>) Q<sup>(&alpha;)<\/sup><sub>mn<\/sub><\/span>\n    <\/li>\n  <\/ol>\n  <div class=\"note\">\n    <span data-lang=\"zh\"><strong>\u6570\u503c\u7a33\u5b9a\u6027\u63d0\u793a\uff1a<\/strong>\u5f53&tau;<sup>(&alpha;)<\/sup>\u63a5\u8fd1\u96f6\u65f6\uff0c\u6d41\u52a8\u5f8b\u5bfc\u6570n&middot;&gamma;&#775;\/&tau;\u53ef\u80fd\u5947\u5f02\u3002\u7a0b\u5e8f\u4e2d\u901a\u5e38\u5f15\u5165\u5fae\u5c0f\u6b63\u5219\u5316\u9879\u6216\u9650\u5236\u6700\u5927\u5bfc\u6570\u503c\u3002\u6b64\u5916\uff0c\u5f53\u67d0\u4e9b\u6ed1\u79fb\u7cfb\u7684\u526a\u5207\u7387\u6781\u5c0f\u65f6\uff0c\u53ef\u5c06\u5176\u89c6\u4e3a&#8221;\u975e\u6fc0\u6d3b&#8221;\u5e76\u4ece\u5c40\u90e8Newton\u7cfb\u7edf\u4e2d\u6d88\u53bb\uff0c\u4ee5\u51cf\u5c0f\u77e9\u9635\u89c4\u6a21\u5e76\u6539\u5584\u6761\u4ef6\u6570\u3002<\/span>\n    <span data-lang=\"en\" style=\"display:none\"><strong>Numerical stability tip:<\/strong> When &tau;<sup>(&alpha;)<\/sup> is near zero, the flow rule derivative n&middot;&gamma;&#775;\/&tau; may become singular. The program usually introduces a small regularization term or caps the maximum derivative value. Additionally, slip systems with extremely small shear rates can be treated as &#8220;inactive&#8221; and eliminated from the local Newton system to reduce matrix size and improve conditioning.<\/span>\n  <\/div>\n  <p data-lang=\"zh\">\n    DDSDDE\u7684\u5bf9\u79f0\u6027\u53d6\u51b3\u4e8e\u6750\u6599\u6a21\u578b\u3002\u5bf9\u4e8e\u5173\u8054\u6d41\u52a8\u6cd5\u5219\u548c\u4e8c\u6b21\u5c48\u670d\u9762\uff0cDDSDDE\u662f\u5bf9\u79f0\u7684\uff1b\n    \u4f46\u5bf9\u4e8e\u975e\u5173\u8054\u6d41\u52a8\u6216\u590d\u6742\u7684\u6f5c\u786c\u5316\u6a21\u578b\uff0cDDSDDE\u53ef\u80fd\u975e\u5bf9\u79f0\u3002\n    \u9ec4\u6c38\u521aUMAT\u4e2d\uff0c\u7531\u4e8e\u91c7\u7528\u5173\u8054\u7684Schmid\u6d41\u52a8\u548c\u4e00\u822c\u7684\u786c\u5316\u77e9\u9635\uff0cDDSDDE\u901a\u5e38\u5177\u6709\u8f7b\u5fae\u7684\u975e\u5bf9\u79f0\u6027\u3002\n    Abaqus\/Standard\u652f\u6301\u975e\u5bf9\u79f0\u521a\u5ea6\u77e9\u9635\uff0c\u4f46\u9700\u8981\u5728*Step\u4e2d\u8bbe\u7f6eUNSymm=YES\u4ee5\u4f18\u5316\u6c42\u89e3\u5668\u6027\u80fd\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The symmetry of DDSDDE depends on the material model. For associative flow rules and quadratic yield surfaces, DDSDDE is symmetric;\n    but for non-associative flow or complex latent hardening models, DDSDDE may be non-symmetric.\n    In Huang&#8217;s UMAT, due to the associative Schmid flow and general hardening matrix, DDSDDE usually has slight non-symmetry.\n    Abaqus\/Standard supports non-symmetric stiffness matrices, but *Step requires UNSymm=YES to optimize solver performance.\n  <\/p>\n<\/div>\n\n<!-- ==================== CHAPTER 7 ==================== -->\n<div class=\"chapter\" id=\"ch7\">\n  <h2>\n    <span class=\"chap-num\">07<\/span>\n    <span data-lang=\"zh\">\u5178\u578b\u53c2\u6570 \/ Typical Parameters<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Typical Parameters<\/span>\n  <\/h2>\n\n  <p data-lang=\"zh\">\n    \u4ee5\u4e0b\u7ed9\u51faFCC\u94dc\u548cFCC\u94dd\u7684\u5178\u578b\u6750\u6599\u53c2\u6570\uff0c\u8fd9\u4e9b\u53c2\u6570\u5728\u5927\u91cf\u6587\u732e\u548c\u5b9e\u9a8c\u4e2d\u5f97\u5230\u9a8c\u8bc1\uff0c\n    \u53ef\u4f5c\u4e3aAbaqus\u8f93\u5165\u6587\u4ef6\u4e2d*User Material\u7684\u53c2\u8003\u503c\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The following presents typical material parameters for FCC copper and FCC aluminum,\n    validated in numerous literature and experiments, serving as reference values for the *User Material keyword in Abaqus input files.\n  <\/p>\n\n  <h3 id=\"ch7-1\">\n    <span data-lang=\"zh\">7.1 FCC\u94dc \/ FCC Copper<\/span>\n    <span data-lang=\"en\" style=\"display:none\">7.1 FCC Copper<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u94dc\u662fFCC\u7ed3\u6784\uff0c\u5ba4\u6e29\u4e0b\u4e3b\u8981\u6ed1\u79fb\u7cfb\u4e3a{111}&lt;110&gt;\uff0c\u517112\u4e2a\u3002\n    \u4ee5\u4e0b\u53c2\u6570\u9002\u7528\u4e8e\u5ba4\u6e29\uff08\u7ea6300K\uff09\u4e0b\u7684\u5355\u6676\u94dc\u6a21\u62df\uff0c\u91c7\u7528\u53cc\u66f2\u6b63\u5272\u786c\u5316\u6a21\u578b\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Copper has an FCC structure. At room temperature, the primary slip systems are {111}&lt;110&gt;, totaling 12.\n    The following parameters are suitable for single-crystal copper simulation at room temperature (~300K) using the hyperbolic secant hardening model.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>\n          <span data-lang=\"zh\">\u53c2\u6570 \/ Parameter<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Parameter<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u6570\u503c \/ Value<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Value<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u5355\u4f4d \/ Unit<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Unit<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr>\n        <td>C11<\/td>\n        <td>168,400<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C11<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C11<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>C12<\/td>\n        <td>121,400<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C12<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C12<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>C44<\/td>\n        <td>75,400<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C44<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C44<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&gamma;&#775;<sub>0<\/sub><\/td>\n        <td>0.001<\/td>\n        <td>s<sup>-1<\/sup><\/td>\n        <td><span data-lang=\"zh\">\u53c2\u8003\u526a\u5207\u7387<\/span><span data-lang=\"en\" style=\"display:none\">Reference shear rate<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>n<\/td>\n        <td>50<\/td>\n        <td>&#8211;<\/td>\n        <td><span data-lang=\"zh\">\u7387\u654f\u611f\u6307\u6570<\/span><span data-lang=\"en\" style=\"display:none\">Rate sensitivity exponent<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&tau;<sub>0<\/sub><\/td>\n        <td>16.0<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u521d\u59cbCRSS<\/span><span data-lang=\"en\" style=\"display:none\">Initial CRSS<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&tau;<sub>s<\/sub><\/td>\n        <td>148.0<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u9971\u548cCRSS<\/span><span data-lang=\"en\" style=\"display:none\">Saturation CRSS<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>H<sub>0<\/sub><\/td>\n        <td>180.0<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u521d\u59cb\u786c\u5316\u6a21\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Initial hardening modulus<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>q<\/td>\n        <td>1.0<\/td>\n        <td>&#8211;<\/td>\n        <td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570\uff08\u5404\u5411\u540c\u6027\u8fd1\u4f3c\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient (isotropic approximation)<\/span><\/td>\n      <\/tr>\n    <\/tbody>\n  <\/table>\n\n  <h3 id=\"ch7-2\">\n    <span data-lang=\"zh\">7.2 FCC\u94dd \/ FCC Aluminum<\/span>\n    <span data-lang=\"en\" style=\"display:none\">7.2 FCC Aluminum<\/span>\n  <\/h3>\n  <p data-lang=\"zh\">\n    \u94dd\u540c\u6837\u662fFCC\u7ed3\u6784\uff0c\u4f46\u5176\u5c42\u9519\u80fd\u8f83\u9ad8\uff0c\u4f4d\u9519\u4ea4\u6ed1\u79fb\u66f4\u5bb9\u6613\u53d1\u751f\uff0c\u56e0\u6b64\u786c\u5316\u884c\u4e3a\u548c\u94dc\u6709\u660e\u663e\u5dee\u5f02\u3002\n    \u4ee5\u4e0b\u53c2\u6570\u9002\u7528\u4e8e\u9ad8\u7eaf\u94dd\u5355\u6676\u5728\u5ba4\u6e29\u4e0b\u7684\u6a21\u62df\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    Aluminum is also FCC, but with higher stacking fault energy, making dislocation cross-slip easier,\n    so its hardening behavior differs significantly from copper.\n    The following parameters are suitable for high-purity single-crystal aluminum at room temperature.\n  <\/p>\n\n  <table>\n    <thead>\n      <tr>\n        <th>\n          <span data-lang=\"zh\">\u53c2\u6570 \/ Parameter<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Parameter<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u6570\u503c \/ Value<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Value<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u5355\u4f4d \/ Unit<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Unit<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u8bf4\u660e \/ Description<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Description<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr>\n        <td>C11<\/td>\n        <td>106,750<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C11<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C11<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>C12<\/td>\n        <td>60,410<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C12<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C12<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>C44<\/td>\n        <td>28,340<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u5f39\u6027\u5e38\u6570C44<\/span><span data-lang=\"en\" style=\"display:none\">Elastic constant C44<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&gamma;&#775;<sub>0<\/sub><\/td>\n        <td>0.001<\/td>\n        <td>s<sup>-1<\/sup><\/td>\n        <td><span data-lang=\"zh\">\u53c2\u8003\u526a\u5207\u7387<\/span><span data-lang=\"en\" style=\"display:none\">Reference shear rate<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>n<\/td>\n        <td>60<\/td>\n        <td>&#8211;<\/td>\n        <td><span data-lang=\"zh\">\u7387\u654f\u611f\u6307\u6570<\/span><span data-lang=\"en\" style=\"display:none\">Rate sensitivity exponent<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&tau;<sub>0<\/sub><\/td>\n        <td>0.75<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u521d\u59cbCRSS<\/span><span data-lang=\"en\" style=\"display:none\">Initial CRSS<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>&tau;<sub>s<\/sub><\/td>\n        <td>25.0<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u9971\u548cCRSS<\/span><span data-lang=\"en\" style=\"display:none\">Saturation CRSS<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>H<sub>0<\/sub><\/td>\n        <td>30.0<\/td>\n        <td>MPa<\/td>\n        <td><span data-lang=\"zh\">\u521d\u59cb\u786c\u5316\u6a21\u91cf<\/span><span data-lang=\"en\" style=\"display:none\">Initial hardening modulus<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td>q<\/td>\n        <td>1.0<\/td>\n        <td>&#8211;<\/td>\n        <td><span data-lang=\"zh\">\u6f5c\u786c\u5316\u7cfb\u6570\uff08\u5404\u5411\u540c\u6027\u8fd1\u4f3c\uff09<\/span><span data-lang=\"en\" style=\"display:none\">Latent hardening coefficient (isotropic approximation)<\/span><\/td>\n      <\/tr>\n    <\/tbody>\n  <\/table>\n\n  <div class=\"note\">\n    <span data-lang=\"zh\"><strong>\u53c2\u6570\u6765\u6e90\u8bf4\u660e\uff1a<\/strong>\u4e0a\u8ff0\u53c2\u6570\u7efc\u5408\u4e86\u591a\u9879\u7ecf\u5178\u5b9e\u9a8c\u548c\u6587\u732e\u6570\u636e\u3002\u5b9e\u9645\u5e94\u7528\u4e2d\uff0c\u5e94\u6839\u636e\u5177\u4f53\u5408\u91d1\u6210\u5206\u3001\u6e29\u5ea6\u3001\u5e94\u53d8\u901f\u7387\u7b49\u6761\u4ef6\u901a\u8fc7\u5355\u6676\u62c9\u4f38\u5b9e\u9a8c\u6821\u51c6\u3002\u7279\u522b\u5730\uff0c\u7387\u654f\u611f\u6307\u6570n\u548c\u53c2\u8003\u526a\u5207\u7387&gamma;&#775;<sub>0<\/sub>\u5e94\u5339\u914d\u5b9e\u9a8c\u5e94\u53d8\u901f\u7387\uff0c\u5426\u5219\u53ef\u80fd\u5bfc\u81f4\u5e94\u529b\u6c34\u5e73\u9884\u6d4b\u504f\u5dee\u3002<\/span>\n    <span data-lang=\"en\" style=\"display:none\"><strong>Parameter source note:<\/strong> The above parameters are synthesized from multiple classic experiments and literature data. In practical applications, they should be calibrated through single-crystal tensile experiments according to specific alloy composition, temperature, strain rate, etc. In particular, the rate sensitivity exponent n and reference shear rate &gamma;&#775;<sub>0<\/sub> should match the experimental strain rate; otherwise, stress level prediction errors may occur.<\/span>\n  <\/div>\n<\/div>\n<!-- ==================== CHAPTER 8 ==================== -->\n<div class=\"chapter\" id=\"ch8\">\n  <h2>\n    <span class=\"chap-num\">08<\/span>\n    <span data-lang=\"zh\">\u4f7f\u7528\u6ce8\u610f\u4e8b\u9879 \/ Usage Notes<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Usage Notes<\/span>\n  <\/h2>\n\n  <p data-lang=\"zh\">\n    \u4ee5\u4e0b\u5217\u51fa\u4f7f\u7528\u9ec4\u6c38\u521aUMAT\u5b50\u7a0b\u5e8f\u65f6\u5e38\u89c1\u7684\u6ce8\u610f\u4e8b\u9879\u548c\u6700\u4f73\u5b9e\u8df5\u5efa\u8bae\uff0c\n    \u6db5\u76d6\u8f93\u5165\u6587\u4ef6\u8bbe\u7f6e\u3001\u6750\u6599\u53c2\u6570\u6821\u51c6\u3001\u6536\u655b\u6027\u8c03\u8bd5\u548c\u540e\u5904\u7406\u7b49\u65b9\u9762\u3002\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The following lists common considerations and best-practice recommendations when using Huang&#8217;s UMAT subroutine,\n    covering input file setup, material parameter calibration, convergence debugging, and post-processing.\n  <\/p>\n\n  <h3 id=\"ch8-1\">\n    <span data-lang=\"zh\">8.1 Abaqus\u8f93\u5165\u6587\u4ef6\u8bbe\u7f6e \/ Input File Setup<\/span>\n    <span data-lang=\"en\" style=\"display:none\">8.1 Input File Setup<\/span>\n  <\/h3>\n  <ul>\n    <li>\n      <span data-lang=\"zh\"><strong>*User Material\u5173\u952e\u5b57\uff1a<\/strong>\u5fc5\u987b\u5728\u8f93\u5165\u6587\u4ef6\u4e2d\u901a\u8fc7*User Material, constants=NPROPS\u5b9a\u4e49PROPS\u6570\u7ec4\u957f\u5ea6\uff0c\u5e76\u786e\u4fdd*Depvar\u7684NSTATV\u8db3\u591f\u5927\u4ee5\u5bb9\u7eb3\u6240\u6709\u72b6\u6001\u53d8\u91cf\u3002\u5bf9\u4e8eFCC 12\u6ed1\u79fb\u7cfb\uff0c\u5efa\u8baeNSTATV >= 70\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>*User Material keyword:<\/strong> The PROPS array length must be defined via *User Material, constants=NPROPS, and *Depvar&#8217;s NSTATV must be large enough for all state variables. For FCC 12 slip systems, NSTATV >= 70 is recommended.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>UNSymm\u9009\u9879\uff1a<\/strong>\u7531\u4e8e\u6676\u4f53\u5851\u6027\u7684\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635\u901a\u5e38\u8f7b\u5fae\u975e\u5bf9\u79f0\uff0c\u5efa\u8bae\u5728*Step\u4e2d\u4f7f\u7528UNSymm=YES\u4ee5\u542f\u7528\u975e\u5bf9\u79f0\u6c42\u89e3\u5668\uff0c\u63d0\u9ad8\u6536\u655b\u901f\u5ea6\u548c\u7a33\u5b9a\u6027\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>UNSymm option:<\/strong> Because the consistent tangent stiffness matrix of crystal plasticity is usually slightly non-symmetric, it is recommended to use UNSymm=YES in *Step to enable the non-symmetric solver, improving convergence speed and stability.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u65f6\u95f4\u6b65\u957f\u63a7\u5236\uff1a<\/strong>\u6676\u4f53\u5851\u6027UMAT\u5bf9\u65f6\u95f4\u6b65\u957f\u654f\u611f\u3002\u5efa\u8bae\u5728\u521d\u59cb\u52a0\u8f7d\u9636\u6bb5\u4f7f\u7528\u8f83\u5c0f\u7684\u65f6\u95f4\u6b65\uff0c\u5f85\u5851\u6027\u53d8\u5f62\u7a33\u5b9a\u540e\u53ef\u9002\u5ea6\u589e\u5927\u3002\u53ef\u901a\u8fc7*Static\u4e2d\u7684INITIAL\u548cMAXIMUM\u53c2\u6570\u63a7\u5236\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Time step control:<\/strong> Crystal plasticity UMAT is sensitive to time step size. It is recommended to use small time steps in the initial loading stage, and moderately increase after plastic deformation stabilizes. Control via INITIAL and MAXIMUM in *Static.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u5355\u5143\u7c7b\u578b\u9009\u62e9\uff1a<\/strong>\u63a8\u8350\u4f7f\u7528C3D8R\uff088\u8282\u70b9\u516d\u9762\u4f53\u51cf\u7f29\u79ef\u5206\u5355\u5143\uff09\u6216C3D8\uff08\u5b8c\u5168\u79ef\u5206\uff09\u3002\u907f\u514d\u4f7f\u7528\u7ebf\u6027\u4e09\u89d2\u5f62\/\u56db\u9762\u4f53\u5355\u5143\u6a21\u62df\u5927\u5e94\u53d8\u5c40\u90e8\u5316\uff0c\u56e0\u4e3a\u5b83\u4eec\u5728\u4e0d\u53ef\u538b\u7f29\u5851\u6027\u53d8\u5f62\u4e2d\u4f1a\u4ea7\u751f\u4f53\u79ef\u9501\u5b9a\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Element type selection:<\/strong> C3D8R (8-node hexahedron with reduced integration) or C3D8 (full integration) are recommended. Avoid linear triangular\/tetrahedral elements for large-strain localization, as they suffer volumetric locking in incompressible plastic deformation.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u6c99\u6f0f\u63a7\u5236\uff1a<\/strong>\u4f7f\u7528C3D8R\u5355\u5143\u65f6\uff0c\u5fc5\u987b\u914d\u5408*Hourglass Stiffness\u6216*Section Controls\u8bbe\u7f6e\u5408\u7406\u7684\u6c99\u6f0f\u521a\u5ea6\uff0c\u4ee5\u6291\u5236\u96f6\u80fd\u6a21\u5f0f\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Hourglass control:<\/strong> When using C3D8R elements, *Hourglass Stiffness or *Section Controls must be set with reasonable hourglass stiffness to suppress zero-energy modes.<\/span>\n    <\/li>\n  <\/ul>\n\n  <h3 id=\"ch8-2\">\n    <span data-lang=\"zh\">8.2 \u6750\u6599\u53c2\u6570\u6821\u51c6 \/ Parameter Calibration<\/span>\n    <span data-lang=\"en\" style=\"display:none\">8.2 Parameter Calibration<\/span>\n  <\/h3>\n  <ul>\n    <li>\n      <span data-lang=\"zh\"><strong>\u5f39\u6027\u5e38\u6570\uff1a<\/strong>C11\u3001C12\u3001C44\u5e94\u4ece\u8d85\u58f0\u5171\u632f\u6216\u5355\u6676\u5f39\u6027\u5b9e\u9a8c\u83b7\u53d6\u3002\u6ce8\u610fAbaqus\u4e2d\u5e94\u529b\u5355\u4f4d\u4e3aMPa\uff0c\u5f39\u6027\u5e38\u6570\u4e5f\u5fc5\u987b\u4ee5MPa\u8f93\u5165\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Elastic constants:<\/strong> C11, C12, C44 should be obtained from ultrasonic resonance or single-crystal elastic experiments. Note that Abaqus uses MPa for stress, so elastic constants must also be entered in MPa.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u521d\u59cbCRSS tau0\uff1a<\/strong>\u53ef\u901a\u8fc7\u5355\u6676\u521d\u59cb\u5c48\u670d\u5e94\u529b\u7ed3\u5408Schmid\u56e0\u5b50\u4f30\u7b97\u3002\u5bf9\u4e8e\u4e0d\u540c\u53d6\u5411\u7684\u5355\u6676\u8bd5\u6837\uff0ctau0\u5e94\u552f\u4e00\uff0c\u800c\u5b8f\u89c2\u5c48\u670d\u5e94\u529b\u968f\u53d6\u5411\u53d8\u5316\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Initial CRSS tau0:<\/strong> Can be estimated from single-crystal initial yield stress combined with the Schmid factor. For single-crystal specimens of different orientations, tau0 should be unique, while macroscopic yield stress varies with orientation.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u7387\u654f\u611f\u6307\u6570n\uff1a<\/strong>\u5e94\u901a\u8fc7\u4e0d\u540c\u5e94\u53d8\u901f\u7387\u4e0b\u7684\u5355\u6676\u62c9\u4f38\u5b9e\u9a8c\u786e\u5b9a\u3002n\u503c\u5bf9\u6d41\u52a8\u5e94\u529b\u6c34\u5e73\u6709\u663e\u8457\u5f71\u54cd\uff0c\u4f46\u5bf9\u786c\u5316\u6a21\u91cf\u5f71\u54cd\u8f83\u5c0f\u3002\u82e5\u5b9e\u9a8c\u6570\u636e\u4e0d\u8db3\uff0c\u53ef\u53c2\u8003\u6587\u732e\u4e2d\u540c\u79cd\u6750\u6599\u7684\u63a8\u8350\u503c\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Rate sensitivity exponent n:<\/strong> Should be determined from single-crystal tensile experiments at different strain rates. n significantly affects flow stress level but has minor effect on hardening modulus. If experimental data is insufficient, refer to recommended values in the literature for the same material.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u786c\u5316\u53c2\u6570\uff1a<\/strong>H0\u3001taus\u548ctau0\u5e94\u901a\u8fc7\u5355\u6676\u5e94\u529b-\u5e94\u53d8\u66f2\u7ebf\u7684\u5851\u6027\u6bb5\u62df\u5408\u786e\u5b9a\u3002\u5efa\u8bae\u4f7f\u7528\u975e\u7ebf\u6027\u6700\u5c0f\u4e8c\u4e58\u6cd5\u5bf9\u5b9e\u9a8c\u6570\u636e\u8fdb\u884c\u4f18\u5316\u62df\u5408\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Hardening parameters:<\/strong> H0, taus, and tau0 should be determined by fitting the plastic portion of single-crystal stress-strain curves. Nonlinear least-squares optimization is recommended for fitting experimental data.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u6f5c\u786c\u5316\u7cfb\u6570q\uff1a<\/strong>\u5404\u5411\u540c\u6027\u8fd1\u4f3c\u4e0b\u53ef\u53d6q=1.0\u3002\u82e5\u8981\u66f4\u7cbe\u786e\u63cf\u8ff0\uff0c\u9700\u8981\u901a\u8fc7\u53cc\u6ed1\u79fb\u6216\u591a\u6ed1\u79fb\u5b9e\u9a8c\u6807\u5b9aBassani\u6a21\u578b\u4e2d\u76847\u4e2a\u6f5c\u786c\u5316\u7cfb\u6570\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Latent hardening coefficient q:<\/strong> In isotropic approximation, q=1.0 can be used. For more accurate description, the 7 latent hardening coefficients in the Bassani model need to be calibrated through double-slip or multi-slip experiments.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u6676\u4f53\u53d6\u5411\uff1a<\/strong>\u5355\u6676\u6a21\u62df\u4e2d\u53d6\u5411\u53c2\u6570\u81f3\u5173\u91cd\u8981\u3002\u5efa\u8bae\u901a\u8fc7EBSD\u5b9e\u9a8c\u6d4b\u91cf\u5b9e\u9645\u8bd5\u6837\u7684Euler\u89d2\uff0c\u800c\u975e\u4ec5\u4f9d\u8d56\u540d\u4e49\u53d6\u5411\u3002\u5373\u4f7f\u5fae\u5c0f\u7684\u53d6\u5411\u504f\u5dee\u4e5f\u53ef\u80fd\u5bfc\u81f4\u5b8f\u89c2\u5e94\u529b-\u5e94\u53d8\u54cd\u5e94\u7684\u663e\u8457\u5dee\u5f02\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Crystal orientation:<\/strong> Orientation parameters are crucial in single-crystal simulations. It is recommended to measure the Euler angles of actual specimens via EBSD, rather than relying solely on nominal orientations. Even small orientation deviations can cause significant differences in macroscopic stress-strain response.<\/span>\n    <\/li>\n  <\/ul>\n\n  <h3 id=\"ch8-3\">\n    <span data-lang=\"zh\">8.3 \u6536\u655b\u6027\u8c03\u8bd5 \/ Convergence Debugging<\/span>\n    <span data-lang=\"en\" style=\"display:none\">8.3 Convergence Debugging<\/span>\n  <\/h3>\n  <ul>\n    <li>\n      <span data-lang=\"zh\"><strong>\u65f6\u95f4\u6b65\u8fc7\u5927\uff1a<\/strong>\u82e5\u51fa\u73b0&#8221;Too many attempts made for this increment&#8221;\u9519\u8bef\uff0c\u6700\u5e38\u89c1\u7684\u539f\u56e0\u662f\u65f6\u95f4\u6b65\u8fc7\u5927\u5bfc\u81f4\u5c40\u90e8Newton\u8fed\u4ee3\u4e0d\u6536\u655b\u3002\u53ef\u901a\u8fc7\u51cf\u5c0f\u6700\u5927\u65f6\u95f4\u6b65\u6216\u542f\u7528\u81ea\u52a8\u6b65\u957f\u63a7\u5236\u89e3\u51b3\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Time step too large:<\/strong> If &#8220;Too many attempts made for this increment&#8221; appears, the most common cause is excessive time step size causing local Newton iteration divergence. Reduce maximum time step or enable automatic step control.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u5947\u5f02Jacobian\uff1a<\/strong>\u5f53\u591a\u4e2a\u6ed1\u79fb\u7cfb\u540c\u65f6\u6fc0\u6d3b\u4e14\u786c\u5316\u77e9\u9635\u63a5\u8fd1\u5947\u5f02\u65f6\uff0cLU\u5206\u89e3\u53ef\u80fd\u5931\u8d25\u3002\u53ef\u901a\u8fc7\u589e\u52a0\u6536\u655b\u5bb9\u5deeTOL\u6216\u5728\u6d41\u52a8\u5f8b\u4e2d\u5f15\u5165\u6b63\u5219\u5316\u6765\u6539\u5584\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Singular Jacobian:<\/strong> When multiple slip systems activate simultaneously and the hardening matrix is nearly singular, LU decomposition may fail. Increase convergence tolerance TOL or introduce regularization in the flow rule.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u8d1f\u526a\u5207\u7387\uff1a<\/strong>\u82e5\u67d0\u4e9b\u6ed1\u79fb\u7cfb\u7684\u5206\u89e3\u526a\u5e94\u529b\u4e3a\u8d1f\u4e14\u7edd\u5bf9\u503c\u5f88\u5927\uff0c\u6d41\u52a8\u5f8b\u53ef\u80fd\u4ea7\u751f\u4e0d\u5408\u7406\u7684\u8d1f\u526a\u5207\u7387\u3002\u7a0b\u5e8f\u4e2d\u5e94\u6709\u7b26\u53f7\u5904\u7406\u903b\u8f91\uff0c\u786e\u4fdd\u7269\u7406\u5408\u7406\u6027\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Negative shear rate:<\/strong> If resolved shear stress on some slip systems is negative with large magnitude, the flow rule may produce unreasonable negative shear rates. The program should have sign-handling logic to ensure physical consistency.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u68c0\u67e5STATEV\u521d\u59cb\u5316\uff1a<\/strong>STATEV\u7684\u521d\u59cb\u503c\u5fc5\u987b\u6b63\u786e\u8bbe\u7f6e\uff0c\u7279\u522b\u662f\u5f39\u6027\u53d8\u5f62\u68af\u5ea6Fe\u5e94\u521d\u59cb\u5316\u4e3a\u5355\u4f4d\u77e9\u9635\uff0cFp\u4e5f\u521d\u59cb\u5316\u4e3a\u5355\u4f4d\u77e9\u9635\uff0cg\u521d\u59cb\u5316\u4e3atau0\u3002\u4efb\u4f55\u9519\u8bef\u7684\u521d\u503c\u90fd\u53ef\u80fd\u5bfc\u81f4\u7b2c\u4e00\u6b65\u5373\u53d1\u6563\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Check STATEV initialization:<\/strong> Initial STATEV values must be set correctly, especially Fe and Fp initialized to identity matrices, and g initialized to tau0. Any incorrect initial value may cause divergence in the first step.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>PNEWDT\u7684\u4f7f\u7528\uff1a<\/strong>\u5f53\u5c40\u90e8\u8fed\u4ee3\u8d85\u8fc7ITMAX\u4ecd\u672a\u6536\u655b\u65f6\uff0c\u5e94\u5728UMAT\u4e2d\u8bbe\u7f6ePNEWDT < 1\uff08\u59820.5\uff09\uff0c\u544a\u77e5Abaqus\u51cf\u5c0f\u5f53\u524d\u589e\u91cf\u6b65\u5e76\u91cd\u8bd5\u3002\u8fd9\u662f\u4e00\u79cd\u6709\u6548\u7684\u81ea\u9002\u5e94\u6b65\u957f\u7b56\u7565\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Use of PNEWDT:<\/strong> When local iteration exceeds ITMAX without convergence, set PNEWDT < 1 (e.g., 0.5) in UMAT to instruct Abaqus to reduce the current increment and retry. This is an effective adaptive step strategy.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u8f93\u51fa\u8c03\u8bd5\u4fe1\u606f\uff1a<\/strong>\u8bbe\u7f6eOUTPUT=2\u53ef\u542f\u7528\u8be6\u7ec6\u8c03\u8bd5\u8f93\u51fa\uff0c\u5305\u62ec\u6bcf\u6b21Newton\u8fed\u4ee3\u7684\u6b8b\u5dee\u8303\u6570\u3001\u526a\u5207\u589e\u91cf\u548c\u5e94\u529b\u72b6\u6001\u3002\u8fd9\u5bf9\u4e8e\u5b9a\u4f4d\u53d1\u6563\u539f\u56e0\u975e\u5e38\u6709\u5e2e\u52a9\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Output debug information:<\/strong> Setting OUTPUT=2 enables verbose debug output, including residual norm, shear increments, and stress state for each Newton iteration. This is very helpful for locating divergence causes.<\/span>\n    <\/li>\n  <\/ul>\n\n  <h3 id=\"ch8-4\">\n    <span data-lang=\"zh\">8.4 \u540e\u5904\u7406\u4e0e\u7ed3\u679c\u9a8c\u8bc1 \/ Post-Processing and Verification<\/span>\n    <span data-lang=\"en\" style=\"display:none\">8.4 Post-Processing and Verification<\/span>\n  <\/h3>\n  <ul>\n    <li>\n      <span data-lang=\"zh\"><strong>\u5e94\u529b-\u5e94\u53d8\u66f2\u7ebf\uff1a<\/strong>\u5c06Abaqus\u8f93\u51fa\u7684\u5b8f\u89c2\u5e94\u529b-\u5e94\u53d8\u66f2\u7ebf\u4e0e\u5355\u6676\u5b9e\u9a8c\u6570\u636e\u8fdb\u884c\u5bf9\u6bd4\uff0c\u9a8c\u8bc1\u6a21\u578b\u9884\u6d4b\u7684\u57fa\u672c\u529b\u5b66\u884c\u4e3a\u662f\u5426\u6b63\u786e\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Stress-strain curve:<\/strong> Compare the macroscopic stress-strain curve from Abaqus with single-crystal experimental data to verify whether the basic mechanical behavior predicted by the model is correct.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u6ed1\u79fb\u7cfb\u6fc0\u6d3b\u5206\u6790\uff1a<\/strong>\u901a\u8fc7\u540e\u5904\u7406STATEV\u4e2d\u5b58\u50a8\u7684\u7d2f\u79ef\u526a\u5207\u5e94\u53d8\uff0c\u53ef\u4ee5\u8bc6\u522b\u54ea\u4e9b\u6ed1\u79fb\u7cfb\u88ab\u6fc0\u6d3b\u53ca\u5176\u76f8\u5bf9\u6d3b\u52a8\u5f3a\u5ea6\u3002\u8fd9\u5bf9\u4e8e\u7406\u89e3\u53d8\u5f62\u673a\u5236\u548c\u7ec7\u6784\u6f14\u5316\u81f3\u5173\u91cd\u8981\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Slip system activation analysis:<\/strong> By post-processing accumulated shear strains stored in STATEV, one can identify which slip systems are activated and their relative activity levels. This is crucial for understanding deformation mechanisms and texture evolution.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u7ec7\u6784\u9a8c\u8bc1\uff1a<\/strong>\u5982\u679c\u542f\u7528\u4e86\u53d6\u5411\u66f4\u65b0\uff08PROPS(70)=1\uff09\uff0c\u5e94\u63d0\u53d6\u53d8\u5f62\u540e\u7684\u6676\u4f53\u53d6\u5411\u5e76\u4e0e\u5b9e\u9a8cEBSD\u7ec7\u6784\u6570\u636e\u8fdb\u884c\u5bf9\u6bd4\uff0c\u8bc4\u4f30\u7ec7\u6784\u6f14\u5316\u9884\u6d4b\u7684\u51c6\u786e\u6027\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Texture verification:<\/strong> If orientation update is enabled (PROPS(70)=1), extract deformed crystal orientations and compare with experimental EBSD texture data to assess the accuracy of texture evolution prediction.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u7f51\u683c\u6536\u655b\u6027\uff1a<\/strong>\u8fdb\u884c\u7f51\u683c\u654f\u611f\u6027\u5206\u6790\uff0c\u786e\u4fdd\u7ed3\u679c\u4e0d\u4f9d\u8d56\u4e8e\u79bb\u6563\u5316\u7a0b\u5ea6\u3002\u6676\u4f53\u5851\u6027\u6a21\u62df\u4e2d\uff0c\u5e94\u53d8\u5c40\u90e8\u5316\u5e26\uff08shear band\uff09\u7684\u5bbd\u5ea6\u901a\u5e38\u53d7\u7f51\u683c\u5c3a\u5bf8\u5f71\u54cd\uff0c\u5e94\u4f7f\u7528\u8db3\u591f\u7ec6\u5bc6\u7684\u7f51\u683c\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Mesh convergence:<\/strong> Perform mesh sensitivity analysis to ensure results do not depend on discretization level. In crystal plasticity simulations, the width of strain localization bands (shear bands) is usually affected by mesh size; sufficiently fine meshes should be used.<\/span>\n    <\/li>\n    <li>\n      <span data-lang=\"zh\"><strong>\u80fd\u91cf\u5b88\u6052\u68c0\u67e5\uff1a<\/strong>\u9a8c\u8bc1\u5851\u6027\u8017\u6563\u529f\u4e0e\u5916\u90e8\u529f\u7684\u5173\u7cfb\uff0c\u786e\u4fdd\u6570\u503c\u5b9e\u73b0\u6ca1\u6709\u80fd\u91cf\u9519\u8bef\u3002Abaqus\u8f93\u51fa\u7684\u80fd\u91cf\u5386\u53f2\uff08ALLIE\u3001ALLKE\u3001ALLPD\u7b49\uff09\u53ef\u7528\u4e8e\u6b64\u76ee\u7684\u3002<\/span>\n      <span data-lang=\"en\" style=\"display:none\"><strong>Energy conservation check:<\/strong> Verify the relationship between plastic dissipation work and external work to ensure the numerical implementation has no energy errors. Abaqus energy history outputs (ALLIE, ALLKE, ALLPD, etc.) can be used for this purpose.<\/span>\n    <\/li>\n  <\/ul>\n\n  <h3 id=\"ch8-5\">\n    <span data-lang=\"zh\">8.5 \u5e38\u89c1\u9519\u8bef\u4e0e\u89e3\u51b3\u65b9\u6848 \/ Common Errors and Solutions<\/span>\n    <span data-lang=\"en\" style=\"display:none\">8.5 Common Errors and Solutions<\/span>\n  <\/h3>\n  <table>\n    <thead>\n      <tr>\n        <th>\n          <span data-lang=\"zh\">\u9519\u8bef\u73b0\u8c61 \/ Error<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Error Symptom<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u53ef\u80fd\u539f\u56e0 \/ Cause<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Possible Cause<\/span>\n        <\/th>\n        <th>\n          <span data-lang=\"zh\">\u89e3\u51b3\u65b9\u6848 \/ Solution<\/span>\n          <span data-lang=\"en\" style=\"display:none\">Solution<\/span>\n        <\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr>\n        <td><span data-lang=\"zh\">\u589e\u91cf\u6b65\u591a\u6b21\u5c1d\u8bd5\u540e\u4ecd\u5931\u8d25<\/span><span data-lang=\"en\" style=\"display:none\">Increment fails after many attempts<\/span><\/td>\n        <td><span data-lang=\"zh\">\u65f6\u95f4\u6b65\u8fc7\u5927\u6216\u6750\u6599\u53c2\u6570\u4e0d\u5408\u7406<\/span><span data-lang=\"en\" style=\"display:none\">Time step too large or material parameters unreasonable<\/span><\/td>\n        <td><span data-lang=\"zh\">\u51cf\u5c0f\u6700\u5927\u65f6\u95f4\u6b65\uff0c\u68c0\u67e5\u6750\u6599\u53c2\u6570\u6570\u91cf\u7ea7<\/span><span data-lang=\"en\" style=\"display:none\">Reduce max time step, check parameter magnitudes<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td><span data-lang=\"zh\">\u5e94\u529b\u8f93\u51fa\u4e3aNaN\u6216Inf<\/span><span data-lang=\"en\" style=\"display:none\">Stress output is NaN or Inf<\/span><\/td>\n        <td><span data-lang=\"zh\">\u9664\u96f6\u9519\u8bef\u6216\u6307\u6570\u6ea2\u51fa<\/span><span data-lang=\"en\" style=\"display:none\">Division by zero or exponent overflow<\/span><\/td>\n        <td><span data-lang=\"zh\">\u68c0\u67e5tau\u548cg\u662f\u5426\u4e3a\u6b63\uff0c\u9650\u5236\u5e42\u5f8b\u6307\u6570\u8303\u56f4<\/span><span data-lang=\"en\" style=\"display:none\">Ensure tau and g are positive, limit power-law exponent range<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td><span data-lang=\"zh\">STATEV\u503c\u5f02\u5e38\u589e\u5927<\/span><span data-lang=\"en\" style=\"display:none\">STATEV values abnormally large<\/span><\/td>\n        <td><span data-lang=\"zh\">\u72b6\u6001\u53d8\u91cf\u672a\u6b63\u786e\u521d\u59cb\u5316<\/span><span data-lang=\"en\" style=\"display:none\">State variables not properly initialized<\/span><\/td>\n        <td><span data-lang=\"zh\">\u5728*Initial Conditions\u4e2d\u6b63\u786e\u8bbe\u7f6eSTATEV\u521d\u503c<\/span><span data-lang=\"en\" style=\"display:none\">Set correct STATEV initial values in *Initial Conditions<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td><span data-lang=\"zh\">\u7ed3\u679c\u4e0e\u5b9e\u9a8c\u5dee\u5f02\u5de8\u5927<\/span><span data-lang=\"en\" style=\"display:none\">Results differ greatly from experiment<\/span><\/td>\n        <td><span data-lang=\"zh\">\u53d6\u5411\u53c2\u6570\u9519\u8bef\u6216\u5355\u4f4d\u4e0d\u4e00\u81f4<\/span><span data-lang=\"en\" style=\"display:none\">Orientation parameters wrong or unit inconsistency<\/span><\/td>\n        <td><span data-lang=\"zh\">\u6838\u5bf9Euler\u89d2\u5355\u4f4d\u548c\u6ed1\u79fb\u7cfb\u5b9a\u4e49\uff0c\u7edf\u4e00\u5e94\u529b\u5355\u4f4d\u4e3aMPa<\/span><span data-lang=\"en\" style=\"display:none\">Check Euler angle units and slip system definitions, unify stress units to MPa<\/span><\/td>\n      <\/tr>\n      <tr>\n        <td><span data-lang=\"zh\">\u8fd0\u884c\u901f\u5ea6\u6781\u6162<\/span><span data-lang=\"en\" style=\"display:none\">Extremely slow execution<\/span><\/td>\n        <td><span data-lang=\"zh\">\u65f6\u95f4\u6b65\u8fc7\u5c0f\u6216Newton\u8fed\u4ee3\u6b21\u6570\u8fc7\u591a<\/span><span data-lang=\"en\" style=\"display:none\">Time step too small or too many Newton iterations<\/span><\/td>\n        <td><span data-lang=\"zh\">\u589e\u5927\u6536\u655b\u5bb9\u5deeTOL\uff0c\u964d\u4f4eITMAX\uff0c\u6216\u68c0\u67e5\u786c\u5316\u53c2\u6570\u662f\u5426\u5bfc\u81f4\u75c5\u6001\u95ee\u9898<\/span><span data-lang=\"en\" style=\"display:none\">Increase TOL, reduce ITMAX, or check if hardening parameters cause ill-conditioning<\/span><\/td>\n      <\/tr>\n    <\/tbody>\n  <\/table>\n<\/div>\n\n<!-- ==================== REFERENCES ==================== -->\n<div class=\"chapter\" id=\"refs\">\n  <h2>\n    <span class=\"chap-num\">Ref<\/span>\n    <span data-lang=\"zh\">\u53c2\u8003\u6587\u732e \/ References<\/span>\n    <span data-lang=\"en\" style=\"display:none\">References<\/span>\n  <\/h2>\n  <p data-lang=\"zh\">\n    \u4ee5\u4e0b\u5217\u51fa\u9ec4\u6c38\u521aUMAT\u5b50\u7a0b\u5e8f\u76f8\u5173\u7684\u6838\u5fc3\u7406\u8bba\u6587\u732e\uff0c\u6309\u53d1\u8868\u65f6\u95f4\u6392\u5e8f\uff1a\n  <\/p>\n  <p data-lang=\"en\" style=\"display:none\">\n    The following lists core theoretical references related to Huang&#8217;s UMAT subroutine, in chronological order:\n  <\/p>\n  <div class=\"ref-item\">\n    [1] Peirce D, Asaro RJ, Needleman A. Material rate dependence and localized deformation in crystalline solids.\n    <em>Acta Metallurgica<\/em>, 1983, 31(12): 1951-1976.\n  <\/div>\n  <div class=\"ref-item\">\n    [2] Peirce D, Shih CF, Needleman A. A tangent modulus method for rate dependent solids.\n    <em>Computers &amp; Structures<\/em>, 1984, 18(5): 875-887.\n  <\/div>\n  <div class=\"ref-item\">\n    [3] Asaro RJ, Needleman A. Texture development and strain hardening in rate dependent polycrystals.\n    <em>Acta Metallurgica<\/em>, 1985, 33(6): 923-953.\n  <\/div>\n  <div class=\"ref-item\">\n    [4] Harren SV, Deve HE, Asaro RJ. Shear band formation in plane strain compression.\n    <em>Acta Metallurgica<\/em>, 1988, 36(9): 2435-2480.\n  <\/div>\n  <div class=\"ref-item\">\n    [5] Bassani JL, Wu TY. Latent hardening in single crystals. II. Analytical characterization and predictions.\n    <em>Proceedings of the Royal Society A<\/em>, 1991, 435(1893): 21-41.\n  <\/div>\n  <div class=\"ref-item\">\n    [6] Huang Y. A user-material subroutine incorporating single crystal plasticity in the ABAQUS finite element program.\n    <em>Harvard University Report MECH-178<\/em>, 1991.\n  <\/div>\n  <div class=\"ref-item\">\n    [7] Kalidindi SR, Bronkhorst CA, Anand L. Crystallographic texture evolution in bulk deformation processing of FCC metals.\n    <em>Journal of the Mechanics and Physics of Solids<\/em>, 1992, 40(3): 537-569.\n  <\/div>\n  <div class=\"ref-item\">\n    [8] Kocks UF, Tome CN, Wenk HR. <em>Texture and Anisotropy: Preferred Orientations in Polycrystals and Their Effect on Materials Properties<\/em>.\n    Cambridge University Press, 1998.\n  <\/div>\n  <div class=\"ref-item\">\n    [9] Roters F, Eisenlohr P, Hantcherli L, et al. Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling: Theory, experiments, applications.\n    <em>Acta Materialia<\/em>, 2010, 58(4): 1152-1211.\n  <\/div>\n  <div class=\"ref-item\">\n    [10] Dunne FPE, Petrinic N. <em>Introduction to Computational Plasticity<\/em>. Oxford University Press, 2005.\n  <\/div>\n<\/div>\n\n<!-- ==================== APPENDIX ==================== -->\n<div class=\"chapter\" id=\"appendix\">\n  <h2>\n    <span class=\"chap-num\">App<\/span>\n    <span data-lang=\"zh\">\u7b26\u53f7\u8868 \/ Notation<\/span>\n    <span data-lang=\"en\" style=\"display:none\">Notation<\/span>\n  <\/h2>\n  <p data-lang=\"zh\">\u672c\u6587\u6863\u4e2d\u4f7f\u7528\u7684\u5173\u952e\u6570\u5b66\u7b26\u53f7\u53ca\u5176\u542b\u4e49\u5982\u4e0b\u8868\u6240\u793a\uff1a<\/p>\n  <p data-lang=\"en\" style=\"display:none\">The key mathematical symbols used in this document and their meanings are listed in the following table:<\/p>\n  <table>\n    <thead>\n      <tr>\n        <th><span data-lang=\"zh\">\u7b26\u53f7 \/ Symbol<\/span><span data-lang=\"en\" style=\"display:none\">Symbol<\/span><\/th>\n        <th><span data-lang=\"zh\">\u542b\u4e49 \/ Meaning<\/span><span data-lang=\"en\" style=\"display:none\">Meaning<\/span><\/th>\n      <\/tr>\n    <\/thead>\n    <tbody>\n      <tr><td><b>F<\/b><\/td><td><span data-lang=\"zh\">\u603b\u53d8\u5f62\u68af\u5ea6 \/ Total deformation gradient<\/span><span data-lang=\"en\" style=\"display:none\">Total deformation gradient<\/span><\/td><\/tr>\n      <tr><td><b>F<sup>e<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u5f39\u6027\u53d8\u5f62\u68af\u5ea6 \/ Elastic deformation gradient<\/span><span data-lang=\"en\" style=\"display:none\">Elastic deformation gradient<\/span><\/td><\/tr>\n      <tr><td><b>F<sup>p<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u5851\u6027\u53d8\u5f62\u68af\u5ea6 \/ Plastic deformation gradient<\/span><span data-lang=\"en\" style=\"display:none\">Plastic deformation gradient<\/span><\/td><\/tr>\n      <tr><td><b>&sigma;<\/b><\/td><td><span data-lang=\"zh\">Cauchy\u5e94\u529b\u5f20\u91cf \/ Cauchy stress tensor<\/span><span data-lang=\"en\" style=\"display:none\">Cauchy stress tensor<\/span><\/td><\/tr>\n      <tr><td><b>&tau;<\/b><\/td><td><span data-lang=\"zh\">Kirchhoff\u5e94\u529b\u5f20\u91cf \/ Kirchhoff stress tensor<\/span><span data-lang=\"en\" style=\"display:none\">Kirchhoff stress tensor<\/span><\/td><\/tr>\n      <tr><td><b>S<sup>e<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u7b2c\u4e8c\u7c7bPiola-Kirchhoff\u5e94\u529b \/ 2nd P-K stress<\/span><span data-lang=\"en\" style=\"display:none\">2nd Piola-Kirchhoff stress<\/span><\/td><\/tr>\n      <tr><td><b>E<sup>e<\/sup><\/b><\/td><td><span data-lang=\"zh\">Green\u5f39\u6027\u5e94\u53d8 \/ Green elastic strain<\/span><span data-lang=\"en\" style=\"display:none\">Green elastic strain<\/span><\/td><\/tr>\n      <tr><td><b>&gamma;&#775;<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u7b2c&alpha;\u6ed1\u79fb\u7cfb\u526a\u5207\u5e94\u53d8\u7387 \/ Shear strain rate on slip system &alpha;<\/span><span data-lang=\"en\" style=\"display:none\">Shear strain rate on slip system &alpha;<\/span><\/td><\/tr>\n      <tr><td><b>&tau;<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u5206\u89e3\u526a\u5e94\u529b \/ Resolved shear stress<\/span><span data-lang=\"en\" style=\"display:none\">Resolved shear stress<\/span><\/td><\/tr>\n      <tr><td><b>g<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u4e34\u754c\u5206\u5207\u5e94\u529b\uff08\u6ed1\u79fb\u5f3a\u5ea6\uff09\/ CRSS (slip strength)<\/span><span data-lang=\"en\" style=\"display:none\">Critical resolved shear stress (slip strength)<\/span><\/td><\/tr>\n      <tr><td><b>P<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">Schmid\u5f20\u91cf \/ Schmid tensor<\/span><span data-lang=\"en\" style=\"display:none\">Schmid tensor<\/span><\/td><\/tr>\n      <tr><td><b>s<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u6ed1\u79fb\u65b9\u5411 \/ Slip direction<\/span><span data-lang=\"en\" style=\"display:none\">Slip direction<\/span><\/td><\/tr>\n      <tr><td><b>m<sup>(&alpha;)<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u6ed1\u79fb\u9762\u6cd5\u5411 \/ Slip plane normal<\/span><span data-lang=\"en\" style=\"display:none\">Slip plane normal<\/span><\/td><\/tr>\n      <tr><td><b>h<sub>&alpha;&beta;<\/sub><\/b><\/td><td><span data-lang=\"zh\">\u786c\u5316\u77e9\u9635 \/ Hardening matrix<\/span><span data-lang=\"en\" style=\"display:none\">Hardening matrix<\/span><\/td><\/tr>\n      <tr><td><b>C<sup>e<\/sup><\/b><\/td><td><span data-lang=\"zh\">\u5f39\u6027\u521a\u5ea6\u5f20\u91cf \/ Elastic stiffness tensor<\/span><span data-lang=\"en\" style=\"display:none\">Elastic stiffness tensor<\/span><\/td><\/tr>\n      <tr><td><b>DDSDDE<\/b><\/td><td><span data-lang=\"zh\">\u4e00\u81f4\u5207\u7ebf\u521a\u5ea6\u77e9\u9635 \/ Consistent tangent modulus<\/span><span data-lang=\"en\" style=\"display:none\">Consistent tangent modulus<\/span><\/td><\/tr>\n      <tr><td><b>STATEV<\/b><\/td><td><span data-lang=\"zh\">\u72b6\u6001\u53d8\u91cf\u6570\u7ec4 \/ State variable array<\/span><span data-lang=\"en\" style=\"display:none\">State variable array<\/span><\/td><\/tr>\n      <tr><td><b>PROPS<\/b><\/td><td><span data-lang=\"zh\">\u6750\u6599\u53c2\u6570\u6570\u7ec4 \/ Material property array<\/span><span data-lang=\"en\" style=\"display:none\">Material property array<\/span><\/td><\/tr>\n    <\/tbody>\n  <\/table>\n<\/div>\n\n<\/div>\n\n<footer style=\"text-align:center; padding: 30px; color: var(--muted); font-size: 0.9rem; border-top: 1px solid var(--rule); margin-top: 20px;\">\n  <p data-lang=\"zh\">\u672c\u6587\u6863\u4ec5\u4f9b\u5b66\u672f\u7814\u7a76\u4e0e\u6559\u5b66\u53c2\u8003\u3002<\/p>\n  <p data-lang=\"en\" style=\"display:none\">This document was for academic research and teaching reference only.<\/p>\n  <p data-lang=\"zh\">\u6700\u540e\u66f4\u65b0\uff1a2026\u5e746\u670829\u65e5<\/p>\n  <p data-lang=\"en\" style=\"display:none\">Last updated: June 29, 2026<\/p>\n<\/footer>\n\n<script src=\".\/_shared\/js\/mermaid.min.js\"><\/script>\n<script>\n  if (typeof mermaid !== 'undefined') {\n    mermaid.initialize({ startOnLoad: true, theme: 'neutral', securityLevel: 'loose' });\n  }\n<\/script>\n<\/body>\n<\/html>\n\n","protected":false},"excerpt":{"rendered":"<p>HuangUMAT\u5355\u6676\u5851\u6027\u5b50\u7a0b\u5e8f\u8be6\u89e3 Huang&#8217;s UMAT Crystal Plasticit [&hellip;]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[7,8],"tags":[],"class_list":["post-261","post","type-post","status-publish","format-standard","hentry","category-abaqus","category-abaqus-subroutine"],"blocksy_meta":[],"_links":{"self":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/261","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=261"}],"version-history":[{"count":2,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/261\/revisions"}],"predecessor-version":[{"id":266,"href":"https:\/\/numsimlab.com\/index.php?rest_route=\/wp\/v2\/posts\/261\/revisions\/266"}],"wp:attachment":[{"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=261"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=261"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/numsimlab.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=261"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}