ADINA-第5章-求解控制

第五章求解控制1.总体说明2.求解控制主要内容总体说明求解控制总体说明：1.不同的分析涉及不同的求解控制，本章主要介绍各个模块求解通用的求解控制选项，分析类型相关的求解控制在《分析类型》一章中详细介绍；2.软件对多数的求解控制设置缺省选项，如果没有明确的目的，可使用这些缺省项；定义分析标题：菜单：Control/Heading说明：分析标题的定义是可选项。分析标题求解控制定义模型自由度：菜单：Control/DegreeOfFreedom…定义求解的自由度范围，点中的自由度为求解自由度；对于2D分析，求解自由度为Y-Translation和Z-Translation；在3D实体单元的模型中，求解自由度为Y-Translation，Y-Translation和Z-Translation，即三个平移自由度；如果求解自由度包括不存在的自由度，则生成求解文件*.dat时，给出警告信息（不影响计算过程和结果）：ADINAdatainputfileE:\result\UMAT\1.datsuccessfullyupdated.***WARNING:Noelementconnectionfornode1directions4,5,6Noelementconnectionfornode2directions4,5,6…定义模型自由度求解控制定义时间函数：菜单：Control/TimeFunction•无论是在静力分析或瞬态分析中，ADINA中任何载荷都必须定义其随时间的变化关系；•对于瞬态问题，时间为真实时间，基于当前时间（步）计算加速度、速度、应变率等物理量；•在静态分析中，时间为“伪时间”，仅作为计数器使用，因此在静态分析中，“时间”可认为其没有单位。载荷时间ADINA中的载荷与时间相关ADINA中的载荷和时间：定义时间函数求解控制时间函数定义了载荷比例因子随时间的变化过程，因此时间函数是没有单位的；缺省的时间函数（标号1）是载荷比例因子为1并随时间没有变化，即从时间0～1e20，载荷比例因子恒定为1；用于计算的载荷随时间变化曲线是定义载荷时输入的单一数值（Magnitude）与时间函数中各个时刻的比例因子乘积的结果；1.0载荷比例因子时间缺省的时间函数定义时间函数求解控制选择时间函数类型，共提供四种函数形式；输入窗口的Value值即为f*(t)；第一种函数时间是通过输入各个时刻的f*(t)值直接定义比例因子，因此用户可以通过控制离散点的多少，来模拟任何函数；第二种为周期函数，窗口中输入的Value如果恒定为1，则是标准周期函数，如果Value值变化，则此时间函数在发展过程中周期不变，幅值发生变化；函数三、四常用于电学分析中；四种时间函数形式定义时间函数求解控制定义时间步长：菜单：Control/TimeStep对于多数分析（包括所有的非线性分析），其求解过程需要分成很多步。由于ADINA规定所有载荷随时间的变化关系，因此通过定义求解过程的时间步大小就决定了每个时间步的载荷大小。窗口中的第一列数值为时间步个数（整数），对应行输入数值为时间步长（实数）；总的求解时间为每行两数值积之累加和。总时间步为26步，总求解物理时间为2.1时间单位；定义时间步长求解控制载荷时间Dtn＝0.4Dtn+4＝0.05任意时间步中的载荷大小在载荷-时间曲线中确定；在右图中，第n个时间步长为0.4，在前一时间步结果收敛时其载荷为Ln0进入此时间步求解时，载荷增加到Ln1；相应于时间步n的载荷增量为Ln1与Ln0之差；Ln0Ln1定义时间步长求解控制流固耦合开关：菜单：Control/AnalysisAssumptions/FluidStructureInteraction（FSI）设定模型是否进行流固耦合求解：如果选择如右图所示选择FluidStructureInteraction则后续菜单中出现定义流固耦合边界的菜单；如果选择IncludesPorousCoupling，并在ElementGroup中选择ElementOption为Porousmedia，则后续菜单中出现定义多孔介质材料特性的菜单；定义分析假设－流固耦合开关求解控制几何非线性开关：菜单：Control/AnalysisAssumptions/Kinematics确定几何非线性形式：1.SmallDisp/Rot和SmallStrains时，模型不存在几何非线性；2.LargeDisp/Rot和SmallStrains时，模型为大位移几何非线性问题；3.LargeDisp/Rot和LargeStrains时，模型为大位移及大应变几何非线性问题；[A]:壳单元的压力更新选项，主要用于频域动力学分析时；[B]:非协调单元选项（适用于2D和3Dsolid），非协调元引入节点无关的多余自由度，有利于解决变形锁定问题，但单元边界位移不再满足协调条件。[A][B]定义分析假设－几何非线性开关求解控制质量矩阵形式选择：菜单：Control/AnalysisAssumptions/MassMatrixLumped：集中质量矩阵（质量矩阵对角化），在进行中心差分方法求解瞬态动力问题时自动采用集中质量矩阵；Consistent：一致质量矩阵，为稀疏矩阵；定义分析假设－质量矩阵形式求解控制定义模型阻尼：菜单：Control/AnalysisAssumptions/DefineRayleighDampingFactors定义整个模型或各个单元组的阻尼，采用标准Rayleigh质量和刚度阻尼形式；右图中对整个模型指定5％的质量阻尼。定义分析假设－模型阻尼求解控制在响应谱分析中定义各阶模态的阻尼比；定义阻尼比：菜单：Control/AnalysisAssumptions/DefineModalDampingFactors定义分析假设－模型阻尼比求解控制温度条件设定：菜单：Control/AnalysisAssumptions/DefineTemperatureSettings用于指定整个模型的初始温度和初始温度梯度（仅对于壳单元）；如果需要指定不同的初始温度或初始温度梯度，在Model/InitialCondition中定义；如果需要指定温度载荷，一般在Model/Loading/Temperature中定义，而不要使用此窗口的Prescribed…选项；定义分析假设－温度条件设定求解控制方程求解设置：菜单：Control/SolutionProcess[A]:求解开始时间（如为重启动分析，择输入前一次求解的终止时间。）[B]:重启动分析可以从静力分析转到动力分析（如预应力模态计算），或从动力分析转到静力分析；ADINA重启动分析具有丰富的功能，一般地可以在重启动模型进行如下修改：自由度（Fix）类边界条件；时间步定义；时间函数定义；约束方程和刚性联接；Rayleigh阻尼；收敛准则，自动时间步长设置等等；定义分析假设－方程求解设置求解控制[C]:开关：出现刚度矩阵特征值为零是否继续计算；对于线性问题，除非模型中定义势流体介质，否则不应出现为零的特征值；对于非线性问题，除非模型中定义了：自动时间步长控制；单元生死；载荷位移控制；接触；势流体；否则不应出现为零的特征值；定义分析假设－方程求解设置求解控制[D]:单元分组当模型规模很大或进行并行求解时使用。一般D1＝12，D2=10*D1,D3=999999；[E]:迭代方法选项（见后面）[F]:收敛准则选项（见后面）[G]:求解器选项（见后面）[E][F][G]定义分析假设－方程求解设置求解控制[E]:迭代方法选项弹出如下选项:[a]:选择迭代方法共有四种：FullNewtonMethod（缺省）ModifiedNewtonMethodBGFSMatrixUpdateMethod[b]:是否使用Line搜索方式[c]:最大迭代次数一般为8～15次；[d]:是否输出每个迭代的结果[a][b][c][d]定义分析假设－方程求解设置-迭代方法求解控制[F]:收敛控制[a]:迭代收敛准则类型Energy(0.001)EnergyandForceEnergyandDisplacementForce（0.01）Displacement(0.01)最常使用的为Energy和EnergyandForce准则；[b]：能量收敛准则[C]：接触力收敛准则[d]：最小参考接触力大小[e]：Line搜索方式收敛准则[f]：Line搜索方式最小参考能量值[g]：力收敛准则[h]：参考力[i]：参考力矩[j]：位移收敛准则设置（同力准则）[b][c][d][e][g][h][i][f][a][j]定义分析假设－方程求解设置-收敛控制求解控制[G]:求解器选择求解器共有四种：DirectIterativeSparseMultigrid缺省求解器为Sparse求解器；如果选用Iteratice迭代求解器，则弹出右侧窗口，一般使用缺省设置即可；定义分析假设－方程求解设置-求解器选择求解控制说明：ADINA-F模块包括显式时间积分求解器Explicit，与自动时间步长中的CFL数相应设置进行显式求解。其它求解控制选项：菜单：Control/MiscellaneousOptions[A][C][B][D][E][F][A]:初始应变作为：InitialStrainsInitialStressInitialStressthatcauseDeformation初始条件输入只有应变没有应力，可将应力作为应变从Model/InitialCondition中输入，然后在此选择InitialStress，则输入数据当作初始应力参与计算；[B]:计算质量特性计算模型的质量和质心坐标；[C]:计算支反力/力矩[D]:在重启动文件中写入载荷向量[E]:指定壳单元刚度为零节点刚度[F]:单元死的时间周期其它求解控制选项求解控制写出重启动文件：菜单：Control/RestartFile（.Res）[A]:写出重启动文件的时间间隔，缺省为0，则在计算结束时写出重启动文件；[B]:缺省为选择Overwrite，即后面写出的重启动文件覆盖前面的重启动文件，即只保留一个重启动状态；写出重启动文件控制求解控制[A][B]OUT文件信息控制：菜单：Control/Printout（.out）/VolumeADINA在求解时生成.out文件，从计算初始化到计算结果，以及如果出现错误的出错信息都写道这个文件中。此文件可以在求解过程中打开查看;此菜单中的三个功能(Volume/ResultatNodes/TimeStep)是确定输出到.out文件中的内容;通常使用缺省设置，否则可能导致.out文件很大。OUT文件信息控制求解控制写出节点结果到OUT文件：菜单：Control/PrintoutControlofNodalResults定义节点的Blocks（组），在求解过程中输出相应的结果到OUT文件，便于查看并监视求解过程。节点结果通常为位移、速度、加速度、压力等等；与下页的定义的结果输出TimeStep联合使用；定义了节点1201，1202，1203，1001，1002共5个节点的结果写出到.out文件中；OUT文件信息控制－写出节点结果到OUT文件求解控制写出节点结果到OUT文件的时间步：菜单：Control/SpecifyPrintoutTimeSteps缺省情况下，每个时间步的计算结果都写到OUT文件中；此功能是将指定时间步的结果输出，主要用于控制OUT文件的大小；与上页的NodalResults联合使用输出结果的时间步为1、2、5、8、11、…OUT文件信息控制－写出节点结果到OUT文件的时间步求解控制结果文件信息控制：菜单：Control/SpecifyPrintoutTimeStepsADINA的结果文件名为.por(t),此组菜单控制输出到.por文件的内容。[A][B][C][D][A]:结果文件.por格式通常采用Binary二进制格式，TEXT格式主要用于对ADINA进行开发的用户使用；[B]:最大