Porohyperelastic analysis of single chondrocyte using AFM and inverse FEA

The aim of this paper is to determine the creep and relaxation responses of single chondrocytes in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of single chondrocytes at the strain-rate of 7.05 s-1. This result was then employed in inverse finite element analysis (FEA) using porohyperelastic (PHE) idealization of the cells to determine their mechanical properties. The PHE model results agreed well with AFM experimental data. This PHE model was then utilized to study chondrocyte’s creep and relaxation behaviors. The results revealed that the effect of fluid was predominant for cell’s mechanical behaviors and that the PHE is a good model for biomechanics studies of chondrocytes.

Tailoring the resonance of bilayer graphene sheets by interlayer sp3 bonds

Graphene-based resonators are envisioned to build the ultimate limit of two-dimensional nanoelectromechanical system due to their ultrasensitive detection of mass, force, pressure and charge. However, such application has been greatly impeded by their extremely low quality factor. In the present work, we explore, using the large-scale molecular dynamics simulation, the possibility of tailoring the resonance properties of a bilayer graphene sheet (GS) with interlayer sp3 bonds. For the bilayer GS resonator with interlayer sp3 bonds, we discovered that the sp3 bonds can either degrade or enhance the resonance properties of the resonator depending on their density and location. It is found that the distribution of sp3 bonds only along the edges of either pristine or hydrogenated bilayer GS, leads to a greatly enhanced quality factor. A quality factor of ~1.18×105 is observed for a 3.07×15.31 nm2 bilayer GS resonator with sp3 bonds, which is more than 30 times larger comparing with that of a pristine bilayer GS. The present study demonstrates that the resonance properties of a bilayer GS resonator can be tuned by introducing sp3 bonds. This finding provides a useful guideline for the synthesis of the bilayer GS for its application as a resonator component.

Correlating flow induced shear stress and chondrocytes activity in micro-porous scaffold using computational fluid dynamics and rapid prototyping

Flow induced shear stress plays an important role in regulating cell growth and distribution in scaffolds. This study sought to correlate wall shear stress and chondrocytes activity for engineering design of micro-porous osteochondral grafts based on the hypothesis that it is possible to capture and discriminate between the transmitted force and cell response at the inner irregularities. Unlike common tissue engineering therapies with perfusion bioreactors in which flow-mediated stress is the controlling parameter, this work assigned the associated stress as a function of porosity to influence in vitro proliferation of chondrocytes. D-optimality criterion was used to accommodate three pore characteristics for appraisal in a mixed level fractional design of experiment (DOE); namely, pore size (4 levels), distribution pattern (2 levels) and density (3 levels). Micro-porous scaffolds (n=12) were fabricated according to the DOE using rapid prototyping of an acrylic-based bio-photopolymer. Computational fluid dynamics (CFD) models were created correspondingly and used on an idealized boundary condition with a Newtonian fluid domain to simulate the dynamic microenvironment inside the pores. In vitro condition was reproduced for the 3D printed constructs seeded by high pellet densities of human chondrocytes and cultured for 72 hours. The results showed that cell proliferation was significantly different in the constructs (p<0.05). Inlet fluid velocity of 3×10-2mms-1 and average shear stress of 5.65×10-2 Pa corresponded with increased cell proliferation for scaffolds with smaller pores in hexagonal pattern and lower densities. Although the analytical solution of a Poiseuille flow inside the pores was found insufficient for the description of the flow profile probably due to the outside flow induced turbulence, it showed that the shear stress would increase with cell growth and decrease with pore size. This correlation demonstrated the basis for determining the relation between the induced stress and chondrocyte activity to optimize microfabrication of engineered cartilaginous constructs.

Atomistic simulation of interfacial behaviour in graphene-polymer nanocomposites

Graphene–polymer nanocomposites have promising properties as new structural and functional materials. The remarkable mechanical property enhancement in these nanocomposites is generally attributed to exceptional mechanical property of graphene and possible load transfer between graphene and polymer matrix. However, the underlying strengthening and toughening mechanisms have not been well understood. In this work, the interfacial behavior of graphene-polyethylene (PE) was investigated using molecular dynamics (MD) method. The interfacial shear force (ISF) and interfacial shear stress (ISS) between graphene and PE matrix were evaluated, taking into account graphene size, the number of graphene layers and the structural defects in graphene. MD results show that the ISS at graphene-PE interface mainly distributes at each end of the graphene nanofiller within the range of 1 nm, and much larger than that at carbon nanotube (CNT)-PE interface. Moreover, it was found that the ISS at graphene-PE interface is sensitive to the layer number.

