有限元法与ANSYS技术-刚度矩阵
有限元分析ANSYS理论与应用（第4版）.例3.1_MATLAB理论求解
有限元分析ANSYS理论与应⽤（第4版）.例3.1_MATLAB理论求解相关帖⼦：有限元分析 ANSYS理论与应⽤(第4版).例3.1_ANSYS.Workbench求解题⽬描述：如图所⽰阳台桁架及其尺⼨。
假设所有杆件均为⽊质材料（道格拉斯红杉），弹性模量E=1.9×106lb/in2，且且⾯积为8in2。
确定每个接头的挠度，以及每个杆件的平均应⼒。
下⾯将MATLAB求解这个问题。
1、将问题结构离散为节点和单元：桁架的每个杆件作为单元，每个杆件的连接点作为节点。
因此，给定的桁架可以⽤5个节点和6个单元进⾏建模。
其中：1ft=12in.Element Node i Node j Length（in.）A（in2）E（lb/in2）θ（°)11236.08 1.9E+06022350.98 1.9E+0613533436.08 1.9E+06042436.08 1.9E+069052550.98 1.9E+064564536.08 1.9E+0602、计算各个单元的刚度矩阵，建⽴整体矩阵，边界条件处理，刚度⽅程及未知位移求解，求解⽀反⼒%% 定义输⼊条件A = 8; %杆件截⾯积E = 1.9E6; % 杆件材料弹性模量L1 = 36; % 1、3、4、6号杆件的长度L2 = 50.9; % 2、5号杆件的长度%% 计算单元刚度矩阵k1 = Bar2D2Node_Stiffness(E, A, L1, 0); %计算单元1刚度矩阵k2 = Bar2D2Node_Stiffness(E, A, L2, 135); %计算单元1刚度矩阵k3 = Bar2D2Node_Stiffness(E, A, L1, 0); %计算单元1刚度矩阵k4 = Bar2D2Node_Stiffness(E, A, L1, 90); %计算单元1刚度矩阵k5 = Bar2D2Node_Stiffness(E, A, L2, 45); %计算单元1刚度矩阵k6 = Bar2D2Node_Stiffness(E, A, L1, 0); %计算单元1刚度矩阵%% 建⽴整体刚度矩阵kk = zeros(10, 10);kk = Bar2D2Node_Assembly(kk, k1, 1, 2);kk = Bar2D2Node_Assembly(kk, k2, 2, 3);kk = Bar2D2Node_Assembly(kk, k3, 3, 4);kk = Bar2D2Node_Assembly(kk, k4, 2, 4);kk = Bar2D2Node_Assembly(kk, k5, 2, 5);kk = Bar2D2Node_Assembly(kk, k6, 4, 5) % 输出整体刚度矩阵%% 边界条件处理k = kk([3478910], [3478910]);%添加位移约束。
ansys提出刚度和质量矩阵
一、连续钢梁的刚度和质最矩阵弹性模量：EX=2.0x 10n Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。
AN SYS命令流：/FILNAM.BEAM MASS AND KNIFF ANALYSIS/TITLE, BEAM MASS AND KNIFF ANALYSIS/PREP7EI1BEAM3MBEX470E11MRNUXYX0.3MRDENS 丄7800SECTYPE,1,BEAM」,,3SECO 阡SET.CENTSECDATAO150・：15,030・02O02,0・(HO0O0 K10O0K210O0一、框架的刚度和质呈:矩阵的提取模型采用右图的集中质屋模型，Ml=2762Kg ,M2=2760Kg, M3=2300Kg,层间刚度分别为：248500、192100. 152200命令流如下:finish/clear/prep7 et4z combinl4keyopt,1,2,1 !—维弹簧单元（Ux平动）et2mass21 keyopt,2,3/2 !3・D mass没有转动惯最0,1041-1）*3,0 *enddo!质最常数r z2,2762r,3,276054,2300type,2*do,i,2,4,lrealje,i*enddo!层间刚度^12,248500G1*******G14J52200 type」•doj,134 realjl+i eJJ+1 *enddo 哟束nselsnode,24JL d，all,uz，0d,all,uy,0 allsel d,l,all,0 /solu antype z7 seoptsubmatB m,all z all solve selistsubmat.B使用该命令流町以得到结构的刚度矩阵和质最矩阵结果: 刚度矩阵：440600. 00 -192100. 000. 0000000 质量矩阵：2762. 00000. 00000000. 0000000•192100.00344300. 00-152200.000. 0000000-152200. 00152200. 000.00000002760. 00000.00000000.00000000.00000002300. 0000imIT12mik-ki附件1：（运行ANSYS命令的输出结果，最后面是刚度和质最矩阵）PRINT CONTENTS OF SUPERELEMENT submatPRINT OPTION = 3HEADER =8 3 2 3 34 0 4 3 01 1 0 0 10 0 1 12 11550 1101 166 103 submat 1078 0922 928 946 953 976982 1027 1072 0 934940 106004406 0 0 30 0419 0 0 0 00 0 0 0HEADER SUMMARY:NUMBER OF ROWS = 3NUMBER OF MATRICES = 2NUMBER OF EDGE PLOT LINES = 3NUMBER OF DEGREES OF FREEDOM PER NODE = 3NUMBER OE DEGREE OF FREEDOM INDICES = 4NUMBER OF NODES = 3NUMBER OF LOAD VECTORS = 1NUMBER OE TRANSFORMATIONS = 0BASE FILE NAME二submatMAXIMUM STIFFNESS二0. 44060E+06DEGREES OF FREEDOM PER NODE =DEGREES OF FREEDOM PER NODE LABELS =UX UY UZDEGREES OF FREEDOM AS GENERATED =4 7 10REORDERED DOF POSITIONSORIGINAL DOF ORDER =DEGREE OF FREEDOM INDICES =1 2 3 4 TITLE =NODES =23 4NODE X Y THXZ2 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 00004 0. 0000 9. 0000 0. 0000EDGE PLOT DATA 二EDGE XI Y1Z21 0. 0000 0. 0000 0. 00002 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 0000GLOBAL DOF SET 二33 65 97GLOBAL DOF SET NODES, LABELS2 UX3 UX4 7 10MASS INFORMATION:Z THXY THYZ 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000Z1 X2 Y20. 0000 0. 0000 3. 0000 0. 0000 0. 0000 6. 0000 0. 0000 0. 0000 9. 00004 UXTOTAL MASS = 7822.0CENTROID (X,Y,Z) = 0. 0000 5. 8228 MOMENT OE INTERTIA ABOUT ORIGIN:IXX = 0.31052E+06 IYY = 0. 0000IXY = 0.0000 IYZ = 0.0000 0. 00001ZZ = 0.31052E+06 IZX = 0. 0000RON 1 MATRIX440600. 00 -192100. 00 0. 0000000ROW 1 MATRIX 22762. 0000 0. 0000000 0. 0000000ROW 2 MATRIX 1-192100. 00 344300. 00 -152200. 00ROW 2 MATRIX 20. 0000000 2760. 0000 0. 0000000ROW 3 MATRIX 10. 0000000 -152200. 00 152200. 00ROW 3 MATRIX 20. 0000000 0. 0000000 2300. 0000LOAD VECTOR 10. 0000000 0. 0000000 0. 0000000。
ansys有限元分析原理
ansys有限元分析原理
ANSYS有限元分析原理是一种数值分析方法，广泛应用于工
程领域。
其核心思想是将复杂的物体或结构划分为许多小的几何单元，称为有限元。
每个有限元由节点和单元组成，其中节点为有限元的角点或自由度，而单元则定义了节点之间的连接关系。
在有限元分析中，首先需要建立物体或结构的有限元模型。
这涉及到将物体或结构的几何形状进行离散化，并定义节点和单元。
通常情况下，物体或结构的复杂性越高，所需要的有限元模型就越精细，节点和单元数量也就越多。
接下来，需要定义物体或结构的边界条件和加载条件。
边界条件包括约束条件和固定边界条件，用于限制节点的位移和旋转。
加载条件包括力、热源、压力等外部作用力，用于模拟实际工程中的加载情况。
有限元分析通过求解有限元模型的全局刚度矩阵和加载向量来计算系统的响应。
根据有限元模型的节点和单元之间的连接关系，全局刚度矩阵可以通过将每个单元的刚度矩阵组合而成。
加载向量则是由加载条件决定的。
最后，通过求解线性方程组，即全局刚度矩阵乘以位移向量等于加载向量的形式，可以得到有限元分析的结果。
位移向量记录了每个节点在加载后的位移情况，从而可以计算各个节点的应力、应变等响应参数。
总之，ANSYS有限元分析原理是将复杂的物体或结构划分为小的几何单元，通过离散化、边界条件和加载条件的定义，以及全局刚度矩阵和加载向量的计算，求解线性方程组，最终得到系统的响应结果。
这个方法在解决工程问题中具有广泛的应用。
ANSYS-STRUCTURAL基础培训-
实体几何模型载荷
培训手册
CFD Analysis with ANSYS/FLOTRAN
优点 改变网格不影响载荷 涉及到的加载实体少 缺点 生成的单元在当前激活的单元座标下，节点为总 体直角座标，因此实体与有限元模型可能有不同 座标系统和载荷方向 实体载荷在凝集分析中不方便，因载荷加在主自 由度上施加关键点约束较繁锁
简例〔续〕
培训手册
CFD Analysis with ANSYS/FLOTRAN
下面以小变形弹性静力问题为例，加以具体介绍。 几何方程：eij=1/2(ui,j+uj,i) 物理方程：sij=aijklekl 平衡方程：sij,j+fi=0 边界条件：
位移已知边界条件 ui=ui 〔在边界Гu上位移已知〕 外力已知边界条件 sij,j+pi=0〔在边界Гp上外力已知〕
2. 弹簧阻尼单元，包括COMBIN系列： COMBIN7,COMBIN14,COMBIN37,COMBIN40等。
3. 质量元，MASS21。
CFD Analysis with ANSYS/FLOTRAN
ANSYS单元分类
初级培训
4 . 梁单元，分为二维和三维单元，弹性和塑性；可定 义各种截面。主要包括： BEAM3，BEAM4，BEAM23，BEAM24，BEAM44， BEAM188，BEAM189。
〔1-1〕
简例〔续〕
培训手册
CFD Analysis with ANSYS/FLOTRAN
即 {F} 1=[K] 1{δ}1 在〔1-1〕中令U11=1，V11=U21=V21=0 则F1x1=K11 ，F1y1=K21 ，F1x2=K31 ，F1y2=K41
上式说明，当节点1沿X方向产生一单位位移， 而单元1的其余节点位移为零时，各节点施于单元1 上的力将组成一平衡力系，表示单元1抵抗位移 U11的刚度。
ansys提出刚度和质量矩阵
2762.0000 0.0000000 0.0000000
附件 1：（运行 ANSYS 命令的输出结果，最后面是刚度和质量矩阵） PRINT CONTENTS OF SUPERELEMENT submat PRINT OPTION = 3 HEADER = 8 4 1 0 0 at 922 982 940 0 419 0 3 0 1 0 1101 1078 928 1027 106004406 0 0 2 4 0 1 166 0 946 1072 0 0 0 3 3 0 12 103 953 0 0 0 0 3 0 1 1155 subm 976 934 3 0 0
一、连续钢梁的刚度和质量矩阵 弹性模量：EX=2.0× 1011 Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。 ANSYS 命令流： /FILNAM,BEAM MASS AND KNIFF ANALYSIS /TITLE, BEAM MASS AND KNIFF ANALYSIS /PREP7 ET,1,BEAM3 MP,EX,1,2.0E11 MP,NUXY,1,0.3 MP,DENS,1,7800 SECTYPE,1,BEAM,I,,3 SECOFFSET,CENT SECDATA,0.15,0.15,0.3,0.02,0.02,0.01,0,0,0,0 K,1,0,0,0 K,2,10,0,0
GLOBAL DOF SET = 33 65 97 GLOBAL DOF SET NODES, LABELS = 2 UX 3 UX MASS INFORMATION: TOTAL MASS = 7822.0 CENTROID (X,Y,Z) = 0.0000 5.8228 MOMENT OF INTERTIA ABOUT ORIGIN: IXX = 0.31052E+06 IYY = 0.0000 IXY = 0.0000 IYZ = 0.0000 ROW 1 MATRIX 1
有限元方法与ANSYS应用第7讲有限元的基础理论与方法 有限元案例分析 动力分析
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析
完全法谐响应分析----加载并求解
步骤：
2 定义分析类型和分析选项
· 选项： Mass Matrix Formulation[LUMPM]
此选项用于指定是采用缺省的分布质量矩阵（取决 于单元类型）还是集中质量矩阵。建议在大多数应用中 采用缺省的分布质量矩阵。但对于某些包含“薄膜”结 构的问题，集中质量近似矩阵经常能产生较好的结果。
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析
三种求解方法----完全法
优点：
· 用单一处理过程计算出所有的位移和应 力。 · 允许定义各种类型的载荷：节点力、外 加的（非零）位移、单元载荷（压力和温 度）。 · 允许在实体模型上定义载荷。
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析
步骤：
9 观察结果
2.派生数据 · 节点和单元应力 · 节点和单元应变 · 单元力 · 节点反作用力，等等。
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析
缩减法谐响应分析
缩减法的分析过程由五个主要步骤组成： 1.建模； 2.加载并求得缩减解； 3.观察缩减解结果； 4.扩展解（扩展过程）； 5.观察已扩展的解结果。 在这些步骤中，第1步的工作与完全法的相同。
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析
任何持续的周期载荷作用在结构系统中 所产生的持续性周期响应（谐响应）。
有限元法分析的基本理论与方法
★ 有限元案例分析—谐响应分析 谐响应分析寻求对已知幅值载荷的
响应振幅。 该载荷随时间以已知频率呈正弦形
式变化。
ansys单元刚度矩阵
ansys单元刚度矩阵
ANSYS单元刚度矩阵是一种表示有限元模型中单元刚度的矩阵形式。
在有限元分析中，刚度矩阵是一个关键的概念，它描述了单元的刚度和其对整个结构的贡献。
刚度矩阵的形式是一个对称矩阵，其中对角线上的元素表示单元的刚度，而非对角线上的元素表示单元在不同方向上的耦合刚度。
ANSYS使用有限元法对结构进行建模和分析。
在建立有限元模型时，将结构分割成许多小的单元，并计算每个单元的刚度矩阵。
这些单元刚度矩阵在组装时组合成整个结构的总刚度矩阵。
然后，通过求解总刚度矩阵的特征值和特征向量，可以确定结构的自由振动频率和模式，并且可以计算结构的应力和应变分布。
ANSYS提供了各种单元类型和材料模型，以便用户可以构建适合其特定应用的模型。
对于不同类型的单元，ANSYS使用不同的数学公式计算单元刚度矩阵。
用户还可以输入自定义的单元刚度矩阵，以便在模型中使用。
总之，ANSYS单元刚度矩阵是有限元模型中的重要概念，可以帮助工程师分析和设计各种结构和设备。
- 1 -。
ansys质量矩阵刚度矩阵提取说课材料
a n s y s质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow:-------------------/aux2fileaux2,file,ematform,longdump,all-------------------2.可以使用/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,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵3.其原理很简单，即使用ansys的超单元即可解决问题。
ansys质量矩阵刚度矩阵提取教学提纲
a n s y s质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow:-------------------/aux2fileaux2,file,ematform,longdump,all-------------------2.可以使用/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,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵3.其原理很简单，即使用ansys的超单元即可解决问题。
有限元单元刚度矩阵计算方法
有限元单元刚度矩阵计算方法
有限元单元刚度矩阵是有限元分析中的一个关键组成部分，它描述了结构中每个元素在承受载荷时的刚度响应。
以下是一个计算有限元单元刚度矩阵的基本步骤：
1. 确定元素类型和参数：首先需要确定所使用的元素类型（例如，杆、梁、板、壳等），以及这些元素的参数，如横截面面积、惯性矩、厚度等。
2. 建立局部坐标系：为每个元素建立一个局部坐标系。
在局部坐标系中，可以方便地描述元素内部的应力和应变。
3. 计算应变矩阵：根据有限元理论，计算元素两端的节点坐标差值，并由此得到应变矩阵。
4. 计算应力矩阵：根据材料的物理性质和胡克定律（Hooke's law），将应变矩阵转换为应力矩阵。
5. 形成刚度矩阵：将应力矩阵乘以相应的刚度系数，得到该元素的刚度矩阵。
6. 组装整体刚度矩阵：将所有元素的局部刚度矩阵组合起来，形成整体结构的刚度矩阵。
7. 施加边界条件和载荷：根据实际问题的边界条件和载荷，对整体刚度矩阵进行修正。
8. 求解线性方程组：通过求解修正后的线性方程组，得到结构中每个节点的位移。
以上步骤仅为有限元分析中的一种基本方法，实际应用中可能还需要考虑更多的因素，如非线性行为、材料失效等。
此外，有限元分析软件（如ANSYS、SolidWorks等）通常已经内置了这些计算过程，用户可以直接调用相应的功能进行有限元分析，而无需手动编写代码。
ansys提取单元刚度矩阵
ansys提取单元刚度矩阵
在使用ANSYS进行有限元分析时，我们通常需要提取单元的刚度矩阵。
刚度矩阵是描述单元刚度的一个重要参数，它可以用来计算单元的应力、应变和位移等参数。
提取单元刚度矩阵的方法如下：
1. 首先，在ANSYS中建立模型并进行网格划分。
确定需要提取刚度矩阵的单元，例如在ANSYS中选择PLANE42单元。
2. 在ANSYS命令窗口中输入“/POST1”命令，进入后处理模式。
3. 输入“ET,LIST”命令，显示所有单元类型的列表。
根据需要选择需要提取刚度矩阵的单元类型。
4. 输入“SET,LIST”命令，选择需要提取刚度矩阵的单元集合。
5. 输入“MAT,LIST”命令，选择需要提取刚度矩阵的材料属性。
6. 输入“SGLSTIF,ALL”命令，提取所有单元的刚度矩阵。
也可以使用“SGLSTIF,LIST”命令，提取指定单元的刚度矩阵。
7. 输入“PRRSOL,U”命令，显示单元的未知位移。
8. 输入“PRRSOL,F”命令，显示单元的未知力。
通过以上步骤，我们就可以成功提取单元的刚度矩阵，并使用它来计算单元的应力、应变和位移等重要参数。
- 1 -。
ANSYS动力学分析中的矩阵
ANSYS动力学分析中的矩阵在ANSYS动力学分析中，矩阵是非常重要的工具。
矩阵作为数据结构的一种形式，可以方便地表示和处理线性方程组和向量的关系。
下面我们将介绍一些在ANSYS动力学分析中常用的矩阵。
1. 刚度矩阵（Stiffness Matrix）：刚度矩阵表示了结构中各个节点之间的刚度关系。
它的行和列对应于结构的自由度，每个元素表示了对应节点之间的刚度。
在动力学分析中，刚度矩阵被用于描述结构在受到外力作用时的刚度响应。
2. 质量矩阵（Mass Matrix）：质量矩阵表示了结构中各个节点之间的质量关系。
它的行和列对应于结构的自由度，每个元素表示了对应节点的质量。
在动力学分析中，质量矩阵被用于描述结构的惯性特性，以及结构在外力作用下的加速度响应。
3. 全局刚度矩阵（Global Stiffness Matrix）：当结构由多个单元组成时，每个单元的刚度矩阵可以通过组装操作得到全局刚度矩阵。
全局刚度矩阵描述了整个结构的刚度关系，可以用于求解结构在受到外力作用时的位移响应。
4. 结构方程矩阵（Structural Equation Matrix）：结构方程矩阵是动力学分析中的一种常用的矩阵形式，它将结构的质量矩阵、刚度矩阵和外力向量组合在一起，用于求解结构在受到外力作用时的位移响应。
结构方程矩阵可以通过求解线性方程组得到结构的位移响应。
5. 阻尼矩阵（Damping Matrix）：阻尼矩阵描述了结构在受到外力作用时的阻尼特性。
阻尼矩阵在动力学分析中被用来模拟结构的阻尼效应，以及计算结构在受到外力作用时的振动响应。
6. 模态矩阵（Modal Matrix）：模态矩阵描述了结构在自由振动状态下的振型特性。
模态矩阵包含了结构的特征向量，它可以用来计算结构的特征值和特征向量，以及结构在外力作用下的振动响应。
7. 自由度矩阵（Degree of Freedom Matrix）：自由度矩阵用于表示结构的自由度信息。
ansys提取质量和刚度矩阵
ansys提取质量和刚度矩阵
在ANSYS中，可以使用命令`*MATRIX, MASS, STIFF`来提取质量和刚度矩阵。
要提取质量矩阵，可以使用以下命令：
```
*MATRIX, MASS
/mesh_needed ! 控制只提取有限元网格中的质量矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后，ANSYS将会在指定的输出文件中生成质量矩阵。
要提取刚度矩阵，可以使用以下命令：
```
*MATRIX, STIFF
/mesh_needed ! 控制只提取有限元网格中的刚度矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后，ANSYS将会在指定的输出文件中生成刚度矩阵。
请注意，以上命令需要在ANSYS命令窗口中执行，并且需要
已经建立合适的有限元模型。
您还需要根据需要修改`/output`命令以指定输出文件的目录和文件名。
ansys提取刚度矩阵案例（命令流及矩阵文件）
ansys提取刚度矩阵案例（命令流及矩阵⽂件）ansys提取刚度矩阵的三种⽅法注：本案例借鉴了王新敏⽼师编著的《ansys⼯程结构数值分析》部分内容结构刚度矩阵K(1)单元刚阵的提取-----详见page356/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLU/OUTPUT,ELEMSTIFF,TXT/DEBUG,-1,,,1SOLVE/OUTPUTFINISH⽤任⼀⽂本编辑器打开ELEMSTIFF.TXT⽂件可得到单元刚度矩阵ELEMSTIFF.TXT(2)⽤HBMAT提取原始刚度矩阵和节点载荷FINISH/CLEAR/FILNAME，HBFILE/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,HBFILE,FULLHBMAT,HBFILE,TXT,ASCII,STIFF,YESFINISH⽤⽂本编辑器打开HBFILE.TXT可看到⽤Harwell-Boeing格式记录的⽂件（仅含⾮零项的上三⾓阵，刚度矩阵为对称矩阵）。
矩阵有12⾏12列，33个⾮零元素HBFILE.TXT（3）⽤HBMAT提取结构刚度矩阵与（2）相同，但施加约束条件即可。
所⽣成的HBMAT.TXT的前5⾏为/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUWRFULL,1SOLVEFILE,⽂件名缺省,FULLHBMAT,⽂件名缺省,TXT,ASCII,STIFF,YESFINISH其意义同上，但数值有变化。
Ansysworkbench17.0提取刚度矩阵的方法
设置文件名称和刚度矩阵类型刚度质量阻尼等nsel3756
Ansysworkbench17.0提 取 刚 度 矩 阵 的 方 法
第一步：导入模型，划分网格 第二步：定义remote point，位置设置在想要提取刚度矩阵的位置。
第三步：添加remote displacement约束，给第二步定义的点添加约束。（本例为全约束）
第四步：在static structural处插入command，内容如下： /solu antype,7 !substructuring分析类型 seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等) nsel,,,,3756 ！此处节点数（3756）由mesh中statistics 查到，nodes 值为3755，则remote point为3756、3757……m,all,all !选择所有节点 allsel !定义所有节点自由度为主自由度 solve !求解 selist,matname,3 !列出整体刚度矩阵 finish 第五步：solve 第六步：在solution中solution information查看结果（结果按照行分布）
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,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵二、整体刚度矩阵的提取（有三种方法：用户程序法、超单元法、HBMAT命令法）1、用户程序法：需要二次开发（略）2、超单元法/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵3、HBMAT命令法提取整体矩阵命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
k
e
BiT BjT
D
Bi
Bj
Bm
t
kii k ji
kij k jj
kim k jm
(3-35)
BmT
kmi kmj kmm
其中 krs Br T DBs t
Et
4 1 2
bcrbrbss1122cbrrccss
br cs
1
2
cr bs
cr cs
1
2
brb s
Kij
Kim
Kin
K
K j1
K ji
K jj
K jm
K jn
考虑到[k]扩充以后，除了对应的i, j, m 双行和双列 上的九个子矩阵之外，其余元素均为零，故（3-33）式
中的单元位移列阵{}e2n×1 便可用整体的位移列阵 {}2n×1 来替代。这样，（3-33）式可改写为
k
Re
2n2n
2n1
2n1
把上式对N个单元进行求和叠加，得
N e1
e
ui
vi
u j
v j
um
T
vm
且假设单元内各点的虚位移为{f *}，并具有与真实位移 相同的位移模式。
故有
f N e
(c)
参照（3-13）式，单元内的虚应变{ *}为
B e
(d)
于是，作用在单元体上的外力在虚位移上所做的功可写
为
({ }e )T Re
(f)
而单元内的应力在虚应变上所做的功为
码的排序一致。各单元的节点力列阵经过这样的扩充之
后就可以进行相加，把全部单元的节点力列阵叠加在一
起，便可得到 (l)式所表示的弹性体的载荷列阵，即
N
R
R e R1T
R2T
RnT T
(p)
e1
这是由于相邻单元公共边内力引起的等效节点力，在叠
加过程中必然会全部相互抵消，所以只剩下载荷所引起
Re BT DBtdxdye
记
ke BTRe ke e
(3-33)
上式就是表征单元的节点力和节点位移之间关系的刚
度方程，[k]e就是单元刚度矩阵。如果单元的材料是均质的， 那么矩阵[D] 中的元素就是常量，并且对于三角形常应变单 元，[B]矩阵中的元素也是常量。当单元的厚度也是常量时， 因，所以（3-28）式可以简化为
第四节
刚度矩阵
一. 单元刚度矩阵
为了推导单元的节点力和节点位移之间的关系，可 应用虚位移原理对图3-2中的单元e进行分析。单元e是在 等效节点力的作用下处于平衡的，而这种节点力可采用 列阵表示为
R e RiT
R
T j
RmT T Ui
Vi
Uj
Vj
T
Um Vm
(a)
假设在单元e中发生有虚位移，则相应的三个节点i、j、 m 的虚位移为
k
N
Re
e1
(r)
上式左边就是弹性体所有单元刚度矩阵的总和，
称为弹性体的整体刚度矩阵（或简称为总刚），记为
[K]。注意到（3-28）式，有
N
N
K k BT DBtdxdy (3-38)
e1
e1
若写成分块矩阵的形式，则
K11 K1i K1 j K1m K1n
Ki1
Kii
j
(
R
e j
)
T
m
(Rme )T
n
T
(n)
其中，子矩阵
Ri
U
e i
Vie T
(i, j, m 轮换)
(o)
是单元节点i上的等效节点力。
(n)式中的省略号处的元素均为零，矩阵号上面的i, j, m 表示在分块矩阵意义下Ri 所占的列的位置。此处 假定了i, j, m 的次序也是从小到大排列的、并且与节点 号
T tdxdy
(g)
这里我们假定单元的厚度t为常量。把（d ）式及（3-16）式 代入上式，并将提到积分号的前面，则有
({ }e )T BT DBetdxdy
根据虚位移原理，由（f）和（h）式可得到单元的虚功方程， 即
({ }e )T Re ({ }e )T BT DBe tdxdy
注意到虚位移是任意的，所以等式两边与相乘的项应该相等， 即得
cr cs 21 brbs
bs
( r = i、j、m；s = i、j、m )
(3-37)
二 整体刚度矩阵
讨论了单元的力学特性之后，就可转入结构的整体 分析。假设弹性体被划分为N个单元和n个节点，对每个 单元按前述方法进行分析计算，便可得到N组形如（3-33） 式的方程。将这些方程集合起来，就可得到表征整个弹 性体的平衡关系式。为此，我们先引入整个弹性体的节
点位移列阵 {}2n×1 ，它是由各节点位移按节点号码以从
小到大的顺序排列组成，即
2n1
T 1
T 2
T T n
(j)
其中子矩阵
i ui vi T
是节点i的位移分量。
(i =1,2, …, n ) (k)
继而再引入整个弹性体的载荷列阵{R}2n×1 ，它是移置 到节点上的等效节点载荷依节点号码从小到大的顺序排列
( r = i、j、m；s = i、j、m )
(3-36)
对于平面应变问题，只要将上式中的E、分别换成E / 1- 2 和 / 1- 即
可。于是
krs
E1 t 41 1 2
1
br bs
1 2
21
cr cs
1 2
crbs 21 br
cs
1
br cs
1 2
21
cr
1 2
的等效节点力。
同样，将（3-35）式的六阶方阵[k]加以扩充，使之
成为2n阶的方阵
1
i
j
m
n
1
kii
kij
kim
i
k
2n2n
k ji
k jj
k jm
j
(q)
kmi
kmj
kmm
m
n
不难看出，（3-35）式中的2×2阶子矩阵[ki j ]将处于上 式中的第i双行、第j双列中。
[k]e =[B]T [D][B]t
(3-34)
与前面讨论过的情况类似，单元刚度矩阵[k]中任一 列的元素分别等于该单元的某个节点沿坐标方向发生单 位位移时，在各节点上所引起的节点力。单元的刚度取 决于单元的大小、方向和弹性常数，而与单元的位置无 关，即不随单元或坐标轴的平行移动而改变。
将（3-30）式写成分块形式，即可得到平面应力问 题中三角形单元的刚度矩阵
组成，即
R
2n1
R1T
R2T
RnT T
(l)
其中子矩阵
Ri Xi
Yi
T
N
U
e i
N
T
Vi
e
e1
e1
(i =1,2, …, n )
(m)
是节点i上的等效节点载荷。
现将各单元的节点力列阵{R}e6×1 加以扩充，使之 成为2n×1阶列阵
R
e 2n1
1
i
( Rie )T