Ansys中单元刚度矩阵、总刚矩阵的提取

合集下载

有限元法与ANSYS技术-刚度矩阵

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)式所表示的弹性体的载荷列阵，即

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中单元刚度矩阵、总刚矩阵的提取一、单元刚度矩阵的提取/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。

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提出刚度和质量矩阵

filnambeammasskniffanalysistitlebeammasskniffanalysisprep7et1beam3mpex120e11mpnuxy103mpdens17800sectype1beami3secoffsetcentsecdata015015030020020010000k21000一框架的刚度和质量矩阵的提取模型采用右图的集中质量模型m12762kgm22760kgm32300kg层间刚度分别为
EDGE PLOT DATA = EDGE Z2 1 0.0000 2 0.0000 3 0.0000 0.0000 0.0000 3.0000 6.0000 0.0000 0.0000 0.0000 0.0000 6.0000 9.0000 0.0000 0.0000 0.0000 0.0000 3.0000 X1 Y1 Z1 X2 Y2
r,12,248500 r,13,192100 r,14,152200 type,1 *do,i,1,3,1 real,11+i e,i,i+1 *enddo !约束 nsel,s,node,,2,4,1 d,all,uz,0 d,all,uy,0 allsel d,1,all,0 /solu antype,7 seopt,submat,3 m,all,all solve selist,submat,3 /
-192100.00 2 0.0000000 1 344300.00 2 2760.0000 1 -152200.00 2 0.0000000
0.0000000
0.0000000
-152200.00
0.0000000
152200.00
2300.0000
0.0000000
0.0000000

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中，我们可以使用命令来计算局部刚度矩阵。

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命令来得到。

ansys质量矩阵刚度矩阵提取教学提纲

ANSYS中几种提取模态方法

ANSYS学说模态提取方法在ANSYS 中有以下几种提取模态的方法：– (1) Block Lanczos 法– (2) 子空间法– (3) PowerDynamics 法– (4) 减缩法– (5) 不对称法– (6) 阻尼法使用何种模态提取方法主要取决于模型大小（相对于计算机的计算能力而言）和具体的应用场合。

(1) Block Lanczos 法Block Lanczos 法可以在大多数场合中使用：- 是一种功能强大的方法，当提取中型到大型模型（50,000 ~ 100,000 个自由度）的大量振型时（40+），这种方法很有效；- 经常应用在具有实体单元或壳单元的模型中；- 在具有或没有初始截断点时同样有效。

（允许提取高于某个给定频率的振型）；- 可以很好地处理刚体振型；- 需要较高的内存。

(2) 子空间法子空间法比较适合于提取类似中型到大型模型的较少的振型（<40）- 需要相对较少的内存；- 实体单元和壳单元应当具有较好的单元形状，要对任何关于单元形状的警告信息予以注意；- 在具有刚体振型时可能会出现收敛问题；- 建议在具有约束方程时不要用此方法。

(3) PowerDynamics 法PowerDynamics 法适用于提取很大的模型（100.000个自由度以上）的较少振型（< 20）。

这种方法明显比Block Lanczos 法或子空间法快，但是：- 需要很大的内存；- 当单元形状不好或出现病态矩阵时，用这种方法可能不收敛；- 建议只将这种方法作为对大模型的一种备用方法。

注: PowerDynamics 方法- 子空间技术使用Power 求解器(PCG) 和一致质量矩阵；- 不执行Sturm 序列检查(对于遗漏模态); 它可能影响多个重复频率的模型；- 一个包含刚体模态的模型, 如果你使用PowerDynamics 方法，必须执行RIGID 命令(或者在分析设置对话框中指定RIGID 设置)。

2016新编ANSYS中导出刚度矩阵

[转贴]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文件，可以看到按单元写出的质量、刚度等矩阵干部教育培训工作总结[干部教育培训工作总结] 年干部教育培训工作，在县委的正确领导下，根据市委组织部提出的任务和要求，结合我县实际，以兴起学习贯彻“三个代表”重要思想新高潮为重点，全面启动“大教育、大培训”工作，取得了一定的成效，干部教育培训工作总结。

