ansys命令流学习笔记5-圆柱形shell单元的复合材料分析
ANSYS命令流学习笔记14-shell单元的铺层复合材料分析
! ANSYS命令流学习笔记14-shell单元的铺层复合材料分析!学习重点：!1、熟悉复合材料的材料特点工程应用中典型的复合材料为纤维增强复合材料。
玻璃纤维增强塑料（玻璃钢）、碳纤维、石墨纤维、硼纤维等高强度和高模量纤维。
复合材料各层为正交各向异性材料（Orthotropic）或者横向各向异性材料（Transversal Isotropic），材料的性能与材料主轴的取向有关。
各向异性Anisotropic，一般的各项同性材料需要两个材料参数弹性模量E和泊松比v。
而各向异性在XYZ有着不同的材料属性，而且拉伸行为和剪切行为互相关联。
定义其几何方程需要21个参数。
正交各向异性orthotropic，在XYZ有着不同的材料属性，而且拉伸行为和剪切行为无关，定义材料需要9个参数：Ex，Ey，Ez，Vxy，Vyz，Vxz，Gxy，Gyz，Gxz。
横向各向异性Transversal Isotropic，属于各向异性材料，但是在某个平面上表现出二维上的各向同性。
!2、熟悉复合材料分析所用的ANSYS单元复合材料单元关键在于能够实现铺层。
不同截面属性的梁单元（beam188, beam189, elbow290），2D对称壳单元（shell208, shell209），3D铺层壳单元（shell181, shell281, shell131, shell132），3D铺层实体单元（solid185, solid186, solsh190, solid278, solid279），均能实现复合材料的搭建。
其中Beam单元和2D对称壳单元很少使用。
SHELL91、SHELL99、SOLID46、SOLID191用于一些以前的分析教程中，但是现在这些单元已经被淘汰，最好选择下列单元区替代他们。
用越来越少的单元做越来越多的事情也是趋势。
Shell208和shell209，2D对称壳单元前者为2节点3自由度单元，后者为3节点3自由度单元，均能用于薄板和中厚板结构（L/h > 5-8）。
ANSYS笔记
1、适用PLANE182和SHELL181建立单元时，旋转成体以后要手动将两个单元删除。
而是用MESH200时，系统会在运算时自动删除或关闭该单元状态。
2、EXPOPT体扫掠相关信息的定义，可用于旋转产生体，定义产生体的单元属性、单元尺寸、是否删除源面上的网格、是否自动选择源面和目标面等。
3、施加周向和径向约束需要在柱坐标系下进行。
4、柱坐标系的平面为XOY平面，在切换到柱坐标系时，要注意工作平面是哪个平面，如果图形不满足XOY平面，应建立局部坐标系，调整平面。
5.局部坐标系中蓝色代表Z轴，白色代表X轴，黄色代表Y轴。
6、转速的施加一般是按照弧度进行施加。
7、考虑预应力模态分析的求解，需要先进性预应力效应打开时静力分析，求解一次，然后点击FINISH，之后再进入求解，选择模态分析进行求解。
一、二1.1、采用三维单元分析平面问题时需要约束其Z向的全部位移。
三1、梁柱铰接可以通过两种形式设置（350）：（1）梁柱连接处共用一个节点，使用梁单元弯矩释放功能，将转动自由度和平动自由度释放，可近似实现铰接。
（2）在建立模型是，梁柱不共用节点，通过约束方程，耦合节点位移实现铰接。
使用CP 命令。
2、BEAM188单元画内力图时，设置KEYOPT,1,3,3设置三次形函数可消除弯矩图出现锯齿状。
3、施加重力加速度时，方向与实际方向相反。
4、单元表中坐标系以单元坐标系为准。
四1、求解塑性极限荷载时，荷载大小未知：（1）可以通过理论求解，将大概的理论之求解出来，通过加载距离理论值相近的值进行求解，选取求解不收敛失败前最后一个荷载值作为极限荷载。
（2）通过不断试错，加一个很大的荷载，加到知道计算不收敛，将最大的荷载乘以不收敛的时刻点作为屈服极限荷载。
2、非线性求解中，时间点（TIME）可以写10也可以写1，最终数值要乘以时间点，最好填写1，方便计算。
求解完成后，打开Results Summary中TIME/FREQ列乘以施加力，可得到每一荷载子步施加的力的大小。
ansys_复合材料分析介绍
SOLID95 是 20 节点的结构实体单元，在 KEYOPT(1)=1 时，其作用与单
层的 SOLID191 单元类似，包括应用方位角和失效准则，还允许非线性材料和大
◆ 料。
5.2.2
BEAM188 和 BEAM189 为三维有限应变梁单元，其截面可以包含多种材
定义材料的叠层结构
复合材料最重要的特征就是其叠层结构。每层材料都有可能由不同的正交各 向异性材料构成，并且其主方向也可能各不相同。对于叠层复合材料，纤维的方 向决定了层的主方向。 有两种方法可用来定义材料层的配置： 通过定义各层材料的性质； 通过定义表示宏观力、力矩与宏观应变、曲率之间相互关系的本构矩阵(只 适合于 SOLID46 和 SHELL99)。
5.2.2.1 定义各层材料的性质
这种方法由下到上一层一层定义材料层的配置。底层为第一层，后续的层沿 单元坐标系的 Z 轴正方向自底向上叠加。如果叠层结构是对称的，可以只定义一 半的材料层。 有时，某个物理层可能只延伸到模型的一部分。为了建立连续的层，可以把 这些中断的层的厚度设置为零，图 5-1 显示了一个四层模型，其中第二层在某处 中断了。
1
及一个特殊的“三明治”选项， 而 SHELL99 则不能。另外 SHELL91 更适用于大 变形的情况。 3、SHELL181—有限应变壳单元 SHELL181 是四节点三维壳单元，每个节点有六个自由度。该单元支持所有 的非线性功能（包括大应变），允许有多达 250 层材料层。应该通过截面命令， 而不是实常数来定义层的信息，可以通过 FC 命令来指定失效准则。 4、SOLID46—三维层状结构体单元 SOLID46 是八节点三维实体单元 SOLID45 的一种叠层形式，其每个节点有 三个自由度(UX, UY, UZ)。它可用来建立叠层壳或实体的有限元模型，每个单元 允许有多达 250 层的等厚材料层， 或者 125 层的厚度在单元面内呈现双线性变 化的不等厚材料层。 该单元的另一个优点是可以用叠加几个单元的方式来对多于 250 层的复合材料建立模型，并允许沿厚度方向的变形斜率连续。用户也可输入 自己的本构矩阵。SOLID46 调整横向的材料特性，以允许在横向上为常应力。与 八节点壳单元相比较，SOLID46 的阶次要低些，因此，如在壳结构应用中要得到 与 SHELL91 或 SHELL99 相同的求解精度，需要更密的网格。 5、SOLID191--层状结构体单元 SOLID191 是 20 节点三维实体单元 SOLID95 的一种叠层形式，其每个节点 有三个自由度(UX, UY, UZ)。它可用以建立厚的叠层壳或实体的有限元模型，每 个单元允许有多达 100 层的材料层。与 SOLID46 类似，SOLID191 可以模拟厚度 上的不连续。SOLID46 可以调整横向的材料特性，以允许在横向上为常应力。这 个单元不支持非线性材料或大挠度。 6、其他 除上述层单元外，还有其它的一些具有层功能的单元： ◆ 挠度。 ◆ SHELL63 是四节点壳单元，可用于对“三明治”壳结构作粗糙、近似 的计算。 象两块金属片之间夹有一层聚合物的问题就很典型，此时聚合物的弯曲 刚度相对于金属片的弯曲刚度来说是一个小量。用户可以用实常数 RMI 来修正 单元的弯曲刚度， 使其等效于由金属片引起的弯曲刚度。从中面到外层纤维的距 离(实常数 CTOP 和 CBOT)可用来获得“三明治”壳的表层输出应力。这种单元 不如 SHELL91 、SHELL99 和 SHELL181 那样用得频繁，故后面不再论述。 ◆ SOLID65 是三维钢筋混凝土实体单元，可以模拟在三个用户指定方向 配筋的各向同性介质。
ANSYS命令流总结(全)
ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。
s就是select，选择。
DIM 是定义数组的意思。
array 数组。
MP命令用来定义材料参数。
K是建立关键点命令。
K,关键点编号,x坐标,y坐标，z坐标。
K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。
NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。
这个命令对于重复的线面都可以用。
这个很简单，压缩关键。
Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50 ：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。
LNSRCH 激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。
SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type： array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX 各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type 为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, VALUELab为参数名称 value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108，超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81，SOLID98，FLUID129，INFIN110，111，FLUID116，130界面单元:INTER192，193，194，195显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK183D有限应EPCDFG另可考虑粘0 变杆 B 弹塑性E-弹性(Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent)等。
(完整版)ANSYS命令流总结(全)
ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。
s就是select，选择。
DIM是定义数组的意思。
array 数组。
MP命令用来定义材料参数。
K是建立关键点命令。
K,关键点编号,x坐标,y坐标，z坐标。
K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。
NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。
这个命令对于重复的线面都可以用。
这个很简单，压缩关键。
Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。
LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。
SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108， 超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81， SOLID98，FLUID129，INFIN110，111，FLUID116，130 界面单元:INTER192，193，194，195 显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent )等。
ansys-SHELL单元
Q:用板壳元shell63单元建模时,如下图示,两个平面互相垂直,如何使这两个平面保持为一个整体来受力.因为实际结构中,比如一个由钢板焊成的箱梁,内有加劲隔板,如何模拟.由于这两个面没有公共边,用了粘贴及搭接都不行,试问怎样简单实现整体性这一目的.A:分网时控制单元长度，使面与面之间有重合节点或相邻节点，然后使用Merge Items或Adjacent Regions节点既可把面与面给“焊起来”。
Adjacent Regions命令时，必须先将一个区域的节点和相邻区域的单元选择好，再自动生成约束方程，是这样的吗，对整体计算的精度应有保证吧Adjacent Regions命令的使用确如楼上所说，但对整体计算精度能否保证就不好妄加评论了。
其实我认为SHELL63单元连接的方法是多种多样的，如何选择关键在于快捷，准确。
当同位置的结点难以实现，Adjacent Regions命令又较困难时，那么还有相邻节点重合（Coincident Nodes），只要知道要求重合节点的距离容差就行了Q:请教ansys中梁壳两种单元一起建模时遇到的问题。
我分析了一个地下结构，楼板用shell63单元。
立柱用BEAM4单元。
结果出现的问题好奇怪，说单元1没有定义材料类型。
其实不是的。
后来我建了一个简单的模型，还是同样的问题。
我非常不解。
我把模型发了上来，哪位好心人能否帮我诊断一下。
谢谢了。
A:/prep7et,1,shell63mp,ex,1,3e10mp,prxy,1,0.2mp,dens,1,2500keyopt,1,3,2!梁壳一起建模时必须要考虑shell中面的转动刚度r,1,1et,2,beam4mp,ex,2,3e10mp,prxy,2,0.2mp,dens,2,2500r,2,0.785,0.049,0.049,1,1,,,0.098,,,,k,1,0,0,10k,2,0,0,0k,3,10,0,0k,4,10,0,10a,1,2,3,4k,9,0,-5,10k,10,0,-5,0k,11,10,-5,0k,12,10,-5,10a,9,10,11,12!a,1,9,12,4wplane,100,0,0,0,10,0,0,0,10,0 csys,4*do,i,1,4,1wpave,0,0,2asbw,all*enddowplane,100,0,0,0,0,0,-10,0,10,0 csys,wp*do,i,1,4,1wpave,0,0,2asbw,all*enddonummrg,allnumcmp,allcsys,0!恢复原始的坐标系l,64,65lsel,s,,,27lesize,all,1latt,2,2,2lmesh,allallselaatt,1,1,1mshkey,0mshape,0esize,2amesh,allallselnummrg,allnumcmp,allallselnsel,s,loc,y,-5d,all,allallselasel,s,loc,y,0sfa,all,,pres,57e3allsel/soluantype,staticsolvefinish我把模型发了上来，出现的错误真是令人匪夷所思。
Ansys复合材料结构分析操作指导书
Ansys10.0 复合材料结构分析操作指导书第一章概述复合材料是两种或两种以上物理或化学性质不同的材料复合在一起而形成的一种多相固体材料，具有很高的比刚度和比强度（刚度和强度与密度的比值），因而应用相当广泛，其应用即涉及航空、航天等高科技领域，也包括游艇、风电叶片等诸多民用领域。
由于复合材料结构复杂，材料性质特殊，对其结构进行分析需要借助数值模拟的方法，众多数值模拟软件中Ansys是个不错的选择。
Ansys软件由美国ANSYS公司开发，是目前世界上唯一一款通过ISO9001质量体系认证的分析设计软件，有着近40年的发展历史，经过多次升级和收购其它CAE(Computer Aided Engineering )软件，目前已经发展成集结构力学、流体力学、电磁学、声学和热学分析于一体的大型通用有限元分析软件，是一款不可多得的工程分析软件。
Ansys在做复合材料结构分析方面也有不俗的表现，此书将介绍如何使用该款软件进行复合材料结构分析。
在开始之前有以下几点需要说明，希望大家能对有限元法有大体的认识，以及Ansys软件有哪些改进，最后给出一些学习Ansys软件的建议。
1、有限元分析方法应用简介有限元法(Finite Element Method，简称FEM)是建立在严格数学分析理论上的一种数值分析方法。
该方法的基本思想是离散化模型，将求解目标离散成有限个单元(Element)，并在每个单元上指定有限个节点(Node)，单元通过节点相连构成整个有限元模型，用该模型代替实际结构进行结构分析。
在对结构离散后，要求解的基本未知量就转变为各个节点位移(Ansys中称之为DOF(Degree Of Freedom)，试想一下，节点的位移包括沿x,y,z轴的平动和转动，也就是节点的自由度)，节点位移通过求解一系列代数方程组得到，在求得节点位移后，利用节点位移和应力、应变之间的关系矩阵就可以求出各个节点上的应力、应变，应用线性插值便可以获得单元内任意位置的位移、应力、应变等信息。
Ansys复合材料结构分析操作指导书---副本
第四章复合材料计算实例在有了前几章知识做铺垫，这一章我们来学习两个复合材料分析的例子，加深复合材料分析的理解，也希望读者能从中收获一些经验。
在这里将第二章的流程图再次拿出来，进一步熟悉ANSYS有限元分析的基本过程。
图7 Ansys 结构分析流程图4.1 层合板受压分析4.1.1 问题描述层合板指的是仅仅由FRP层叠而成的复合板材，中间不包含芯材，板材的性能不仅与纤维的弹性模量、剪切模量有关，还与纤维的铺层方向有着密切关系。
本例中的板材有4层厚度为0.025m的单元板复合而成，单元板的铺层方向为0°、90°、90°、0°，见图13所示。
单元板的材料属性见表4.1。
表 4.1 单元板材料属性EX/MP EY/MP EZ/MP GXY/MP GYZ/MP GXZ/MP PRXY PRYZ PRXZ12.5 300 300 50 20 50 0.25 0.25 0.01图13 复合材料板4.1.2 求解步骤根据问题描述，所要分析的问题为壳体结构的复合材料板，可以采用SOLID46单元建立3D有限元模型进行分析。
结合图7的一般步骤进行分析。
步骤一：选取单元类型，设置单元实常数⑴、在开始一个新分析前，需要指定文件保存路径和文件名。
文件保存路径GUI：【Utility Menu】|【File】|【Change Directory】见图14指定新的文件名GUI：【Utility Menu】|【File】|【Change Jobname】见图15所示图14 指定文件保存路径图15 修改文件名⑵、选取单元类型1）选取单元类型的GUI操作：【Main Menu】|【Preprocessor】|【Element Type】|【Add/Edit/Delete】，执行后弹出Element Types对话框。
2）在Element Types对话框点击Add定义新的单元类型，弹出“Library of Element Types”对话框，见图16所示，按图中所示选择，单元类型参考号输入框中输入数字1。
Ansys命令流大全剖析
Ansys彭彭(沈阳化工大学机械工程学院，辽宁沈阳110142)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。
点要依次序输入，输入的顺序会决定面的法线方向。
如果超过四个点，则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。
8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。
[整理]Ansys复合材料结构分析总结.
![[整理]Ansys复合材料结构分析总结.](https://img.taocdn.com/s3/m/00857e6c16fc700aba68fc2b.png)
Ansys复合材料结构分析总结说明：整理自Simwe论坛，复合材料版块，原创fea_stud，大家要感谢他呀目录1# 复合材料结构分析总结（一）——概述篇5# 复合材料结构分析总结（二）——建模篇10# 复合材料结构分析总结（三）——分析篇13# 复合材料结构分析总结（四）——优化篇做了一年多的复合材料压力容器的分析工作，也积累了一些分析经验，到了总结的时候了，回想起来，总最初采用I-deas，到MSC.Patran、Nastran，到最后选定Ansys为自己的分析工具，确实有一些东西值得和大家分享，与从事复合材料结构分析的朋友门共同探讨。
（一）概述篇复合材料是由一种以上具有不同性质的材料构成，其主要优点是具有优异的材料性能，在工程应用中典型的一种复合材料为纤维增强复合材料，这种材料的特性表现为正交各向异性，对于这种材料的模拟，很多的程序都提供了一些处理方法，在I-Deas、Nastran、Ansys中都有相应的处理方法。
笔者最初是用I-Deas下建立各项异性材料结合三维实体结构单元来模拟（由于研究对象是厚壁容器，不宜采用壳单元），分析结果还是非常好的，而且I-Deas强大的建模功能，但由于课题要求要进行压力容器的优化分析，而且必须要自己写优化程序，I-Deas的二次开发功能开放性不是很强，所以改为MSC.Patran，Patran 提供了一种非常好的二次开发编程语言PCL（以后在MSC的版中专门给大家贴出这部分内容），采用Patran结合Nastran的分析环境，建立了基于正交各项异性和各项异性两种分析模型，但最终发现，在得到的最后结果中，复合材料层之间的应力结果始终不合理，而模型是没有问题的（因为在I-Deas中，相同的模型结果是合理的），于是最后转向Ansys，刚开始接触Ansys，真有相见恨晚的感觉，丰富的单元库，开放的二次开发环境（APDL 语言），下面就重点写Ansys的内容。
在ANSYS程序中，可以通过各项异性单元（Solid 64）来模拟，另外还专门提供了一类层合单元（Layer Elements）来模拟层合结构（Shell 99, Shell 91, Shell 181, Solid 46 和Solid 191）的复合材料。
最新ansys_复合材料分析汇总
a n s y s_复合材料分析第五章复合材料5.1 复合材料的相关概念复合材料作为结构应用已有相当长的历史。
在现代，复合材料构件已被大量应用于飞行器结构、汽车、体育器材及许多消费产品中。
复合材料由一种以上具有不同结构性质的材料构成，它的主要优点是具有很高的比刚度(刚度与重量之比)。
在工程应用中，典型复合材料有纤维和叠层型材料，如玻璃纤维、玻璃环氧树脂、石墨环氧树脂、硼环氧树脂等。
ANSYS程序中提供一种特殊单元--层单元来模拟复合材料。
利用这些单元就可以作任意的结构分析了(包括非线性如大挠度和应力刚化等问题)。
对于热、磁、电场分析，目前尚未提供层单元。
5.2 建立复合材料模型与铁或钢等各向同性材料相比，建立复合材料的模型要复杂一些。
由于各层材料性能为任意正交各向异性，材料性能与材料主轴取向有关，在定义各层材料的材料性能和方向时要特别注意。
本节主要探讨如下问题：选择合适的单元类型；定义材料层；确定失效准则；应遵循的建模和后处理规则。
5.2.1 选择合适的单元类型用于建立复合材料模型的单元类型有SHELL99、SHELL91、SHELL181、SOLID46和SOLID191 五种单元。
但 ANSYS/Professional 只能使用SHELL99 和 SHELL46 单元。
具体应选择哪一类单元要根据具体应用和所需计算结果类型等来确定。
所有的层单元允许失效准则计算。
1、SHELL99--线性层状结构壳单元SHELL99 是一种八节点三维壳单元，每个节点有六个自由度。
该单元主要适用于薄到中等厚度的板和壳结构，一般要求宽厚比应大于10。
对于宽厚比小于10的结构，则应考虑选用 SOLID46 来建立模型。
SHELL99 允许有多达 250 层的等厚材料层，或者 125 层厚度在单元面内呈现双线性变化的不等材料层。
如果材料层大于 250，用户可通过输入自己的材料矩阵形式来建立模型。
还可以通过一个选项将单元节点偏置到结构的表层或底层。
ANSYS分析实例与工程应用命令流学习笔记
ANSYS分析实例与工程应用命令流学习笔记1大纲静力分析：2杆、3梁、5薄膜和板壳、4实体单元梁单元：简化计算，结构总体受力情况实体单元：较复杂的结构，局部细节的受力情况稳定性分析：6振动、模态分析：7简单振动和梁的振动、8膜板和实体振动2杆系结构的静力分析2.1铰接杆在外力作用下的变形二维杆单元LINK1*AFUN,DEG：三角函数默认为弧度，改为角度后处理：结构变形图、显示节点位移和杆件应力2.2人字形屋架的静力分析后处理：杆单元的轴力、轴向应力、轴向应变2.3超静定拉压杆的反力计算后处理：节点反力2.4平行杆件与刚性梁连接的热应力问题定义3点的UY为耦合自由度，即三者的UY位移相等温度（增量）后处理：寻找特定位置的节点和单元，并从单元表中提取它们的内力2.5端部有间隙的杆的热膨胀二维带厚度的平面应力单元PLANE42、二维接触单元CONTACT26温度（始、末）后处理：定义水平应力和铅直应力单元表，并提取3号单元的应力结果*Status,ParmFINISH定义数组变量，将计算结果通过数组变量输出到文件3梁的弯曲静力分析3.1单跨等截面超静定梁的平面弯曲二维弹性梁单元BEAM3后处理：定义以两端弯矩和剪力的单元表，并列出单元表数据并用单元表数据绘制剪力图和弯矩图更细的节点划分方案，更精细3.2四跨连续梁的内力计算体素建模：keypoint, line, area, volume便于细分单元3.3七层框架结构计算3.4工字形截面外伸梁的平面弯曲3.5矩形截面梁的纵横弯曲分析考虑应力强化效应后处理：迭代过程3.6空间刚架静力分析三维梁单元BEAM43.7悬臂梁的双向弯曲三维8节点耦合场实体单元SOLID5三维20节点固体单元SOLID92三维10节点耦合场实体单元SOLID98三维结构实体自适应单元SOLID147定义宏程序，对应四种工况，各种结果差别不大3.8圆形截面悬臂杆的弯扭组合变形三维直管单元PIPE16（只定义外直径，不定义内直径）3.9悬臂等强度梁的弯曲四边形壳单元SHELL63（这里用退化的三角形单元，并使用节点耦合自由度保证模型的对称变形）三维非对称锥形梁单元BEAM44（定义横截面主轴，单元宽度线性变化）计算结果都很好，但壳体单元更能模拟出等强度梁的实际几何形状，更直观，截面定义更简单。
ansys 复合材料分析
第五章复合材料5.1 复合材料的相关概念复合材料作为结构应用已有相当长的历史。
在现代，复合材料构件已被大量应用于飞行器结构、汽车、体育器材及许多消费产品中。
复合材料由一种以上具有不同结构性质的材料构成，它的主要优点是具有很高的比刚度(刚度与重量之比)。
在工程应用中，典型复合材料有纤维和叠层型材料，如玻璃纤维、玻璃环氧树脂、石墨环氧树脂、硼环氧树脂等。
ANSYS程序中提供一种特殊单元--层单元来模拟复合材料。
利用这些单元就可以作任意的结构分析了(包括非线性如大挠度和应力刚化等问题)。
对于热、磁、电场分析，目前尚未提供层单元。
5.2 建立复合材料模型与铁或钢等各向同性材料相比，建立复合材料的模型要复杂一些。
由于各层材料性能为任意正交各向异性，材料性能与材料主轴取向有关，在定义各层材料的材料性能和方向时要特别注意。
本节主要探讨如下问题：选择合适的单元类型；定义材料层；确定失效准则；应遵循的建模和后处理规则。
5.2.1 选择合适的单元类型用于建立复合材料模型的单元类型有SHELL99、SHELL91、SHELL181、SOLID46和SOLID191 五种单元。
但 ANSYS/Professional 只能使用 SHELL99 和SHELL46 单元。
具体应选择哪一类单元要根据具体应用和所需计算结果类型等来确定。
所有的层单元允许失效准则计算。
1、SHELL99--线性层状结构壳单元SHELL99 是一种八节点三维壳单元，每个节点有六个自由度。
该单元主要适用于薄到中等厚度的板和壳结构，一般要求宽厚比应大于10。
对于宽厚比小于10的结构，则应考虑选用 SOLID46 来建立模型。
SHELL99 允许有多达 250 层的等厚材料层，或者 125 层厚度在单元面内呈现双线性变化的不等材料层。
如果材料层大于 250，用户可通过输入自己的材料矩阵形式来建立模型。
还可以通过一个选项将单元节点偏置到结构的表层或底层。
2、SHELL91--非线性层状结构壳单元SHELL91 与 SHELL99 有些类似，只是它允许复合材料最多只有 100 层，而且用户不能输入自己的材料性能矩阵。
ansys命令流语法
ansys命令流语法ANSYS命令流语法是使用ANSYS软件进行仿真分析的关键部分。
它是一种将命令以特定顺序组合在一起的方式，以实现特定的分析目标。
本文将介绍ANSYS命令流语法的基本语法规则和常用命令，以及如何使用它们进行仿真分析。
一、ANSYS命令流语法的基本语法规则1. 命令的基本格式：命令[选项] [参数1, 参数2, ...]2. 命令的执行顺序：ANSYS命令流是按照命令的顺序逐条执行的。
如果需要改变执行顺序，可以使用条件语句、循环语句等控制结构。
3. 注释：可以在命令流中添加注释，以"!"开头。
注释部分不会被执行，可以用于解释命令的用途或添加说明。
4. 变量和参数：可以使用变量和参数来存储和传递数据。
变量以"$"开头，参数以"%"开头。
二、常用命令1. Preprocessor命令：用于定义和准备分析模型的预处理操作。
- /PREP7：进入预处理器界面。
- ET，MP，REAL等：定义单元类型、材料属性、实数等。
- K，L，A等：创建节点、单元、区域等。
2. Solution命令：用于设置和运行分析求解器。
- /SOLU：进入求解器界面。
- SOLVE，ANTYPE等：设置分析类型、求解选项等。
- D，S等：定义边界条件、加载条件等。
3. Postprocessor命令：用于后处理和分析结果的可视化。
- /POST1：进入后处理器界面。
- PLOT，PDEF等：绘制图形、定义图形属性等。
- PRINT，*VWRITE等：输出结果数据。
三、使用ANSYS命令流语法进行仿真分析使用ANSYS命令流语法进行仿真分析的一般步骤如下：1. 导入几何模型：使用CAD软件创建几何模型，并将其导入ANSYS 中。
2. 定义材料属性：根据实际材料的物理特性，使用MP命令定义材料属性。
3. 网格划分：使用网格划分命令划分几何模型，生成有限元网格。
ansys 复合材料添加分析流程
ansys 复合材料添加分析流程下载温馨提示:该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
文档下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!并且，本店铺为大家提供各种各样类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，如想了解不同资料格式和写法，敬请关注!Download tips: This document is carefully compiled by theeditor.I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!ANSYS复合材料分析流程详解在现代工程设计中，复合材料因其独特的性能组合而被广泛应用于航空、航天、汽车和能源等领域。
[整理]shell单元厚度复合材料.
![[整理]shell单元厚度复合材料.](https://img.taocdn.com/s3/m/155a4d326fdb6f1aff00bed5b9f3f90f76c64dd9.png)
[整理]shell单元厚度复合材料.6.1 复合材料结构分析基本过程6.1.1 概述复合材料是由两种或两种以上物理或化学性质不同的材料复合在一起而形成的一种多相固体材料，其主要优点是具有很高的比刚度(刚度与密度之比)。
复合材料作为结构材料应用已有很长的历史。
目前，复合材料的应用已非常普遍，其应用范围涉及航空、航天、军事、民用等诸多领域。
ANSYS程序提供了一种特殊的单元、层单元来模拟复合材料，利用这些单元就可以进行任意的复合材料结构分析。
复合材料结构分析也包括建模、加载求解及后处理3个基本步骤，其中加载求解及后处理基本同于一般的结构分析过程，建模部分具有其特殊性，下面主要对其建模部分进行详细讨论。
6.1.2 建立复合材料模型与一般的各向同性材料相比，复合材料的建模过程要相对复杂。
由于各层材料性能为任意正交各向异性，材料性能与材料主轴取向有关，所以在定义各层材料性能和方向时要特别注意。
在本节中，主要探讨以下4个问题。
1.选择适当的单元类型用于建立复合材料模型的单元有SHELL99、SHELL91、SHELL181、SOLSH190、SOLID46、SOLID186和SOLID191 7种单元。
单元类型的选择主要根据具体的应用和所需计算的结果类型来确定。
(1)SHELL99单元SHELL99是一种8节点3D壳单元，每个节点有6个自由度。
该类单元主要适用于薄到中等厚度的板和壳结构，一般要求结构宽厚比大于10。
对于宽厚比小于10的结构，则应考虑选择SOLID46单元建模。
SHELL99允许有多达250层的等厚度材料层，或者是125层的厚度在单元面内成双线性变化的不等厚度材料层。
如果材料层大于250，用户可通过输入自定义的材料矩阵来建立模型。
SHELL99单元可进行失效分析。
另外，该类型单元可以将单元节点偏置到结构的表层或底层。
(2)SHELL91单元SHELL91和SHELL99相类似，只是它允许的复合材料最多有100层，而且用户不能输入自定义的材料性能矩阵。
ansys命令流学习笔记5-圆柱形shell单元的复合材料分析
ansys命令流学习笔记5-圆柱形shell单元的复合材料分析! ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析!学习重点：!1、熟悉单元坐标系下的铺层当零件形状为规则圆筒时，如何进行铺层？建立局部的柱坐标系，将需要铺层单元坐标设置为局部坐标系，进行铺层即可。
譬如圆筒铺层的单元坐标系要建立局部圆柱坐标系。
如果还使用笛卡尔坐标系，铺层也能进行，但是铺层方向有很大不同，求解结果也会异常。
所以划分网格时，确认单元坐标系选择，划分网格之后，检查单元坐标系情况。
确认铺层方向符合预期要求。
本例中要特别注意横向（即Y向）是否符合要求。
!2、熟悉圆面的建模和局部坐标系建立不解释!3、熟悉利用MPC施加扭矩APDL如何对一个圆周施加扭矩？在圆心处建立一个节点，然后用MPC单元连接圆心节点和圆周节点，然后在圆心节点上施加一个扭矩即可。
注意将MPC单元的属性改为刚性梁。
注意这里MPC单元的利用，也是自己的一些理解。
很多细节也不知道如何在APDL实现。
!问题描述! 传动轴长度为1m，壁厚0.003m，直径0.08m，铺层共十层，角度为-45/45/-45/45/-45/45/-45/45/-45/45。
一端固定，一端圆周施加扭矩M=2000N·m。
复合材料为横向正交各向异性Ex，Ey，Ez，Vxy，Vyz，Vxz，Gxy，Gyz，Gxz分别为195e9Pa, 35e9Pa, 35e9Pa,0.28, 0.3, 0.3, 15e9Pa, 3.78e9Pa, 15e9Pa。
应力失效参数：+X:767E6Pa; -X:392E6Pa; +Y:20E6Pa; -Y:70E6Pa; +Z:30E6Pa; -Z:55E6Pa; Sxy: 41E6Pa; Syz: 30E6Pa; Sxz:41E6Pa。
应变失效参数：+X:0.05; -X:0.045; +Y:0.08; -Y:0.06; +Z:0.04; -Z:0.045; Sxy: 0.035; Syz: 0.042; Sxz:0.025。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
! ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析
!学习重点：
!1、熟悉单元坐标系下的铺层
当零件形状为规则圆筒时，如何进行铺层？建立局部的柱坐标系，将需要铺层单元坐标设置为局部坐标系，进行铺层即可。
譬如圆筒铺层的单元坐标系要建立局部圆柱坐标系。
如果还使用笛卡尔坐标系，铺层也能进行，但是铺层方向有很大不同，求解结果也会异常。
所以划分网格时，确认单元坐标系选择，划分网格之后，检查单元坐标系情况。
确认铺层方向符合预期要求。
本例中要特别注意横向（即Y向）是否符合要求。
!2、熟悉圆面的建模和局部坐标系建立
不解释
!3、熟悉利用MPC施加扭矩
APDL如何对一个圆周施加扭矩？在圆心处建立一个节点，然后用MPC单元连接圆心节点和圆周节点，然后在圆心节点上施加一个扭矩即可。
注意将MPC单元的属性改为刚性梁。
注意这里MPC单元的利用，也是自己的一些理解。
很多细节也不知道如何在APDL实现。
!问题描述
! 传动轴长度为1m，壁厚0.003m，直径0.08m，铺层共十层，角度为-45/45/-45/45/-45/45/-45/45/-45/45。
一端固定，一端圆周施加扭矩M=2000N·m。
复合材料为横向正交各向异性Ex，Ey，Ez，Vxy，Vyz，Vxz，Gxy，Gyz，Gxz分别为195e9Pa, 35e9Pa, 35e9Pa,0.28, 0.3, 0.3, 15e9Pa, 3.78e9Pa, 15e9Pa。
应力失效参数：+X:767E6Pa; -X:392E6Pa; +Y:20E6Pa; -Y:70E6Pa; +Z:30E6Pa; -Z:55E6Pa; Sxy: 41E6Pa; Syz: 30E6Pa; Sxz: 41E6Pa。
应变失效参数：+X:0.05; -X:0.045; +Y:0.08; -Y:0.06; +Z:0.04; -Z:0.045; Sxy: 0.035; Syz: 0.042; Sxz:0.025。
!APDL命令：
finish
/clear
/title, composite shaft
/prep7
et,1,shell181 !选择单元181
keyopt,1,8,1 !保存每一层的数据
et,2,184
keyopt,2,1,1 !定义MPC184单元，利用其施加扭矩。
将其属性定义为刚性梁
mptemp,1,0
mpdata,ex,1,,195e9
mpdata,ey,1,,35e9
mpdata,ez,1,,35e9
mpdata,prxy,1,,0.28
mpdata,pryz,1,,0.3
mpdata,prxz,1,,0.3
mpdata,gxy,1,,15e9
mpdata,gyz,1,,3.78e9
mpdata,gxz,1,,15e9 !定义各向同性材料
fc,1,s,xten,767e6
fc,1,s,yten,20e6
fc,1,s,zten,30e6
fc,1,s,xcmp,-392e6
fc,1,s,ycmp,-70e6
fc,1,s,zcmp,-55e6
fc,1,s,xy,41e6
fc,1,s,yz,30e6
fc,1,s,xz,41e6 !定义应力失效准则
fc,1,epel,xten,0.05
fc,1,epel,yten,0.08
fc,1,epel,zten,0.04
fc,1,epel,xcmp,-0.045
fc,1,epel,ycmp,-0.06
fc,1,epel,zcmp,-0.045
fc,1,epel,xy,0.035
fc,1,epel,yz,0.042
fc,1,epel,xz,0.025 !定义应变失效准则
!用到Tasi-Wu失效准则还需定义应力耦合系数，默认为-1，-1，-1。
实际值需要通过双轴试验测定，较难得到。
sectype,1,shell,, !sectype, secID,type,subtype,name,refinekey,定义截面类型
secdata, 3e-4,1,-45,3 !积分点为3个
secdata, 3e-4,1,45,3
secdata, 3e-4,1,-45,3
secdata, 3e-4,1,45,3
secdata, 3e-4,1,-45,3
secdata, 3e-4,1,45,3
secdata, 3e-4,1,-45,3
secdata, 3e-4,1,45,3
secdata, 3e-4,1,-45,3
secdata, 3e-4,1,45,3 !定义铺层
secplot,1 !查看ID为1的section
n,1,0,0,1 !建立一个node，为MPC做准备
cylind,0.04, ,0,1,0,360 !在坐标原点，创建一个圆柱体
vdele,all, , ,0
asel,s,loc,z,0,1
asel,u,loc,z,0.5
adele,all, , ,1
allsel !通过加减乘除，得到最终的面
local,11,1,0,0,0 !定义局部圆柱坐标系，原点000，ID=11
type, 1
mat, 1
esys, 11
secnum, 1
esize,0.007
amesh,all !网格划分属性设置，主要选择单元坐标系。
csys,0
type,2
mat,2
e,1,21
*do,i,1,18
e,1,163+i
*enddo
*do,i,1,17
e,1,2754+i
*enddo
!不知道怎么设置MPC连接。
出此下策。
!如果直接GUI操作，不需要指定MPC单元，直接用contact manage中的contact wizard生成约束方程，约束自由度即可。
f,1,mz,2000
nsel,s,loc,z,0
d,all,all
allsel
finish
!##################------------------------
/solu
solve !会提示两个警告，一个说是没开大变形，一个说是没定义MPC单元材料属性，可忽略
finish !求解完成
!##################------------------------
/post1
asel,s,loc,z,0.5, , ,1 !仅仅查看面部分，不看MPC单元。
layer,0
plnsol,u,sum
plnsol,s,eqv !查看整体结构
layer,2 !查看第2层结果
plnsol,s,eqv
layer,0 !切换回查看整体结果
plnsol,fail,emax
plnsol,fail,smax
plnsol,fail,twsr !校核三种准则下的危险系数，均小于1，则合格。
finish。