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。
- 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。