现总结报告如下：一、基本情况全县共有干部**人，其中中共党员**人，大学本科以上学历**人，大专学历**人，中专学历**人，高中及以下学历**人。

**年，以县委党校、县行政学校为主阵地，举办各类培训**期，培训在职干部**人，占在职干部总数的**.*%，培训农村党员、干部**人，其中：举办科级领导干部轮训班*期，培训**人；举办科级领导干部“三个代表”重要思想专题学习班*期，培训**人；举办科级以下公务员培训班*期，培训**人；举办企业经营管理者培训班*期，培训**人；举办专业技术人员培训班*期，培训**人；举办非中共党员干部培训班*期，培训**人；举办理论骨干培训班*期，培训**人；举办妇女干部培训班*期，培训**人；举办基层团干培训班*期，培训**人；举办农村党支部书记、村主任培训班各*期，培训**人，达到了每年培训在职干部五分之一的要求，超额完成了培训任务。

ansys常见问题解答

ANSYS常见问题及应用技巧本篇开始讲述ANSYS在使用过程中常见的问题和在使用时一些常用的使用技巧，对与初学者来说，理解和弄清楚这些问题的根源和掌握这些使用技巧，能够更好的理解ANSYS这个软件本身。

1.ANSYS中的等效应力是什么物理含义？1）ANSYS中等效应力最大应力s1有什么区别，平常讨论应力分布，应该用等效应力还是最大应力s1呢？2）计算等效应力时是否需要输入等效泊松比呢？3）在实际的应用中，例如在讨论平板上的圆孔应力集中的应力分布问题时，应该用等效应力来描述应力集中的现象，还是采用主应力s1来反应集中的程度呢？还是采用一个单方向的sx来说明问题呢？答：1）这个等效应力应该就是弹塑性力学里的VonMises应力，他主要考察的是材料在各个方向上的应力差值，因为在实验室里获得材料强度都是单向载荷作用下的强度（当然现在也有三轴应力实验仪），所以有时候材料所受的单向载荷可能很大，但并没有造成破坏，这是就是看他的等效应力，具体计算公式是: σ等效=sqrt{0.5[(σ1-σ2)^2+(σ2-σ3)^2+(σ3-σ1)^2]}2）等效应力是三项主应力的组合如s，int即为max(si-sj),si,sj为三项主向应力。

i，j=1，2，3 i≠j即tresca型s,eqv为sqrt(0.5*∑(si-sj)**2)，i，j=1，2，3 i≠j即mises型3）个人认为应该采用等小应力来描述应力集中的现象，因为在实际中很难找到真正的单轴拉压的情况，一般结构的受力都没有这么简单，所以在分析的时候需要用等效应力来将各主应力进行转化，因此应该用等效应力来描述应力集中的现象。

4）等效泊松比就是泊松比，等效应力计算时不会用到泊松比，不过在计算mises 等效应变时会用到。

对于泊松比的取值原则应遵循以下两条：a:对于elastic & thermal strains 泊松比取为材料的泊松比；b:对于plastic creep hyperelastic strains 泊松比取为0.5。

ANSYS中整体、单元刚度和质量矩阵的提取

ANSYS中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。

1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。

基于ansys的一个用户开发程序例子编了一个程序（附件中）。

开发环境：compaqfortran6.5运行环境：win2000。

一个主文件self.for；另一个文件matrixout.f90用于矩阵输出；binlib.lib 为ansys提供的库文件，将其引入项目中（也可直接扔进debug目录）；.full文件由子空间迭代模态分析获得。

运行编译后的可执行文件.exe其他文件见/f/EE24A2F87F524606.html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。

定义超单元，然后列出超单元的刚度矩阵即可。

下面是一个小例题，自可明白。

/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,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵1.3、HBMAT命令法提取整体矩阵13.1、命令说明命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。

ext---输出矩阵文件的扩展名，缺省为.matrix。

form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。

