有限元基础教学-3杆系有限元.pptVIP

2006－7－8 哈尔滨工业大学 王焕定 杆系结构单元分析子程序 pi=ACOS(-one); gamrad=gamma(iel)*pi/d180; cg=COS(gamrad) sg=SIN(gamrad); den=ell*SQRT(xl*xl+zl*zl) IF(den/=zero)THEN r0(1,1)=xl/ell; r0(1,2)=yl/ell; r0(1,3)=zl/ell r0(2,1)=(-xl*yl*cg-ell*zl*sg)/den; r0(2,2)=den*cg/(ell*ell) r0(2,3)=(-yl*zl*cg+ell*xl*sg)/den; r0(3,1)=(xl*yl*sg-ell*zl*cg)/den r0(3,2)=-den*sg/(ell*ell); r0(3,3)=(yl*zl*sg+ell*xl*cg)/den ELSE r0(1,1)=zero; r0(1,3)=zero; r0(2,2)=zero; r0(3,2)=zero r0(1,2)=one; r0(2,1)=-cg; r0(3,3)=cg; r0(2,3)=sg; r0(3,1)=sg END IF DO i=1,3; DO j=1,3 ; x=r0(i,j); DO k=0,9,3 t(i+k,j+k)=x; tt(j+k,i+k)=x END DO; END DO; END DO DO i=1,12; DO j=1,12; sum=zero; DO k=1,12 sum=sum+km(i,k)*t(k,j) END DO; cc(i,j)=sum; END DO; END DO DO i=1,12; DO j=1,12; sum=zero; DO k=1,12 sum=sum+tt(i,k)*cc(k,j) END DO; km(i,j)=sum; END DO; END DO END SELECT; RETURN; END SUBROUTINE rigid_jointed 关于三维杆定位向量的说明 γ的定义——转动坐标。 1）“竖直”单元：单元长轴平行于整体坐标y轴. γ定义为整体坐标轴z和局部坐标轴z′之间的夹角（以指向整体坐标轴x为正）,如下图所示； X Y, x′ Z z′ y′ γ 关于三维杆定位向量的说明 γ的定义——转动坐标。 2）非“竖直”单元：单元绕局部坐标轴x′转动的角度（从x′的正方向看）。 Y Z z′ y′ γ * * 杆系结构有限元分析 最基本的概念都在这里，必须下功夫学好 杆系结构单元分析 单元位移模式的建立方法 单元刚度方程的建立方法 单元分析的总结及物理实质 杆系结构单元分析子程序 单元位移模式的建立方法 以等直拉压杆为例 结构中拆出的单元如图所示。 广义坐标法 设任意点位移为 u=a1+ a2x 广义坐标，边界条件只两个 幂级数简单 右手系 i j x y 1 2 u1,F1 u2,F2 p EA，l 利用边界条件可得 a1＝ u1； a2＝(u2- u1)／l 将广义坐标代回 u=a1+ a2x，整理后可得 u=(1-x／l)u1+ u2x／l 形函数及性质 形函数 自然坐标 本点处为1；它点处为0；?处总和为1 任意点的位移可用形函数表为 u=(1-x／l)u1+ u2x／l=Ｎ1u1＋Ｎ２u２ 单元位移模式的建立方法 以等直杆扭转为例 结构中拆出的单元如图所示。 试凑法 设任意点自然坐标为 ? ,为满足“本1，它0” 可设 N1=1-? ，N2= ? 。 由性质试凑得到 右手系 i j x y 1 2 ?1,M1 ?2,M2 m GJ，l 单元位移模式的建立方法 ? = N1 ? 1+ N2 ? 2 任意点的位移可用形函数表为 以等直弯曲杆单元为例 设杆端位移（不计轴向变形） 则杆中任意点挠度可表示为 由此可见， N1是图示v1=1的挠曲线 单元位移模式的建立方法 因此可由如下微分方程积分获得 同理，可得 单元位移模式的建立方法 积分二次，利用边界条件