资源描述
,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第7章,材料性质,MSC/NASTRAN 可处理多种材料性质,NASTRAN 可处理的,适于,线性静力分析,的材料类型,：,各向同性材料（MAT1）,二维各向异性材料（MAT2）,轴对称体正交异性材料（MAT3）,二维正交异性材料（MAT8）,三维各向异性材料（MAT9）,层复合材料PCOMP,各向同性材料（MAT1）,各向同性材料在各方向都具有同样的材料性质,典型应力应变曲线,当应力超过弹性极限，材料进入非线性，需用非线性分析方法,材料常数E、G、NU满足 。需提供E、G、NU中两个,质量密度RHO用于计算重力载荷及动力分析,热膨胀系数A和参考温度TREF仅用于热分析,结构阻尼GE不用于静力分析,MSC/NASTRAN中用MAT1卡描述，格式如下,名,称,内,容,MID,材料标识号（整数,0,）。,E,扬氏模量（实数,0.0,或空白）。,G,剪切模量（实数,0.0,或空白）。,MU,泊桑比（,-1.00,）。,Gij,材料性质矩阵,(,实数,),。,RHO,质量密度,(,实数,),。,Ai,热膨胀系数向量,(,实数,),。,TREF,参考温度,(,实数,),。,GE,结构单元阻尼系数,(,实数,),。,ST、,SC,、SS,分别为拉伸、压缩和剪切应力极限,(,实数，用于计算安全裕度,),。,MCSID,材料坐标系标识号（整数,0,，或空白）。,用PCOMP卡进行复合材料分析时，MAT2卡自动生成,轴对称体正交异性材料（MAT3）,轴对称体正交异性材料，应力应变关系,其中，为轴对称体横剖面坐标系,MAT3仅适用于CTRIAX6单元,为保证对称性，必须满足如下关系,MAT3,卡格式,名,称,内,容,MID,材料标识号（整数,0,）。,EX、,ETH,、,EZ,X、和,Z,方向扬氏模量（实数,0,）,MUXTH、,MUTHZ,、,MUZX,泊桑比（实数）。,RHO,质量密度（实数）。,GZX,剪切模量（实数,0.0,）。,AX、,ATH,、,AZ,热膨胀系数（实数）。,TREF,参考温度（实数）。,GE,阻尼系数（实数）。,二维正交异性材料（MAT8）,二维正交异性材料：,平面应力应变关系,横向应力横向应变关系,MAT8卡只适用于,板（壳）,单元，格式如下,名,称,内,容,MID,材料标识号（整数,0,）。,E1,纵向弹性模量,(,实数,),。,E2,横向弹性模量,(,实数,),。,MU12,泊桑比（,1,方向单轴载荷作用下的）,(,实数,),。,G12,面内剪切模量,(,实数,0.0),。,G1Z,1-Z平面受剪的剪切模量,(,实数,0.0,或空白,),。,G2Z,2-Z平面受剪的剪切模量,(,实数,0.0,或空白,),。,RHO,质量密度,(,实数,),。,Ai,i方向的热膨胀系数,(,实数,),。,TREF,参考温度,(,实数,),。,X,t,、X,c,分别为纵向拉伸、压缩时的允许应力或允许应变,(实数,0.0),。,Y,t,、Y,c,分别为横向拉伸、压缩时的允许应力或允许应变,(实数,0.0),。,S,允许的面内剪切应力或应变,(,实数,0.0),。,GE,结构阻尼系数,(,实数,),。,F12,TsaiWu理论中的交互项,(,实数,),。,STRN,最大应变理论要求的应力（或应变）识别码，若,STRN=1.0,，则X,t,、Y,t,、X,c,、Y,c,及各项输入应变值；若,STRN,为空白，则为应力值。,三维各向异性材料（MAT9）,三维各向异性材料，应力应变关系,MAT9卡格式如下,名,称,内,容,MID,材料标识号（整数,0,）。,Gij,材料坐标系中的对称材料性质矩阵元素,(,实数,),。,RHO,质量密度,(,实数,),。,Ai,热膨胀系数,(,实数,),。,TREF,参考温度,(,实数,),。,GE,结构阻尼系数,(,实数,),。,MAT9卡适用于体元CHEXA、CPENTA和CTETRA,层复合材料（PCOMP）,层复合材料，NASTRAN提供材料性质卡PCOMP，格式如下,名,称,内,容,PID,性质标识号（整数,0,）。,Z0,参考面至底面之距,(,实数,),。,NSM,单位面积非结构质量,(,实数,),。,SB,胶接材料允许剪应力,(,实数,0.0),。,FT,破坏准则识别码（,BCD,值）：,FT =,HILL,，,HILL,准则；FT =,HOFF,，,Hoffman,准则；FT =,TSNI,，,Tsai-Wu,准则；FT =,STRN,，最大应变破坏准,则。,TREF,参考温度,(,实数,),。,LAM,叠层排列方式识别码（,BCD,值）：,LAM =,SYM,，对称铺层，仅需输入半铺,层；LAM为空白，需输入全部层数据。,MIDi,不同层的标识号,ID,，各层是以底层为,1,号依次定义的。,Ti,各铺层的厚度,(,实数,),。,THETAi,每层纵向与单元材料轴的夹角,(,实数,),。,SOUTi,应力或应变输出请求（,YES,或,NO,）。,第8章,静力载荷,Nastran中，每一类载荷可以单独或以任何线性组合形式施加给结构,。,集中力和力矩,集中力和力矩直接施加给结点,集中力用FORCE、FORCE1和FORCE2卡定义,FORCE卡格式,概 述,名,称,内,容,SID,载荷集标识号（整数,0,）。,G,结点标识号（整数,0,）。,CID,坐标系标识号（整数,0,，缺省值,= 0,）。,F,比例系数（实数）。,Ni,集中力向量的分量，在,CID,坐标系中定义（实数，只少有一个,Ni,0,）。,例子：集中力F作用于悬臂梁自由端,自由域格式为：,或,FORCE , 100 , 2 , , 10. , 0. , -1. , 0.,或,FORCE , 100 , 2 , , 1. , 0. , -10. , 0,FORCE1和FORCE2卡定义集中力,方向：FORCE1用两个结点的连线,FORCE2用四个结点组成两个向量（G1-G2，G3-G4）的向量积，,格式 （下面force1、force2位置有错，刚好相反）,名,称,内,容,SID,载荷集标识号（整数,0,）。,G,结点标识号（整数,0,）。,F,力的系数（实数）。,Gi,定义力方向的结点标识号（整数,0,）。,作用于结点的集中力矩，用MOMENT卡定义，格式,SID,载荷集标识号（整数,0,）。,G,为力矩作用的结点标识号（整数,0,）。,CID,坐标系标识号（整数,0,，或空白）。,M,比例系数（实数）。,Ni,力矩向量分量，在由,CID,定义的坐标系中量度（实数，只少有一个,Ni,0,）。,MOMENT1和MOMENT2也可定义集中力矩，同FORCE1和FORCE2。,分布载荷,作,用于一维单元上的分布载荷,（,PLOAD1,）,用PLOAD1卡对一维单元（CBAR、CBEAM和CBEND）施加分布载荷；,对CBAR和CBEAM单元，分布载荷可沿单元全长或部分长度来施加；,对CBEND单元，分布载荷限沿单元全长线性变化,PLOAD1卡的格式如下,名,称,内,容,SID,载荷集标识号（整数,0,）。,EID,CBAR、,CBEAM,或,CBEND,单元的识别号（整数,0,）。,TYPE,载荷类型。用如下字符表示：,“,FX,”、“,FY,”或“,FZ,”：分别为基本坐标系中,X,、,Y,或,Z,方向的力；,“,MX,”、“,MY,”或“,MZ,”：,分别为单元坐标系中,X,、,Y,或,Z,方向的力矩；,“,MXE,”、“,MYE,”或“,MZE,”：,分别为单元坐标系中,X,、,Y,或,Z,方向的力矩；,SCALE,为,X1,、,X2,确定比例系数，用如下字符表示：,“,LE,”（实际长度）,，,Xi,值是沿单元轴的实际距离，若,X1,X2,，,Pi,值是单元每单位长度的载荷密度；,“,FR,”（比例长度）,，,Xi,值是沿单元轴距离与单元总长度之比，如果,X1,X2,，,Pi,值是单元每单位长度的载荷密度；,“,LEPR,”（投影长度）,，,Xi,值是沿单元轴的实际距离；,“,FRPR,”,（比例投影长度,），,X,i,值为沿单元轴实际距离与总长之比，,Pi,值为单元每单位投影长度载荷密度；,X1，,X2,从梁（,CBAR,、,CBEAM,或,CBEND,）端,A,起算的沿梁轴线之距（实数，,X2,可为空白，）。,P1，,P2,分别为,X1,和,X2,处的载荷系数（实数，或空白）。若,X1,X2,，则在,X1,和,X2,之间载荷呈线性分布，而,X1,，,X2,处之单位长度载荷分别为,P1,及,P2,；若,X1=X2,，或,X2,为空白，表示在,X1,处作用一大小为,P1,的集中力。,例1 均布载荷,用比例长度“FR”来确定X1和X2，有X1 = 0.0, X2 = 1.0，P1= P2=12.6磅/英寸,用“,LE,”来确定,X1,和,X2,例2，线性分布载荷,对沿梁轴线线性分布载荷,PLOAD1,卡为：,例3 ，集中载荷,对作用于梁上的集中载荷,PLOAD1卡为,作用于二维单元上的均布压力（PLOAD2）,对有相同均布压力的多个二维单元（CQUAD4或CTRIA3），PLOAD2卡十分方便PLOAD2卡格式：,或,名称,内 容,SID,载荷集识别号（整数0），由情况控制指令LOAD = SID选取。,P,压力值（实数），压力方向按单元连接结点顺序，以右手定则定义。,EIDi,单元识别号（整数0，或空白，对于替换形式，EID2EID1,并且要求其间的所有单元号是实际存在的）。,PLOAD,名称,内,容,SID,载荷集标识号（整数,0,）。,Gi,结点号,P,压力（实数），方向按单元结点连接顺序，以右手定则确定。,作,用于二维或三维单元面上的分布压力,（,PLOAD4,）,PLOAD4定义多种二维或三维单元面上的压力载荷,这种压力载荷可垂直或不垂直于单元面,可在单元面各角点输入不同压力值,PLOAD4,卡的格式如下,替换格式（仅适于面单元）,名称,内 容,SID,载荷集标识号（整数0）。,单元标识号（整数0，对于“THRU”选择，EID10，或为空白）。,G3,与G1结点同一面上的对角结点编号，仅对于体单元CHEXA和CPENTA才要求的。对于CPENTA单元的三角面是不要求的。,G4,四面体元CTETRA不受压力的角结点编号（整数0）。,CID,坐标系标识号（整数0，缺省值= 0）。,N1，N2，N3,在坐标系CID中定义的向量分量，用于定义载荷密度向量方向（实数）。,例,1,：,作,用于曲板上的均布压力,（,PLOAD4,）,用,PLOAD4,卡替换形式定义,例,2,：,作,用于六面体,CHEXA,单元上的均布压力载荷,重力和离心力（GRAV，RFORCE）,在NASTRAN中，重力载荷用GRAV卡来施加，格式,名称,内,容,SID,重力载荷集标识号（整数,0,）。,CID,坐标系标识号（整数,0,，缺省值为零）。,G,重力加速度向量比例系数（实数）。,Ni,在坐标系,CID,中定义的加速度向量分量（实数，只少有一个,Ni,0,）。,旋转引起的离心惯性力，用,RFORCE,卡定义。格式,名称,内,容,SID,离心载荷集标识号（整数,0,，为零时代表基本坐标系原点）。,G,结点标识号（整数,0,）。,CID,定义选转向量的坐标系标识号（整数,0,，缺省值为,0,）。,A,R1,R2,R3,单位时间旋转角速度标量系数（实数）。,旋转向量的分量，该向量过,G,点（实数）。,METHOD,计算离心力所用的方法（整数,= 1,，,2,，或空白）。,METHOD=1（或空白），耦合质量矩阵；,METHOD=2，集中质量矩阵,强迫位移,静力分析中，用DEFORM卡定义,强迫变形,DEFORM卡在情况控制集中用指令DEFORM = SID来选取,DEFORM,卡格式,名称,内,容,SID,强迫变形集标识号（整数,0,）。,EIDi,单元标识号（整数,0,）。,Di,变形（实数，正值表示伸长）。,对静力分析中的,强迫位移,，用,SPCD,施加,SPCD卡由情况控制指令LOAD = SID来选取,DEFORM格式,名称,内 容,SID,载荷集识别号（整数0）。,Gi,结点或标量点标识号（整数0）。,Ci,位移分量号（0整数6，1至6个的唯一整数，无嵌入空白）。,Di,与Gi和Ci对应的强迫位移值（实数）。,热 载,温度场用结点温度和单元温度数据定义,（ TEMP和TEMPD模型数据卡定义结点温度,TEMP卡,格式,名称,内 容,SID,温度集标识号（整数0）。,G,结点标识号（整数0）。,T,温度值（实数）。,一张TEMP卡最多可定义三个结点的温度，由情况控制TEMP = SID选取,TEMPD,卡,定义结点温度场，格式,名称,内,容,SID,载荷集标识号（整数,0,）。,T,结点温度原设定值（实数）。,该卡最多可定义四组温度场，它由情况控制指令,TEMP = SID,选取。,对于,一维单元,：,ROD,、,BAR,、,BEAM,、,BEND,和,CONROD,或,TUBE,，用,TEMPRB,卡来定义温度场，格式,：,第二继续卡的替换形式,名称,内 容,SID,温度集标识号（整数0）。,EIDn,单元标识号。,TA、TB,分别为在端点A和B剖面之平均温度（实数）。,在端点j处i方向的等效线性梯度（只用于CBAR、CBEAM和CBEND，实数）。,端点j剖面上i位置的温度（i位置由PBAR、PBEAM或PBEND卡所定义，只用于应力恢复，实数）。,二维单元,，用,TEMPP1,和,TEMPP3,卡来定义温度场，,TEMPP1,卡,定义二维单元面上温度场，格式,继续卡的替换格式为,名称,内 容,SID,温度集标识号（整数0）。,EIDn,单元标识号（整数0）。,T,平均温度（实数）。,等效温度梯度（实数，不用于膜单元）。,T1、T2,单元性质卡PSHELL上规定的计算应力点的温度（实数）。,TEMPP3,卡,定义,二维单元,剖面,温度梯度,，格式,第三继续卡的替换格式为,名称,内 容,SID,温度集标识号（整数0）。,EIDn,单元标识号（整数0）。,Z0,底面至参考面之距（实数）。,Zi,平面i至参考面之距（实数）。,T0,底面温度（实数）。,Ti,平面i的温度（实数）。,组合载荷,LOAD,模型数据卡可用来进行静力载荷的线性组合（叠加）,所组合的载荷是用,FORCE,、,MOMENT,、,FORCE1,、,MOMENT1,、,FORCE2,、,MOMENT2,、,PLOAD,、,PLOAD1,、,PLOAD2,、,PLOAD3,、,PLOAD4,、,PLOADX,、,SLOAD,、,SPCD,、,RFORCE,或,GRAV,卡来定义,LOAD,卡格式：,名称,内 容,SID,载荷集标识号（整数0）。,S,总比例系数（实数）。,Si,分比例系数（实数）。,Li,可组合载荷卡上定义的载荷集标识号（整数0）。,由,LOAD,定义的组合载荷为,例子,:,(1),Y,向,15.2,磅的集中力,F,作用于结点,12,上,(2),6.4,英寸,-,磅关于,X-,轴的集中弯矩,M,作用于结点,127,上。,要求,:,组合载荷为：,2,倍集中力,F,和,3,倍的集中力矩,M,。即,式中,30,和,40,分别为集中力,F,和集中力矩,M,的载荷集标识号。,在,情况控制集中,：,LOAD = 2 2,在,模型数据集,中,注意：,NASTRAN,可用情况控制,SUBCASE,在一次程序运行中分析多种载荷工况（每一种子情况定义一种唯一的载荷工况）,还可用情况控制指令,SUBCOM,和,SUBSEQ,定义子情况（不同载荷条件）的线性组合。,例子,：说明,SUBCOM,和,SUBSEQ,应用,F1,2,4,3,1,5,Z,F2,X,第一种载荷情况：同时加上F1和F2；,第二种载荷情况：同时加上F1和半倍的反向F2。,在情况控制集中：,SUBCASE 1,（选取垂直载荷F1）,LABEL = VERTLCAL LOAD 2000.0 LBS,LOAD = 10,SUBCASE 2,（选取水平载荷F2）,LABEL = HORIZONTAL LOAD 1000.0 LBS,LOAD = 20,SUBCASE 100,（组合F1和F2）,SUBTITLE = COMBINATION OF THE PREVIOUS SUBCASES,LABEL = FIST COMBINATION,SUBSEQ = 1. , 1.,SUBCASE 200,（组合F1和,-0.5F2,）,LABEL = SECOND COMBINATION,SUBSEQ = 1. , -0.5,在模型数据集中：,FORCE , 10 , 3 , , -2000. , 0. , 0. , 1.,FORCE , 20 , 5 , , 1000. , 1. , 0. , 1.,第9章,约束处理,MSC/NASTRAN两种基本约束类型：,单点约束（SPC）和多点约束（MPC）,单点约束,单点约束是约束单个自由度的一种约束,用途：,支持一个结构，或将一结构“固定”于地面；,处理对称或反对称边界条件；,消除结构分析中未用自由度；,消除非常弱的耦合自由度,；,给结点施加零或非零强迫位移,永久性固定约束,对,永久性固定约束，用GRID模型数据卡中第8域PS来定义;,如果分析模型中仅有几个结点需施加零位移约束，则采用GRID卡。,单点约束（,SPC,）,用,SPC,模型数据卡施加一组单点约束或,强迫位移,SPC,的格式,名 称,内 容,SID,单点约束集标识号（整数0）。,Gi,结点标识号。,Ci,位移分量号（1整数6，可填由1至6组成的多至6个唯一的整数，其间不允许堪入空格）。,Di,由Gi 和Ci指定的自由度强迫位移值（实数，缺省值为0.0）。,SID是SPC集识别号，由情况控制指令集中SPC=SID 指令来选取,一个SPC卡可定义两组（G、C、D）值,单点约束（SPC1）,SPC1模型数据卡用于定义,零位移单点约束,，格式,替换格式,名 称,内 容,SID,单点约束集标识号（整数0）。,C,分量号（任何唯一的由1至6的整数组合，不允许堪入空格。对于结点必须为整数；标量点则为空白）。,Gi,结点或标量点标识号（整数0或为替代格式中的“THRU”，要求G1G2）。,SID是由情况控制指令集中SPC = SID 指令来选取，,任何数量的SPC1卡可用来定义一个约束集。,AUTOSPC,在MSC/NASTRAN中，利用在模型数据集中加入参数卡：,PARAM，AUTOSPC，YES,可自动识别与消除刚度矩阵的奇异性,消除未用自由度和非常弱的耦合自由度,边界条件例题,1,一根梁，边界条件为一端固支、另一端铰支。具有5个结点，4个梁元,对于结点1需约束全部6个自由度，三个移动（1，2，3）和三个转动（4，5，6）。对于结点5，需约束三个移动自由度（1，2，3）和两个转动自由度（4， 6）。,用SPC卡，假定情况控制约束集中已选定SPC=100。强迫位移的值为零，取缺省值。,SPC卡为,用自由域格式，则为：,SPC，100，1，123456， ，5，12346,边界条件例题,2,一四边固支板，四个板元，,9,个结点。,边界结点1，2，3，4，6，7，8和9必须约束所有6个自由度（1，2，3，4，5，6）。假定在情况控制约束集已选定SPC=100。,采用SPC1卡，则有,多点约束,NASTRAN,中，多点约束（,MPC,）描述两个或多个自由度间的线性关系：,式中 u,j,为节点或标量点的任何自由度；,R,j,为用户定义的比例系数。,多点约束用途：,(1) 将两个结点间的相关运动定义为一个自由度；,(2) 将几个运动平均值定义为一个自由度；,(3) 在结构部件间提供铰链或滑动连接；,(4) 连接不同种类的单元，将具有转动自由度的单元与仅有移动自由度的单元相连接（壳元与体元的连接就是一例）；,(5) 获得作用于结构或结构部件上的合力；,(6) 将一个力分配到结构中几个点。当一个力是未知的时候，这是实际有用的，例如，要求压缩流体的力；,(7) 连接具有不一致结点的单元，例如，在结构内部改变网格尺寸便会出现这种情况；,(8) 替换具有刚性连接的极刚硬结构元件。这只当需改善刚度矩阵的数值状态时才用，且是不推荐的，最好采用刚性元来处理；,(9)定义结点的运动分量，该分量方向不是与结点的局部坐标系轴向一致,多点约束用,MPC,卡,描述，其格式如下,域,内 容,SID,多点约束集标识号，由情况控制指令选定。,Gi,结点或标量点识别号。,Ci,自由度分量号。,Ai,比例系数。,第10章,线性静力分析,基本有限元方程,静力分析的基本有限元方程为：,Ku=P,式中 K为结构的弹性刚度矩阵；,u为结点广义位移向量；,P为结点载荷向量。,刚度矩阵按结点自由度装配而成的，对应,位移集称之为g-集；,该矩阵往往是可能奇异的，不能直接被分解的。一般需进行如下运算：,多点约束减缩（MPC），消去线性相关自由度（可选的）；,单点约束减缩（SPC），消去刚体运动自由度；,静力减缩（OMIT），减小求解问题的规模（可选）,自由体支持（SUPORT），分析惯性卸载问题,上面运算由一个或多个模型数据卡所控制：,运,算,模,型,数,据,卡,多点约束,MPC，MPCADD，RROD，RBAR，RTRPLT，RBE1，RBE2，RBE3，RSPLINE,单点约束,SPC，,SPC1,，,SPCADD,，,GRID,，,GRDSET,，,PARAM,，,AUTOSPC,，,YES,静力减缩,OMIT，,OMIT1,，,ASET,，,ASET1,自由体支持,SUPORT,解题过程,确定物理模型,简单物理模型：承受集中载荷作用简支梁。,材料特性：弹性模量 E=30E6,泊桑比 v=0.3,载荷： 承受集中力F，见图,边界条件：两端简支。,求解：确定在载荷F作用点处梁的挠度和应力，不考虑横向剪切的影响。,建立有限元模型的输入数据文件,按MSC/NASTRAN的输入数据格式，建立相应有限元模形的数据文件,第1步,确定分析类型，形成输入文件的执行控制语句集,ID MPM，CH 12 EXAMPLE,SOL 101,TIME 100,CEND,第,2,步,离散物理模型以形成有限元模型,对承受集中力梁弯曲问题。因杆单元CROD不能承弯，故不能选用它。,CBAR是常用的等剖面直梁元，满足本问题的要求。,最后形成具有四个结点、三个CBAR单元、在结点3处作用集中力F的两端铰支梁模型,第,3,步,形成模型数据集,建立有限元模型的“几何”、“单元”、“材料”、“载荷”与“约束”五类数据,几何数据,GRID,卡如下：,若采用自由域格式,则为,:,GRID , 1 , 0. , 0. , 0.,GRID , 2 ,10. , 0. , 0.,GRID , 3 , 20. , 0. , 0.,GRID , 4 , 30. , 0. , 0.,单元数据,CBAR单元卡中第6、7、8项为梁端点GA处的定位向量分量（X1，X2，X3）；,该向量确定梁元坐标系和主平面1的位置，据此计算梁剖面性质。,尽管的选择有点任意，但其方向必须与梁剖面的惯性主平面相一致。,定义平面1和Y-单元轴,单元性质卡标识号为101，三个梁单元的CBAR卡数据,CBAR单元的性质PBAR卡：,a) PBAR卡第5项和第6项。,两者不能颠倒。,应力恢复系数Ci、Di、Ei和Fi的选取，由用户在梁剖面（Y-Z平面）上任选的四点局部坐标值。,考虑横向剪切影响，需剪力面积系数K1和K2填入适当的值，对矩形剖面，K1=K2=5/6。若不填，其缺省值为无穷大，意味着剪切对挠度无影响。,对于无剪切影响，PBAR卡为,材料性质,梁的材料为钢材，弹性模量，泊桑比。采用,MAT1,卡来定义，材料标识号选为,201,载荷数据,结点3处有100磅集中力，方向向下。用FORCE卡定义。选用载荷标识号选为10，,该FORCE卡如下,边界条件,（,1）可用GRID卡第8项（PS）填入固定约束自由度分量。,（2）考虑到1和4结点约束相同，采用SPC1卡描述。SPC1卡,第4步,提出分析输出要求，建立情况控制指令集,位移和单元应力输出，如下的情况控制指令是必须的,DISP = ALL（输出全部结点位移）,STRESS = ALL（打引所有单元的应力）,为了,检查计算结果,，采用如下指令是合适的：,FORCE = ALL（打引所有单元力）,SPCF = ALL（打引所有约束反力）,输入数据文件的输出采用如下指令：,FCHO= BOTH （输入文件以分类和未分类两种形式输出）,输出结果页面的标题和子标题可用,TITLE和SUBTITLE指令：,TITLE = HINGED BEAM,SUBTITLE = WITH CONCENTRATED FORCE,对载荷和约束的选取如下：,标识号应与所选模型数据卡的标号相对应,最后，形成了如下情况控制指令集：,CEND,ECHO = BOTH,DISP = ALL,STRESS = ALL,FORCE = ALL,SPCF = ALL,SPC = 100,LOAD = 10,TITLE = HINGED BEAM,SUBTITLE = WITH CONCENTRATED FORCE,上述情况控知指令要求放在,CEND,之后，相互的顺序是任意的。,第5步,完成输入文件,运行MSC/NASTRAN并输出结果,MSC/NASTRAN,的计算结果放在,.F06,文件中,。,检验结果,错误信息,MSC/NASTRAN,对有限元模型、运算过程错误信息，记录在文件,.F06,中。,残余力向量误差,残余力向量误差,EPSILON,是十分小的（,10,-6,），这表明稳定的数字特性，不存在大的舍入误差，所求解的刚度矩阵特性好（不存在奇异性、病态条件）。,位移量级判断,梁的挠度量级为英寸，这同梁的长度与剖面尺寸相比均是很小的量，符合小位移假设，结果是合理的。如果得到大于几英寸的挠度结果，,载荷太大，或是材料特性数据（弹性模量）不正确，或该问题根本不是线性问题。,平衡检查,为检验静力平衡，需计算支持点的反力,与经典理论结果的比较,
展开阅读全文