CICC导体缩径过程的有限元仿真.pdfVIP

维普资讯 · 44· 焊管 ．第30卷第4期．2007年7月 ●应 用 与 开 发 CICC导体缩径过程的有限元仿真 刘志宏，郁 杰，吴杰峰 (中国科学院等离子体物理研究所，合肥230031) 摘 要：利用显式动力学有限元分析软件ANSYS／LS—DYNA仿真 了CICC(CableinConduit Conductor)导体的缩径过程，研究了不同轧制速度和摩擦系数条件下轧制力的变化规律。模 拟了导体缩径过程中卷入和稳定轧制时导体 内应力、应变的变化过程，分析 了缩径过程中钢管 的弹复行为。另外对轧制仿真过程中网格密度对计算精度的影响也进行了计算和比较。其结 果为 CICC导体的研制提供 了依据。 关键词：CICC；缩径；有限元法 中图分类号：TG335．71 文献标志码：A 文章编号：1001—3938(2007)04—0044—03 1 CICC导体简介 CICC(CableinConduitConductor)导体是将 由超导丝和铜丝绞制成的超导电缆用不锈钢外壳 包裹，然后将整个导体压制成需要的形状和内部 空隙率而得到的。它常用在超导核聚变实验装置 中做成超导线圈以产生高能量的磁场。在2006 年启动 的由中国、美国、俄罗斯等 7国参与 的 图2 双辊缩径模型 rrER(国际热核试验反应堆)项 目中，中国承担了 很大部分CICC导体的研制任务。 2 显式动力学有限元理论 CICC导体的缩径过程具有复杂的几何非线 性、物理非线性、边界非线性和接触非线性等特征， 在动力学有限元分析中，系统的求解方程可 描述为 并具有缩径量小、单根长度长、尺寸精度要求高等 特点，所以传统的解析或半解析方法无法合理求解 ga(t)+Ca(t)+Ka(t)=Q(t)， (1) 式中：d(t)，d(t)分别是系统的节点加速度向量 其三维变形问题。本研究采用显式动力学有限元 和节点速度向量； ，C，K和 Q(t)分别是系统的 分析软件 ANSYS／LS—DYNA对其缩径成形工艺 质量矩阵、阻尼矩阵、刚度矩阵和节点载荷向量， 进行全流程数值仿真，模拟了不同轧制速度、摩擦 分别由各 自的单元矩阵和向量集成。 系数和网格密度条件下的轧制过程。ITER项 目中 与静力分析相比较，在动力分析中，由于惯性 所使用的一种TF(纵向场)线圈的CICC导体截面 力和阻尼力出现在平衡方程中，因此引入了质量 结构见图1，其双辊缩径模型见图2。 矩阵和阻尼矩阵，最后得到的求解方程不是代数 方程组，而是常微分方程组。在显式动力学有限 元分析软件 ANSYS／LS—DYNA 中，用直接积分 法中的中心差分格式对运动方程进行积分。在中 心差分法中，加速度和速度可以用位移表示为 1 d= (口f_△一2a+口… )，