提取刚度矩阵
有限元刚度矩阵和质量矩阵提取
有限元刚度矩阵和质量矩阵提取一、概述有限元方法是一种常用的数值计算方法，它将复杂的物理问题离散化为简单的几何体元素，并在每个元素内部进行近似计算。
在有限元分析中，刚度矩阵和质量矩阵是两个重要的矩阵，它们提供了系统的结构信息和物理特性。
本文将介绍有限元刚度矩阵和质量矩阵提取的方法。
二、有限元刚度矩阵提取1. 刚度矩阵定义刚度矩阵是描述结构物体在受到外力作用下所产生的应变能与外力之间关系的一个重要参数。
对于一个n自由度系统，其刚度矩阵K为n*n的实对称正定矩阵。
2. 刚度矩阵推导假设一个二维平面三角形单元，其节点数为3个，分别为1、2、3号节点，其自由度数为6个（每个节点有2个自由度）。
则该单元刚度矩阵K可以表示为：K = [k11 k12 k13 k14 k15 k16;k21 k22 k23 k24 k25 k26;k31 k32 k33 k34 k35 k36;k41 k42 k43 k44 k45 k46;k51 k52 k53 k54 k55 k56;k61 k62 k63 k64 k65 k66]其中，kij表示单元局部坐标系中第i个自由度受到第j个自由度作用时的刚度系数。
对于三角形单元，其刚度矩阵可以通过以下公式推导得到：kij = ∫∫B^TDBdΩ其中，B为单元形函数的梯度矩阵，D为材料弹性模量与泊松比的组合参数，Ω为单元面积。
3. 刚度矩阵组装在有限元分析中，通常需要将多个单元组装成一个整体系统。
这时需要将各个单元的局部刚度矩阵按照节点编号和自由度顺序组装成全局刚度矩阵。
三、有限元质量矩阵提取1. 质量矩阵定义质量矩阵是描述结构物体在振动或运动过程中所具有的惯性特性的一个重要参数。
对于一个n自由度系统，其质量矩阵M为n*n的实对称正定矩阵。
2. 质量矩阵推导假设一个二维平面三角形单元，其节点数为3个，分别为1、2、3号节点，其自由度数为6个（每个节点有2个自由度）。
则该单元质量矩阵M可以表示为：M = [m11 m12 m13 m14 m15 m16;m21 m22 m23 m24 m25 m26;m31 m32 m33 m34 m35 m36;m41 m42 m43 m44 m45 m46;m51 m52 m53 m54 m55 m56;62 63 64 65 66]其中，mij表示单元局部坐标系中第i个自由度的质量。
abaqus提取整体刚度矩阵
abaqus提取整体刚度矩阵abaqus是一款常用的有限元分析软件，主要用于模拟和分析各种工程结构的力学行为。
在abaqus中，整体刚度矩阵是一个重要的概念，它能够描述结构在受力作用下的刚度特性。
本文将深入探讨abaqus 如何提取整体刚度矩阵，并分享对该概念的观点和理解。
一、整体刚度矩阵的概念整体刚度矩阵是指在有限元分析中，将结构划分成若干个离散的单元后，通过单元刚度矩阵的叠加得到的描述结构整体刚度特性的矩阵。
整体刚度矩阵反映了结构在受力作用下的刚度响应，是进行结构力学分析的重要工具。
二、abaqus提取整体刚度矩阵的方法在abaqus中，提取整体刚度矩阵的方法主要有以下步骤：1. 创建有限元模型：需要在abaqus中创建一个准确表达所研究结构的有限元模型。
这包括定义结构的几何形状、材料性质以及边界条件等。
2. 定义材料属性：在有限元分析中，材料的力学性质对整体刚度矩阵具有重要影响。
在abaqus中需要明确定义结构中所使用的材料的力学性质，包括弹性模量、泊松比等。
3. 定义加载条件：接下来，需要定义结构在受力作用下的加载条件。
这可以是施加在结构上的力或约束条件等。
4. 进行力学分析：有了有限元模型、材料属性和加载条件后，就可以进行力学分析。
在abaqus中，通常使用有限元方法求解结构的响应，得到结构的位移和应力等。
5. 提取整体刚度矩阵：通过分析结果，abaqus提供了方便的工具来提取整体刚度矩阵。
用户可以在abaqus的后处理模块中选择相应的输出选项来得到整体刚度矩阵的结果。
三、对整体刚度矩阵的理解整体刚度矩阵是结构力学分析中的一个关键概念，对于研究和理解结构的强度和刚度特性具有重要意义。
整体刚度矩阵可以用来计算结构在受力作用下的位移、应力和应变等响应，进而评估结构的安全性和可靠性。
从数学角度看，整体刚度矩阵是由单元刚度矩阵叠加得到的。
单元刚度矩阵描述了单个有限元单元在特定边界条件下的刚度特性。
提取单元刚度矩阵
单元刚度矩阵的提取刚度矩阵在有限元求解过程中扮演者非常重要的角色，以最小位能原理求解过程为例最终越是转换为含有结构刚度矩阵的能量泛函的取值问题。
有限元过程中涉及到三类刚度：单元刚度矩阵，组合结构刚度矩阵和最终求解刚度矩阵。
其中单元刚度矩阵：仅与单元的自身自由度有关，同一编号的单元矩阵的维数是固定。
组合结构刚度：矩阵根据求解的初始变量个数决定刚度矩阵的维数，属于单元组装后的初始刚度，维数和整个单元初始变量个数相等。
最终求解刚度矩阵：代入边界条件简化后的刚度。
以《Finite Element Analysis-Theory and Application With ANSYS》中的梁单元例子为例，解释刚度提取过程：此模型的单元刚度矩阵：（学则beam3梁单元后，该单元包含两节点，每个节点具有三个自由度，因此对应单元刚度矩阵为6*6的方阵）组合结构刚度矩阵：（该结构含有三个节点，每个节点具有三个原始自由度，因此组合结构刚度矩阵具有9*9阶的形式）最终求解刚度矩阵：（由于边界条件的存在，该结构中，1,3点的自由度不存在，求解参数中有六个参数已知，因此对最终求解刚度矩阵为三阶方阵）通过最终的刚度矩阵组成的方程，求解出2节点的位移解，再以这些原始解得出应力，应变，支反力的其他的解。
ansys实现过程：提取思路如下：通过/debug提取单元刚度矩阵，通过filname.full文件提取后两者的矩阵ansys实现过程如下：finish/clear/filname,k,1/prep7N,1N,2,120N,3,120,-108et,1,beam3mp,ex,1,3.0e7mp,prxy,1,0.3R,1,7.65,204,10E,1,2E,2,3/debug,-1,,,1,,,,,/OUTPUT,k,MATRix,finish/soluallsel,allsolve!生产单个刚度矩阵,查看k.MATRix文件看出单元1,2各次刚度矩阵如下/AUX2FILE,'k','full','HBMAT,'k','txt',' ',ASCII,STIFF,NOFINISH!读取k.txt文件，得出初始结构刚度刚度矩阵如下：是个 9*9阶矩阵，含有非零参数23个，（以harwell-boeing文件格式生成的，只要该格式图区方式详见站内）结果如下：!为加入边界条件的刚度矩阵finish/soluSFBEAM,1,1,PRES,800/12nsel,s,,,1,3,2/GOD,all, , , , , ,ALLallsel,all!/debug,-3,,1,,,,,,!/OUTPUT,k4,MATRix,solve/AUX2FILE,'k','full',' 'HBMAT,'k','txt',' ',ASCII,STIFF,NOFINISH!读取k.txt文件，得出求解结构刚度刚度矩阵如下：。
apdl提取刚度矩阵
APDL提取刚度矩阵1. 引言在有限元分析中，刚度矩阵是一个重要的概念。
刚度矩阵描述了结构体系的刚度特性，是进行结构分析和求解的基础。
在ANSYS的APDL（ANSYS Parametric Design Language）中，我们可以通过一系列命令和操作来提取刚度矩阵。
本文将详细介绍如何在APDL中提取刚度矩阵，并探讨相关的技巧和注意事项。
2. APDL提取刚度矩阵的基本步骤提取刚度矩阵的基本步骤如下： 1. 创建几何模型：在APDL中，首先需要创建一个几何模型，可以使用命令或者图形界面进行建模。
2. 定义材料属性：根据实际情况，定义材料的弹性模量、泊松比等力学性质。
3. 定义单元类型：选择适当的单元类型，如梁单元、壳单元或体单元，并设置相应的单元属性。
4. 定义边界条件：定义结构的边界条件，如约束和加载。
5. 生成刚度矩阵：使用APDL提供的命令，生成刚度矩阵。
6. 导出刚度矩阵：将刚度矩阵导出到外部文件进行后续分析。
3. 创建几何模型在APDL中，可以使用命令行或者图形界面来创建几何模型。
命令行方式更加灵活，可以通过输入命令来精确控制几何体的生成。
图形界面方式则更加直观，通过鼠标操作可以创建各种几何体。
4. 定义材料属性在APDL中，可以使用*MAT命令来定义材料属性。
根据实际情况，需要输入材料的弹性模量、泊松比等力学性质。
这些参数将用于计算刚度矩阵。
5. 定义单元类型在APDL中，可以使用*ELEMENT命令来定义单元类型。
根据实际情况选择适当的单元类型，如梁单元、壳单元或体单元，并设置相应的单元属性。
不同类型的单元具有不同的自由度和刚度特性，选择合适的单元类型对于提取准确的刚度矩阵至关重要。
6. 定义边界条件在APDL中，可以使用*BOUNDARY命令来定义结构的边界条件。
边界条件包括约束和加载，用于模拟实际结构的受力和支撑情况。
正确定义边界条件对于提取准确的刚度矩阵非常重要。
hypermesh刚度矩阵提取
hypermesh刚度矩阵提取标题：Hypermesh刚度矩阵提取与应用引言：在工程领域中，结构分析是一个关键的步骤，用于评估和优化设计。
而刚度矩阵是结构分析中的重要概念之一，它能够描述结构在应力和应变下的响应。
Hypermesh是一款常用于有限元分析（FEA）的前处理软件，具有强大的刚度矩阵提取功能，本文将详细讨论Hypermesh中刚度矩阵提取的方法及其应用。
一、Hypermesh简介1.1 Hypermesh的概述Hypermesh是由Altair Engineering公司开发的一款有限元前处理软件，它具有友好的用户界面和丰富的功能，被广泛应用于航天航空、汽车、机械和建筑工程等领域。
1.2 Hypermesh的刚度矩阵提取功能Hypermesh具有强大的刚度矩阵提取功能，可以根据给定的几何模型、材料属性和约束条件，自动生成结构的刚度矩阵。
刚度矩阵是描述结构在外力作用下的刚性行为的矩阵，可以用于解析分析和优化设计。
二、刚度矩阵的基本概念2.1 刚度矩阵的定义刚度矩阵是一个N×N矩阵，其中N是结构中自由度的数目。
刚度矩阵的元素表示结构中各自由度之间的相互关系和相互作用。
2.2 刚度矩阵的性质刚度矩阵是对称的，正定的，且具有零在对角线的性质。
这些性质使得刚度矩阵可以表示结构系统的稳定性和刚性。
三、Hypermesh中刚度矩阵提取的方法3.1 几何建模在Hypermesh中，需要先进行几何建模，包括创建结构模型、定义节点和单元等。
可以选择直接绘制几何模型，或导入CAD模型进行后续编辑。
3.2 材料属性定义根据结构的材料特性，在Hypermesh中定义材料属性，比如弹性模量、泊松比等。
这些材料属性将用于刚度矩阵的计算。
3.3 约束条件的设定在刚度矩阵的提取过程中，需要设定结构的边界条件，即约束条件。
Hypermesh提供了丰富的边界条件选择，可以固定节点的位移或设定节点的约束力等。
3.4 刚度矩阵的提取通过上述步骤的设定，可以直接在Hypermesh中提取结构的刚度矩阵。
apdl提取刚度矩阵
APDL提取刚度矩阵1. 引言在工程领域中，刚度矩阵是一个非常重要的概念，它描述了结构体系中各个节点之间的刚度关系。
通过提取刚度矩阵，我们可以得到结构体系的整体刚度特性，从而进行结构分析和设计。
ANSYS Parametric Design Language (APDL) 是一种用于有限元分析的编程语言，可以方便地进行结构分析和后处理。
本文将介绍如何使用APDL来提取刚度矩阵。
2. 刚度矩阵的定义在开始讲解如何提取刚度矩阵之前，我们首先需要了解什么是刚度矩阵。
在有限元分析中，我们通常将连续介质离散化为有限个单元，并将每个单元看作是一个简单的力学系统。
对于每个单元，我们可以定义一个局部坐标系，并通过节点坐标和材料参数来计算该单元的刚度。
多个单元组成整个结构体系时，我们可以将每个单元的局部坐标系转换为全局坐标系，并组装得到整个结构体系的刚度矩阵。
3. APDL提取刚度矩阵的步骤下面将介绍使用APDL提取刚度矩阵的具体步骤：3.1 定义节点和单元首先，我们需要在APDL中定义节点和单元。
节点是结构体系中的重要组成部分，而单元则是连接节点的基本构件。
通过定义节点和单元，我们可以建立结构体系的几何形状，并为后续的刚度矩阵提取做准备。
3.2 定义材料参数和边界条件在APDL中，我们可以通过命令来定义材料参数和边界条件。
材料参数包括弹性模量、泊松比等与材料性质相关的参数；而边界条件包括固支、荷载等与结构约束相关的条件。
这些参数和条件将直接影响到刚度矩阵的计算。
3.3 划分单元在APDL中，我们可以使用命令来划分单元。
单元划分是将整个结构体系分割成多个小块的过程，每个小块称为一个单元。
通过合理地划分单元，我们可以更加精确地计算刚度矩阵。
3.4 计算局部刚度矩阵对于每个单元，我们需要计算其局部刚度矩阵。
局部刚度矩阵描述了单元内部的刚度特性，可以通过单元的材料参数和几何形状来计算得到。
在APDL中，我们可以使用命令来计算局部刚度矩阵。
abaqus 提取刚度矩阵
abaqus 提取刚度矩阵【原创实用版】目录1.Abaqus 简介2.刚度矩阵的定义和作用3.Abaqus 中提取刚度矩阵的方法4.注意事项和常见问题5.结论正文1.Abaqus 简介Abaqus 是一款广泛应用于工程领域的有限元分析软件，它可以帮助工程师在设计过程中预测和优化产品的性能。
在 Abaqus 中，刚度矩阵是一个关键的组成部分，它在分析结构的强度、刚度和稳定性等方面起着重要作用。
2.刚度矩阵的定义和作用刚度矩阵是一个描述结构刚度特性的矩阵，它表示了结构在受到外力作用时，各节点之间的相对位移与外力之间的比例关系。
刚度矩阵可以帮助我们了解结构的刚度分布和变形情况，从而指导我们进行合理的结构设计和优化。
3.Abaqus 中提取刚度矩阵的方法在 Abaqus 中，提取刚度矩阵的方法相对简单。
以下是具体操作步骤：（1）打开 Abaqus 软件，导入或创建一个有限元模型。
（2）在模型浏览器中，选择“Assembly”→“Instances”，在实例列表中选择需要提取刚度矩阵的模型。
（3）在“Instances”对话框中，选择“Visualize”→“Mesh”，查看模型的网格划分情况。
（4）在“Instances”对话框中，选择“K”→“Global”，打开“Global”对话框。
（5）在“Global”对话框中，选择“Stiffness”选项卡，勾选“Stiffness matrix”选项，然后单击“Change”按钮。
（6）在弹出的“Stiffness Matrix”对话框中，选择所需的刚度矩阵类型（如：总刚度矩阵、局部刚度矩阵等），并设置相应的输出选项。
（7）单击“OK”按钮，刚度矩阵将会被提取并显示在 Abaqus 的“Results”窗口中。
4.注意事项和常见问题在提取刚度矩阵时，需要注意以下几点：（1）确保模型已经正确划分网格，否则提取的刚度矩阵可能不准确。
（2）在选取刚度矩阵类型时，需要根据实际需求进行选择。
ansys提取刚度,质量矩阵
ansys提取刚度矩阵和质量矩阵单元刚度矩阵和单元质量矩阵提取：直接用/debug命令获得（建模－加约束－分析），具体参看下例finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,filename,out,, ! 将输出信息送到filename.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中---------------------华丽的分割线－－－－－－－－－－－－－－－整体刚度矩阵和质量矩阵提取有两种方法：（1）二次开发，用自空间法进行模态分析，从.full 文件中提取，具体见/forum/viewthread.php?tid=10860&extra=&page=1（2）利用超单元获得。
建模－网格划分－分析类型－分析选项－定义主自由度－加约束条件－分析具体参看下面一个小例题，自可明白。
/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加荷载/soluantype,7 !substructuring分析类型seopt,filename,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,filename,3 !列出整体刚度矩阵。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
============为什么折腾这个文档========我有一个计算线性动力学方程组的瞬态、谐响应和静力学的python程序，现希望开发一个将ANSYS组集好的总体矩阵导入该PYTHON程序中的接口。
该问题可分解为：[STEP1] [ANSYS]->[包含矩阵信息的文件][STEP2] [包含矩阵信息的文件]->[python通用数据对象][STEP3] [python通用数据对象]->[程序特定数据对象]->[进行计算]因此检索了一些帖子，基本上完成了这项工作，本文是对[STEP1]和[STEP2]的整理，并且利用[STEP3]对结果进行了验证============主要内容==================1，了解从ANSYS中提取总体矩阵和载荷向量的方法；2，了解提取出来的矩阵是怎样表示的；3，说明在Python中，如何读取这样的矩阵；4，构造一个简单的算例，说明整个【建模】-【提取】-【读取】过程及其正确性；=========站内检索综述====================检索词：提取矩阵得到21个结果，代表性的帖子有下面这9个：编号[1]标题：ansys中怎样提取质量，刚度，阻尼矩阵？地址：/forum-vi ... fromuid-159019.html 要点：pengweicai给出了一段网上最常见的提取代码,该程序以fortran 写成，可以利用.full文件以及一些列约定将ANSYS中的总体矩阵读入FORTRAN中。
编号[2]标题：如何得知HBMAT命令提取的质量、刚度矩阵对应的自由度？地址：/forum-vi ... fromuid-159019.html 要点：提出了使用HBMAT命令提取稀疏矩阵时常见的问题：我们如何知道提取出来的信息是怎么储存的呢？编号[3]标题：[分享]ANSYS中整体、单元刚度和质量矩阵的提取地址：/forum-vi ... fromuid-159019.html 要点：在该帖子的7楼，其实已经给出了帖子[2]中问题的解答，即HBMAT 中提取出来的矩阵是Harwell-Boeing格式的，并且给出了该格式的细节，可惜是英文的，没引起多少关注。
编号[4]标题：帮我看看提取的刚度与质量矩阵地址：/forum-vi ... fromuid-159019.html 要点：这个帖子所示的矩阵并非是使用HBMAT命令提出出来的，而应该是SELIST命令列举出来的未压缩的矩阵，后续楼层的回帖给了大家一个提示，即有可能提取出来的矩阵是引入了边界条件的（即删除了被约束的行和列的）。
编号[5] 标题：提取刚度矩阵的问题地址：/forum-vi ... fromuid-159019.html 要点：本帖作者的工作是基于单元刚度矩阵的，因此ANSYS中提取的单元刚度矩阵是否处于总体坐标系就成为问题。
该问题并非本文内容，但仍值得关注。
编号[6]标题：提取刚度矩阵丢失节点的问题地址：/forum-vi ... fromuid-159019.html 要点：帖子[5]作者的又一帖，在这里帖子[5]的问题得到了欧阳中华老师的回答。
编号[7]标题：提取刚度矩阵的ANSYS操作过程地址：/forum-vi ... fromuid-159019.html 要点：实际上这就是使用HBMAT从ANSYS中提取总体矩阵的全过程！只是还有一些细节待确定。
编号[8]标题：提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法地址：/forum-vi ... fromuid-159019.html 要点：给出了利用“不减缩的”子结构方法来得到总体矩阵的方法（这也是网络上常见的代码之一）编号[9]标题：质量矩阵、刚度矩阵如何提取？地址：/forum-vi ... fromuid-159019.html 要点：16443在5楼的回帖中给出了提取刚度矩阵的三种方法=======站外检索略述========================百度检索：提取矩阵比较好的帖子有：编号[10]来源:百度文库标题:怎样从ansys中提取单元刚度矩阵与质量矩阵地址:/view/3cf5e567f5335a8102d220d9.html 要点:这应该就是16443在帖子[9]中回复的内容了，全面的总结了在帖子[3,4,5,9]中涉及的问题。
编号[11]来源:中华钢结构标题:ansys刚度矩阵Harwell-Boeing格式的具体含义讨论地址:/forum/viewthread.php?tid=184007要点:如题，后续楼层给出了一些将矩阵读入ANSYS的APDL（好不容易读出来，又读进去干嘛呢……）编号[12]来源:simwe标题:关于ANSYS（质量、刚度、阻尼）矩阵Harwell-boeing格式数据的说明地址:/archiver/tid-924778.html要点:比[11]更透彻的HB格式说明!====================================================================1.从ANSYS中提取总体矩阵的方法================================= =============================================================1,用/DEBUG命令2,子结构法3,HBMAT详见帖子[10]PS.个人感觉HBMAT方法最靠谱，一是它的格式（Harwell-boeing）在很多场合都是通用的，二是BHMAT命令是文档化的、功能就是用来提取总体刚度矩阵的命令。
因此，相比于子结构法的剑走偏锋，/DEBUG命令的繁复，HBMAT命令方法更“标准”一些，因此在后文只关注此方法。
=======2.BH格式的矩阵是如何表示的===================================HBMAT命令并不是很复杂的命令，稍复杂的地方是采用该命令提取出来的矩阵是经过压缩的，称为Harwell-boeing格式，也叫Compressed Sparse Column格式。
其具体压缩和还原方式见帖子[3]（English）或[11][12](中文)=======3.如何在Python中读入BH格式的矩阵===============================上文说过，Harwell-boeing格式，也叫Compressed Sparse Column格式，而Python.scipy中就有这样的稀疏矩阵：1.class scipy.sparse.csc_matrix(arg1, shape=None, dtype=None, copy=False, dims=None, nzmax=None)可以通过HB文件中直接读取的行标指针，行标和数据创建，例如：1.>>> indptr = array([0,2,3,6])2.>>> indices = array([0,2,2,0,1,2])3.>>> data = array([1,2,3,4,5,6])4.>>> csc_matrix( (data,indices,indptr), shape=(3,3) ).todense()5.matrix([[1, 0, 4],6.[0, 0, 5],7.[2, 3, 6]])对应的HB文件应为(*号部分表示并非本例关注的数据)：1.Rainyboy Testing Matrix in BH format2.*** 4 6 63.RRA ** ** **4.(I14) (I14) (d25.15) (d25.15)5.06. 27. 38. 69.010.211.212.013.114.215.116.217.318.419.520.6由文件头可知，indptr的长度为4，因此0,2,3,6就是indotr的内容indices的长度为6，因此后续的0,2,2,0,1,2就是indices的内容data的长度为6，因此后续的1,2,3,4,5,6就是data的内容=======4.一个【建模】-【提取】-【读取】-【计算】的例子===============【建模APDL】1.FINISH2./CLEAR3./TITLE,CASE STUDY _BEAM _BEAM3 BY RAINYBOY4./PREP75./ESHAPE,1 !显示壳单元厚度6.!**********************7.!几何参数表8.!**********************9.*SET,L_HORI,0.1 !横梁的长度10.*SET,TA,0.005 !正方形截面的边长11.*SET,MESHCOUNT,2 !每段的分网数12.*SET,IZZ,TA*TA*TA*TA/12 !转动惯量13.*SET,IYY,TA*TA*TA*TA/12 !转动惯量14.!**********************15.!材料参数表16.!**********************17.*SET,MEX,1.78E11 !弹性模量18.*SET,MPRXY,0.3 !泊松比19.*SET,MDENS,7850 !密度20.!**********************21.!相关设置22.!**********************23.MP,EX,1,MEX !设置材料弹性模量24.MP,PRXY,1,MPRXY !设置材料泊松比25.MP,DENS,1,MDENS !设置材料密度26.BETAD,1E-5!BETA阻尼27.ET,1,BEAM3 !设置平面梁单元28.R,1,TA*TA,IZZ,TA !设置截面参数29.!DMPRAT,0.10000 !阻尼比30.!**********************31.!几何->分网32.!**********************33.TYPE,1 !指定分网类型34.MAT,1!指定材料类型35.REAL,1 !指定实参数36.K,1,0,0,0 !建立三个关键点37.K,2,L_HORI,0,038.L,1,2!建立几何体39.ALLSEL,ALL40.LESIZE,ALL,,,MESHCOUNT !设置线段分网数41.LMESH,ALL !分网42.!**********************43.!几何约束44.!**********************45.ALLSEL,ALL46.NSEL,S,LOC,X,0 !选择固定端节点47.D,ALL,ALL !设置为约束所有自由度48.ALLSEL,ALL49.NSEL,S,LOC,X,L_HORI50.F,ALL,FY,10 !力载荷51.ALLSEL,ALL52.save【提取APDL】1.!进行一次QRDAMP分析，以生成包含K、M、C和RHS的FULL文件2./SOLU3.ANTYPE,MODAL4.MODOPT,QRDAMP,2,25.SOLVE6.!将对应的矩阵提取到文件中7./AUX28.FILE,re,FULL9.HBMAT,K_RHS,txt,ASCII,,STIFF,YES10.HBMAT,M,txt,ASCII,,MASS,YES11.HBMAT,C,txt,ASCII,,DAMP,YES12.FINISH【ANSYS谐响应分析】(计算完毕后，手动把受力点的频响结果存在ree.txt中)1./SOLU2.ANTYPE,HARM3.HARFRQ,0,5024.NSUBST,2515.KBC,16.HROPT,FULL7.HROUT,OFF8.LUMPM,09.EQSLV, ,1e-008,10.SOLVE【读取&计算】(运行APP_HP_From_ANSYS.py之前在当前目录准备刚才ANSYS计算目录下的K_RHS.txt,M.txt,K.txt,ree.txt)1.# -*- coding: cp936 -*-2.#2011.4 重构动力学计算程序3.#使之可计算线性问题的瞬态和谐响应4.#使之可从文件读入矩阵5.#范雨****************6.7.#本文件包含：一个应用8.#从ANSYS中导出一个模型的M,K,C和F信息，进行谐响应分析。