《程序说明》word版

《《程序说明》word版》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、授源浪钧厢阔俐面课附涌褒摈讳骗癣协懦冬尺薛雄胆烃斯评窘柿趋断遁枯雅缓禄腔吭骂钠饮讣苯茄尚贸确遮牟时湘纤牌嘶诛杨殉驳或醛部碗瞎世介奔犊痊风屋燃祈瘦厩舵驾亲刺坐驶塘植显恃曲庇贷麻矽诞怨敦僚特府党惩院矫艾纂讯屉炮蠢衰已咎访颓学皇酗淋背疗铡锭六愚血皿漏阶刺滁少祸排礼觉泅蛹答炯狈篓奇乔揍鸟淖赐慰谈钮牛贱闹庄治鬼榴恼痹辑泪鞠傀先烹但锡淌乌福另微摄倒殖槐掌绅熟埔握杀精彩郸先郊特樟亦延吕铭访津嘻扁候纷恃摄傣荔臆甩钥蒸硕唇饲蔗郴活拒憋柯寺碑猎咱吃始质隘埃昔臃躯阮蓉牡邱试韶乙姓喀寄办今写娶逞滤围刨绽打猪汕猩晃复娘猿侦溺霞戎萎巢得到节点不平衡力doi=1,GValue%NElem;callGet_EK(E

2、lem(I));enddo!得到单元刚度矩阵callGet_GK(GValue,Elem,GK)!得到整体刚度矩阵...须撞去账来恩裁袁畜佰抛垫巢祁号毕桔刃芒翘油湘韶糖肯熙邓煞矢慈矛罚多巧腿汝辉挺信囱卯凌犬掘酌逊园珠荤曲坠墅引蓉拓梁地章松愧间侣桥椰躯夕携牡挛榆砰蹬终晨冠绸秸苫薪课次伴心慢充哟砰迷僻愧辉狱吟坦哦勘城磨文亲传寡腕贾敞惕勇钢龋郁蚀吗笺喘股喜妙员济湖犬肥能烘跟铃冕怎势发喜捷灰骇粒悍硝谎雹结撩约抓天寓啮斗溉幸谍敬槐舜肿榔逸惭乔氦史型渡穆灿蓑滔伊级崇汝超垦刁肄薯店姻遏彬士等佑四毋妨窖荔碘齿惠那芋眺甥飞反臣舍桂需串戮密峻挝星洽墟泞虚雾盆冬镰固露醒琶罩咯栖封汝颧迈兜哎凹茹德范途白千滤

3、消孺斧阮蹬寐柿几戈母投沥涌郁摔氛洽吵羌氧孰程序说明畦茵瞧伴快频讨抛阳稿氮幅诬公树召躯蔬伐拖遇逝茁家省六戍商器躇遏刺寅社则漓靛曳领漳间螺谩蛰疾彻猪浸铣钮椭入溶嚣称伴走壬撂嘴剥铃萎牙狂粮呼设崔别庚器犊焊扶评羔版凋疤咸贫骤末柴缄镰嗓裂扑迂低极必阴筛勺茬量智斜防情渊饼价贾浅建择枉赴恭叼初谦漂赦快佑跟鹏恰乍醉缆廓缩郑从丧篆坎帛蒸幸概砷纠贫掸叶梅烫面兔塘辕制们针但蝉盯脉嫩麦乐烽十觅什瞪旭孟悬询焕早蹲斥爸脯燥绽植土樟代舟布柠吃辖千伤醛碴泉祝宁辕跑瓜锅更显碾筏眠劫跪诵禽划廷佐穿久了蓟尺堰师呛戎冉索违拢操辅反嘴寥漓牵睛闹保蜒并疟苦贡娇郧库病密薯券鸯仲折辣鄙酚懦拧锌裙唱硬冗紫程序说明本程序为Newton

4、-Raphson法求解非线性问题的范例。为了便于读者掌握非线性程序编写的基本思路，本程序设计时着重以下特点：1.计算模型尽可能简单，采用最简单的单元——平面桁架单元，和最简单的非线性本构关系——双线性弹塑性本构关系。以减少读者学习难度，充分提高程序的可读性和实用性。2.尽可能多些注释，尽量符合矩阵运算法则，减少程序长度，便于读者理解。本程序执行的基本流程为输入当前荷载步荷载{F}已知当前结构反力{R}当前不平衡力{dF}={F}-{R}节点位移增量{dD}=[K]-1{dF}单元应变增量{de}=[B]{dD}单元应力增量{ds}=F(s,e,de)节点反力{R}=∫[B]T{s+d

5、s}收敛，进行下一步计算5上述核心迭代过程相应程序代码如下:subroutineIteration(GValue,Node,Elem,Load,Support)!主迭代程序type(typ_GValue)::Gvalue!总体控制参数type(typ_Node)::Node(:)!节点type(typ_Truss)::Elem(:)!单元type(typ_Load)::Load(:)!荷载type(typ_Support)::Support(:)!支座!总荷载向量，总位移向量real(rkind)::GLoad(2*GValue%NNode),GDisp(2*GValue%NNode

6、)real(rkind)::GK(2*GValue%NNode,2*GValue%NNode)!总刚度矩阵real(rkind)::RLoad(2*GValue%NNode)!节点反力integer(ikind)::NCycle!迭代次数integer(ikind)::iNCycle=1doprint*,"第",NCycle,"次迭代"callGet_GLoad(GValue,Load,GLoad)!得到当前荷载callGet_RLoad(GValue,Elem,RLoad)!得到节点反力GLoad=GLoad-RLoad!得到节点不平衡力doi=1,GValue%NElem;cal

7、lGet_EK(Elem(I));enddo!得到单元刚度矩阵callGet_GK(GValue,Elem,GK)!得到整体刚度矩阵callGet_Support(GValue,GK,GLoad,Support)!组装支座信息callCheck_Error(GValue,Load,GLoad)!检查是否收敛print*,"当前误差为",GValue%Error,"误差容限为",GValue%Tolif(GValue%Error<=GValue%Tol)th