摘要:
检查post processor是否正确project 布局每次移动大小和的0.8倍(步长=0.8)1.0倍可见重叠现象仍然很严重把最大迭代次数设为100后不解算重叠左图迭代100次,步长1.0,可以正确解算(为什么commands半径不一样???因为这些元素自身大小也是依赖于布局决定)右图迭代1000次,interfaces位置发生变化步长改为0.2,也可以正确解算,但布局差异较大to... 阅读全文
摘要:
确定较好的边绑定参数没有对轮廓进行重采样平滑参数为0.1,2没有对轮廓进行重采样平滑参数为0.5,8对轮廓做2倍重采样,平滑参数0.1,2出现了一些“主干道”2倍重采样平滑参数0.5,2“主干道”较为平滑2倍重采样平滑参数0.1,104倍重采样平滑参数0.5,10得出一些经验,平滑迭代次数较低的时候,曲线不容易分散采样倍数太高时,计算时间增多,曲线不光滑为了把一些边更好地绑定在一起,有... 阅读全文
摘要:
搞清楚平面布局类的各个参数orthoLayout此类的getBoundingBox()函数可以返回包围盒separationcOverhangoptionsmargin结果200.40010.400可见separation控制线之间的距离此值不能是0201000cOverhang似乎没什么影响,但是它的值也不能是0200.40100margin对结果没有什么影响当结果有多个连通分量时,m... 阅读全文
摘要:
检查不同平均距离计算方法的结果1.所有点对距离的调和平均2.每个点与其他点最短距离的算术平均总半径 =5622.2026没有明显变化。用第二种方法似乎布局反而更稀疏。3. 采用2的方案,再把sparseFactor从0.2改为1.0总的坐标范围反而变小,总半径 =3889.0559按道理,avgDist *scaleRate = 2 * avgSize * m_sparseF... 阅读全文
摘要:
检查不同布局算法备注CircularLayout可以非连通FastMultipoleMultilevelEmbedderFMMMLayout可以非连通StressMajorization未成功必须连通GEMLayout未成功BalloonLayout未成功需要为树?PlanarizationLayout来自为知笔记(Wiz) 阅读全文
摘要:
当前一些不合理现象来自为知笔记(Wiz) 阅读全文
摘要:
11月17日c++调用fortran用intel fortran 编译器编译fortran 代码,设置以下地方fortran->compatibility->enable f77 runtime compatibility truefortran->compatibility->use f77 integer constants truefortran->... 阅读全文
摘要:
11月22日当前任务stiffness matrix 究竟怎么算?什么是tangent stiffnesmatrix?在Exact Corotational Linear FEM Stiffness Matrix之中linear stiffness matrix Krest 是指小变形的情形下,节点的力f = Krest * 节点位移xtangent stiffness matri... 阅读全文
摘要:
11月21日了解tangent stiffness matrix相关论文了解动力学方程相关论文来自为知笔记(Wiz) 阅读全文
摘要:
11月20日ImplicitBackwardEulerSparse::DoTimestep的总步骤:1通过当前的位移,计算出内力和刚度矩阵2 用刚度矩阵和质量矩阵构造阻尼矩阵3 用求解器解出qdelta,就是下一步的速度增量(为什么?求解器的矩阵在哪里更新?) 求解器实际上求解tangentStiffnessMatrix * qdelta =bufferConstrained 于... 阅读全文
摘要:
11月19日成功运行了Vega测试代码发现Maya 的obj支持体网格结构obj的面格式可以只指定顶点索引,不知道发现和贴图坐标,例如f 1 2 3VolumetricMesh规定每个元素(element)的顶点数是一样的,例如四面体网格点数必须是4,不允许同时存在点数不一样的元素想知道质量矩阵的计算方法发现质量矩阵的大小与元素的数目不相等,而与顶点数目有关质量矩阵是一个方阵,行数=列数=点数*3... 阅读全文
摘要:
11月26日utilities里面几个程序displayObj 显示obj文件generateInterpolantgenerateInterpolationMatrixgenerateSurfaceMesh veg体网格转obj面网格,内部点保留objMergeFilesgenerateMassMatrix 计算质量矩阵largeModalDeformationFactory 可以计算模态系数i... 阅读全文
摘要:
11月25日发现几种术语tangent linear vibration modelinear modal analysis(LMA) modelinear modal analysis(LMA)估计是同线性弹力模型结合的,此时力是位移的线性函数,以此得到刚度矩阵tangent linear vibration mode估计是根据某点的力对位置导数得到tangent stiffness matri... 阅读全文
摘要:
11月24日彻底明白一致质量矩阵、刚度矩阵的计算方法明白了中心差商求解器的计算流程下一步是要了解降维积分的做法来自为知笔记(Wiz) 阅读全文
摘要:
当前任务读入面网格,完成体网格化读入面网格以及控制它的参数,能够构造一个体网格黑盒tetgen支持的文件格式tetgen命令行参数-p Tetrahedralizes a piecewise linear complex (PLC). 体网格化,网格没有洞-Y Preserves the input surface mesh (does not modify it).-r Reconstructs... 阅读全文
摘要:
11月27日任务:研究reducedDynamicSolver-rt,了解降维后的积分方法代码结构:函数initScene()准备必要的资源,按顺序为:调用了ReadMatrixFromDisk_()从磁盘读取低维基向量,放入ModalMatrix 对象renderingModalMatrix之中,这个类可以由高维力算出低维力。同时分配低维坐标q 低维力 fq fqBase。其中fqBase表示在... 阅读全文
摘要:
12月8日开始准备数值实现tangentStiffnessMatrix在ForceModel类中得到,同时还可以得到力的值首先构造TetMesh对象,再构造CorotationalLinearFEM对象,最后构造ForceModel对象。即可从任意变形得到力和tangentStiffnessMatrix发现一点,ForceModel计算出来的力是每个节点受到的内力的反向来自为知笔记(Wiz) 阅读全文
摘要:
12月9日构造各个矩阵内力f1 = f0 + K * dxK为tangent stiffnes matrixdx为当前变形下的一个微小位移对coratational模型调用testTangentStiffnessMatrix函数,发现原来的wrapped stiffness matrix并不准确。只有使用准确的tangent stiffness matrix ,误差才为位移的高阶无穷小。(与步长长... 阅读全文