ANSYS中整体、单元刚度和质量矩阵的提取

ANSYS中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。

1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。

基于ansys的一个用户开发程序例子编了一个程序（附件中）。

开发环境：compaqfortran6.5运行环境：win2000。

运行编译后的可执行文件.exe其他文件见/f/EE24A2F87F524606.html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。

定义超单元，然后列出超单元的刚度矩阵即可。

下面是一个小例题，自可明白。

ext---输出矩阵文件的扩展名，缺省为.matrix。

form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。

ANSYS问题汇总

1：ansys中的等效应力是什么物理含义？ (2)2．ansys后处理中负值的应力是压应力还是拉应力？ (3)3解决非线性分析不收敛的技巧！ (3)4非线性计算完的收敛图线，如何看他的收敛性呢，每条颜色的线代表什么意思呢？ (4)5求教accat及lccat命令 (5)6有关分块后的merge问题。

(5)7请教如何用APDL命令程序提取单元相关信息。

(5)8．ansys的刚度矩阵是在那一步骤生成的，如何读出，其格式如何 (5)9在混凝土的计算中，如何选择裂缝模型 (6)10．请问TB命令怎么用？？TB命令是用在非线性材料里吗？那么mp和TB有什么区别啊？一般什么情况下可以用TB命令？ (7)11．ansys 如何输入yield stress. 712、将Ansys计算结果输出到Tecplot的完整命令流！ (7)13、【原创】将数组中数据导入表中命令流，然后用曲线画出 (8)14、有一个问题，就是我想看我在建模时用的是什么单位，本来是mm的，我可能用成cm 了，怎么查看？1115、ansys是否可以自动搜索实体边界 (11)16、问一下疲劳计算得出的结果都是什么意思？？补充资料 (11)17、温度荷载问题 (12)18、如何把指定位置的节点的节点号提取出来？用什么命令？ (12)19、划分网格后修改单元属性问题 (12)20、快速去掉窗口的ansys标志 (13)21、merge节点与glue-mesh的区别 (13)22、glue不glue的区别？？？ (13)23、画等应力线大全，呵呵 (13)24、荷载步数太多了，该怎么办？ (14)25、我能否在建立好模型后再定义我的元素属性 (14)26、怎样划分不平行于全球坐标系xy的面 (15)27、如何查看层间应力值？ (15)28、在ANSYS中用表面效应单元加任意方向的荷载 (15)29、撰写论文中图像的灰度处理 (16)30、提取材料号2的所有体积命令是什么呢?. 1731、请教：solve后的warning怎样才能不用手动点OK?. 1732、循环计算生成的转矩结果数据如何保存成TXT文件 (17)33、请问不在同一平面的两个面如何合成一个面啊？？？ (17)34、关于hardpoint。

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中两种方式给材料添加材料属性1 第一种在划分网格之前指定main menu/preprocessor/meshing/mesh attributes/default attribs 显现meshing attributes对话框，在【mat】material number下拉框当选择你需要的材料序号。

单击ok然后划分网格，那么这次划分的网格的材料属性为选择的材料序号的属性。

2 第二种在划分网格以后指定先划分好网格点击select/entities/ 第一项选择areas ，第二项选择by num/pick,然后点击ok ，弹出面积选择框，选定面积，点击ok，完成面积选择点击select/entities,第一项选择elements，第二项选择attached to ，第三项选择areas，表示所要选择的单元为已选定面积中的单元，点击ok，选中面中的所有单元。

点击plot/replot,将只显示已选定的单元和面积。

点击main menu/preprocessor/material pros/change mat num,在new material number 文本框中输入你需要的材料序号，在elements No. to modefied 输入all 表示所选定的所有单元对应的材料属性转为此材料属性。

ansys多种材料如何设置材料属性呀,用什么命令GUI方式楼上正解，或Proprecessor->Meshing->Mesh Attributes->Picked Volume s命令为VATT本人喜爱在划分单元前先选好材料、实常数等再划分，命令流如下：type,1mat,1real,11vmesh,all关于其他不同材料，方式相同：ansys中的等效应力是什么物理含义它与最大应力s1有什么区别，平常讨论应力散布，应该用等效应力仍是最大应力s1呢1）计算等效应力时是不是需要输入等效泊松比呢仿佛有效泊松比的默许值是。

