汽车仿真分析基础-2010指南.pptVIP

接触类型分类 接触类型可以按下述方法分类 按接触点搜索方法分类：connectivity, bucket sort 接触反力计算方法：penalty, constraint 接触面的定义方法：面-面接触(symmetry)，点-面接触(non-symmetry) 将接触点搜索方法，反力计算和接触面定义组合起来，形成接触类型： SURFACE_TO_SURFACE（C, P,S） NODE_TO_SURFACE（C, P,N） AUTOMATIC_SURFACE_TO_SURFACE（B, P, S） CONSTRAINT_SURFACE_TO_SURFACE（B,C, S） AUTOMATIC_SINGLE_SURFACE（B, P, N） 使用Soft constraint，penalty接触类型使用constraint方法计算反力。 1.接触面搜索 接触处理是计算中最耗时的工作，而其中最主要的部分是接触面搜索。 有接触可能的零件，segment在输入文件中已经定义，接触面搜索的工作在于确定哪个部分与哪个部分相接触。在程序内部，调查某一特定的slave节点与哪个master segment接触。 为了选出相接触的segment，需要找到与slave节点最近的master节点（构成master segment的节点）。这个工作是单纯地计算slave节点和各master节点间的距离，选出最近的节点，但这需要花费许多时间，因此，需设法节约计算时间。 接触面搜索方法特征 1.接触面搜索 方法之一是利用单元形状的连续性。将上一个时间步接触的master segment事先记忆下来，首先确认该segment是否接触。如果没有可能接触，就从与slave节点最靠近的、构成master segment的节点开始顺序搜索。 方法之二是利用bucket sort方法。这个方法是将节点分成组，在每个组内计算节点距离，达到减少计算工作量的目的。 无论使用哪种方法，第一步的工作都是要找到最近的master节点。 接触面搜索方法特征 2.利用单元形状的连续性 利用master segment形状的连续性计算接触点。在接触范围限定的情况下，效率高，但当网格不连续时，可能带来接触不安定。 接触面搜索方法特征 3.利用三维节点的空间布置的方法 特征：接触范围广时，搜索高效，不需考虑接触面的方向。 说明：将包含网格的空间分为许多立方体（MAX5000）；搜索包含master，slave的立方体。 分析错误：立方体最大为5000，立方体大小依存于网格。如果有单元发散而产生大变形时，一个立方体内包含segment将超过2000，程序会给出错误信息（more than 2000 nodes in bucket)。Bucket sort并非每个时间步都做，single_surface（25），automatic_single_surface(100)。这种搜索方法与segment指向无关，壳单元必须考虑板厚。 接触面搜索方法特征 3.接触面搜索方法不同的影响 接触面搜索方法特征 3.接触面搜索方法不同的影响 接触面搜索方法特征 　　由接触面搜索，确定与slave节点接触的segment后，下面计算slave节点与该segment的接触位置。 　　 接触反力计算方法特征 　　确定接触位置后，给slave节点和master节点加上接触反力。接触反力的计算方法有两种。 接触面反力计算方法特征 penalty方法 (2) Constraint 方法 k:接触弹簧刚度，l:接触面穿透深度，n:接触面法向矢量，fsi: penalty系数， K:体积弹性模量，A:接触面积，V:接触单元体积，dt:时间增量，ms:slave节点集中质量 (1) Penalty方法的特征： 该方法的关键在于在接触节点和接触面上设定有垂直弹簧。弹簧刚度以接触面上的单元为单位计算，而且用户可以进行刚性参数的调节。 　　最初，确认slave节点是否穿透主面，如果没穿透，当然不进行接触反力计算。如有穿透，接触反力被加载在从节点和主面节点之间，接触反力的大小与穿透量成反比。 接触反力计算方法特征 (1) Penalty方法的特征： 接触反力计算方法特征 If l0 If l0 立体单元 壳单元 li:壳单元对角线长，fsi: penalty系数（通常取0.1）， Ki:体积弹性模量，Ai:接触segment面积，Vi:接触单元体积 接触面定义 面-面接触(symmetry)，点-面接触(non-symmetry) SURFACE_TO_SURFA