Delineating breast cancer cell interactions with engineered bone microenvironments

The mechanisms leading to colonization of metastatic breast cancer cells (BCa) in the skeleton are still not fully understood. Here, we demonstrate that mineralized extracellular matrices secreted by primary human osteoblasts (hOBM) modulate cellular processes associated with BCa colonization of bone. A panel of four BCa cell lines of different bone-metastatic potential (T47D, SUM1315, MDA-MB-231, and the bone-seeking subline MDA-MB-231BO) was cultured on hOBM. After 3 days, the metastatic BCa cells had undergone morphological changes on hOBM and were aligned along the hOBM's collagen type I fibrils that were decorated with bone-specific proteins. In contrast, nonmetastatic BCa cells showed a random orientation on hOBM. Atomic force microscopy-based single-cell force spectroscopy revealed that the metastatic cell lines adhered more strongly to hOBM compared with nonmetastatic cells. Function-blocking experiments indicated that β1-integrins mediated cell adhesion to hOBM. In addition, metastatic BCa cells migrated directionally and invaded hOBM, which was accompanied by enhanced MMP-2 and -9 secretion. Furthermore, we observed gene expression changes associated with osteomimickry in BCa cultured on hOBM. As such, osteopontin mRNA levels were significantly increased in SUM1315 and MDA-MB-231BO cells in a β1-integrin-dependent manner after growing for 3 days on hOBM compared with tissue culture plastic. In conclusion, our results show that extracellular matrices derived from human osteoblasts represent a powerful experimental platform to dissect mechanisms underlying critical steps in the development of bone metastases.

Horror

The monstrous landscape and the revenge of nature are recurring motifs in Australian cinema. In the horror genre, the idea of the monstrous landscape emerges from, and builds upon, an established tradition in Australian cinema in which landscape functions not just as a setting for action, but also as a character in its own right. Rather than a picturesque wilderness or countryside, or a serene natural world untainted by civilisation – representations common in landscape cinema celebrating positive aspects of the Australian ‘outback’ – the monstrous landscape is a dangerous, malevolent and threatening force. Drawing upon themes also common in Australian Gothic narratives such as entrapment in a hostile environment, isolation and fear of the unknown (Turcotte, 1988, see also Jonathan Rayner’s essay in this volume), the monstrous landscape acts according to its own logic indecipherable to non-Indigenous Australians and is represented in terms of its alien-ness and inhuman horror.

Stiffness-damping matching method of an ECAS system based on LQG control

A novel method of matching stiffness and continuous variable damping of an ECAS (electronically controlled air suspension) based on LQG (linear quadratic Gaussian) control was proposed to simultaneously improve the road-friendliness and ride comfort of a two-axle school bus. Taking account of the suspension nonlinearities and target-height-dependent variation in suspension characteristics, a stiffness model of the ECAS mounted on the drive axle of the bus was developed based on thermodynamics and the key parameters were obtained through field tests. By determining the proper range of the target height for the ECAS of the fully-loaded bus based on the design requirements of vehicle body bounce frequency, the control algorithm of the target suspension height (i.e., stiffness) was derived according to driving speed and road roughness. Taking account of the nonlinearities of a continuous variable semi-active damper, the damping force was obtained through the subtraction of the air spring force from the optimum integrated suspension force, which was calculated based on LQG control. Finally, a GA (genetic algorithm)-based matching method between stepped variable damping and stiffness was employed as a benchmark to evaluate the effectiveness of the LQG-based matching method. Simulation results indicate that compared with the GA-based matching method, both dynamic tire force and vehicle body vertical acceleration responses are markedly reduced around the vehicle body bounce frequency employing the LQG-based matching method, with peak values of the dynamic tire force PSD (power spectral density) decreased by 73.6%, 60.8% and 71.9% in the three cases, and corresponding reduction are 71.3%, 59.4% and 68.2% for the vehicle body vertical acceleration. A strong robustness to variation of driving speed and road roughness is also observed for the LQG-based matching method.

Second-order inelastic analysis of composite framed structures based on the refined plastic hinge method

