ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析

合集下载

（完整版）ANSYS最常用命令流+中文注释（超级大全）

（完整版）ANSYS最常用命令流+中文注释（超级大全）ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。

keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。

如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。

同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。

ansys_复合材料分析介绍

2
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有限元分析的基本过程。

图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的一般步骤进行分析。

步骤一：选取单元类型，设置单元实常数⑴、在开始一个新分析前，需要指定文件保存路径和文件名。

2）在Element Types对话框点击Add定义新的单元类型，弹出“Library of Element Types”对话框，见图16所示，按图中所示选择，单元类型参考号输入框中输入数字1。

(完整版)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中的坐标系及相应的命令

ANSYS中的坐标系及相应的命令2009-06-23 22:01ANSYS中的坐标系坐标系用于定义几何结构的空间位置，规定节点的自由度，定义材料的线性方向，以及改变图形显示和列表。

ANSYS中的坐标系有：总体坐标系，局部坐标系，节点坐标系，单元坐标系，显示坐标系，结果坐标系。

同一时刻只能有一个坐标系被激活。

总体坐标系：用于确定几何结构的空间位置，是绝对参考系。

如：笛卡尔坐标系(CSYS,0)，柱坐标系(CSYS,1)，球坐标系(CSYS,2)。

局部坐标系：由用户自己创建的(坐标系编号从11开始)，原点相对于总体坐标系的原点偏离了一定的距离或各轴相对于总体坐标系偏转了一定的角度。

定义的方法有：在特定位置(笛卡尔坐标系)定义(LOCAL)；通过已有节点定义(CS)；通过已有关键点定义(CSKP)；以当前定义的工作平面的原点为中心定义(CSWPLA)；通过已激活的坐标系定义(CLOCAL)。

删除局部坐标系(CSDELE)。

查看局部坐标系(CSLIST)。

节点坐标系：用于定义节点自由度的方向，需要在不同于总体坐标系的方向施加位移约束时用到。

每个节点都有自己的节点坐标系，默认为平行于总体笛卡尔坐标系。

定义的方法有：定义节点时直接设定(N)；将节点坐标系旋转到当前激活的坐标系的方向(NROTAT,可以批量操作)；按照给定的旋转角度旋转(NMODIF)；通过新坐标系各轴的方向余弦旋转(NANG)。

显示节点坐标系(NLIST)。

此外节点复制(NGEN)时，节点坐标系也一并复制。

单元坐标系：用于规定正交材料特性的方向和面力结果的输出方向。

每个单元均有各自的单元坐标系，默认为：线单元X轴正方向由该单元的I节点指向J节点；壳单元X轴正方向由该单元的I节点指向J节点，Z轴与壳面垂直并且通过I点，其正方向有单元的I、J、K节点按右手准则确定，Y轴垂直于X轴和Z轴；2D实体和3D实体单元的单元坐标系总是平行于总体笛卡尔坐标系。

(完整版)史上最全的ANSYS命令流大全

《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1，/PREP7 ! 加载前处理模块2，/CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/FILENAME, EX10.5 ! 定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题4，F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力6，FINISH ! 退出模块命令7，/POST1 ! 加载后处理模块8，PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9，ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRSETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORXETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_STETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_COETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSXETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY*GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果，存入变量STRSS_ST;*GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果，存入变量STRSS_CO10 FINISH !退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI !申明采用国际单位制14 /NUMBER, 2 !只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析(STA TIC或者0)OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL !指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第1个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容OUTRES,ALL,0 !设置将所有数据不记录到数据库。

Ansys的复合材料分析

– No material property matrix input – Supports plasticity, largestrain behavior; laminated and sandwich structures – Failure Modeling through TB, FAIL
© 2005 ANSYS, Inc.
8
ANSYS, Inc. Proprietary
SHELL99 Linear Layered Structural Shell Element
• Element Definition
– 8node, 3D shell element with six degrees of freedom at each node – Thin to moderately thick plate and shell structures with a sideto thickness ratio of roughly 10 or greater
© 2005 ANSYS, Inc.
rietary
Benefits of Composites
• Stronger and stiffer than metals on a density basis • Capable of high continuous operating temperatures • Highly corrosion resistant • Electrically insulating/conducting/selectively conducting properties • Tailorable thermal expansion properties • Exceptional formability • Outstanding durability

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。