空间桁架程序设计
第六章 空间桁架程序设计
第一节 概述
一、计算模型及计算
方法
快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载
1( 计算模型及整体坐标系下单元刚度矩阵的形成
选取计算模型时,应以杆简联结点和支座结点作为计算结点,各结点均为光滑的理想铰结点;
以任意两结点间的杆简为计算单元,各单元支承受轴力;非结点荷载要转化为等效结点荷载,
各单元抗拉压刚度相同。
由于结构为空间桁架,所以,结构各结点的位移分量为
T ,,,,,,,xyz
在局部坐标系下,单元的杆端位移列阵和杆端力列阵分别为
eeux,,,,ii,,,,
,,,,vyii,,,,ee,,,,,wz,,,,Fiiiiee,,,,,,,,,,,F,, (6-1) ,,,,,,,,ux,Fjjjj,,,,,,,,
,,,,vyjj,,,,
,,,,wzjj,,,,
式中u、v、w和u、v、w分别为结点i、j沿局部坐标系方向的线位移,x、y、ziiijjj
见图6-1。杆端力应与杆端位移一一对应,图中没再绘出。
当单元的杆端位移分量为任意值时,可写出空间桁架单元刚度方程。以矩阵
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示为 u,,X100,100,,,,eee,,,,,,ii,,,,vY,,000000
,,,,,,ii,,,,,,wZ000000EA,,,,,, (6-2) ,ii,,,,l,,uX,100100,,,,,,jj,,,,,,vY000000,,,,,,j,,j,,,,wZ000000,,,,,,
jj
eee简写成 (6-3) F,k,
式中 100,100,,e
,,
,,000000
,,
,,000000EA (6-4) ,,k,el,,,100100,,
,,000000,,
,,000000,,e称为单元?在局部坐标系下的刚度矩阵。
在空间桁架中,各杆方向不尽相同。每根杆件采用各自的局部坐标系,这对于单元
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
十分方便。而在结构整体分析中,所有的杆件都应该采用统一的坐标系,即整体坐标系。因
此在进行整体分析前,应将局部坐标系下的单元刚度矩阵转换为整体坐标系下的单元刚度矩
阵。
在整体坐标系中,单元杆件位移列阵和杆端力列阵分别为
eeuX,,,,ii,,,,
,,,,vYii,,,,ee,,,,,wFZ,,,,iiiiee,,,,,,,,,,,F,, ,,,,,,,,,uFXjjjj,,,,,,,,
,,,,vYjj,,,,
,,,,wZjj,,,,由坐标转换可知两种坐标系的位移分量存在如下关系:
ulll000u,,,,,,ixxxyxzi
,,,,,,
,,,,,,vlll000viyxyyyzi,,,,,,
,,,,,,wlll000wizxzyzzi,,,,,,, ,,,,,,u000lllujxxxyxzj,,,,,,
,,,,,,v000lllvjyxyyyzj,,,,,,
,,,,,,w000lllwjzxzyzzj,,,,,,在两种坐标系中,杆端力的转换关系为
eeeF,TF (6-11)
将式(6-9)与式(6-11)带入局部坐标系下单元刚度方程(6-3)并化简,可得
eeTeeeF,TkT,
artefactes (duresa en HRC30 dins) de fitxer clau 1 i cuc cuc engranar superfície no a un segon 2 i mandíbula graus no paral•leles pinça de cable 19 més tres cinquenes parts i nas agulla pinça, 1 i elèctric amb de vora portar liquidació més de 0,2 m/m, dent boca liquidació és superior a 0,4 m/m 2 i altres treballs del despatx vora sobre 0,4 m/m despatx de boca de dent és més gran que la pinça de canonada de 0.5-0.8M/m 20 1 i pinça mare greus desgast 2 i vareta de cargol, i Eix de rosca per lligar la liquidació és superior a 1-1.5M / m 21 Taiwan cargol 1, cargol canya i liquidació de rosca i la liquidació de distribució corredissa és superior a 1,5 a 2 m/m 2 i ontologia ha cruixit 3 i liquidació de la mandíbula és més gran que 0.3M / m 22 en directe els millors 1 i els millors radial batec és més gran que 0,07 m/m 2 i la part superior de la cua con part greu desgast 3 i que s'executen fracàs 23 tornavís 1 i vora duresa prou o alta 2 i nanses Departament esquerdes 24 Serra de metalls Marc 1 i abans i després responsable galeria trenat, és no un avió dins eina elèctrica 25 1, i Part elèctrica fracàs 2 i mecànica part danyat no pot córrer 26 roda 1 i restants diàmetre diàmetre original de gruix 2 i restant cèntim de gruix original de dos terços 3, segons tramitació situació roda consum força complir no processament requisits 27 1 de diamant i per la direcció de diamant acostar implementació 28 diversos mesura 1 i els plans de sector verificació de mesura d'aplicació individual de ferralla 29 altra eina 1 i per la situació específica mitjançant l'ús que va proposar fer supervisor de tecnologia eina validació 18 article màquina i proporciona eina de treball equipades
eeTee令 (6-12) k,TkT
eee则 (6-13) F,k,
这就是整体坐标系下单元刚度方程。将式(6-4)、式(6-10)带入式(6-12)化简,并
将分别用表示,则整体坐标系下空间桁架单元刚度矩l、l、lcos,、cos,、cos,xxxyxz
阵为
,,,,,,,,,,22,,coscoscoscoscos,cos,coscos,coscos,,22,,,,,,,,,,,,coscoscoscoscos,coscos,cos,coscos(6-14) ,,22,,,,,,,,,,,,coscoscoscoscos,coscos,coscos,cosEAe,,k,,,22,,l,,,,,,,,,,,cos,coscos,coscoscoscoscoscoscos,,22,,,,,,,,,,,,,coscos,cos,coscoscoscoscoscoscos,,,,22,cos,cos,cos,cos,,cos,cos,cos,cos,cos,cos,,,
2( 支承条件的引入及整体刚度矩阵得组集
整个结构的刚度矩阵式按照整体坐标系下单元刚度矩阵各元素的下标,采用直接刚度法组集成的。整体坐标系下结构刚度矩阵中得各元素,是由该坐标系下单元刚度矩阵中各元素按单元定位数组,“对号入座、同号相加”组集而成的。
与平面桁架的计算过程类似,由于在形成整体刚度矩阵之前,已经引入了支承条件,故此时的刚度矩阵不再是奇异矩阵。整个结构刚度方程
(6-15) K,,PFFF
表示的是自由结点位移与结点荷载的自由结点分量之间的关系。解方程组式(6-15)可直接的出资由结点位移,进而求出整个结构各结点的总位移列阵:
T (6-16) ,,,,,,FR
e?从结构的结点位移向量Δ中取出各单元两端的杆端位移分量δ,进而可求出局部坐标系下的单元杆端力
eeeeeeF,TF,Tk, (6-17)
第二节 空间桁架内力和位移计算的框图与程序
一、程序标识符说明
空间桁架静力分析程序名为STSAP(Space Truss Structural Analysis Program),其主要标
识符说明如下:
TITLE(20)——算力标题。整型数组,输入参数。
NJ——结点总数。整型变量,输入参数。
N——结构的自由度。即整体刚度矩阵的阶数。整型变量,输入参数。
NNE——单元总数。整型变量,输入参数。
NMT——单元类型数。同类型单元E、A相同。整型变量,输入参数。
NPJ——结点载荷数。整型变量,输入参数。
JN(3,100)——结点位移号组数。整型变量,输入参数。
X(100),Y(100),Z(100)——结点坐标数组,X(I)、Y(I)、Z(I)分别为I号结点的x坐
标、y坐标、z坐标。实型数组,输入参数。
JE(2,100)——单元两端结点号数组。整型数组,输入参数。
JEA(50)——单元类型数组。JEA(I)为I单元类型号。同类型的单元弹性模量、截面积相同。实型数组,输入参数。
EA(2,50)——各类型单元的物理、几何性质数组,EAI(1,I)、EAI(2,I)分别表示I号类型单元的弹性模量、截面积。实型数组,输入参数。
JPJ(100)——结点荷载位移数组,JPJ(I)为第I个结点相应位移分量的位移号。整型数组,输入参数。
PJ(100)——结点荷载数值数组。PJ(I)为第I各结点荷载的数值。实型数组,输入参数。
M(6)——单元定位数组,即单元两端位移号数组。整型数组,输入参数。 AK(200,200)——结构整体刚度矩阵,实型数组。
AKE(6,6)——整体坐标系下单元刚度矩阵。实型数组。
T(6,6)——存坐标转换矩阵或其转置矩阵。
P(200)——结点荷载,解方程后存结点位移。
FE(6)——存整体坐标系下的单元刚度矩阵与单元杆端位移的乘积。 FF(100)——单元轴力数组。
READ——子程序,读取数据。
CSL——子程序,计算单元常数。
STIFFN——子程序,计算整体坐标系下的单元刚度矩阵。
CALM——子程序,形成单元定位数组。
SLOV——子程序,解方程并求自由结点位移。
MULV6——子程序6 阶矩阵与6阶矩阵元素列阵相乘。
TRANS——子程序,计算坐标转换矩阵。
SQRT——标准函数,计算非负式数的平方根。
MAX0、MIN0——FORTRAN语言标准函数,分别为求整型量的最大值、最小值。
二、框图
空间桁架内力及位移计算总框图如图6-2所示。
artefactes (duresa en HRC30 dins) de fitxer clau 1 i cuc cuc engranar superfície no a un segon 2 i mandíbula graus no paral•leles pinça de cable 19 més tres cinquenes parts i nas agulla pinça, 1 i elèctric amb de vora portar liquidació més de 0,2 m/m, dent boca liquidació és superior a 0,4 m/m 2 i altres treballs del despatx vora sobre 0,4 m/m despatx de boca de dent és més gran que la pinça de canonada de 0.5-0.8M/m 20 1 i pinça mare greus desgast 2 i vareta de cargol, i Eix de rosca per lligar la liquidació és superior a 1-1.5M / m 21 Taiwan cargol 1, cargol canya i liquidació de rosca i la liquidació de distribució corredissa és superior a 1,5 a 2 m/m 2 i ontologia ha cruixit 3 i liquidació de la mandíbula és més gran que 0.3M / m 22 en directe els millors 1 i els millors radial batec és més gran que 0,07 m/m 2 i la part superior de la cua con part greu desgast 3 i que s'executen fracàs 23 tornavís 1 i vora duresa prou o alta 2 i nanses Departament esquerdes 24 Serra de metalls Marc 1 i abans i després responsable galeria trenat, és no un avió dins eina elèctrica 25 1, i Part elèctrica fracàs 2 i mecànica part danyat no pot córrer 26 roda 1 i restants diàmetre diàmetre original de gruix 2 i restant cèntim de gruix original de dos terços 3, segons tramitació situació roda consum força complir no processament requisits 27 1 de diamant i per la direcció de diamant acostar implementació 28 diversos mesura 1 i els plans de sector verificació de mesura d'aplicació individual de ferralla 29 altra eina 1 i per la situació específica mitjançant l'ús que va proposar fer supervisor de tecnologia eina validació 18 article màquina i proporciona eina de treball equipades
开始
(1)标题计数组说明 (7)子程序READ
(2)输入并打印原始数据 (8)子程序CSL
(3)组集整体刚度矩阵 (9)子程序STIFFN
(4)形成结点荷载列阵 (10)子程序CALM
(5)解方程并打印杆端位移 (11)子程序SLOV
(6)计算并打印杆端力 (12)子程序MULV6
(13)子程序TRANS 结束