Composite steel-concrete structures experience non-linear effects which arise from both instability-related geometric non-linearity and from material non-linearity in all of their component members. Because of this, conventional design procedures cannot capture the true behaviour of a composite frame throughout its full loading range, and so a procedure to account for those non-linearities is much needed. This paper therefore presents a numerical procedure capable of addressing geometric and material non-linearities at the strength limit state based on the refined plastic hinge method. Different material non-linearity for different composite structural components such as T-beams, concrete-filled tubular (CFT) and steel-encased reinforced concrete (SRC) sections can be treated using a routine numerical procedure for their section properties in this plastic hinge approach. Simple and conservative initial and full yield surfaces for general composite sections are proposed in this paper. The refined plastic hinge approach models springs at the ends of the element which are activated when the surface defining the interaction of bending and axial force at first yield is reached; a transition from the first yield interaction surface to the fully plastic interaction surface is postulated based on a proposed refined spring stiffness, which formulates the load-displacement relation for material non-linearity under the interaction of bending and axial actions. This produces a benign method for a beam-column composite element under general loading cases. Another main feature of this paper is that, for members containing a point of contraflexure, its location is determined with a simple application of the method herein and a node is then located at this position to reproduce the real flexural behaviour and associated material non-linearity of the member. Recourse is made to an updated Lagrangian formulation to consider geometric non-linear behaviour and to develop a non-linear solution strategy. The formulation with the refined plastic hinge approach is efficacious and robust, and so a full frame analysis incorporating geometric and material non-linearity is tractable. By way of contrast, the plastic zone approach possesses the drawback of strain-based procedures which rely on determining plastic zones within a cross-section and which require lengthwise integration. Following development of the theory, its application is illustrated with a number of varied examples.

Higher-order non-linear analysis of steel structures. Part II : refined plastic hinge formulation

In the companion paper, a fourth-order element formulation in an updated Lagrangian formulation was presented to handle geometric non-linearities. The formulation of the present paper extends this to include material non-linearity by proposing a refined plastic hinge approach to analyse large steel framed structures with many members, for which contemporary algorithms based on the plastic zone approach can be problematic computationally. This concept is an advancement of conventional plastic hinge approaches, as the refined plastic hinge technique allows for gradual yielding, being recognized as distributed plasticity across the element section, a condition of full plasticity, as well as including strain hardening. It is founded on interaction yield surfaces specified analytically in terms of force resultants, and achieves accurate and rapid convergence for large frames for which geometric and material non-linearity are significant. The solutions are shown to be efficacious in terms of a balance of accuracy and computational expediency. In addition to the numerical efficiency, the present versatile approach is able to capture different kinds of material and geometric non-linearities on general applications of steel structures, and thereby it offers an efficacious and accurate means of assessing non-linear behaviour of the structures for engineering practice.

Second-order elastic finite element analysis of steel structures using a single element per member

Finite element frame analysis programs targeted for design office application necessitate algorithms which can deliver reliable numerical convergence in a practical timeframe with comparable degrees of accuracy, and a highly desirable attribute is the use of a single element per member to reduce computational storage, as well as data preparation and the interpretation of the results. To this end, a higher-order finite element method including geometric non-linearity is addressed in the paper for the analysis of elastic frames for which a single element is used to model each member. The geometric non-linearity in the structure is handled using an updated Lagrangian formulation, which takes the effects of the large translations and rotations that occur at the joints into consideration by accumulating their nodal coordinates. Rigid body movements are eliminated from the local member load-displacement relationship for which the total secant stiffness is formulated for evaluating the large member deformations of an element. The influences of the axial force on the member stiffness and the changes in the member chord length are taken into account using a modified bowing function which is formulated in the total secant stiffness relationship, for which the coupling of the axial strain and flexural bowing is included. The accuracy and efficiency of the technique is verified by comparisons with a number of plane and spatial structures, whose structural response has been reported in independent studies.

On the advanced analysis of steel frames allowing for flexural, local and lateral-torsional buckling

Detailed procedure for second-order analysis has been coded in the newest Eurocode 3 and the Hong Kong steel code (2005). The effective length method has been noted to be inapplicable to analysis of shallow domes of imperfect members exhibiting snap-through buckling, to portals with leaning columns and others. On the other hand, the advanced analysis is not limited to buckling design of these structures. This paper demonstrates its application to the design of a simple plane sway portal and a three diminsional non-sway steel building. The results by the advanced analysis and the first-order linear analysis are compared and the technique for practical second-order analysis steel structures is described. It is observed that the use of a straight element by itself cannot model the buckling resistance of columns governed by different buckling curves for hot-rolled and cold-formed sections of various shapes like I, H, hollow etc. Also the curvature of the conventional cubic Hermite element is not varied by the external axial force and thus it cannot simulate the response of a buckling column. Thus its use for second-order analysis is basically unacceptable. A technique for additional checking of beams undergoing lateral-torsional buckling is also suggested making the advanced analysis a complete design tool for conventional steel frames.

