Abaqus学习资料摘取.docx

材料的塑性性能可以用它的屈服点和屈服后硬化特性来描述。从弹性到塑性的转变发生在材料应力应变曲线上的某个确定点，即通常所说的弹性极限或屈服点(见图8-2)。屈服点的应力叫做屈服应力。大部分金属的屈服应力为材料弹性模量的0.05%-0.1%。金属在到达屈服点之前的变形只产生弹性应变，在卸载后可以完全恢复。然而，一旦应力超过屈服应力就会开始产生永久(塑性)变形。与这种永久变形相关的应变称为塑性应变。在屈服后的区间里，弹性和塑性应变共同组成了金属的变形。金属的刚度在材料屈服后会显著下降(见图8-2)。已屈服的延性金属在卸载后将恢复它的初始刚度(见图8-2)。通常，塑性变形会提高材料重新加载时的屈服应力：这一特性称为硬化。金属塑性的另一个重要特点就是非弹性变形与材料几乎不可压缩的特性相关，这一效应给用于模拟弹塑性的单元类型的选用带来很大的限制。承受拉力的金属在塑性变形时，可能会在材料失效时经历局部的高度伸长与变细，称为颈缩(见图8-2) 。金属的工程应力(力除以变形前的截面面积)称为名义应力，与之对应的为名义应变(伸长量除以原长)。金属在发生颈缩时的名义应力远低于材料的极限强度。这种材料特性是由试件几何形状、实验本身特点以及应力应变测量方法引起的。例如，由相同材料的压缩实验所得到的应力应变曲线就不会出现颈缩区域，因为试件在受压变形时不会变细。因此，描述金属塑性的数学模型应该能够考虑拉伸和压缩的不同特性，并与结构几何形状和加载特性无关。为了实现这一目的，应当把已十分熟悉的名义应力的定义, F A0 , 和名义应变, l l0（这里用下标0 代表材料未变形状态下的值），替换为能考虑在有限变形中面积的改变的应力与应变的新度量方法。有限变形中的应力应变度量只有考虑在极限l dl 0 的情况下，拉伸和压缩应变是相同的，例如：其中是l 当前长度， 0 l 是原始长度， 为真实应变或对数应变。与真实应变对应的应力是真实应力，定义为：其中F 是材料承受的力，A 是当前面积。承受有限变形的延性材料的真实应力真实应变曲线对于拉伸和压缩是一致的。在ABAQUS 中定义塑性在ABAQUS 中必须用真实应力和真实应变定义塑性。ABAQUS 需要这些值并相应地在输入文件中解释这些数据。大多数实验数据常常是用名义应力和名义应变的值给出的。这时，必须应用下列公式把塑性材料的数据从名义应力/应变的值转换为真实应力/应变的值。先把名义应力表达为如下公式：表达式两边同时加上1，并取自然对数可以得到真实应力和名义应力之间的关系为：考虑塑性变形的不可压缩性，并假定弹性变形也是不可压缩的，由体积不变得：当前面积和原面积的关系为：将A 的定义代入到真实应力的定义式中，得到：其中ll0也可以写为1+nom ，代入上式便得到真实应力和名义应力、名义应变之间的关系：ABAQUS 中经典金属塑性模型定义了大部分金属的屈服后特性。ABAQUS用连接给定数据点的一系列直线的方法来逼近材料的光滑应力应变曲线。因为可以用任意多个点来逼近实际的材料性质；所以，有可能非常逼真地模拟材料的真实性质。在塑性数据中将材料的真实屈服应力定义为真实塑性应变的函数。第一组数据定义材料的初始屈服应力，因此，该组数据中的塑性应变值应该为零。在塑性性能的材料试验数据中所提供的应变不可能仅仅是材料的塑性应变，而应当是材料的总应变。所以必须将总应变分解为弹性和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值，从总应变中减去弹性应变，就得到了塑性应变(见图8-3)。其关系式为：其中 pl ：真实塑性应变，t：总体真实应变，elel ：真实弹性应变， ：真实应力，E ：杨氏模量。试验数据转换为ABAQUS 输入数据的例题图8-4 中的名义应力应变曲线可以作为一个例子用来示范如何将材料塑性特性的试验数据转换为ABAQUS 适用的输入格式。名义应力应变曲线上的六个点用来确定塑性数据。第一步是用前述公式将名义应力和名义应变转化为真实应力和应变。一旦得到这些值，就可以用上述塑性应变与总应变和弹性应变之间的关系式来确定与每个屈服应力值相关的塑性应变。表8.1 中给出了转换后的数据。当应变很小时，真实值和名义值之间差别很小，而当应变很大时，二者之间就会有明显的差别；因此，如果模拟的应变比较大，就一定要向ABAQUS 提供合适的应力应变数据，这是极为重要的。生成等值线图的算法需要的是节点数据，而ABAQUS/Standard计算的却是积分点的单元变量。ABAQUS/CAE 将积分点的数据外推从而得到节点的数据。外推算法的阶次由单元类型决定： 对二次减缩积分单元，ABAQUS/CAE 采用线性外推算法计算单元变量的节点数值。为了显示VonMises 应力的等值线图，ABAQUS/CAE 在每个单元内将应力分量从积分点外推到节点位置并计算Von Mises 应力。如果由节点周围各个单元计算出来的Mises应力之间的差值落在了所定义的平均阀值（specified averaging threshold）内，则由各个单元的应力不变量计算节点平均Mises 应力。外推过程会导致不变量超出弹性极限。若绘制应力张量的每一个分量的等值线图(变量S11，S22，S33，S12，S23 和S13)，可以发现在固定端截面上的单元应力有明显的变化。这将导致外推的节点应力高于积分点应力，从而由此计算出的von Mises 应力也会偏高。积分点的Meses 应力不会超出单元材料的当前屈服应力，但应力等值线中的外推节点应力可能会超出。另外，独立的应力分量可以超出当前的屈服应力值；因为只要求Mises 应力小于等于当前的屈服应力值。绘制等效塑性应变图材料的等效塑性应变(PEEQ)是用来表示材料非弹性变形的标量。如果该变量大于零，则材料屈服。可以在PEEQ 等值线图中确定连接环已屈服的部分。从主菜单中选择Result Field Output，在弹出的对话框中的输出变量列表中选择PEEQ可以绘出等效应变的等值线图。激活Contour Plot Options 对话框，设置等效塑性应变等值线的最小值为一个很小的数值（例如：1.E-4），那么，ABAQUS/CAE 绘制的模型中的所有深兰色区域将仍保持弹性特性(见图8-11)。从中可以清楚的发现，连接环与母体相连的部分有明显的屈服。等值线图例中给出的最大塑性应变为0.1082。当然，这里面可能包括外推过程的误差，但这种对应变的外推误差要比对应力的小。利用可视化查询工具检查塑性应变最大的单元积分点处的PEEQ 值。可以发现模型中积分点处最大等效塑性应变约为0.067，这并不一定表明外推误差很大，因为在峰值塑性变形附近的应变梯度很大。应变势能对于超弹性材料，ABAQUS 不用杨氏模量和泊松比，而用应变势能(U)来表达应力应变关系。有多种应变势能可利用：多项式模型、Ogden 模型、Arruda-Boyce 模型、van der Waals 模型。也可以使用较简单的多项式模型，包括Mooney-Rivlin 模型、 neo-Hookean 模型、简约的多项式模型和Yeoh 模型。多项式形式的应变势能是最常用的一种形式。可以表达为：其中U 是应变势能； J el 是弹性体积比； I 1 和I 2 是材料的变形度量；N ， C ij 和D i 是与温度有关的材料参数。参数C ij 描述材料的剪切特性，参数D i 引入可压缩性。如果D i 设置为0，ABAQUS 将认为材料为完全不可压缩， 并忽略上述公式的第二项。如果项数N 为1，初始剪切模量0 ，体积模量K 0 的公式如下：如果材料是不可压缩的，则应变能密度为：该表达式即通常所说的Moongy-Rovlin材料模型。如果C 01 为0，则称neo-Hookean 材料模型。图817 用于定义材料超弹性的各种试验的变形模式接触分析逻辑流程图ABAQUS 在每个增量步开始之前检察所有接触相互作用状态，以判断从属节点是脱开还是闭合。在图11.7 中p 表示从属节点上的接触压力，h 表示从属节点对主面的侵入距离。如果一个节点是闭合的，ABAQUS 确定它是在滑动还是粘结。ABAQUS 对每个闭合节点加以约束，而对那些接触状态从闭合到脱离变化的节点撤除约束。然后ABAQUS 再次进行迭代并用计算修正值来改变模型。在检验力或力矩的平衡前，ABAQUS 先检查从属节点上接触状态的变化。若节点在迭代后间隙变为负的或零，则它的状态由脱离变为闭合。若节点在迭代后接触压力变为负的，则它的状态则由闭合变为脱开。如果检测到当前迭代步的接触状态有变化，ABAQUS 将它标识为严重不连续迭代（severe discontinuityiteration），且不进行平衡检验。在第一次迭代结束后，ABAQUS 通过改变接触约束来反映接触状态的改变，然后进行第二次迭代。ABAQUS 重复这个过程，直到接触状态不再变化才结束迭代。接着的迭代为第一次平衡迭代，并且ABAQUS 进行正常的平衡收敛检查。如果收敛检查失败，ABAQUS 将进行另一次迭代。每当一个严重不连续迭代发生时，ABAQUS 将内部平衡迭代计数器重新置零。这个平衡迭代的计数用于确定是否因收敛慢而放弃这个增量步。ABAQUS 重复整个过程直至获得收敛的结果，如图11.7 所示。在信息和状态文件中，每完成一个增量步就会总结显示有多少次严重不连续迭代，和多少次平衡迭代。增量步的总迭代数是这两者之和。通过区分这两类迭代，可以看到ABAQUS 非常适合处理接触计算和很恰当地完成平衡迭代。如果严重不连续迭代数很多，而只有很少的平衡迭代，那么ABAQUS 对确定合适的接触状态就会出现困难。在默认情况下，ABAQUS 会地放弃那些超过12 个严重不连续迭代的增量步，而改用更小的增量步。如果没有严重不连续迭代，接触状态从一个增量步到另一个增量步之间没有改变。可能的刚体运动取决于模型的维数。施加约束的规则如下： 边界条件不应施加在从面上，它会与接触约束发生冲突。 边界条件不应施加在有线性约束方程的自由度上。如果在有约束方程的自由度上施加边界条件，ABAQUS 将在数据文件中给出错误消息。 为了容易检查约束处的反力，边界条件不应施加在载荷作用处。*OUTPUT：定义输出请求该选项用来写入接触、单元、能量、节点或诊断输出到输出数据库中。是一个统领的关键词，下面跟随的子关键字有*CONTACT OUTPUT、*ELEMENT OUTPUT、*NODE OUTPUT等。Standard中的输出必需的相互排斥的参数：DIAGNOSTICS：设置DIAGNOSTICS=YES(默认)表明详细的诊断信息将被写入输出数据库。NO则抑制。FIELD：包含该参数表明与*OUTPUT联合使用的输出请求将以场的形式写入输出数据库。HISTORY：包含该参数表明与*OUTPUT联合使用的输出请求将以历程的形式写入输出数据库。可选参数：FREQUENCY：设置该参数等于分期的输出频率。输出将总是在每个分析步的最后被写入输出数据库。设置FREQUENCY=0则抑制输出。如果该参数和NUMBER INTERVAL、TIME INTERVAL和TIME POINT被省略，输出将在每个增量步写入，除了*DYNAMIC和*MODAL DYNAMIC，输出会每隔10个增量步写入。FREQUENCY、NUMBER INTERVAL、TIME INTERVAL和TIME POINT是相互排斥的。MODE LIST：包含该参数表明期望输出的一系列模态会在数据行列表显示。该参数仅在*FREQUENCY、*COMPLEX FREQUENCY或*BUCKLE有用，并且包含FIELD。NUMBER INTERVAL：设置该参数等于输出数据写入时间间隔的数量。如果该参数和FREQUENCY、TIME INTERVAL和TIME POINT被省略，输出将在每个增量步写入，除了*DYNAMIC和*MODAL DYNAMIC，输出会每隔10个增量步写入。FREQUENCY、NUMBER INTERVAL、TIME INTERVAL和TIME POINT是相互排斥的。TIME MARKS：设置TIME MARKS=YES(默认)将在由NUMBER INTERVAL、TIME INTERVAL或TIME POINT参数控制的准确时间写入。=NO则将在由NUMBER INTERVAL、TIME INTERVAL或TIME POINT参数控制的准确时间后的增量步结束时写入。TIME POINTS：设置该参数等于*TIME POINTS的名字来定义输出写入的时间点。下面的参数是可选的，而且仅当包含FIELD和HISTORY参数时有效：OP：设置OP=NEW(默认)表明前一个分析步定义的所有输出数据请求都将被移除。并定义新的输出请求。设置OP=ADD表明所定义的输出请求被添加到签名分析步的输出请求中。 设置OP=REPLACE则当前输出请求会代替相同类型(比如field)的输出请求和频率。如果没有匹配的请求，输出请求将被看作是ADD。TIME INTERVAL：设置该参数等于要写入输出状态的时间间隔。VARIABLE：设置VARIBALE=ALL表明所有应用到该过程的变量和材料类型将被写入输出数据库。设置VARIABLE=PRESELECT表明当前过程的默认输出变量被写入输出数据库。额外的输出请求可以与*OPTION其他选项联合使用得到。如果省略该参数，则只有那些各自输出选项的输出变量请求写入输出数据库。包含MODE LIST参数的数据行，列表期望的特征模态第一行1、执行一系列期望的特征模态Explicit分析中的输出必需的、相互排斥的参数：FIELD：包含该参数表明输出场类型HISTORY：表明输出历程类型下列参数参数是可选的并只用于FIELD：NUMBER INTERVAL：设置该参数等于输出数据写入时间间隔的数量。Explicit总是写分析步开始时的结果，例如，如果NUMBER INTERVAL=10，explicit会写入11个输出数据包含分析步开始时的值。该参数默认值必须是正整数或0。0表示抑制所有输出。如果该参数被省略，值是20。TIME MARKS：设置TIME MARKS=NO(默认)在增量步结束立刻写出结果到输出数据库。=YES则在准确时间写入结果。TIME POINTS：未完设置该参数等于*TIME POINTS的名字来定义输出写入的时间点。基本形式：*OUTPUT，FIELD，VARIABLE=PRESELECT含义：将ABAQUS默认的场变量写入ODB文件。*OUTPUT，HISTORY，VARIABLE=PRESELECT含义：将ABAQUS默认的历史变量写入ODB文件。*PREPRINT：设置DAT文件中记录的内容。可选参数：CONTACT：(该参数只用于standard分析)设置contact=YES，则打印由接触对定义数据产生的接触约束的详细信息；NO(默认)，抑制该输出。ECHO：设置ECHO=YES打印输入数据的响应；NO(默认)抑制该输出。History：设置History=YES打印历史数据；NO(默认)抑制该输出。Model：设置Model=YES打印模型定义数据；NO(默认)抑制该输出。Parsubstitution：设置parsubstitution=YES打印输入模型参数化自由原始数据的更改版本；NO(默认)抑制该输出。Parvalues：设置parvalues=YES打印原始输入文件的更改版本并先死用来模型参数化的参数和相应数值；NO(默认)则抑制该输出。该选项没有数据行。Abaqus的默认设置：*Preprint, echo=NO, model=NO, history=NO, contact=NO。表示：在DAT文件中不记录对INP文件的处理过程，以及详细的模型和历史数据。*PRE-TENSION SECTION：关联一个预紧节点到一个预紧截面必需参数：NODE：设置该参数等于预紧节点号码或包含预紧节点的节点集(只包含一个节点)名称。必需的、相互排斥的参数：ELEMENT：设置该参数等于用来定义预紧截面的杆或梁单元号码，或单元集。仅一个单元。SURFACE：设置该参数等于*SURFACE定义的面的名称，用来定义预紧截面。数据行，定义截面法向，可选的第一行(仅一行)1、法向的第一个分量2、法向的第二个分量3、法向的第三个分量如果该数据行被省略，则standard会为连续单元计算一个平均法向；而杆或梁单元的法向是两个端点方向。*RESTART：保存和重用数据和分析结果警告：该选项可能产生非常大的数据使用在standard中以下参数至少一个是必需的READ：包含该参数表示该分析是以前分析的重启动。基本模型定义数据(单元、材料、节点)不能在这样的重启动被改变。但是可以添加单元集、节点集和幅值表，而且与该零件有关的、已经分析的历史数据也能被更改。WRITE：包含该参数表示将写入重启动数据。使用READ参数时可选的其它参数：END STEP：该参数表明用户希望结束设置重启动的分析中的当前分析步。当用户希望重定义载荷历史、输出选项等，可使用该参数；如果包含该参数，数据必须包含更多的分析步以定义分析如果继续。如果省略该参数，软件会继续分析直到完成当前定义重启动的分析步。INC：设置该参数等于由STEP参数指定的分析步内的增量数，其后面的分析将会恢复。如果省略该参数，重启动将在指定的分析步的结尾开始。ITERATION：如果新的分析是由前一个直接循环分析重启动，设置该参数等于指定分析步内的积分数。由于新分析只能从前一个直接循环分析的载荷循环结尾重启动，所以如果INERATION参数已经设置，则INC参数可忽略。STEP：设置该参数等于重启动的分析步数。如果省略，分析将重启动最后一个可用的分析步。使用WRITE参数时可选的其它参数：FREQUENCY：该参数指定重启动将写入的增量。例如，FREQUQNCY=2，则将在增量步2、4、6写入重启动信息。没有数据行基本形式：*RESTART，WRITE，FREQUENCY=0(含义：不输出用于重启动分析的数据。)*RIGID BODY：定义刚体及其属性必需参数：REF NODE：设置该参数等于刚体参考节点的节点号或节点集。可选参数：ANALYTICAL SURFACE：要分配到刚性体的解析面名称ELSET：要分配给刚性体的单元集合名称，一个单元不能分配给多个刚体ISOTHERMAL：只用于热力耦合分析。设置=YES指定等温线刚体，默认是NO。PIN NSET：设置该参数等于包含销钉类型的节点名来分配给刚性体。该参数可用来添加节点到刚体或重定义节点类型。销钉类型的节点只有移动自由度与刚体关联。POSITION：如果参考节点由用户定义，设置POSITION=INPUT(默认)。如果参考点位于刚体的质心，则设置POSITION=CENTERTIE NSET：设置该参数等于包含绑定类型的节点名来分配给刚体。该参数可用来添加节点到刚体或重定义节点类型。绑定类型的节点的移动和转动自由度都与刚体关联。仅当指定单元集包含刚性单元时的可选参数DENSITY：仅用于explicit分析。单元集内所有刚性单元的密度NODAL THICKNESS：仅用于explicit分析。表明刚性单元的厚度不从数据行读入，而是由*NODAL THICKNESS选项指定。OFFSET：仅用于explicit分析。定义从单元中面到参考面的距离。由于对刚性单元来说，没有单元级别的计算，所以指定的偏移值只影响由刚性单元构成的刚性面的接触对的操作。该参数接受正值或负值或标记SPOS和SNEG。正值代表单元法向方向。当OFFSET=0.5(或SPOS)，刚性单元的顶面是参考面。当当OFFSET=-0.5(或SNEG)，刚性单元的底面是参考面。默认是OFFSET=0，表明刚性单元的中面是参考面。Standard分析中无数据行对于explicit中的R2D2、RB2D2、RB3D2的数据行第一行，仅一行1、单元的横截面面积，默认是0对于explicit中的RAX2、RB2D2、R3D4的数据行第一行，仅一行1、单元的厚度，默认是0基本格式：*RIGID BODY，ref node=参考点集名称，analytical surface=解析面集名称*RIGID SURFACE：定义解析刚性面*ROTARY INERTIA：定义刚体转动惯量该选项用来定义与ROTARYI单元相关的刚体转动惯量的数值。也用在standard分析中，定义质量比例阻尼(对于直接积分动力学分析)和复合阻尼(对于模态动力学分析)。必需参数：ELSET：包含ROTARYI单元的单元集可选参数：ALPHA：只用于standard分析。当用于直接积分动力学，设置该参数等于 因子来为ROTARYI单元创建质量比例阻尼。在模态分析时忽略该参数。默认是0COMPOSITE：只用于standard分析。当用于模态动力学计算复合阻尼因子时，设置该参数等于临界阻尼。在直接积分动力学时忽略该参数。默认是0定义转动惯量的数据行第一行，仅一行1、关于局部坐标轴1的转动惯量 2、关于局部坐标轴2的转动惯量 3、关于局部坐标轴3的转动惯量 4、 5、 6、*SECTION CONTROLS：指定截面控制警告：对于沙漏控制，使用大于默认值会产生额外的刚度响应，甚至当值太大时有时导致不稳定。默认沙漏控制参数下出现沙漏问题表明网格太粗糙，因此，更好的解决办法是细化网格而不是施加更大的沙漏控制。该选项用来为减缩积分单元选择非默认的沙漏控制方法，和standard中的修正的四面体或三角形单元或缩放沙漏控制的默认系数；在explicit中，也为8节点块体单元选择非默认的运动方程：为实体和壳选择二阶方程、为实体单元激活扭曲控制、缩放线性和二次体积粘度、设置当单元破损时是否删除他们、或为上述完全破损的单元指定一标量退化参数。等必需参数：NAME：名字可选参数：DISTORTION CONTROL：只用于explicit分析。=YES激活约束防止负体积单元出现或其他可压缩材料的过度变形，这对超弹材料是默认的。DISTORTION CONTROL参数对线性动力学不可用并且不能防止单元由于时间不稳定、沙漏不稳定或不切实际的物理变形造成的扭曲。 =NO不激活约束，对于除了超弹材料等单元都是默认值。HOURGLASS：设置HOURGLASS=COMBINED定义沙漏控制的单元粘性-刚度形式；HOURGLASS=ENHANCED基于假定的增强应变方法来控制沙漏；HOURGLASS=RELAX STIFFNESS使用整合的粘弹性形式控制沙漏；HOURGLASS=STIRRNESS对于standard分析除了超弹材料和修正的四面体和三角形外的单元默认，为所有减缩积分单元定义沙漏控制是严格的弹性；HOURGLASS=VISCOUS为缩减积分单元定义沙漏阻尼。数据行：第一行，仅一行1、对于沙漏刚度的比例因子，影响移动自由度。如果为空，默认值是1.0。建议范围是0.232、对于沙漏刚度的比例因子，影响转动自由度。如果为空，默认值是1.0。建议范围是0.233、对于沙漏刚度的比例因子，在explicit中，影响小应变壳单元的超出平面的移动自由度，如果为空，默认值是1.0，建议范围是0.23。不用于standard中。4、explicit中的线性体积粘度的比例系数。如果为空，默认值是1.0，建议范围是0.01.0。不用于standard中。5、explicit中的二次体积粘度的比例系数。如果为空，默认值是1.0，建议范围是0.01.0。不用于standard中。 接触分析1、 塑性材料和接触面上都不能用C3D20R和C3D20单元，这可能是你收敛问题的主要原因。如果需要得到应力，可以使用C3D8I (在所关心的部位要让单元角度尽量接近90度)，如果只关心应变和位移，可以使用C3D8R, 几何形状复杂时，可以使用C3D10M。2、 接触对中的slave surface应该是材料较软，网格较细的面。3、 接触面之间有微小的距离，定义接触时要设定“Adjust=位置误差限度”，此误差限度要大于接触面之间的距离，否则ABAQUS会认为两个面没有接触：*Contact Pair, interaction=SOIL PILE SIDE CONTACT, small sliding, adjust=0.2.4、 定义tie时也应该设定类似的position tolerance:*Tie, name=ShaftBottom, adjust=yes, position tolerance=0.15、 msg文件中出现zero pivot说明ABAQUS无法自动解决过约束问题，例如在桩底部的最外一圈节点上即定义了tie，又定义了contact, 出现过约束。解决方法是在选择tie或contact的slave surface时，将类型设为node region, 然后选择区域时不要包含这一圈节点（我附上的文件中没有做这样的修改）。6、 接触定义在哪个分析步取决于你模型的实际物理背景，如果从一开始两个面就是相接触的，就定义在initial或你的第一个分析步中；如果是后来才开始接触的，就定义在后面的分析步中。边界条件也是这样。7、 我在前面上传的文件里用*CONTROL设了允许的迭代次数18，意思是18次迭代不收敛时，才减小时间增量步（ABAQUS默认的值是12）。一般情况下不必设置此参数，如果在msg文件中看到opening和closure的数目不断减小（即迭代的趋势是收敛的），但12次迭代仍不足以完全达到收敛，就可以用*CONTROL来增大允许的迭代次数。8、 桩头掉在了地表下，说明接触定义得不正确。可能接触面的距离还是大于*contact pair 中的adjust=0.02, 可改为adjust=0.29、 原则上，90度的圆弧应该划分10个单元，适当少一些可能也行。*contact pair中的adjust=0.005，还是太小，在后处理时可以看到，接触面之间的距离大于0.005。把adjust设置大一些没关系，比如adjust=0.1。10、 网格不好也可能产生过约束问题，不要只考虑边界条件啊！ NUMERICAL SINGULARITY WHEN PROCESSING NODE 15294 D.O.F. 2 RATIO =2.48305E+11, 说明NODE 15294 所在的实体在方向2上出现无限大的刚体位移。可以在此实体上的任意一点和地面之间定义一个很软的弹簧，以消除刚体位移。方法是：interaction模块，菜单special / springs-dashpots / create, 选connect points to