上海市力学学会成立于1959年9月20日，到今年正好是一个甲子．60个春秋真实写照了上海市力学学会的建立、成长、发展与壮大的过程．60年来，在上海市科协的领导和中国力学学会的指导下，在各团体成员单位的大力支持下，上海市力学学会一代代力学工作者薪火相传、开拓耕耘，历届理事会及下属工作委员会和专业委员会团结奋进、不断创新，广大会员踊跃参与，为上海力学大力发展、力学人才精心培养、力学知识科学普及、力学服务国家发展战略和区域经济建设中做出了重要贡献，学会持续发展壮大．在庆祝力学学会成立六十周年之际，让我们共同抒怀上海力学六十年薪火相传、回顾近10年重要发展，展望上海力学新时代．
针对植入Z-pin后碳纤维增强复合材料的微观结构，通过施加Z-pin周期性边界约束条件，建立了Z-pin复合材料单层板单胞细观模型．考虑固化过程中树脂体积收缩、弹性模量随固化度变化和纤维因Z-pin进入偏转因素，运用有限单元法计算了单胞结构在固化成型工艺过程中树脂和纤维应力发展和分布，并研究了Z-pin直径和分布密度对单层板面内残余应力的影响．结果表明：凝胶点之前，树脂模量和残余应力很小，凝胶点之后，树脂模量和残余应力增加较快；残余应力分布与纤维偏转有关；Z-pin直径和分布密度增加会使固化残余应力增大．
本文通过局部展开算法研究了多孔介质中排列结构对非线性渗流规律的影响．计算结果表明：规则排列结构的渗流标度律和渗流方向相关，然而非规则排列结构的渗流标度律不依赖于渗流方向而只依赖于多孔介质的几何结构．伴随着排列结构的非规则程度增强，不同方向的渗流阻力变化规律的差异性减小．当流场输运区域较为复杂且不存在较宽直通道时，渗流阻力在较宽雷诺数范围内体现出较为明显的标度律特性．周期性强非规则排列结构中不同区域大小的模型却具有相似的渗流阻力的标度律特性，同时标度律特性与孔隙度的大小相关性不大．阻力非线性标度律由多孔介质结构的不同而取值为2至3之间．
金属材料疲劳寿命由裂纹萌生和裂纹扩展寿命两部分组成，其中对于萌生寿命中的小裂纹分析是精确描述裂纹萌生寿命的关键．而小裂纹在扩展过程中由于尺寸相对较小，导致传统线弹性断裂力学预测方法失效，需要对其进行改进，考虑裂纹尖端塑性区引起的残余压应力对小裂纹扩展速度的影响．本文针对此问题进行了初步分析，通过对塑性区引起的残余应力的量化，结合小裂纹门槛值特性，提出了一种经验型修正的小裂纹扩展模型，用于定量预测裂纹的萌生寿命．使用铝合金6082-T6缺口试样进行了疲劳实验，并与理论结果进行了对比，验证了所提模型的有效性．
斜拉桥中拉索承受着多种端部激励，可激发大幅空间振动．以斜拉索为对象，探究不同端部激励间相位差对其非线性振动的影响．首先，推导斜拉索无量纲离散控制方程，引入考虑相位的三向端部激励得到一般化模型；然后，针对拉索下端存在的纵桥向、竖向和横桥向激励的两两组合，受大幅或小幅激励，及其在主共振区或主参数共振区几组因素，共计12种工况，采用数值分析法分别研究了各工况下不同激励相位差时的斜拉索稳态响应．研究发现：激励相位差能加剧与激励频率相近的面内、外模态振动；在任意端部激励组合下，激励相位差不仅可使斜拉索非线性振动出现定量变化，还可改变内共振的表现形式．面内、外激励组合下，相位差对拉索响应幅值的影响以π为周期变化，且当相位差趋于π/2 + kπ (k = 0, 1, 2…)时影响最为突出；而面内激励组合下，以2π为变化周期，当相位差为π + 2kπ (k = 0, 1, 2, …)时其对稳态幅值的影响最显著．其原因是：面外激励关于拉索所在的竖直面对称，故其本质上以π为周期；而面内激励无此对称性，仍以2π为周期．因此，有无面外激励参与决定了激励间相位差对斜拉索响应的影响规律．
当周期激励频率远小于系统固有频率时，会存在快慢耦合效应，与单项激励不同，参外联合激励不仅会导致快子系统平衡曲线和分岔行为的复杂化，也会产生一些特殊的非线性现象，为此，本文以两耦合Hodgkin-Huxley细胞模型为例，引入周期参外联合激励，探讨在频域不同尺度耦合时该系统的簇发振荡的特点及其分岔机制．通过建立相应的快慢子系统，得到慢变参数变化下的快子系统的各种分岔模式以及相应的分岔行为，结合转换相图，揭示耦合系统随激励幅值变化时的动力学行为及其机理．研究表明，在激励幅值较小时，系统表现为概周期振荡，两频率分别近似于快子系统平衡曲线由Hopf分岔引起的两稳定极限环的振荡频率．概周期解随激励幅值的增加进入簇发振荡，导致这些簇发振荡的主要原因是在慢变参数变化的部分区间内，存在唯一稳定的平衡曲线，使得系统的轨迹逐渐趋向该平衡曲线，产生沉寂态，并随着慢变参数的变化，由分岔进入激发态．同时，快子系统中参与簇发振荡的稳定吸引子随激励幅值的变化也会不同，导致不同形式的簇发振荡．另外，与单项激励下的情形不同，联合激励时快子系统的部分稳定吸引子掩埋在其它稳定吸引子内，从而失去对簇发振荡的影响．
利用复变函数方法并结合双准周期Riemann边值问题理论，获得了含双周期分布非均匀相（夹杂/界面层）的复合材料在远场均匀反平面应力下弹性场的全场解答．该解答可用于对纳米夹杂复合材料的应力进行分析，结合平均场理论也用于预测纳米夹杂复合材料的有效性能．计算结果表明：当夹杂尺度在纳米量级时，应力和有效反平面剪切模量具有明显的尺度依赖性，并且随着夹杂尺寸的增加，趋近于不考虑界面效应时的结果；界面层厚度和性能对应力和有效反平面剪切模量明显变化时所对应的夹杂尺度范围和趋近于无界面效应结果的快慢有显著影响；当界面厚度足够薄时，界面层模型可用于模拟零厚度界面情况．
对于一般任意支撑的连续梁结构动力稳定性问题，已有的计算方法求解过程都很复杂，给工程设计带来极大的不便．本文提出了一个简化的分析方法，利用现有的商业软件，只需求得连续梁的自然频率及静力屈曲（失稳）荷载，就可容易得到结构的动力失稳区域，当考虑结构阻尼对不稳定区域的影响时，可将阻尼矩阵表达为Rayleigh阻尼的形式．研究结果表明：采用本文计算方法与已有的理论计算方法得到的连续梁主参数共振的不稳定边界非常吻合，而本文计算方法更为简单，计算结果可靠，计算精度高，可满足工程设计的需要．
为研究地震作用下地铁车站及其上盖结构体系的地震反应特征，本文以一典型的地铁车站及上盖结构体系为背景，采用ABAQUS软件建立土-地铁车站及其上盖结构体系的三维有限元模型，利用等效线性化方法考虑土的非线性，以不同类型的单向地震波和双向地震波（含水平向和竖向）作为输入，对土-地铁车站及其上盖结构体系进行了地震反应分析，比较了不同类型地震波作用下地铁车站及其上盖结构体系地震反应的差异，探讨了竖向地震动效应和转换梁刚度的影响．本文的算例结果表明：与普通基岩波和El-centro反演波输入相比较，在长周期基岩地震波作用下地铁车站及其上盖结构体系的内力响应和位移响应均最大．这一地铁车站及其上盖结构体系的薄弱部位位于地铁车站底层第2列、第3列柱的柱端和侧墙与底板连接处．竖向地震动对地铁车站及其上盖结构体系中柱的轴力响应值有较大影响．对于本文的动力计算工况，适当弱化转换梁刚度可降低上盖结构层间位移角的动力响应值．
建立了轴向压力作用下悬臂裂纹梁边界支承和裂纹损伤程度识别方法．首先，将悬臂梁边界非完整支承等效为竖向和扭转弹簧、梁中开裂纹等效为内部扭转弹簧，利用Laplace变换，得到了边界弹性支承、考虑轴向压力二阶效应、具有任意裂纹数目Euler-Bernoulli悬臂梁弯曲挠度的解析解．其次，提出了边界弹性支承弹簧柔度和裂纹等效扭转弹簧柔度的识别方法．最后，通过数值试验，考察了轴向压力，裂纹深度以及测量误差等对识别结果的影响，说明了本文考虑轴向压力二阶效应的悬臂梁边界支承弹簧柔度及裂纹等效扭转弹簧柔度识别方法的适用性和可靠性，结果表明：相比于应变测量误差，挠度测量误差对裂纹损伤程度识别结果影响更加敏感，且轴向压力对裂纹损伤程度识别影响较小，因此，应严格控制挠度的测量误差．同时，边界支承扭转弹簧柔度的识别误差大于其竖向弹簧柔度识别误差．这些结果为实际工程中边界非完整支承悬臂裂纹梁的参数识别提供了指导．
谐波减速器和力矩传感器等柔性元件因其独特性能而广泛应用在空间机器人关节系统中，以获取高减速比．但同时这些柔性元件的存在为空间机械臂系统引入了关节柔性，使得对其稳定控制变得更为复杂．基于此，文中讨论了基于自适应回归小波神经网络(Self-Recurrent Wavelet Neural Networks, SRWNN)的弹性关节空间机械臂系统动力学建模及级联智能滑模控制．首先，利用级联系统理论及第二类拉格朗日方法推导出了由外环刚性臂子系统和内环关节电机转子子系统组成的系统级联动力学模型；其次，为两个子系统分别设计了内、外环自适应滑模回归小波神经网络控制．外环控制算法以期望轨迹为控制量，而其控制信号作为抑制弹性关节振动的内环控制算法的控制量，整个控制系统由内、外环控制系统叠加而成；而后，基于Lyapunov稳定性理论证明了整个控制系统的稳定性并设计了自适应回归小波神经网络的各权值参数在线学习算法．所提的控制算法有效地消除了模型不确定的影响，避免了复杂的求导计算和角加速度可测的要求，同时，控制系统设计过程中未涉及惯常奇异摄动双时标分解操作，在理论上适合任意大小的关节柔性刚度．最后，系统对比仿真试验证明了所提的级联智能控制算法优于惯常基于奇异摄动法和基于柔性铰补偿奇异摄动法的控制方案．
研究判定广义Birkhoff系统稳定性的三重组合梯度方法．首先，给出4类三重组合梯度系统的定义和微分方程；其次，得到广义Birkhoff系统成为三重组合梯度系统的条件，从而将广义Birkhoff系统化成三重组合梯度系统；最后，利用组合梯度系统的性质来研究系统的稳定性，举例说明结果的应用．
为运用荷载挠度曲线进行损伤识别，通过力法推导了梁结构在均布荷载作用下的挠度曲率理论公式，提出通过损伤前后的挠度曲率差进行损伤定位，针对多跨连续梁均布荷载下挠度曲率指标存在的损伤识别漏判问题，提出逐跨均布荷载挠度曲率指标以避免其影响，并建立了挠度曲率与损伤程度的理论关系式，可对损伤程度进行较精确的定量描述．通过一简支梁和一三跨连续梁算例，考虑多种损伤工况，分析了指标在不同程度噪声水平下的抗干扰能力，验证了挠度曲率损伤指标应用于实际的可行性．
插值型重构核粒子法的形函数具有离散点插值特性和不低于核函数的高阶光滑性，因而不仅可以直接施加本质边界条件，同时也保证了较高的计算精度．本文将弹性动力学方程作拉氏变换后，在变换域内用插值型重构核粒子法求解，最后再借助Durbin数值反演方法求得时间域的解．针对典型的弹性动力学问题，给出了插值型重构核粒子法的数值算例，并验证了本文方法的有效性．
挤出胀大的数值模拟是非牛顿流体研究中具有挑战性的问题．本文运用格子Boltzmann方法(LBM)分析Oldroyd-B和多阶松弛谱PTT粘弹流体的挤出胀大现象，采用颜色模型模拟出口处粘弹流体和空气的两相流动，通过重新标色获得两种流体的界面，并最终获得胀大的形状．Navier-Stokes方程和本构方程的求解采用双分布函数模型．将胀大的结果与解析解、实验解和单相自由面LBM结果进行了比较，发现格子Boltzmann两相模型结果与解析解和实验结果相吻合，相比于单相模型，收敛速度更快，解的稳定性更高．研究了流道尺寸对胀大率的影响，并对挤出胀大的内在机理进行了分析．
喷丸强化是一种有效提高工件抗疲劳性能的表面处理工艺，经喷丸强化处理后的航空薄壁弧形结构件发生变形是航空制造领域的共性难题．本文基于变截面曲梁平面弯曲理论，建立了弧形件喷丸诱导应力与变形方程，将得到的构件最大变形与有限元方法计算结果进行对比，所建立理论模型能有效预测弧形件喷丸强化后的变形．讨论了喷丸压力和喷丸区域对喷丸变形的影响，结果表明：喷丸强化工艺仅作用于弧形件腹板表面时，弧形件变形以径向收缩为主；喷丸强化工艺仅作用于弧形件筋板表面时，弧形件变形以径向扩张为主；与局部喷丸相比，对弧形件所有表面喷丸强化时，弧形件弯曲变形减小，轴向伸长增大．
本文应用空间滤波方法：FST(Filter-space technique)方法，研究二维Rayleigh-Bénard(RB)湍流热对流系统中湍动能、热能和拟涡能的能量输运．研究中Rayleigh数(Ra)选取为1x10^8、1x10^9和1x10^10，Prandtl数(Pr)固定为4.38．我们展示了的结果表明，在二维RB系统中，三个Ra数下全场的平均湍动能和平均拟涡能在不同滤波尺度下的能量输运与Kraichnan在1967年预测的二维湍流中的级串理论有所偏差，而中心区域的能量都是向小尺度输运的．结果还揭示了瞬时能量输运的一些局部特性，包括它们在小尺度上不对称的分布．
本文针对现有的损伤识别方法不能满足部分结构损伤识别精度要求的现状，对结构的小损伤精确识别方法开展研究．以长细结构为研究对象，对具有不同损伤位置和损伤程度的圆柱形的轻阻尼梁结构进行了数值分析和实验研究，应用数值计算方法和实验确定的特征向量和特征频率对长细结构裂缝参数进行识别计算．本文在研究过程中编制了一个创新性的预测程序，通过其一次性生成目标函数图来选择合适的初始参数，从而对识别结果进行分析．研究结果表明，应用本文提出的识别方法，裂缝位置的识别误差可以控制在0.05 %~0.28 %范围内，裂缝深度识别误差低于7 %．
采用统一强度理论并考虑材料拉伸与压缩弹性模量的差异性，建立均匀内压作用下双层厚壁圆筒的应力表达式，获得了其内压相应的弹性极限解答、塑性极限解答，并分析拉压强度比、拉压模量系数、统一强度理论参数、半径比及分层半径对弹性、塑性极限内压的影响规律．研究结果表明：弹性、塑性极限内压随拉压强度比的增加而减小，但随统一强度理论参数、半径比的增加而增大；弹性极限内压随分层半径的增加呈现先增大后减小变化，随拉压模量系数的增加而一直减小；塑性极限内压与拉压模量系数、分层半径无关．应用于实际工程时，可根据所得结果选择合理的壁厚及分层半径，再根据材料特性确定其他参数，以便更加准确地计算结构的受力状况．
本文利用双摄动方法求解缓慢变化管道中Johnson-Segalman(J-S)流体流动的渐近解．将管道的扩张（或收缩）角度和粘弹性参数分别作为双摄动的参数，由流函数和涡量函数的形式，推导出压力和壁面剪切应力的渐近解．在此基础上，分析了管道角度，粘弹性参数和雷诺数等参数对压力以及剪切力影响．主要结论如下：(1) 管道扩张角度增加时，流向同一位置处径向压力以及壁面剪切应力随扩张角度减小；(2) 在同一扩张管道中，径向压力随着流向位移减小，收缩管与之相反；(3) 扩张角度与雷诺数对流场起主导作用，粘弹性系数起次要作用．
使用高仿真精度和高计算效率的有限元模型进行飞机颤振边界预测具有重要意义．本文提出了基于Huth紧固件柔度计算方法的民机气动弹性模型的改进建模方法，计及了固定前缘的刚度及其与机翼盒段之间的弹性效应，建立了全机气动弹性模型并进行了固有特性和颤振特性分析．分析结果与全机地面共振试验、风洞试验结果的偏差均在±1 %以内，分析时长与现有模型基本保持不变，验证了改进建模方法的有效性．
高速电梯运行时，井道内气压的快速变化会对轿厢运行的稳定性及安全性产生极大影响．为解决上述问题，提出了三种在井道壁面上开风口的设计方案，运用计算流体力学分析软件FLUENT对电梯运行过程中的空气流动状态进行数值仿真，同时改变风口形状以及数量，对井道内风速变化、轿厢所受气动阻力及气动力矩进行分析．分析结果表明，在提出的三种开风口方案中，上下圆形风口为最优，与未开风口时相比，轿厢气动阻力减小22.05 %，侧翻力矩减小29.87 %，俯仰力矩减小18.02 %．最后进行试验验证，表明了开风口方案的有效性，研究结果可为井道风口的设计提供理论依据．