Material yielding by both axial and bending spring stiffness at elevated temperature

Material yielding is typically modeled either by plastic zone or plastic hinge methods under the context of geometric and material nonlinear finite element methods. In fire analysis of steel structures, the plastic zone method is widely used, but it requires extensively more computational efforts. The objective of this paper is to develop the nonlinear material model allowing for interaction of both axial force and bending moment, which relies on the plastic hinge method to achieve numerical efficiency and reduce computational effort. The biggest advantage of the plastic-hinge approach is its computational efficiency and easy verification by the design code formulae of the axial force–moment interaction yield criterion for beam–column members. Further, the method is reliable and robust when used in analysis of practical and large structures. In order to allow for the effect of catenary action, axial thermal expansion is considered in the axial restraint equations. The yield function for material yielding incorporated in the stiffness formulation, which allows for both axial force and bending moment effects, is more accurate and rational to predict the behaviour of the frames under fire. In the present fire analysis, the mechanical properties at elevated temperatures follow mainly the Eurocode 3 [Design of steel structures, Part 1.2: Structural fire design. European Committee for Standisation; 2003]. Example of a tension member at a steady state heating condition is modeled to verify the proposed spring formulation and to compare with results by others. The behaviour of a heated member in a highly redundant structure is also studied by the present approach.

Turbulence in an inundated urban environment during a major flood : implications in terms of people evacuation and sediment deposition

Floods through inundated urban environments constitute a hazard to the population and infrastructure. A series of field measurements were performed in an inundated section of the City of Brisbane (Australia) during a major flood in January 2011. Using an acoustic Doppler velocimeter (ADV), detailed velocity and suspended sediment concentration measurements were conducted about the peak of the flood. The results are discussed with a focus on the safety of individuals in floodwaters and the sediment deposition during the flood recession. The force of the floodwaters in Gardens Point Road was deemed unsafe for individual evacuation. A comparison with past laboratory results suggested that previous recommendations could be inappropriate and unsafe in real flood flows.

Experimental verification of the modified spring-mass theory of fiber Bragg grating accelerometers using transverse forces

A fiber Bragg grating (FBG) accelerometer using transverse forces is more sensitive than one using axial forces with the same mass of the inertial object, because a barely stretched FBG fixed at its two ends is much more sensitive to transverse forces than axial ones. The spring-mass theory, with the assumption that the axial force changes little during the vibration, cannot accurately predict its sensitivity and resonant frequency in the gravitational direction because the assumption does not hold due to the fact that the FBG is barely prestretched. It was modified but still required experimental verification due to the limitations in the original experiments, such as the (1) friction between the inertial object and shell; (2) errors involved in estimating the time-domain records; (3) limited data; and (4) large interval ∼5 Hz between the tested frequencies in the frequency-response experiments. The experiments presented here have verified the modified theory by overcoming those limitations. On the frequency responses, it is observed that the optimal condition for simultaneously achieving high sensitivity and resonant frequency is at the infinitesimal prestretch. On the sensitivity at the same frequency, the experimental sensitivities of the FBG accelerometer with a 5.71 gram inertial object at 6 Hz (1.29, 1.19, 0.88, 0.64, and 0.31 nm/g at the 0.03, 0.69, 1.41, 1.93, and 3.16 nm prestretches, respectively) agree with the static sensitivities predicted (1.25, 1.14, 0.83, 0.61, and 0.29 nm/g, correspondingly). On the resonant frequency, (1) its assumption that the resonant frequencies in the forced and free vibrations are similar is experimentally verified; (2) its dependence on the distance between the FBG’s fixed ends is examined, showing it to be independent; (3) the predictions of the spring-mass theory and modified theory are compared with the experimental results, showing that the modified theory predicts more accurately. The modified theory can be used more confidently in guiding its design by predicting its static sensitivity and resonant frequency, and may have applications in other fields for the scenario where the spring-mass theory fails.