ANSYS中几种提取模态的方法

ANSYS中⼏种提取模态的⽅法在 ANSYS 中有以下⼏种提取模态的⽅法：– (1) Block Lanczos 法– (2) ⼦空间法– (3) PowerDynamics 法– (4) 减缩法– (5) 不对称法– (6) 阻尼法使⽤何种模态提取⽅法主要取决于模型⼤⼩（相对于计算机的计算能⼒⽽⾔）和具体的应⽤场合。

(1) Block Lanczos 法Block Lanczos 法可以在⼤多数场合中使⽤：- 是⼀种功能强⼤的⽅法，当提取中型到⼤型模型（50,000 ~ 100,000 个⾃由度）的⼤量振型时（40+），这种⽅法很有效；- 经常应⽤在具有实体单元或壳单元的模型中；- 在具有或没有初始截断点时同样有效。

（允许提取⾼于某个给定频率的振型）；- 可以很好地处理刚体振型；- 需要较⾼的内存。

(2) ⼦空间法⼦空间法⽐较适合于提取类似中型到⼤型模型的较少的振型（<40）- 需要相对较少的内存；- 实体单元和壳单元应当具有较好的单元形状，要对任何关于单元形状的警告信息予以注意；- 在具有刚体振型时可能会出现收敛问题；- 建议在具有约束⽅程时不要⽤此⽅法。

(3) PowerDynamics 法PowerDynamics 法适⽤于提取很⼤的模型（100.000个⾃由度以上）的较少振型（< 20）。

这种⽅法明显⽐Block Lanczos 法或⼦空间法快，但是：- 需要很⼤的内存；- 当单元形状不好或出现病态矩阵时，⽤这种⽅法可能不收敛；- 建议只将这种⽅法作为对⼤模型的⼀种备⽤⽅法。

注: PowerDynamics ⽅法- ⼦空间技术使⽤ Power 求解器 (PCG) 和⼀致质量矩阵；- 不执⾏ Sturm 序列检查 (对于遗漏模态); 它可能影响多个重复频率的模型；- ⼀个包含刚体模态的模型, 如果你使⽤ PowerDynamics ⽅法，必须执⾏ RIGID 命令 (或者在分析设置对话框中指定RIGID 设置)。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

