4杆系结构单元解决方案.ppt

例5-1 一维拉杆 图示阶梯形直杆，各段长度均为，横截面积分别为3A，2A，A，材料重度为γ，弹性模量E。求结点位移和各段杆中内力。 离散化：将单元划分为3个单元，4个结点。 单元刚度矩阵： 等效结点荷载：按静力等效原则，有: 整体结构平衡方程 例5-4 变截面梁 有一变截面梁，一端固定，另一端铰支。梁长为2l，固支端的截面尽寸为b×1.6h，铰支端的截面尺寸为b×h。梁上作用均布载荷p0。求梁端的约束反力。 离散化 将梁划分成2个单元，3个结点。每个单元 长度为,截面取平均截面。 对号入座，组合整体刚度矩阵 荷载等效结点力向量 引入边界条件 例5-2 两根桁架 两根杆件的横截面积为A，弹性模量为E，垂直杆长为，两杆铰接处受到水平方向的外力P。求结点位移和杆中的内力。 ①单元：i=1，j=2，α=45°,单刚为： ②单元： i=2，j=3，α=-90°,单刚为： （4）结构整体平衡方程 （5）引入约束 1、3结点约束，划去1、2、5、6行与列，得 （7）单元内力 整体节点位移变换到单元节点位移。 （7）单元内力 对于①单元，i=1,j=2,杆长 （7）结点力量 （利用整体平衡方程） 将结点位移分量代入整体平衡方程，得结点力量 5.7 桁架结构有限元程序（TRUSS ） 1、程序使用步骤： 1） 输入程序 2）编译调试，无错 4）建立输入数据文件（文件名：T1.DAT） 5）在程序中输入节点力 6）编译执行得计算结果（文件名：OUTT1.，可改变） 注：每个算例都需修改程序，数组说明和节点力输入。 5.7 桁架结构有限元程序（TRUSS ） 2、输入数据格式（文件名：T1.DAT） 1）总信息：NN，NE，ND，NFIX，E 结 点 数；单元数；总自由度数；约束总数；弹性模量 2）单元节点编号i，j，A： LOC（I,1），LOC（I,2），AREA（I）,单元顺序 3）节点坐标x，y：CX（J），CY（J），节点顺序 4）约束自由度号：IFIX（K）,结点顺序(x向在前，y向在后) 5.7 桁架结构有限元程序（TRUSS ） 3、载荷信息 在程序中输入节点力F(i)。 4、输出数据（文件名：OUTT1.） 1）总信息：NN，NE，ND，NFIX，E 2）结点号和结点力：I，F(2*I-1)，F(2*I) 3）结点号和结点位移：I，F(2*I-1),F(2*I) （求解后F存入结点位移） 5.7 桁架结构有限元程序（TRUSS ） 5、 TRUSS主程序 数组大小、等价节点力（程序中输入）、输入输出文件名 5.7 桁架结构有限元程序（TRUSS ） 5.7 桁架结构有限元程序（TRUSS ） 6、桁架算例 6、桁架算例 单位统一 :cm，N 建立数据文件T1.DAT 在程序中输入载荷，在结点1的Y向载荷为-10kN。 5.8 刚架有限元程序（FRAME） 1、程序使用步骤： 1） 输入程序 2）编译调试，无错 4）建立输入数据文件（文件名：F1.DAT） 5）在程序中输入节点力 6）编译执行得计算结果（文件名：OUTF1.，可改变） 注1：根据算例情况需修改数组大小。 注2：集中载荷处、变截面处均设结点，集中载荷按结构坐标系输入； 注3：程序中仅考虑了刚架单元作用均布横向载荷情形。 5.8 刚架有限元程序（FRAME） 2、输入数据格式（文件名：T1.DAT） 1）总信息：NN，NE，ND，NFIX，E 结 点 数；单元数；总自由度数(3N）；约束总数；弹性模量 2）单元i，j，A，I：单元顺序 LOC（I,1），LOC（I,2），AREA（I）,FINT(i) 3）节点坐标x，y：CX（J），CY（J），节点顺序 4）约束自由度号：IFIX（K）,结点顺序(x向在前，y向在后) 5）NP，NVD：集中载荷个数和横向均布力载荷组数。 6）如NP≠0,顺序输入NP组数据： I ， F（3*I-1）， F（3*I-2）， F（3*I） 结点号 ，集中力向分量，集中力向分量 ，弯矩 7）如NVD≠0，顺序输入NVD组数据：IE ，Q 单元号 ，均布载荷值 5.8 刚架有限元程序（FRAME） 3、输出数据（文件名：OUTF1.） 1）总 信 息：NN，NE，ND，NFIX，E 2）单元信息