伽辽金法在悬臂梁大挠度问题中的应用
应用无网格伽辽金方法分析结构大变形问题
5 影响域
影响域是影响最小二乘精度和计算量的另一个重要参数，在移动最小二乘法的收敛性 证明中，当节点的密度趋向于无穷大的时候，影响半径应该趋向于零。也就是说，影响半 径是和节点密度密切相关的。而且，其选取对结果影响非常明显，尤其是对应力的影响。
结构工程师,Structural Engineers,2003. 66. 增刊/Sup,18~22
因此，式（1）可改写为
u h ( x) = ∑ ni ( x)ui*
i =1
n
(6)
这里 ni ( x) 即为 i 节点的形函数在 x 点的值
ni ( x) = ∑ p j ( x)[ A −1 ( x) B( x)] ji
j =1
m
(6a)
形函数关于坐标的偏导数为
1 −1 ni ,k ( x) = ∑ { p j ,k ( x)[ A −1 ( x) B( x)] ji + p j ( x)[ A,− k ( x ) B ( x ) + A ( x ) B, k ( x )] ji } (7)
1 引
言
无网格伽辽金方法（EFGM）是近年来兴起的一种新型数值计算方法，其基本思路是 利用移动最小二乘法，根据积分点附近一定影响范围（称作影响域）内的节点的位移，用 最小二乘插值得到积分点附近的近似位移场函数。它突破了传统有限元分析中单元网格的 限制，极大地简化了前后处理工作。并且，因为不使用单元网格，在结构发生大变形的情 况下，也不会出现网格畸变的问题，因此，近年来得到广泛的重视和迅速发展。 在无网格伽辽金方法中，移动最小二乘法（MLS）是影响计算结果的一个关键问题， 在实际计算过程中，移动最小二乘法有三个关键性的参数需要事先确定：第一个是基函数
j =1 m
《弹性力学》试题参考答案(参考题)
《弹性力学》试题参考答案（答题时间：100分钟）一、填空题（每小题4分）1．最小势能原理等价于弹性力学基本方程中： 平衡微分方程 ， 应力边界条件 。
2．一组可能的应力分量应满足： 平衡微分方程 ，相容方程（变形协调条件） 。
3．等截面直杆扭转问题中， M dxdy D=⎰⎰2ϕ的物理意义是 杆端截面上剪应力对转轴的矩等于杆截面内的扭矩M 。
4．平面问题的应力函数解法中，Airy 应力函数ϕ在边界上值的物理意义为 边界上某一点（基准点）到任一点外力的矩 。
5．弹性力学平衡微分方程、几何方程的张量表示为：0,=+i j ij X σ ，)(21,,i j j i ij u u +=ε。
二、简述题（每小题6分）1．试简述力学中的圣维南原理，并说明它在弹性力学分析中的作用。
圣维南原理：如果物体的一小部分边界上的面力变换为分布不同但静力等效的面力（主矢与主矩相同），则近处的应力分布将有显著的改变，但远处的应力所受影响可以忽略不计。
作用：（1）将次要边界上复杂的面力（集中力、集中力偶等）作分布的面力代替。
（2）将次要的位移边界条件转化为应力边界条件处理。
2．图示两楔形体，试分别用直角坐标和极坐标写出其应力函数ϕ的分离变量形式。
题二（2）图（a ）⎩⎨⎧=++= )(),(),(222θθϕϕf r r cy bxy ax y x （b ）⎩⎨⎧=+++=)(),(),(33223θθϕϕf r r dy cxy y bx ax y x3．图示矩形弹性薄板，沿对角线方向作用一对拉力P ，板的几何尺寸如图，材料的弹性模量E 、泊松比 μ 已知。
试求薄板面积的改变量S ∆。
题二（3）图设当各边界受均布压力q 时，两力作用点的相对位移为l ∆。
由q E)1(1με-=得，)1(2222με-+=+=∆Eb a q b a l设板在力P 作用下的面积改变为S ∆，由功的互等定理有：l P S q ∆⋅=∆⋅将l ∆代入得：221b a P ES +-=∆μ显然，S ∆与板的形状无关，仅与E 、μ、l 有关。
应用弹塑性力学习题解答
应用弹塑性力学习题解答目录第二章习题答案设某点应力张量的分量值已知，求作用在过此点平面上的应力矢量，并求该应力矢量的法向分量。
解该平面的法线方向的方向余弦为而应力矢量的三个分量满足关系而法向分量满足关系最后结果为利用上题结果求应力分量为时，过平面处的应力矢量，及该矢量的法向分量及切向分量。
解求出后，可求出及，再利用关系可求得。
最终的结果为已知应力分量为，其特征方程为三次多项式，求。
如设法作变换，把该方程变为形式，求以及与的关系。
解求主方向的应力特征方程为式中：是三个应力不变量，并有公式代入已知量得为了使方程变为形式，可令代入，正好项被抵消，并可得关系代入数据得，，已知应力分量中，求三个主应力。
解在时容易求得三个应力不变量为，，特征方程变为求出三个根，如记，则三个主应力为记已知应力分量，是材料的屈服极限，求及主应力。
解先求平均应力，再求应力偏张量，，，，，。
由此求得然后求得，，解出然后按大小次序排列得到，，已知应力分量中，求三个主应力，以及每个主应力所对应的方向余弦。
解特征方程为记，则其解为，，。
对应于的方向余弦，，应满足下列关系（a）（b）（c）由（a）,(b)式，得，，代入（c）式，得，由此求得对，，代入得对，，代入得对，，代入得当时，证明成立。
解由，移项之得证得第三章习题答案取为弹性常数，，是用应变不变量表示应力不变量。
解：由，可得，由，得物体内部的位移场由坐标的函数给出，为，，，求点处微单元的应变张量、转动张量和转动矢量。
解：首先求出点的位移梯度张量将它分解成对称张量和反对称张量之和转动矢量的分量为，，该点处微单元体的转动角度为电阻应变计是一种量测物体表面一点沿一定方向相对伸长的装置，同常利用它可以量测得到一点的平面应变状态。
如图所示，在一点的3个方向分别粘贴应变片，若测得这3个应变片的相对伸长为，，，，求该点的主应变和主方向。
解：根据式先求出剪应变。
考察方向线元的线应变，将，，，，，代入其中，可得则主应变有解得主应变，，。
(完整版)《弹性力学》试题参考答案
《弹性力学》试题参考答案（答题时间：100分钟）一、填空题（每小题4分）1．最小势能原理等价于弹性力学基本方程中： 平衡微分方程 ， 应力边界条件 。
2．一组可能的应力分量应满足： 平衡微分方程 ，相容方程（变形协调条件） 。
3．等截面直杆扭转问题中， 的物理意义是 杆端截面上剪应力对转轴的矩等于M dxdy D=⎰⎰2ϕ杆截面内的扭矩M 。
4．平面问题的应力函数解法中，Airy 应力函数在边界上值的物理意义为 边界上某一点（基准ϕ点）到任一点外力的矩 。
5．弹性力学平衡微分方程、几何方程的张量表示为： ，。
0,=+i j ij X σ)(21,,i j j i ij u u +=ε二、简述题（每小题6分）1．试简述力学中的圣维南原理，并说明它在弹性力学分析中的作用。
圣维南原理：如果物体的一小部分边界上的面力变换为分布不同但静力等效的面力（主矢与主矩相同），则近处的应力分布将有显著的改变，但远处的应力所受影响可以忽略不计。
作用：（1）将次要边界上复杂的面力（集中力、集中力偶等）作分布的面力代替。
（2）将次要的位移边界条件转化为应力边界条件处理。
2．图示两楔形体，试分别用直角坐标和极坐标写出其应力函数的分离变量形式。
ϕ题二（2）图（a ） （b ）⎩⎨⎧=++= )(),(),(222θθϕϕf r r cy bxy ax y x ⎩⎨⎧=+++= )(),(),(33223θθϕϕf r r dy cxy y bx ax y x 3．图示矩形弹性薄板，沿对角线方向作用一对拉力P ，板的几何尺寸如图，材料的弹性模量E 、泊松比 μ 已知。
试求薄板面积的改变量。
S∆题二（3）图设当各边界受均布压力q 时，两力作用点的相对位移为。
由得，l ∆q E)1(1με-=)1(2222με-+=+=∆Eb a q b a l 设板在力P 作用下的面积改变为，由功的互等定理有：S ∆lP S q ∆⋅=∆⋅将代入得：l ∆221b a P ES +-=∆μ显然，与板的形状无关，仅与E 、、l 有关。
一种求解多阶梯悬臂梁自由振动问题的新方法与实验验证_谢瑾荣
模态 测 试 实 验 中， 采用冲击激励法测量结构频 率和振型 。 实验时采用多点激励 、 单点响应的方式 。 将试件分成 17 等 份， 即 共 有 17 个 测 点 。 实 验 用 力 锤输入激励信号， 加速度传感器采集某一点的响应 信号， 拾振点固定 放 在 第 5 个 测 点 。 首 先 对 第 1 个
A New Method for Solving Free Vibration of Cantilever Beam with Multiple Steps and Its Experimental Validation
Xie Jin-rong ，Zhou Cui-ying ， Cheng Ye ， Lv Zhong-rong
［收稿日期］ 2011 0817
（ 1）
Earthquake Resistant Engineering and Retrofitting
August
2012
第 34 卷第 4 期
谢瑾荣， 等： 一种求解多阶梯悬臂梁自由振动问题的新方法与实验验证
· 53 ·
EI 1 （ x 1 ） EI （ X ） = 0 N 11 0 S（ x） = 0
… O … N 12 0 …
EI i （ x i ） N 13 N 21 N 14 N 22 0
0
（ 2）
EI i （ x i ） = E 其中，
bh 3 i ； b 为 梁 的 横 截 面 宽 度； E 为 材 12
料的弹性模量 。 将复合单元法单元形函数组装成整体形函数： … N 24 0 N i1 … O N i2 N i3 … … 0 0 N i4 j11 j21 j i1 j12 j22 j i2 … … O … j1 j j2 j j ij
卷积型伽辽金法求解任意边界梁的动力学问题
卷积型伽辽金法求解任意边界梁的动力学问题本文以卷积型伽辽金法为基础，探讨了其在任意边界梁动力学问题中的应用。
首先介绍了卷积型伽辽金法的基本原理和优势，然后详细阐述了其在任意边界梁动力学问题中的具体应用，包括边界条件的处理、频率响应函数的计算和动态响应的求解等。
最后，通过一个实例验证了卷积型伽辽金法在任意边界梁动力学问题中的有效性和实用性。
关键词：卷积型伽辽金法；任意边界梁；动力学问题；频率响应函数；动态响应1. 引言任意边界梁动力学问题是结构工程领域中的一个重要课题，其研究对于提高结构工程的安全性和可靠性具有重要意义。
在这类问题中，我们需要求解结构在外界作用下的动态响应，以评估结构的稳定性和安全性。
然而，由于梁的边界条件和载荷条件通常是复杂的，因此求解任意边界梁的动力学问题是一个具有挑战性的任务。
近年来，卷积型伽辽金法在结构工程领域中得到了广泛应用。
该方法具有高效、精确、通用等优势，能够很好地解决任意边界梁动力学问题。
本文就是基于卷积型伽辽金法，探讨其在任意边界梁动力学问题中的应用。
2. 卷积型伽辽金法的基本原理卷积型伽辽金法是一种基于积分变换的方法，其基本思想是将结构的位移函数表示为伽辽金函数的线性组合形式，然后利用卷积积分的性质将结构的动态响应计算转化为频率响应函数的计算。
卷积型伽辽金法的基本方程组如下所示：$$begin{aligned}&int_{0}^{L} G(x, s) u(s, omega) d s=frac{sin (omega L)}{omega} f(x, omega)&u(x, t)=frac{1}{2 pi} int_{-infty}^{+infty} U(x, omega) e^{i omega t} d omegaend{aligned}$$其中，$G(x,s)$是伽辽金函数，$u(x,t)$是结构在时间$t$时刻$x$处的位移，$f(x,omega)$是结构在频率$omega$时刻$x$处的外力，$U(x,omega)$是结构在频率$omega$时刻$x$处的频率响应函数。
悬臂梁挠度公式推导过程
悬臂梁挠度公式推导过程1. 悬臂梁的基本概念1.1 悬臂梁的定义悬臂梁的定义简单明了，一头固定、另一头自由。
就像是你在玩秋千，一头被绳子固定住，而你在秋千的另一头尽情摇摆。
1.2 挠度的重要性挠度是个技术词，但它其实就是测量梁弯曲程度的指标。
你能想象吗？一个结构如果挠度太大，可能就要“出大事”了，比如变形、开裂，甚至坍塌，真是个不小的麻烦！2. 挠度的推导过程2.1 载荷与反应在推导之前，我们得先了解一下载荷。
想象一下，你在悬臂梁的自由端放了个大西瓜，这个西瓜的重量就是载荷。
这个时候，梁会因为重力而弯曲。
我们需要计算出这弯曲的程度，嘿，这就是我们的目标！。
2.2 力学基本原理这时就得用到力学的基本原理了。
我们通常使用的公式是 ( y = frac{F cdotL^3{3EI )，其中的F是载荷，L是梁的长度，E是材料的弹性模量，而I是截面的惯性矩。
听起来有点复杂，但我们可以想象，F越大、L越长，挠度也会跟着增大，简单粗暴！3. 公式的应用与实际意义3.1 实际应用在实际工程中，这个挠度公式就像是建筑师的“心灵密码”。
无论是设计桥梁还是大楼，都要考虑挠度，确保它们在使用时不会像小鸡一样“摇摇欲坠”。
3.2 安全第一不要小看这个挠度，它关系到我们生活的安全。
比如，你家楼上的阳台，设计时可不能让它往下垂，万一有人站上去，哎呀，那可真是“屋漏偏逢连夜雨”了！总之，悬臂梁挠度的推导虽然看似枯燥，但其实背后蕴含着无数的智慧和安全考虑。
希望下次你看到悬臂梁时，不再只是觉得它是一根普通的梁，而是会想到它的挠度与安全的重要性。
生活中，哪怕是微不足道的小事，也总是有它背后的大道理，不是吗？。
弹性力学试题及答案
《弹性力学》试题参考答案（答题时间：100分钟）一、填空题(每小题4分)1．最小势能原理等价于弹性力学基本方程中：平衡微分方程, 应力边界条件。
2．一组可能的应力分量应满足：平衡微分方程，相容方程（变形协调条件) 。
3．等截面直杆扭转问题中，的物理意义是杆端截面上剪应力对转轴的矩等于杆截面内的扭矩M.4．平面问题的应力函数解法中,Airy应力函数在边界上值的物理意义为边界上某一点(基准点)到任一点外力的矩。
5．弹性力学平衡微分方程、几何方程的张量表示为：，。
二、简述题（每小题6分）1．试简述力学中的圣维南原理，并说明它在弹性力学分析中的作用。
圣维南原理：如果物体的一小部分边界上的面力变换为分布不同但静力等效的面力（主矢与主矩相同），则近处的应力分布将有显著的改变，但远处的应力所受影响可以忽略不计。
作用：（1）将次要边界上复杂的面力（集中力、集中力偶等）作分布的面力代替。
（2）将次要的位移边界条件转化为应力边界条件处理。
2．图示两楔形体，试分别用直角坐标和极坐标写出其应力函数的分离变量形式。
题二（2）图(a）（b）3．图示矩形弹性薄板，沿对角线方向作用一对拉力P，板的几何尺寸如图,材料的弹性模量E、泊松比 已知.试求薄板面积的改变量.题二（3）图设当各边界受均布压力q时，两力作用点的相对位移为。
由得，设板在力P作用下的面积改变为，由功的互等定理有：将代入得：显然，与板的形状无关，仅与E、、l有关。
4．图示曲杆，在边界上作用有均布拉应力q，在自由端作用有水平集中力P.试写出其边界条件（除固定端外）。
题二（4）图（1）;(2）（3）5．试简述拉甫（Love）位移函数法、伽辽金（Galerkin）位移函数法求解空间弹性力学问题的基本思想，并指出各自的适用性Love、Galerkin位移函数法求解空间弹性力学问题的基本思想:（1）变求多个位移函数或为求一些特殊函数，如调和函数、重调和函数。
（2)变求多个函数为求单个函数(特殊函数）。
压杆大挠度问题的里兹_伽辽金近似解
(2) 可
得:y
=-
1 K2
dθ ds
( 14)
将 ( 4) 式代入有 :
y
απ π = K2L sin L s
=
αL π
Pcr sin π s PL
αλ π = π sin L s
( 15 )
最大挠度 δ出现在 x = L /2处 ,由 ( 15) 式得 :
δ/L =α·λ/π
( 16 )
52
低 温 建 筑 技 术 2005年第 5期 (总第 107期 )
3
取 γ = - ( 1 + 9λ) + ( 1 + 9λ) 2 + 8
( 13)
综合 ( 13) 、( 10) 和 ( 8) 式可以求得 α。
2 压杆的位移
将
( 1) 式中以
dy 代替 ds
sinθ, 积分并代入边界条件
李朋主等 :压杆大挠度问题的里兹 - 伽辽金近似解
51
压杆大挠度问题的里兹 - 伽辽金近似解
李朋主 , 赵 锐 , 刘晚成
(东北林业大学土木工程学院 , 哈尔滨 150040 )
【摘 要 】 通过非线性振动中常用的里兹 2伽辽金法 ,推出了两端铰支压杆大挠度问题的一个新近似公式 ,计 算结果表明具有较高的精度 。 【关键词 】 压杆 ;大挠度 ;近似解 ;里兹 2伽辽金法 【中图分类号 】 TU31112 【文献标识码 】 B 【文章编号 】 1001 - 6864 (2005) 05 - 0051 - 02
10 °
20 °
40°
60°
80°
90°
11518 100°
11678 110°
11885 120°
【国家自然科学基金】_大挠度问题_基金支持热词逐年推荐_【万方软件创新助手】_20140802
科研热词 非线性 误差分析 计量学 触觉再现 标定 有限元法 斜板 挠度 悬臂梁 小波 大挠度 大变形 固支 后屈曲 变积方法 变截面薄壁箱梁 压杆稳定性 势能变分原理 加权残值法 剪力滞 两类变量广义变分原理
2010年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
科研热词 驻值余能原理 饱和多孔弹性梁 非线性稳定性 非保守力 软土路基 薄板理论 荷载分担比 等效连续化 热颤振 桩承式加筋路堤 极限环颤振 大挠度 多孔介质理论 壁板颤振 双层网格圆底扁球壳 几何非线性 修正迭代法 galerkin截断法
推荐指数 3 11 1 1 1
2013年 序号 1 2 3 4 5 6 7 8 9 10 11 12
科研热词 蠕变屈曲 精确解 粘弹性板 柔性梁 挠度函数 悬臂梁 小波伽辽金方法 大挠度弯曲 大挠度 多解问题 几何非线性 jacobi椭圆函数
推荐指数 1 1 1 1 1 1 1 1 1 1 1 1
2014年 序号 1 2 3 4 5
2014年 科研热词 钢悬链式立管 触地点 涡激振动 流固耦合 大挠度细长梁 推荐指数 1 1 1 1 1
2012年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
科研热词 大挠度 风力机叶片 非线性 钢梁 载荷 计算模型 薄板大挠度弯曲问题 算例 火载荷 梁 极限温度准则 无网格方法 振动 挥舞振动 挠度 弯曲方程 幅频特性 屈服函数 对偶互易方法 夹层圆板 基本解方法 周期解 变形 参数分析 几何非线性
变分法
变分法综述1.变分法1.1.变分法起源变分法是17世纪末发展起来的一门数学分支，主要是古典变分法，它理论完整，在力学、光学、物理学、摩擦学、经济学、宇航理论、信息论和自动控制论等诸多方面有广泛应用。
20世纪中叶发展起来的有限元法，其数学基础之一就是变分法。
[1]变分法是处理泛函的数学领域，和处理函数的普通微积分相对。
譬如，这样的泛函可以通过未知函数的积分和它的导数来构造。
变分法最终寻求的是极值函数：它们使得泛函取得极大或极小值。
有些曲线上的经典问题采用这种形式表达：一个例子是最速降线，在重力作用下一个粒子沿着该路径可以在最短时间从点A 到达不直接在它底下的一点B 。
在所有从A 到B 的曲线中必须极小化代表下降时间的表达式。
变分法的关键定理是欧拉－拉格朗日方程。
它对应于泛函的临界点。
在寻找函数的极大和极小值时，在一个解附近的微小变化的分析给出一阶的一个近似。
它不能分辨是找到了最大值或者最小值（或者都不是）。
变分法在理论物理中非常重要：在拉格朗日力学中，以及在最小作用量原理在量子力学的应用中。
变分法提供了有限元方法的数学基础，它是求解边界值问题的强力工具。
它们也在材料学中研究材料平衡中大量使用。
而在纯数学中的例子有，黎曼在调和函数中使用狄力克雷原理。
最优控制的理论是变分法的一个推广。
[2]同样的材料可以出现在不同的标题中，例如希尔伯特空间技术，摩尔斯理论，或者辛几何。
变分一词用于所有极值泛函问题。
微分几何中的测地线的研究是很显然的变分性质的领域。
极小曲面（肥皂泡）上也有很多研究工作，称为Plateau 问题。
1.2变分问题类型固定边界的变分问题，可动边界的变分问题，条件极值变分问题和参数形式的变分问题。
[3]（1）古典变分问题举例 例1：最速降线或捷线问题（Brachistorone or curve of Steepest descent ）问题。
这是历史上出的第一个变分法问题，1696年约翰·伯努利提出的。
梁的挠度和转角问题分析
科学技术创新2018.06梁的挠度和转角问题分析王爽焦之森（齐齐哈尔大学建筑与土木工程学院，黑龙江齐齐哈尔161000）对简支梁、外伸梁的变形问题的解析计算方法有很多种，常见的有积分法[1-5]、能量法[1-5]、叠加法[1-5]、奇异函数法[1-5]和共轭梁法[1-5]等，在用积分法求解简支梁、外伸梁的变形问题时须求解多个积分常数，计算繁琐；奇异函数法仍属于积分法，求解过程也须解积分常数；如果仅计算某一截面的位移，能量法较为简单，不过仍须进行积分计算[6]。
本文通过间接叠加法，来介绍简支梁、外伸梁等结构在受载荷作用时挠度及转角问题的简单求解方法，即将简支梁、外伸梁等结构在受载荷作用时挠度及转角问题，转化为有初始转角的悬臂梁受载荷时的变形问题，使简支梁、外伸梁等结构在受载荷作用时挠度及转角问题的求解过程的思维难度得到很大程度的降低，从而问题变得更容易理解。
1原理介绍与例题分析悬臂梁具有一个固定端，当悬臂梁受已经与水平线外荷载作用时，靠近固定端的载面不发生转动，转角为零。
如果有一个悬臂梁，在未荷载时，形成一个小的角度θB ，如图1所示。
图1有初始转角的悬臂梁x 轴为水平方向，梁轴线与x 轴成角θB ，即θB 为初始转角，此梁称为有初始转角的悬臂梁。
在未受荷载时，相对于x 轴，自由端已经有一挠度为θB l 。
根据叠加法，当加一静荷载F 时，自由端的挠度ω=θB l+Fl 33EI 转角为θB +Fl22EI。
应用初始转角悬臂梁概念，只要知道悬臂梁在集中力偶、集中力和均布载荷作用下自由端的挠度和转角公式，就可以通过叠加法，求解简支梁、外伸梁、的变形问题。
跨长l ，刚度EI 的悬臂梁在集中力偶Me ，集中力F ，均布荷载q 作用下，自由端的挠度和转角公式列出如下Mel 22EI ，Mel EI ，Fl 33EI，Fl 23EI ，ql 48EI ，ql 36EI。
下面举几个例子。
例1.如图例2-1所示简支梁端受集中力偶Me 作用，求端截面转角。
悬臂长板的大挠度弯曲实验
3.1 ANSYS分析与计算（SHELL43壳体单元）
图2. 悬臂长板有限元模型加载图
3.1 ANSYS分析与计算（SHELL43壳体单元）
图3. 非线性迭代求解过程收敛情况显示
3.1 ANSYS分析与计算（SHELL43壳体单元）
图4. 变形前后的有限元模型
3.1 ANSYS分析与计算（SHELL43壳体单元）
4. 结果分析与讨论
利用上面程序，分别用SHELL43壳体单元和 SHELL181壳体单元计算了悬臂长板在端部集中力偶 作用下的大挠度弯曲，计算结果如表1所示。和理论 值相比，两者误差都不大，但是有限元变形壳单元 SHELL181精度更高。说明这两种壳体单元对大变形 分析的模拟精度是非常高的。
4. 结果分析与讨论
附录:数值计算程序
CP,1,ROTY,ALL TORQ=7.854 F,2,MY,TORQ NSEL,ALL FINISH ! (11) 进入求解模块,设置求解选项后,开始求解 /SOLU AUTOTS,ON NSUBST,10 LNSRCH,ON OUTPR,BASIC,LAST OUTRES,ALL,ALL SOLVE FINISH ! (12)在求解模块重启动分析 /SOLU ANTYPE,,REST F,2,MY,TORQ*2 SOLVE FINISH ! (13)在后处理模块,显示计算结果 /POST1 RSYS,SOLU SET,2 SHELL,TOP PRNSOL,S,COMP PRNSOL,DOF FINISH
附录:数值计算程序
FINISH /CLEAR, NOSTART ! (1)设置工程参数、分析类型、单元类型、实常数和材料参数 /FILNAM, EX6.2 /PREP7 /TITLE, EX6.2 LARGE DEFLECTION OF A CANTILEVERED PLATE SMRT,OFF ANTYPE,STATIC NLGEOM,ON ! 打开大变形选项 ET,1,SHELL43 R,1,1 MP,EX,1,1800 MP,NUXY,,0 ! (2)定义关键点、线和面,并划分为面单元 K,1 K,2,12 K,3,12,1 K,4,,1 L,1,2 L,3,4 LESIZE,ALL,,,10 ESIZE,,1 A,1,2,3,4 /VIEW,1,0.61,-0.64,0.47 /VUP,1,Z AMESH,1
变截面悬臂梁最大挠度简化计算
Ｆ
Ｍ ｘ＝Ｆ（ ｘ － ａ ）
６ ｂ
：
１
（ ７ ａ） （ ７ ｂ）
均布 荷 载 作 用 ，
况 ， 由式 （ １ Ｏ） 理 论 公式 计 算 得最 大 挠
度 为１ ． ７ ０ ７ ７ ｍｍ ，而 按 式 （ １ １） 得 最 大 挠度 为１ ． ５ ２ ７ Ｘ １ ． １ １ ６ ＝ １ ． ７ ０ ４ １ ｍｍ。 二
鲁 ＝ ｈ 一 ｈ
．
（ ４ ）
、 ＿ ｒ ，
叩
（ ２）
相应地 ，截面位 置Ｘ 处的截面特性
和 内力 如 下 ：
＝
１ ． ２ 最大挠 度计算原理
对 图１ 所 示 的悬 臂 构 件 ，最 大 挠 度
本 文利 用经 典结 构 力学 中 的单位
ｂ ｈ ＝ｂ Ａ ｘ
构 件 的 挠度 计 算 提供 了简 单 实 用 的 计 算 方法。
＋
０
一 ｘ
』
０
“ｘ
ｏ
』
弯矩和剪 力是其主要 内力 ，采用变截 面 则易 于充 分发 挥材 料 强度 ，节 约 了材 １构件参数和最大挠 度计 算原理 料 ，因此 ，在土木工程 中变截面构件较 为 常 见 ，如 挑 梁 、 门式 刚 架 的 梁 和柱 、 １ ． １构件参数
形 产 生 的 挠度 。
算其挠度 时，其显式表达一般较 复杂甚
至 无 法 求 解 。 因 此 ， 变 截 面 构 件 的 挠 度 ， 常 采 用 近 似 方 法 计 算 。 其 中 ，共 轭梁 法…需计算虚梁 内力图的合力及位
图１ 构件 参 数
构 件 端部 和 根 部 的截 面 面 积和 截 面 置 ，相 应 地 ，计 算 也 较 为 复 杂 ；等效 刚 度 法 将 构 件 分 若 干段 后 ，利 用 最 大 挠 模 量 分 别 如 下 ： ２ 最大挠度的理论解 度 相 等 得 到 其 等 效 刚 度 ，再 利用 等截 面 端 部截 面 ， 梁 的挠 度 计 算 公 式 得 到 具体 挠 度 ，该 法 的计 算 工 作 量 较 大 。 而 在工 程 实践 中 ，
大挠度Euler-Bernoulli梁单元的算法研究
大挠度Euler-Bernoulli梁单元的算法研究周建博;冯志强;周洋靖【摘要】通过有限元方法研究Euler-Bernoulli梁变形的求解.形状函数矩阵通过Lagrange插值函数和两节点Hermite单元构造.然后根据形状函数矩阵推导大挠度几何非线性单元刚度矩阵.再利用C++编程,开发出一套可用于求解梁杆结构大挠度问题的算法.最后通过典型算例来验证本算法的计算精度.【期刊名称】《黑龙江科学》【年(卷),期】2018(009)013【总页数】4页(P10-13)【关键词】Euler-Bernoulli梁;有限元方法;大挠度;算法【作者】周建博;冯志强;周洋靖【作者单位】西南交通大学力学与工程学院,成都610031;西南交通大学力学与工程学院,成都610031;西南交通大学力学与工程学院,成都610031【正文语种】中文【中图分类】TH113工程中有大量的结构使用梁杆构件，其安全性和可靠性是设计所应追求的主要目标。
随着科学技术的发展，高强、轻质材料在显著提高了梁杆结构承载能力的同时，又降低自重，越来越受到人们的重视。
在静力学条件下，由于轴向和横向荷载的共同作用，梁杆构件将表现出明显的二阶效应(轴力效应)，即产生附加的横向位移和弯矩(拉力使横向位移和弯矩减小、压力使横向位移和弯矩增大)。
其实质是轴向荷载导致了构件刚度的改变[1]。
这个效应在静力学中表现为几何非线性。
如果此时仍然使用线性分析方法，将导致较大的计算误差。
严格意义上来讲，自然界的问题都是非线性问题。
对于很多非线性问题，可以近似用线性分析，也能达到相应的精度要求。
为了更准确地求解梁杆结构，很多学者开展了深入研究。
18世纪，基于Euler-Bernoulli的两个假设：第一，刚性横截面假定：变形前垂直梁中心线的平剖面，变形后仍然为平面。
第二，变形后的横截面仍与变形后的中性轴垂直。
于是建立了欧拉梁理论，在实际中应用广泛[2]。
Bank[3]给出了弹性体的非线性增量形式的完全拉格朗日描述和更新拉格朗日描述。
悬臂梁集中载荷大挠度弯曲变形的一种解_葛如海
dθ si n(θB - θ)
( 4)
则
L* = 1 F (θB )
( 5)
2a
又
dy = dssi nθ
dx = dsco sθ
将 ( 3) 式代入并积分得:
∫yB dy = 0
∫ 1
θB
2a 0
sinθ dθ si n(θB - θ)
∫l dx = 0
∫ 1
θ B
2a 0
cosθ dθ si n(θB - θ)
1 2
Fx2
(θB
)
[ co sθB
+
sinθB
Fy Fx
(θB ) (θB )
]
令
FT (θB ) =
1 2
Fx2
(θB
)
[ cosθB
+
si nθB
ห้องสมุดไป่ตู้
Fy Fx
(θB ) (θB )
]
( 14)
则
Tw A =
EI l
FT
(θB
)
( 15)
积分函数 FT (θB ) 可根据 Fy (θB )、 Fx (θB ) 求出近似解 ,最后可得
在图三所示工况由 ( 1) 式可知:
ddθs = -
Tw EI
=
PE′I ( yB -
y)
分析过程与前相似 ,令
b2 =
P′ EI
( 16)
则
ddθs = b2 ( yB - y )
对等式两端微分 ,两端同乘 dθ后再积分得
1 2
(
ddθs ) 2
=
b2 co sθ+
c
在梁的自由端 ddθs = -
板的数值解法—伽辽金积分—伽辽金解法
板的数值解法—伽辽金积分—伽辽金解法
李从尊
【期刊名称】《天津纺织工学院学报》
【年(卷),期】1989(000)002
【摘要】本文从伽辽金积分出发,从而导出一个常微分方程,然后再利用伽辽金解法解这个常微分方程,这样就得到了板的挠曲面方程。
【总页数】6页(P61-66)
【作者】李从尊
【作者单位】无
【正文语种】中文
【中图分类】O241.4
【相关文献】
1.基于Matlab的直管流致振动的伽辽金解法 [J], 曹建华
2.Euler-Bernoulli梁的伽辽金解法 [J], 袁晓光
3.无单元伽辽金弱形式积分方法的研究 [J], 王立鹏;王欣彦;战洪仁;曾祥福;张先珍;王翠华
4.基于单位分解积分的伽辽金无网格方法研究 [J], 曾清红;卢德唐
5.稳定节点积分伽辽金无网格法的应力计算方法研究 [J], 王东东;李凌;张灿辉因版权原因，仅展示原文概要，查看原文内容请购买。