/solu
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
(d25.15)
*vread,rhsval(1),hbfile,txt,,,rhscrd,,,ls0+ptrcrd+indcrd+valcrd
(d25.15)
*dim,smatr,,nrow,ncol
*do,icol,1,ncol
stacol=pointr(icol)
endcol=pointr(icol+1)
第4行：格式（2A16,2A20），分别表示列指针格式、行索引格式、系数矩阵数值格式、右边项数值格式。
第5行：格式（A3,11X,2I14），A3各列分别表示右边项格式、应用高斯起始矢量、应用eXact求解矢量；两个整数分别表示右边项列数、行索引数。三个字符中的第1个字符可取：F---全部存贮（如节点荷载向量的全部元素）、M---与系数矩阵相同方法。
3、HBMAT命令法提取整体矩阵
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：
=ASCII：ASCII码格式；
=BIN：二进制格式。
matrix---定义输出矩阵的类型，其值可取：
=STIFF：输出刚度矩阵。可用于写入了.FULL文件的任何类型的分析。
=MASS：输出质量矩阵。可用于特征值屈曲、子结构分析、模态分析。
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,1, , , , , , ,1
MSHAPE,0,3D
MSHKEY,1CM,_Y,VOLU
VSEL, , , ,1
CM,_Y1,VOLU
*dim,rhsval,,rhscrd
*vread,pointr(1),hbfile,txt,,,ptrcrd,,,ls0
(f14.0)
*vread,rowind(1),hbfile,txt,,,indcrd,,,ls0+ptrcrd
(f14.0)
*vread,values(1),hbfile,txt,,,valcrd,,,ls0+ptrcrd+indcrd
根据Harwell-Boeing文件格式，可读取矩阵的任意行列元素的数值，也可编程还原为满矩阵存储，以便它用，很显然这种提取方式比较方便。如当生成.FULL文件后，可采用命令/AUX2$FILE,mywork,full$HBMAT,mystiff,txt,ASCII,STIFF,YES$FINISH将二进制mywork.full文件输出为ASCII码文件mystiff.txt，并输出右边项。
命令流：
/aux2
file,hbfile,full
hbmat,hbfile,txt,,ascii,stiff,yes
finish
!从hbfile.txt读入数据，并还原为满矩阵存储
*dim,contline,,5
*vread,contline(1),hbfile,txt,,,5,,,1
=DAMP：输出阻尼矩阵。仅用于有阻尼的模态分析。
rhs---右边项输出控制（右边项指用矩阵所表示方程的等号右端矢量，这里可为节点荷载向量），如rhs=YES则输出，如rhs=NO则不输出。
模态分析时，因仅LANB和QR法可生成完整的质量矩阵，因此也仅采用这两种方法时才可使用HBMAT命令得到质量矩阵文件。
*enddo
*enddo
pointr=
rowind=
values=
rhsval=
icol=
irow=
ls0=
stacol=
endcol=
truerow=
totcrd=
ptrcrd=
indcrd=
valcrd=
rhscrd=
FINISH
/cle
FINISH
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
第2行：格式（5I14），分别表示该文件的总行数（不包括文件头）、矩阵列指针的总行数、矩阵行索引的总行数、矩阵元素数值的总行数、右边项总行数。
第3行：格式（A3,11X,4I14），分别为矩阵类型、矩阵行数、矩阵列数、矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）、单元元素数（对组装后的矩阵此值为0）。
FITEM,2,34
FITEM,2,47
FITEM,2,57
FITEM,2,68
FITEM,2,80
FITEM,2,90
!*
/GO
D,P51X, , , , , ,ALL, , , , ,
FLST,2,1,1,ORDE,1
FITEM,2,79
第6行后：矩阵元素值（单列）。
矩阵类型用3个字符表示，第1个字符可取：R---实数矩阵、C---复数矩阵、P---仅矩阵结构（无元素数值）；第2个字符可取：S---对称矩阵、U---不对称矩阵、H---Hermitian矩阵、Z---病态对称矩阵；R---带状矩阵；第3个字符可取：A---组装的矩阵、E---单元矩阵（未组装）。对称矩阵只存储下三角元素，如结构刚度矩阵为对称矩阵，Harwell-Boeing格式则仅记录下三角元素。
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
(5f14.0)
totcrd=contline(1)
ptrcrd=contline(2)
indcrd=contline(3)
valcrd=contline(4)
rhscrd=contline(5)
*vread,contline(1),hbfile,txt,,,4,,,2
(a3,11x,4f14.0)
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,2,9,1,ORDE,9
FITEM,2,1
FITEM,2,14
FITEM,2,24
*do,irow,stacol,endcol-1
truerow=rowind(irow)
smatr(truerow,icol)=values(irow)
*enddo
*enddo
*do,irow,1,nrow
*do,icol,1,ncol
smatr(irow,icol)=smatr(icol,irow)
!*
/GO
F,P51X,FY,-1000
/solu
/output,elemstiff,txt ! 进入求解层，定义输出单元刚度矩阵文件名和扩展名
/debug,-1,,,1 ! 设置输出单元刚度矩阵和单元荷载向量
solve
/output
finish ! 求解，设置输出到终端
FINISH
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中
2-----------------------------------------------------
nrow=contline(2)
ncol=contline(3)
strline=
contline=
*if,rhscrd,eq,0,then
ls0=4
*else
ls0=5
*endif
*dim,pointr,,ptrcrd
*dim,rowind,,indcrd
*dim,values,,valcrd
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中